Microsatellite evolution--a reciprocal study of repeat lengths at homologous loci in cattle and sheep

Fine-scale genetic structure in the critically endangered red-fronted macaw in the absence of geographic and ecological barriers

Behavioural and socio-cultural traits are recognized in the restriction of gene flow in species with high cognitive capacity and complex societies. This isolation by social barriers has been generally overlooked in threatened species by assuming disrupted gene flow due to population fragmentation and decline. We examine the genetic structure and ecology of the global population of the Critically Endangered red-fronted macaw (Ara rubrogenys), an endemic species to the inter-Andean valleys of Bolivia. We found a fine-scale genetic structuring in four genetic clusters. Genetic diversity was higher in wild compared to captive-bred macaws, but similar to that of captive wild-caught macaws. We found no clear evidence of severe genetic erosion in the population in recent decades, but it was patent in historic times, overlapping with drastic human habitat transformation and macaw persecution over millennia. We found no evidence of geographical and ecological barriers, owing to the high dispersal ability, nesting and foraging habits between genetic clusters. The lack of genetic intermixing despite long-distance foraging and seasonal movements suggests recruitment in natal colonies and other social factors reinforcing philopatry-related genetic structure. Conservation efforts should be specifically focussed on major threats in each genetic cluster as independent conservation units, and also considered in ex-situ management.

De-novo transcriptome analysis unveils differentially expressed genes regulating drought and salt stress response in Panicum sumatrense

Screening the transcriptome of drought tolerant variety of little millet (Panicum sumatrense), a marginally cultivated, nutritionally rich, susbsistent crop, can identify genes responsible for its hardiness and enable identification of new sources of genetic variation which can be used for crop improvement. RNA-Seq generated ~ 230 million reads from control and treated tissues, which were assembled into 86,614 unigenes. In silico differential gene expression analysis created an overview of patterns of gene expression during exposure to drought and salt stress. Separate gene expression profiles for leaf and root tissue revealed the differences in regulatory mechanisms operating in these tissues during exposure to abiotic stress. Several transcription factors were identified and studied for differential expression. 61 differentially expressed genes were found to be common to both tissues under drought and salinity stress and were further validated using qRT-PCR. Transcriptome of P. sumatrense was also used to mine for genic SSR markers relevant to abiotic stress tolerance. This study is first report on a detailed analysis of molecular mechanisms of drought and salinity stress tolerance in a little millet variety. Resources generated in this study can be used as potential candidates for further characterization and to improve abiotic stress tolerance in food crops.

The genetic structure of the spectacled bear (Tremarctos ornatus; Ursidae, Carnivora) in Colombia by means of mitochondrial and microsatellite markers

The spectacled bear (Ursidae: Tremarctos ornatus) is an emblematic umbrella species and one of the top carnivores in the Andean mountains. It is also listed as vulnerable by IUCN and as endangered by CITES. We analyzed the genetic structure of this species in nine geographical regions representing the three Andean Cordilleras in Colombia. We sequenced six mitochondrial genes in 115 spectacled bears; a subset of these specimens (n = 61) were genotyped at seven nuclear microsatellites. We addressed three objectives: 1) determine the genetic diversity and historical demographic changes of the spectacled bear in Colombia; 2) determine phylogeographic patterns of genetic divergence among spectacled bear populations in Colombia; and 3) estimate the levels of gene flow among different regions of Colombia. Our analyses show evidence of high mitochondrial genetic diversity in spectacled bears, both in Colombia as well as in each of the nine regions, most particularly Norte de Santander, Nariño, and Antioquia-Córdoba. In addition, we detected population expansion in Colombia that occurred around 24,000 years ago, followed by a population decrease during the last 7,000 years, and a sudden expansion in the last 300 years. Phylogenetic analyses showed few well-supported clades, with some haplotypes detected in all the departments and Colombian Andean Cordilleras, and other haplotypes restricted to certain geographical areas (Antioquia, Norte de Santander, Cundinamarca, and Nariño). We detected significant genetic heterogeneity among some departments and among the three Colombian Andean Cordilleras for both mitochondrial and nuclear genes. Nevertheless, the moderate levels of gene flow estimated from FST statistics suggest that geographical barriers have not been definitive obstacles to the dispersion of the spectacled bear throughout Colombia. Despite these gene flow estimates, significant spatial autocorrelation was detected for spectacled bear in Colombia, where two kinds of spatial patterns were discovered: genetic patches of 144 km of diameter, and isolation by distance among bears separated from 578 to 800 km. The two most northern spectacled bear populations of Colombia (Norte de Santander and Antioquia) also were the two most differentiated. Their distinctiveness may qualify them as distinct Management Units (MUs) in the context of conservation policies for the spectacled bear in Colombia.

Genetics of the Andean bear (Tremarctos ornatus; Ursidae, Carnivora) in Ecuador: when the Andean Cordilleras are not an Obstacle

One of the top carnivores in the Andean mountains is the Andean bear (Tremarctos ornatus, Ursidae), the only bear in South America. This is a flagship and key umbrella species in Ecuador because its conservation has a positive impact on the conservation of many other species in the Andes. But to preserve, first one must know the genetic characteristics of a species, among other things. For this, we analyzed six mitochondrial genes and seven nuclear DNA microsatellites of 108 Andean bear specimens sampled throughout Ecuador. We adopted three strategies for analyzing the data: by Province, by Region (north vs south), and by Cordillera. Four main results were obtained. First, the mitochondrial genetic diversity levels were elevated, but there were no differences in genetic diversity by Province or by Cordillera. By Regions, southern Ecuador had higher genetic diversity levels than to northern Ecuador. The genetic diversity for the microsatellites was only medium for the Andean bear at this country. Second, there was clear and significant evidence of female population expansions, for the overall sample, by Province, Region, and Cordillera. This population expansion was determined to have occurred in the time interval of 30,000-20,000 years ago (YA), during the last phase of the Pleistocene. We detected a population decrease to have occurred more recently, within the last 5000 years. It continued until about 300-200 YA when a population increase was again detected. Third, there were, practically, no phylogeographic pattern nor genetic differentiation among Andean bear populations in Ecuador by Province or by Cordillera for either mitochondrial or microsatellite markers. There was a little more genetic differentiation between northern and southern areas. Fourth, there was no trace of significant spatial genetic structure for the Andean bear in Ecuador in agreement with the genetic differentiation analyses. This shows that the Andean Cordilleras in this country did not present an obstacle to the dispersion of this species. Therefore, all of the Andean bear specimens in Ecuador should be treated as a unique Management Unit (MU) for conservation purposes, differently to that determined for other countries as Colombia.

Genomic diversity and selection sweeps identified in Indian swamp buffaloes reveals it's uniqueness with riverine buffaloes

The present investigation was focused to study genomic diversity of Indian swamp buffalo populations through reduced representation approach (ddRAD). The heterozygosity (F_ST) among the swamp buffaloes was 0.11 between Assam and Manipuri; 0.20 between swamp (Manipuri) and riverine buffaloes; 0.30 between swamp (Manipuri) and cattle. The average observed and expected heterozygosity in swamp buffalo populations was 0.254 and 0.221 respectively. The Inbreeding coefficient (F_IS) value was 0.02 among the swamp buffaloes. PCA and structure analysis revealed Manipuri swamp buffalo was genetically distinct and closely related to Nagaland swamp buffalo and least to Assam swamp buffalo. Identification of selective sweeps revealed 1087 regions to have undergone selection related to immune response, adaptation and nervous system. A total of 3451 SSRs were identified in the genome of swamp buffaloes. The study evidenced the genomic diversity in the swamp buffalo populations and its uniqueness in comparison with riverine buffalo and cattle.

Development of microsatellite loci for two New World vultures (Cathartidae)

Objective

Use next-generation sequencing to develop microsatellite loci that will provide the variability necessary for studies of genetic diversity and population connectivity of two New World vulture species.

Results

We characterized 11 microsatellite loci for black vultures (Coragyps atratus) and 14 loci for turkey vultures (Cathartes aura). These microsatellite loci were grouped into 3 multiplex panels for each species. The number of alleles among black vulture samples ranged from 2 to 11, and 3 to 48 among turkey vulture samples.

Electronic supplementary material

The online version of this article (10.1186/s13104-019-4295-z) contains supplementary material, which is available to authorized users.

Geo-referenced population-specific microsatellite data across American continents, the MacroPopGen Database

Population genetic data from nuclear DNA has yet to be synthesized to allow broad scale comparisons of intraspecific diversity versus species diversity. The MacroPopGen database collates and geo-references vertebrate population genetic data across the Americas from 1,308 nuclear microsatellite DNA studies, 897 species, and 9,090 genetically distinct populations where genetic differentiation (FST) was measured. Caribbean populations were particularly distinguished from North, Central, and South American populations, in having higher differentiation (FST = 0.12 vs. 0.07-0.09) and lower mean numbers of alleles (MNA = 4.11 vs. 4.84-5.54). While mammalian populations had lower MNA (4.86) than anadromous fish, reptiles, amphibians, freshwater fish, and birds (5.34-7.81), mean heterozygosity was largely similar across groups (0.57-0.63). Mean FST was consistently lowest in anadromous fishes (0.06) and birds (0.05) relative to all other groups (0.09-0.11). Significant differences in Family/Genera variance among continental regions or taxonomic groups were also observed. MacroPopGen can be used in many future applications including latitudinal analyses, spatial analyses (e.g. central-margin), taxonomic comparisons, regional assessments of anthropogenic impacts on biodiversity, and conservation of wild populations.

Genome-wide characterization of simple sequence repeats in Pyrus bretschneideri and their application in an analysis of genetic diversity in pear

Background

Pear (Pyrus spp.) is an economically important temperate fruit tree worldwide. In the past decade, significant progress has been made in pear molecular genetics based on DNA research, but the number of molecular markers is still quite limited, which hardly satisfies the increasing needs of geneticists and breeders.

Results

In this study, a total of 156,396 simple sequence repeat (SSR) loci were identified from a genome sequence of Pyrus bretschneideri ‘Dangshansuli’. A total of 101,694 pairs of SSR primers were designed from the SSR loci, and 80,415 of the SSR loci were successfully located on 17 linkage groups (LGs). A total of 534 primer pairs were synthesized and preliminarily screened in four pear cultivars, and of these, 332 primer pairs were selected as clear, stable, and polymorphic SSR markers. Eighteen polymorphic SSR markers were randomly selected from the 332 polymorphic SSR markers in order to perform a further analysis of the genetic diversity among 44 pear cultivars. The 14 European pears and their hybrid materials were clustered into one group (European pear group); 29 Asian pear cultivars were clustered into one group (Asian pear group); and the Zangli pear cultivar ‘Deqinli’ from Yunnan Province, China, was grouped in an independent group, which suggested that the cultivar ‘Deqinli’ is a distinct and valuable germplasm resource. The population structure analysis partitioned the 44 cultivars into two populations, Pop 1 and Pop 2. Pop 2 was further divided into two subpopulations. Results from the population structure analysis were generally consistent with the results from the UPGMA cluster analysis.

Conclusions

The results of the present study showed that the use of next-generating sequencing to develop SSR markers is fast and effective, and the developed SSR markers can be utilized by researchers and breeders for future pear improvement.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4822-7) contains supplementary material, which is available to authorized users.

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.

Genetic signature of strong recent positive selection at interleukin-32 gene in goat

OBJECTIVE

Identification of the candidate genes that play key roles in phenotypic variations can provide new information about evolution and positive selection. Interleukin (IL)-32 is involved in many biological processes, however, its role for the immune response against various diseases in mammals is poorly understood. Therefore, the current investigation was performed for the better understanding of the molecular evolution and the positive selection of single nucleotide polymorphisms in IL-32 gene.

METHODS

By using fixation index (F_ST ) based method, IL-32 (9375) gene was found to be outlier and under significant positive selection with the provisional combined allocation of mean heterozygosity and F_ST . Using nucleotide sequences of 11 mammalian species from National Center for Biotechnology Information database, the evolutionary selection of IL-32 gene was determined using Maximum likelihood model method, through four models (M1a, M2a, M7, and M8) in Codeml program of phylogenetic analysis by maximum liklihood.

RESULTS

IL-32 is detected under positive selection using the F_ST simulations method. The phylogenetic tree revealed that goat IL-32 was in close resemblance with sheep IL-32. The coding nucleotide sequences were compared among 11 species and it was found that the goat IL-32 gene shared identity with sheep (96.54%), bison (91.97%), camel (58.39%), cat (56.59%), buffalo (56.50%), human (56.13%), dog (50.97%), horse (54.04%), and rabbit (53.41%) respectively.

CONCLUSION

This study provides evidence for IL-32 gene as under significant positive selection in goat.

How do SNP ascertainment schemes and population demographics affect inferences about population history?

Background

The selection of variable sites for inclusion in genomic analyses can influence results, especially when exemplar populations are used to determine polymorphic sites. We tested the impact of ascertainment bias on the inference of population genetic parameters using empirical and simulated data representing the three major continental groups of cattle: European, African, and Indian. We simulated data under three demographic models. Each simulated data set was subjected to three ascertainment schemes: (I) random selection; (II) geographically biased selection; and (III) selection biased toward loci polymorphic in multiple groups. Empirical data comprised samples of 25 individuals representing each continental group. These cattle were genotyped for 47,506 loci from the bovine 50 K SNP panel. We compared the inference of population histories for the empirical and simulated data sets across different ascertainment conditions using F_ST and principal components analysis (PCA).

Results

Bias toward shared polymorphism across continental groups is apparent in the empirical SNP data. Bias toward uneven levels of within-group polymorphism decreases estimates of F_ST between groups. Subpopulation-biased selection of SNPs changes the weighting of principal component axes and can affect inferences about proportions of admixture and population histories using PCA. PCA-based inferences of population relationships are largely congruent across types of ascertainment bias, even when ascertainment bias is strong.

Conclusions

Analyses of ascertainment bias in genomic data have largely been conducted on human data. As genomic analyses are being applied to non-model organisms, and across taxa with deeper divergences, care must be taken to consider the potential for bias in ascertainment of variation to affect inferences. Estimates of F_ST, time of separation, and population divergence as estimated by principal components analysis can be misleading if this bias is not taken into account.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1469-5) contains supplementary material, which is available to authorized users.

Cross-priming of microsatellite loci in subfamily cyprininae (family Cyprinidae): their utility in finding markers for population genetic analysis in three Indian major carps

This study is aimed to identify polymorphic microsatellite markers and establish their potential for population genetics studies in three carp (family cyprinidae; subfamily cyprininae) species, Labeo rohita, Catla catla and Cirrhinus mrigala through use of cyprinid primers. These species have high commercial value and knowledge of genetic variation is important for management of farmed and wild populations. We tested 108 microsatellite primers from 11 species belonging to three different cyprinid subfamilies, Cyprininae, Barbinae and Leuciscinae out of which 63 primers (58.33%) successfully amplified orthologous loci in three focal species. Forty-two loci generated from 29 primers were polymorphic in these three carp species. Sequencing of amplified product confirmed the presence of SSRs in these 42 loci and orthologous nature of the loci. To validate potential of these 42 polymorphic loci in determining the genetic variation, we analyzed 486 samples of three focal species collected from Indus, Ganges and Brahmaputra river systems. Results indicated significant genetic variation, with mean number of alleles per locus ranging from 6.80 to 14.40 and observed heterozygosity ranging from 0.50 to 0.74 in the three focal species. Highly significant (P < 0.00001) allelic homogeneity values revealed that the identified loci can be efficiently used in population genetics analysis of these carp species. Further, thirty-two loci from 19 primers were useful for genotyping in more than one species. The data from the present study was compiled with cross-species amplification data from previous results on eight species of subfamily cyprininae to compare cross-transferability of microsatellite loci. It was revealed that out of 226 heterologous loci amplified, 152 loci that originated from 77 loci exhibited polymorphism and 45 primers were of multispecies utility, common for 2-7 species.

Anagenetic speciation in Ullung Island, Korea: genetic diversity and structure in the island endemic species, Acer takesimense (Sapindaceae)

Anagenetic speciation is an important mode of speciation in oceanic islands; one-fourth of the endemic plants are estimated to have been derived via this process. Few studies, however, have critically examined the genetic consequences of anagenesis in comparison with cladogenesis (involved with adaptive radiation). We hypothesize that endemic species originating via anagenetic speciation in a relatively uniform environment should accumulate genetic variation with limited populational differentiation. We undertook a population genetic analysis using nine nuclear microsatellite loci of Acer takesimense, an anagenetically derived species endemic to Ullung Island, Korea, and its continental progenitor A. pseudosieboldianum on the Korean Peninsula. Microsatellite data reveal a clear genetic distinction between the two species. A high F value in the cluster of A. takesimense was found by Bayesian clustering analysis, suggesting a strong episode of genetic drift during colonization and speciation. In comparison with A. pseudosieboldianum, A. takesimense has slightly lower genetic diversity and possesses less than half the number of private and rare alleles. Consistent with predictions, weak geographical genetic structure within the island was found in A. takesimense. These results imply that anagenetic speciation leads to a different pattern of specific and genetic diversity than often seen with cladogenesis.

Anagenetic speciation in Ullung Island, Korea: genetic diversity and structure in the island endemic species, Acer takesimense (Sapindaceae)

Anagenetic speciation is an important mode of speciation in oceanic islands; one-fourth of the endemic plants are estimated to have been derived via this process. Few studies, however, have critically examined the genetic consequences of anagenesis in comparison with cladogenesis (involved with adaptive radiation). We hypothesize that endemic species originating via anagenetic speciation in a relatively uniform environment should accumulate genetic variation with limited populational differentiation. We undertook a population genetic analysis using nine nuclear microsatellite loci of Acer takesimense, an anagenetically derived species endemic to Ullung Island, Korea, and its continental progenitor A. pseudosieboldianum on the Korean Peninsula. Microsatellite data reveal a clear genetic distinction between the two species. A high F value in the cluster of A. takesimense was found by Bayesian clustering analysis, suggesting a strong episode of genetic drift during colonization and speciation. In comparison with A. pseudosieboldianum, A. takesimense has slightly lower genetic diversity and possesses less than half the number of private and rare alleles. Consistent with predictions, weak geographical genetic structure within the island was found in A. takesimense. These results imply that anagenetic speciation leads to a different pattern of specific and genetic diversity than often seen with cladogenesis.

Genetic structure and diversity of coffee (Coffea) across Africa and the Indian Ocean islands revealed using microsatellites

BACKGROUND AND AIMS

The coffee genus (Coffea) comprises 124 species, and is indigenous to the Old World Tropics. Due to its immense economic importance, Coffea has been the focus of numerous genetic diversity studies, but despite this effort it remains insufficiently studied. In this study the genetic diversity and genetic structure of Coffea across Africa and the Indian Ocean islands is investigated.

METHODS

Genetic data were produced using 13 polymorphic nuclear microsatellite markers (simple sequence repeats, SSRs), including seven expressed sequence tag-SSRs, and the data were analysed using model- and non-model-based methods. The study includes a total of 728 individuals from 60 species.

KEY RESULTS

Across Africa and the Indian Ocean islands Coffea comprises a closely related group of species with an overall pattern of genotypes running from west to east. Genetic structure was identified in accordance with pre-determined geographical regions and phylogenetic groups. There is a good relationship between morpho-taxonomic species delimitations and genetic units. Genetic diversity in African and Indian Ocean Coffea is high in terms of number of alleles detected, and Madagascar appears to represent a place of significant diversification in terms of allelic richness and species diversity.

CONCLUSIONS

Cross-species SSR transferability in African and Indian Ocean islands Coffea was very efficient. On the basis of the number of private alleles, diversification in East Africa and the Indian Ocean islands appears to be more recent than in West and West-Central Africa, although this general trend is complicated in Africa by the position of species belonging to lineages connecting the main geographical regions. The general pattern of phylogeography is not in agreement with an overall east to west (Mascarene, Madagascar, East Africa, West Africa) increase in genome size, the high proportion of shared alleles between the four regions or the high numbers of exclusive shared alleles between pairs or triplets of regions.

Evidence of stable genetic structure across a remote island archipelago through self-recruitment in a widely dispersed coral reef fish

We used microsatellite markers to assess the population genetic structure of the scribbled rabbitfish Siganus spinus in the western Pacific. This species is a culturally important food fish in the Mariana Archipelago and subject to high fishing pressure. Our primary hypothesis was to test whether the individuals resident in the southern Mariana Island chain were genetically distinct and hence should be managed as discrete stocks. In addition to spatial sampling of adults, newly-settled individuals were sampled on Guam over four recruitment events to assess the temporal stability of the observed spatial patterns, and evidence of self-recruitment. We found significant genetic structure in S. spinus across the western Pacific, with Bayesian analyses revealing three genetically distinct clusters: the southern Mariana Islands, east Micronesia, and the west Pacific; with the southern Mariana Islands being more strongly differentiated from the rest of the region. Analyses of temporal samples from Guam indicated the southern Mariana cluster was stable over time, with no genetic differentiation between adults versus recruits, or between samples collected across four separate recruitment events spanning 11 months. Subsequent assignment tests indicated seven recruits had self-recruited from within the Southern Mariana Islands population. Our results confirm the relative isolation of the southern Mariana Islands population and highlight how local processes can act to isolate populations that, by virtue of their broad-scale distribution, have been subject to traditionally high gene flows. Our results add to a growing consensus that self-recruitment is a highly significant influence on the population dynamics of tropical reef fish.

Genome-wide characterization of simple sequence repeats in cucumber (Cucumis sativus L.)

Background

Cucumber, Cucumis sativus L. is an important vegetable crop worldwide. Until very recently, cucumber genetic and genomic resources, especially molecular markers, have been very limited, impeding progress of cucumber breeding efforts. Microsatellites are short tandemly repeated DNA sequences, which are frequently favored as genetic markers due to their high level of polymorphism and codominant inheritance. Data from previously characterized genomes has shown that these repeats vary in frequency, motif sequence, and genomic location across taxa. During the last year, the genomes of two cucumber genotypes were sequenced including the Chinese fresh market type inbred line '9930' and the North American pickling type inbred line 'Gy14'. These sequences provide a powerful tool for developing markers in a large scale. In this study, we surveyed and characterized the distribution and frequency of perfect microsatellites in 203 Mbp assembled Gy14 DNA sequences, representing 55% of its nuclear genome, and in cucumber EST sequences. Similar analyses were performed in genomic and EST data from seven other plant species, and the results were compared with those of cucumber.

Results

A total of 112,073 perfect repeats were detected in the Gy14 cucumber genome sequence, accounting for 0.9% of the assembled Gy14 genome, with an overall density of 551.9 SSRs/Mbp. While tetranucleotides were the most frequent microsatellites in genomic DNA sequence, dinucleotide repeats, which had more repeat units than any other SSR type, had the highest cumulative sequence length. Coding regions (ESTs) of the cucumber genome had fewer microsatellites compared to its genomic sequence, with trinucleotides predominating in EST sequences. AAG was the most frequent repeat in cucumber ESTs. Overall, AT-rich motifs prevailed in both genomic and EST data. Compared to the other species examined, cucumber genomic sequence had the highest density of SSRs (although comparable to the density of poplar, grapevine and rice), and was richest in AT dinucleotides. Using an electronic PCR strategy, we investigated the polymorphism between 9930 and Gy14 at 1,006 SSR loci, and found unexpectedly high degree of polymorphism (48.3%) between the two genotypes. The level of polymorphism seems to be positively associated with the number of repeat units in the microsatellite. The in silico PCR results were validated empirically in 660 of the 1,006 SSR loci. In addition, primer sequences for more than 83,000 newly-discovered cucumber microsatellites, and their exact positions in the Gy14 genome assembly were made publicly available.

Conclusions

The cucumber genome is rich in microsatellites; AT and AAG are the most abundant repeat motifs in genomic and EST sequences of cucumber, respectively. Considering all the species investigated, some commonalities were noted, especially within the monocot and dicot groups, although the distribution of motifs and the frequency of certain repeats were characteristic of the species examined. The large number of SSR markers developed from this study should be a significant contribution to the cucurbit research community.

Utility of sequenced genomes for microsatellite marker development in non-model organisms: a case study of functionally important genes in nine-spined sticklebacks (Pungitius pungitius)

Background

Identification of genes involved in adaptation and speciation by targeting specific genes of interest has become a plausible strategy also for non-model organisms. We investigated the potential utility of available sequenced fish genomes to develop microsatellite (cf. simple sequence repeat, SSR) markers for functionally important genes in nine-spined sticklebacks (Pungitius pungitius), as well as cross-species transferability of SSR primers from three-spined (Gasterosteus aculeatus) to nine-spined sticklebacks. In addition, we examined the patterns and degree of SSR conservation between these species using their aligned sequences.

Results

Cross-species amplification success was lower for SSR markers located in or around functionally important genes (27 out of 158) than for those randomly derived from genomic (35 out of 101) and cDNA (35 out of 87) libraries. Polymorphism was observed at a large proportion (65%) of the cross-amplified loci independently of SSR type. To develop SSR markers for functionally important genes in nine-spined sticklebacks, SSR locations were surveyed in or around 67 target genes based on the three-spined stickleback genome and these regions were sequenced with primers designed from conserved sequences in sequenced fish genomes. Out of the 81 SSRs identified in the sequenced regions (44,084 bp), 57 exhibited the same motifs at the same locations as in the three-spined stickleback. Di- and trinucleotide SSRs appeared to be highly conserved whereas mononucleotide SSRs were less so. Species-specific primers were designed to amplify 58 SSRs using the sequences of nine-spined sticklebacks.

Conclusions

Our results demonstrated that a large proportion of SSRs are conserved in the species that have diverged more than 10 million years ago. Therefore, the three-spined stickleback genome can be used to predict SSR locations in the nine-spined stickleback genome. While cross-species utility of SSR primers is limited due to low amplification success, SSR markers can be developed for target genes and genomic regions using our approach, which should be also applicable to other non-model organisms. The SSR markers developed in this study should be useful for identification of genes responsible for phenotypic variation and adaptive divergence of nine-spined stickleback populations, as well as for constructing comparative gene maps of nine-spined and three-spined sticklebacks.

First report on the presence of fire blight resistance in linkage group 11 of Pyrus ussuriensis Maxim

Fire blight, caused by the gram-negative bacterium Erwinia amylovora (Burrill) Winslow et al., is a dangerous disease of pome fruits, including pear. A pear breeding program for fire blight resistance was initiated in 2003 at the Department of Pomology, Warsaw University of Life Sciences, Poland. Since several Asian species are considered to be potential sources of resistance to fire blight, the susceptible Pyrus communis 'Doyenne du Comice' was crossed with the resistant P. ussuriensis. The F1 full-sib progeny composed of 155 seedlings was tested for susceptibility to fire blight by artificial shoot inoculation. A framework linkage map of both parents was constructed based on 48 AFLP and 32 SSR markers and covered a length of 595 cM and 680 cM in 'Doyenne du Comice' and P. ussuriensis, respectively. For the first time a putative QTL for fire blight resistance in P. ussuriensis linkage group 11 was identified. Another putative QTL in linkage group 4 of 'Doyenne du Comice' seems to indicate that sources of fire blight resistance can be identified also in the susceptible cultivars.

Development of chickpea EST-SSR markers and analysis of allelic variation across related species

Despite chickpea being the third important grain legume, there is a limited availability of genomic resources, especially of the expressed sequence tag (EST)-based markers. In this study, we generated 822 chickpea ESTs from immature seeds as well as exploited 1,309 ESTs from the chickpea database, thus utilizing a total of 2,131 EST sequences for development of functional EST-SSR markers. Two hundred and forty-six simple sequence repeat (SSR) motifs were identified from which 183 primer pairs were designed and 60 validated as functional markers. Genetic diversity analysis across 30 chickpea accessions revealed ten markers to be polymorphic producing a total of 29 alleles and an observed heterozygosity average of 0.16 thereby exhibiting low levels of intra-specific polymorphism. However, the markers exhibited high cross-species transferability ranging from 68.3 to 96.6% across the six annual Cicer species and from 29.4 to 61.7% across the seven legume genera. Sequence analysis of size variant amplicons from various species revealed that size polymorphism was due to multiple events such as copy number variation, point mutations and insertions/deletions in the microsatellite repeat as well as in the flanking regions. Interestingly, a wide prevalence of crossability-group-specific sequence variations were observed among Cicer species that were phylogenetically informative. The neighbor joining dendrogram clearly separated the chickpea cultivars from the wild Cicer and validated the proximity of C. judaicum with C. pinnatifidum. Hence, this study for the first time provides an insight into the distribution of SSRs in the chickpea transcribed regions and also demonstrates the development and utilization of genic-SSRs. In addition to proving their suitability for genetic diversity analysis, their high rates of transferability also proved their potential for comparative genomic studies and for following gene introgressions and evolution in wild species, which constitute the valuable secondary genepool in chickpea.

Gene flow and population subdivision in a pantropical plant with sea-drifted seeds Hibiscus tiliaceus and its allied species: evidence from microsatellite analyses

The genetic differentiation and structure of Hibiscus tiliaceus, a pantropical plant with sea-drifted seeds, and four allied species were studied using six microsatellite markers. A low level of genetic differentiation was observed among H. tiliaceus populations in the Pacific and Indian Ocean regions, similar to the results of a previous chloroplast DNA (cpDNA) study. Frequent gene flow by long-distance seed dispersal is responsible for species integration of H. tiliaceus in the wide distribution range. On the other hand, highly differentiated populations of H. tiliaceus were detected in West Africa, as well as of Hibiscus pernambucensis in southern Brazil. In the former populations, the African continent may be a geographical barrier that prevents gene flow by sea-drifted seeds. In the latter populations, although there are no known land barriers, the bifurcating South Equatorial Current at the north-eastern horn of Brazil can be a potential barrier to gene flow and may promote the genetic differentiation of these populations. Our results also suggest clear species segregation between H. tiliaceus and H. pernambucensis, which confirms the introgression scenario between these two species that was suggested by a previous cpDNA study. Our results also provide good evidence for recent transatlantic long-distance seed dispersal by sea current. Despite the distinct geographical structure observed in the cpDNA haplotypes, a low level of genetic differentiation was found between Pacific and Atlantic populations of H. pernambucensis, which could be caused by transisthmian gene flow.

Genome-wide analysis of microsatellite polymorphism in chicken circumventing the ascertainment bias

Studies of microsatellites evolution based on marker data almost inherently suffer from an ascertainment bias because there is selection for the most mutable and polymorphic loci during marker development. To circumvent this bias we took advantage of whole-genome shotgun sequence data from three unrelated chicken individuals that, when aligned to the genome reference sequence, give sequence information on two chromosomes from about one-fourth (375,000) of all microsatellite loci containing di- through pentanucleotide repeat motifs in the chicken genome. Polymorphism is seen at loci with as few as five repeat units, and the proportion of dimorphic loci then increases to 50% for sequences with approximately 10 repeat units, to reach a maximum of 75%-80% for sequences with 15 or more repeat units. For any given repeat length, polymorphism increases with decreasing GC content of repeat motifs for dinucleotides, nonhairpin-forming trinucleotides, and tetranucleotides. For trinucleotide repeats which are likely to form hairpin structures, polymorphism increases with increasing GC content, indicating that the relative stability of hairpins affects the rate of replication slippage. For any given repeat length, polymorphism is significantly lower for imperfect compared to perfect repeats and repeat interruptions occur in >15% of loci. However, interruptions are not randomly distributed within repeat arrays but are preferentially located toward the ends. There is negative correlation between microsatellite abundance and single nucleotide polymorphism (SNP) density, providing large-scale genomic support for the hypothesis that equilibrium microsatellite distributions are governed by a balance between rate of replication slippage and rate of point mutation.

Genomic conservation of cattle microsatellite loci in wild gaur (Bos gaurus) and current genetic status of this species in Vietnam

BACKGROUND

The wild gaur (Bos gaurus) is an endangered wild cattle species. In Vietnam, the total number of wild gaurs is estimated at a maximum of 500 individuals. Inbreeding and genetic drift are current relevant threats to this small population size. Therefore, information about the genetic status of the Vietnamese wild gaur population is essential to develop strategies for conservation and effective long-term management for this species. In the present study, we performed cross-species amplification of 130 bovine microsatellite markers, in order to evaluate the applicability and conservation of cattle microsatellite loci in the wild gaur genome. The genetic diversity of Vietnamese wild gaur was also investigated, based on data collected from the 117 successfully amplified loci.

RESULTS

One hundred-thirty cattle microsatellite markers were tested on a panel of 11 animals. Efficient amplifications were observed for 117 markers (90%) with a total of 264 alleles, and of these, 68 (58.1%) gave polymorphic band patterns. The number of alleles per locus among the polymorphic markers ranged from two to six. Thirteen loci (BM1314, BM2304, BM6017, BMC2228, BMS332, BMS911, CSSM023, ETH123, HAUT14, HEL11, HEL5, ILSTS005 and INRA189) distributed on nine different cattle chromosomes failed to amplify wild gaur genomic DNA. Three cattle Y-chromosome specific microsatellite markers (INRA124, INRA126 and BM861) were also highly specific in wild gaur, only displaying an amplification product in the males. Genotype data collected from the 117 successfully amplified microsatellites were used to assess the genetic diversity of this species in Vietnam. Polymorphic Information Content (PIC) values varied between 0.083 and 0.767 with a mean of 0.252 while observed heterozygosities (Ho) ranged from 0.091 to 0.909 (mean of 0.269). Nei's unbiased mean heterozygosity and the mean allele number across loci were 0.298 and 2.2, respectively.

CONCLUSION

Extensive conservation of cattle microsatellite loci in the wild gaur genome, as shown by our results, indicated a high applicability of bovine microsatellites for genetic characterization and population genetic studies of this species. Moreover, the low genetic diversity observed in Vietnamese wild gaur further underlines the necessity of specific strategies and appropriate management plans to preserve this endangered species from extinction.

Historical and ecological determinants of genetic structure in arctic canids

Wolves (Canis lupus) and arctic foxes (Alopex lagopus) are the only canid species found throughout the mainland tundra and arctic islands of North America. Contrasting evolutionary histories, and the contemporary ecology of each species, have combined to produce their divergent population genetic characteristics. Arctic foxes are more variable than wolves, and both island and mainland fox populations possess similarly high microsatellite variation. These differences result from larger effective population sizes in arctic foxes, and the fact that, unlike wolves, foxes were not isolated in discrete refugia during the Pleistocene. Despite the large physical distances and distinct ecotypes represented, a single, panmictic population of arctic foxes was found which spans the Svalbard Archipelago and the North American range of the species. This pattern likely reflects both the absence of historical population bottlenecks and current, high levels of gene flow following frequent long-distance foraging movements. In contrast, genetic structure in wolves correlates strongly to transitions in habitat type, and is probably determined by natal habitat-biased dispersal. Nonrandom dispersal may be cued by relative levels of vegetation cover between tundra and forest habitats, but especially by wolf prey specialization on ungulate species of familiar type and behaviour (sedentary or migratory). Results presented here suggest that, through its influence on sea ice, vegetation, prey dynamics and distribution, continued arctic climate change may have effects as dramatic as those of the Pleistocene on the genetic structure of arctic canid species.

Effects of flowering tree density on the mating system and gene flow in Shorea leprosula (Dipterocarpaceae) in Peninsular Malaysia

Pristine tropical rainforests in Southeast Asia have rich species diversity and are important habitats for many plant species. However, the extent of these forests has declined in recent decades and they have become fragmented due to human activities. These developments may reduce the genetic diversity of species within them and, consequently, the species' ability to adapt to environmental changes. Our objective in the study presented here was to clarify the effect of tree density on the genetic diversity and gene flow patterns of Shorea leprosula Miq. populations in Peninsular Malaysia. For this purpose, we related genetic diversity and pollen flow parameters of seedling populations in study plots to the density of mature trees in their vicinity. The results show that gene diversity and allelic richness were not significantly correlated to the mature tree density. However, the number of rare alleles among the seedlings and the selfing rates of the mother trees were negatively correlated with the density of the adult trees. Furthermore, in a population with high mature tree density pollination distances were frequently <200 m, but in populations with low adult tree density the distances were longer. These findings suggest that the density of flowering trees affects selfing rates, gene flow and, thus, the genetic diversity of S. leprosula populations. We also found an individual S. leprosula tree with a unique reproductive system, probably apomictic, mating system.

Biased distribution of microsatellite motifs in the rice genome

Microsatellites are useful tools to study the extent of divergence between two taxonomic groups that show high sequence similarity. We have compared microsatellite distribution to illustrate genetic variation between the two rice genomes, Oryza sativa L. ssp. indica and Oryza sativa L. ssp. japonica. Microsatellite distribution proved to be non random as certain regions of very high microsatellite density have been identified. Microsatellite density in the subspecies japonica was computed marginally higher than in the subspecies indica in the genomic regions compared between the two subspecies. Unexpectedly high microsatellite densities were observed in 5'-untranslated regions of genes. These regions also displayed a clear motif bias. Some of the longest microsatellite repeats were found in intron sequences. Frequency, as well as motif bias was also noted with respect to the association of microsatellites with transposable elements. Microsatellite mutability values were exemplarily estimated for 90 loci by aligning the microsatellite containing regions between the two genomes. Poor rates of finding an orthologue corresponded with high microsatellite mutability in rice. These insights are likely to play a significant role in selecting microsatellite loci to be used in molecular breeding and studying evolutionary dynamics of the two subspecies.

A predicted microsatellite map of the passerine genome based on chicken-passerine sequence similarity

Abstract We present a predicted passerine genome map consisting of 196 microsatellite markers distributed across 25 chromosomes. The map was constructed by assigning chromosomal locations based on the sequence similarity between 550 publicly available passerine microsatellites and the draft chicken genome sequence published by the International Chicken Genome Sequencing Consortium. We compared this passerine microsatellite map with a recently published great reed warbler (Acrocephalus arundinaceus) linkage map derived from the segregation of marker alleles in a pedigree of a natural population. Twenty-four microsatellite markers were shared between the two maps, distributed across ten chromosomes. Synteny was maintained between the predicted passerine microsatellite map and the great reed warbler linkage map, confirming the validity and accuracy of our approach. Possible applications of the predicted passerine microsatellite map include genome mapping; quantitative trait locus (QTL) discovery; understanding heterozygosity-fitness correlations; investigating avian karyotype evolution; understanding microsatellite mutation processes; and for identifying loci conserved in multiple species, unlinked loci for use in genotyping sets and sex-linked markers.

Microsatellite length differences between humans and chimpanzees at autosomal Loci are not found at equivalent haploid Y chromosomal Loci

When homologous microsatellites are compared between species, significant differences in mean length are often noted. A dominant cause of these length differences is ascertainment bias due to selection for maximum repeat number and repeat purity when the markers are being developed. However, even after ascertainment bias has been allowed for through reciprocal comparisons, significant length differences remain, suggesting that the average microsatellite mutation rate differs between species. Two classes of mechanism have been proposed: rapid evolution of enzymes involved in the generation and repair of slippage products (enzyme evolution model) and heterozygote instability, whereby interchromosomal events at heterozygous sites offer extra opportunities for mutations to occur (heterozygote instability model). To examine which of these hypotheses is most likely, we compared ascertainment bias and species length differences between humans and chimpanzees in autosomal and Y chromosomal microsatellites. We find that levels of ascertainment bias are indistinguishable, but that interspecies length differences are significantly greater for autosomal loci compared with haploid Y chromosomal loci. Such a pattern is consistent with predictions from the heterozygote instability model and is not expected under models of microsatellite evolution that do not include interchromosomal events such as the enzyme evolution model.

Microsatellites for ecologists: a practical guide to using and evaluating microsatellite markers

Recent improvements in genetic analysis and genotyping methods have resulted in a rapid expansion of the power of molecular markers to address ecological questions. Microsatellites have emerged as the most popular and versatile marker type for ecological applications. The rise of commercial services that can isolate microsatellites for new study species and genotype samples at reasonable prices presents ecologists with the unprecedented ability to employ genetic approaches without heavy investment in specialized equipment. Nevertheless, the lack of accessible, synthesized information on the practicalities and pitfalls of using genetic tools impedes ecologists' ability to make informed decisions on using molecular approaches and creates the risk that some will use microsatellites without understanding the steps needed to evaluate the quality of a genetic data set. The first goal of this synthesis is to provide an overview of the strengths and limitations of microsatellite markers and the risks, cost and time requirements of isolating and using microsatellites with the aid of commercial services. The second goal is to encourage the use and consistent reporting of thorough marker screening to ensure high quality data. To that end, we present a multistep screening process to evaluate candidate loci for inclusion in a genetic study that is broadly targeted to both novice and experienced geneticists alike.

Quantifying ascertainment bias and species-specific length differences in human and chimpanzee microsatellites using genome sequences

Surveys of variability of homologous microsatellite loci among species reveal an ascertainment bias for microsatellite length where microsatellite loci isolated in one species tend to be longer than homologous loci in related species. Here, we take advantage of the availability of aligned human and chimpanzee genome sequences to compare length difference of homologous microsatellites for loci identified in humans to length difference for loci identified in chimpanzees. We are able to quantify ascertainment bias for a range of motifs and microsatellite lengths. Because ascertainment bias should not exist if a microsatellite selected in one species is as likely to be longer as it is to be shorter than its homologue, we propose that the nature of ascertainment bias can provide evidence for understanding how microsatellites evolve. We show that bias is greater for longer microsatellites but also that many long microsatellites have short homologues. These results are consistent with the notion that growth of long microsatellites is constrained by an upper length boundary that, when reached, sometimes results in large deletions. By evaluating ascertainment bias separately for interrupted and uninterrupted repeats we also show that long microsatellites tend to become interrupted, thereby contributing a second component of ascertainment bias. Having accounted for ascertainment bias, in agreement with results published elsewhere, we find that microsatellites in humans are longer on average than those in chimpanzees. This length difference is similar among repeat motifs but surprisingly comprises two roughly equal components, one associated with the repeats themselves and one with the flanking sequences. The differences we find can only be explained if microsatellites are both evolving directionally under a biased mutation process and are doing so at different rates in different closely related species.

Quantification of cattle DNA using quantitative competitive PCR with sheep DNA as competitor

A novel method was developed to enable accurate and high-throughput measurement of cattle DNA concentration using quantitative competitive PCR, with sheep DNA as competitor. While quantitative competitive PCR has been used extensively for the quantification of specific RNA or DNA molecules, they have required development of internal standards with matching primer binding sites and similar amplification efficiencies to the target molecule. To develop such as assay can constitute a significant work-up. Instead, by utilizing the tendency of microsatellites developed in one species to amplify homologous loci across closely related species removes the need for internal standard development. Two cattle microsatellite markers were identified that produced distinct sheep specific peaks in an electropherogram. A standard graph was plotted for various dilutions of a cattle standard and a constant amount of sheep competitor. The sheep DNA, which is co-amplified with the cattle template in the PCR reaction served as the internal standard. The cattle DNA concentration of an unknown sample was determined by relating the ratio of sheep to cattle PCR product peaks to the standard curve. The standard deviation between replicate measurements of cattle DNA was 0.52 ng/microl using this method.