Correction method for null alleles in species with variable microsatellite flanking regions, a case study of the black-lipped pearl oyster Pinctada margaritifera

In bivalves, heterozygote deficiencies and departures from Hardy-Weinberg equilibrium (HWE) in microsatellite analysis are common and mainly attributed to inbreeding, genetic patchiness (Walhund effect), or null alleles. We checked for the occurrence of null alleles at 3 microsatellite loci in 3 populations of black-lipped pearl oyster, Pinctada margaritifera, using a step-by-step method to re-amplify homozygotes and null individuals with redesigned primer pair combinations. After amplification with original primer pairs, the 3 populations exhibited null alleles, absence of structure, and significant departure from HWE for all 3 loci due to heterozygote deficiencies. After 3 re-amplification steps, with modified primer sets, all loci were corrected for null alleles. Once corrected, all populations appeared at HWE, demonstrating that null alleles were responsible for the initial disequilibrium of the populations. Furthermore, analysis from corrected genotypes demonstrates significant genetic differentiation for one population from the other 2.     

Isolation and characterization of 19 novel microsatellite loci in Rhododendron aureum and Rhododendron brachycarpum (Ericaceae)

Nineteen polymorphic microsatellite loci are described for Rhododendron aureum Georgi (Ericaceae), an endangered species in Korea, and the closely related Rhododendron brachycarpum. The sequences containing repeat motifs were identified using low-depth next generation sequencing and 148 microsatellite loci were examined for amplification success and the detection of polymorphisms. All 19 loci were polymorphic and the number of alleles for these markers varied from 4 to 25, with an average of 15 alleles per locus. Three R. aureum loci showed departure from Hardy–Weinberg equilibrium (HWE) and one R. brachycarpum locus deviated significantly from HWE. These microsatellite loci will be useful for investigating the genetic diversity and population structure of both species.     

Isolation and characterization of ten polymorphic microsatellite loci in Ixodes arboricola, and cross-amplification in three other Ixodes species

We characterized ten polymorphic microsatellite loci from the tree-hole tick, Ixodes arboricola. Loci were screened in 11–18 individuals from three Belgian populations and five to ten alleles were found at each locus. Seven loci did not show deviations from Hardy–Weinberg equilibrium conditions and there were no indications for null alleles at these loci. The three other loci showed significant heterozygote deficiencies in at least one population, and a high potential for the occurrence of null alleles. We observed no effect of potential host DNA on the scoring of the microsatellites. Cross-amplification of the microsatellites was tested in eight specimens of three congeneric species: I. ricinus, I. hexagonus and I. frontalis. Depending on the species, six or seven of the loci were amplified in ≥4 of the 8 specimens and were polymorphic in each of these species (except for Ixaf 11 in I. frontalis and I. ricinus). These loci thus provide a tool for population genetic analysis of I. arboricola. The suitability of these markers needs to be further investigated in its congeners.     

Sixteen microsatellite loci for the Bufo boreas group

The four species that comprise the Bufo boreas group of toads are critically imperiled in all or portions of their geographical ranges. We present data from 16 microsatellite loci isolated from B. boreas that cross‐amplify in these four species. These markers have proven useful in the analyses of population structure and conservation genetics, and provide a powerful tool for future researchers who seek to understand the conservation genetics of these rare toad species. Polymorphism was assessed for 339 individuals from seven populations representing the four species. All loci were polymorphic (X? = 8, range of four to 19 alleles). Three loci were not in Hardy–Weinberg equilibrium (HWE) in one population, and one of these loci was out of HWE in a second population (P < 0.003 after Bonferroni correction for multiple tests). However, there were no systematic deviations from HWE across all study populations. Small populations in fragmented habitat appear to explain the existing HWE deviations.     

Microsatellite variability and heterozygote excess in Elymus trachycaulus populations from British Columbia in Canada

Microsatellite markers were used to assess the genetic diversity and population structure in four populations of Elymus trachycaulus from British Columbia and one population of Elymus alaskanus from Northwest Territories. Fourteen microsatellite loci were used in this study. Our results indicated that E. trachycaulus is highly polymorphic, with an average percentage of polymorphic loci of 96.5% over the four populations. Average expected heterozygosity values (H_E or gene diversity) varied from 0.418 to 0.585 with a mean of 0.497. Most of the genetic variation was found within populations (85%) and the differentiation among populations was found to be 15% (F_st = 0.15). Interpopulation genetic distances corresponded well with the geographic distance between the population sites of origin, as well as morphological characteristics. Tests for Hardy–Weinberg equilibrium (HWE) for all loci and all populations revealed that all loci significantly differ from HWE. Subsequent analysis indicated that departure from HWE at some loci was due to an excess of heterozygotes. Possible explanations for heterozygote excess are discussed. The most likely reason for observed heterozygote excess could be due to the polyploidy nature of the species.     

Isolation and characterization of 12 polymorphic microsatellite markers for the frog Pelophylax hubeiensis

Twelve polymorphic microsatellite loci were isolated from an (AC)n- and (AG)n-enriched DNA library for the endemic Chinese frog Pelophylax hubeiensis (Ranidae). The number of alleles per locus ranged from two to eight, with a mean of 5.17. The observed and expected heterozygosities ranged from 0.226 to 0.839 and from 0.204 to 0.826, with means of 0.568 and 0.656, respectively. No significant linkage disequilibrium was detected among these loci. However, two significant deviations from HWE were discovered at loci Pehu-11 and Pehu-12 (P<0.05). MICRO-CHECKER tests showed that null alleles could be present at locus Pehu-12. These polymorphic microsatellite loci can be employed for exploring mating mechanisms, population genetic structure and other relevant genetic investigations of P. hubeiensis.     

Effects of salinity on larval stages of the rhizocephalan barnacle Loxothylacus texanus: survival and metamorphosis in response to the host, Callinectes sapidus

It is known that the rhizocephalan barnacle Loxothylacus texanus infects the greater blue crab, Callinectes sapidus, in the Gulf of Mexico and adjacent waters, however, factors that affect the prevalence and distribution of this parasite, particularly the dispersive larval stages of this organism, are not well understood. In the current study, the effects of salinity on larval survival and the metamorphosis of L. texanus in response to postmolt host exoskeleton were examined. Acute and acclimated responses were similar. Larval survival was highest in the 20-35‰ range, with 100% mortality of nauplii at all salinities <20‰ and >50‰. L. texanus cyprids were able to metamorphose over a broad range of salinities (15-60‰). In several cases, metamorphosis was actually greatest at high salinities (40-50‰). These data predict that L. texanus larvae would be concentrated in portions of Gulf of Mexico waters with salinities >20‰ such as the mouths of estuaries and bays. Conversely, upper regions of estuaries may be inhospitable to the dispersive (naupliar) stage of the parasite and may serve as a refuge from infection for host crabs.     

Heterozygote deficits in cyst plant‐parasitic nematodes: possible causes and consequences

Deviations of genotypic frequencies from Hardy–Weinberg equilibrium (HWE) expectations could reveal important aspects of the biology of populations. Deviations from HWE due to heterozygote deficits have been recorded for three plant‐parasitic nematode species. However, it has never been determined whether the observed deficits were due (i) to the presence of null alleles, (ii) to a high level of consanguinity and/or (iii) to a Wahlund effect. The aim of the present work was, while taking into the possible confounding effect of null alleles, to disentangle consanguinity and Wahlund effect in natural populations of those three economically important cyst nematodes using microsatellite markers: Globodera pallida, G. tabacum and Heterodera schachtii, pests of potato, tobacco and sugar beet, respectively. The results show a consistent pattern of heterozygote deficiency in the three nematode species sampled at the spatial scale of the host plant. We demonstrate that the prevalence of null alleles is weak and that heterozygote deficits do not have a single origin. Our results suggested that it is restricted dispersal that leads to heterozygote deficits through both consanguinity and substructure, which effects can be linked to soil movement, cyst density, and the number of generations per year. We discuss potential implications for the durability of plant resistances that are used to protect crops against parasites in which mating between relatives occur. While consanguineous mating leads to homozygosity at all loci, including loci governing avirulence/virulence, which favours the expression of virulence when recessive, the Wahlund effect is expected to have no particular effect on the adaptation of nematodes to resistances.     

Population structure and cryptic genetic variation in the mango fruit fly,Ceratitis cosyra (Diptera,Tephritidae)

The fruit fly Ceratitis cosyra is an important agricultural pest negatively affecting the mango crop production throughout Africa and also feeding on a variety of other wild and cultivated hosts. The occurrence of deeply divergent haplotypes, as well as extensive morphological variability, previously suggested possible cryptic speciation within Ceratitis cosyra. Here we provide the first large-scale characterisation of the population structure of Ceratitis cosyra with the main objective of verifying cryptic genetic variation. A total of 348 specimens from 13 populations were genotyped at 16 polymorphic microsatellite loci. Hardy-Weinberg equilibrium (HWE) deviations were observed in 40.4% of locus-population combinations and suggested the occurrence of genetic substructuring within populations. Discriminant Analysis of Principal Components (DAPC) showed genetic divergence between the vast majority of vouchers from Burundi and Tanzania (plus a few outliers from other African countries) and all other specimens sampled. Individual Bayesian assignments confirmed the existence of two main genotypic groups also occurring in sympatry. These data provided further support to the hypothesis that Ceratitis cosyra might include cryptic species. However, additional integrative taxonomy, possibly combining morphological, ecological and physiological approaches, is required to provide the necessary experimental support to this model.     

Development of eighteen tetranucleotide microsatellite markers in Tibetan Macaque (Macaca thibetana) and genetic diversity analysis of captive population

In this study, eighteen tetranucleotide microsatellite loci were isolated from AAAG-enriched and GATA-enriched libraries of the Tibetan macaque, Macaca thibetana. These loci were tested on 53 individuals of M. thibetana, and most loci were proved to be highly polymorphic. A total of 109 alleles were detected with an average of 6.06 alleles per locus. The PIC values of these loci ranged from 0.192 to 0.879, with an average of 0.624. The observed and expected heterozygosities ranged from 0.170 to 0.800 and from 0.217 to 0.898, with an average of 0.583 and 0.675, respectively. 5 loci significantly deviated from Hardy–Weinberg equilibrium (HWE). Significant linkage disequilibrium (LD) was found between 9 pairs of loci. The newly identified polymorphic markers would facilitate the study of M. thibetana on the population structure and genetic diversity.     

Pleiotropy of FRIGIDA enhances the potential for multivariate adaptation

An evolutionary response to selection requires genetic variation; however, even if it exists, then the genetic details of the variation can constrain adaptation. In the simplest case, unlinked loci and uncorrelated phenotypes respond directly to multivariate selection and permit unrestricted paths to adaptive peaks. By contrast, ‘antagonistic’ pleiotropic loci may constrain adaptation by affecting variation of many traits and limiting the direction of trait correlations to vectors that are not favoured by selection. However, certain pleiotropic configurations may improve the conditions for adaptive evolution. Here, we present evidence that the Arabidopsis thaliana gene FRI (FRIGIDA) exhibits ‘adaptive’ pleiotropy, producing trait correlations along an axis that results in two adaptive strategies. Derived, low expression FRI alleles confer a ‘drought escape’ strategy owing to fast growth, low water use efficiency and early flowering. By contrast, a dehydration avoidance strategy is conferred by the ancestral phenotype of late flowering, slow growth and efficient water use during photosynthesis. The dehydration avoidant phenotype was recovered when genotypes with null FRI alleles were transformed with functional alleles. Our findings indicate that the well-documented effects of FRI on phenology result from differences in physiology, not only a simple developmental switch.     

Genetic diversity analysis by microsatellite markers in four captive populations of the sika deer (Cervus nippon)

The Chinese sika deer (Cervus nippon) is a rare and vulnerable animal in China for medical use. In this study, the genetic diversity and population genetic structure of 113 Chinese sika deer from 4 populations (Linyi Farm, LF; Linyi Park, LP, Yangzhou Farm, YF; Yangzhou Zoo, YZ) were investigated with 14 microsatellite loci. Eighty-three alleles were detected at the 14 loci in all populations. The expected heterozygosity ranged from 0.257 to 0.863 and the observed heterozygosity from 0.226 to 0.821. The polymorphism information content at different loci ranged from 0.217 to 0.825. The results of the HWE (Hardy–Weinberg equilibrium) tests indicated that only four loci (CEH-5, BL42, Mber70, and CEH-2) were in HWE (P > 0.01). The mean number of alleles per population ranged from 3.21 to 5.64, observed and expected heterozygosities ranged from 0.568 to 0.685, respectively. Positive inbreeding coefficient (F_IS) values were found in every population. F_ST values ranged from 0.101 in the LF to 0.155 in the YZ. The genetic identity ranged from 0.1236 to 0.1645. The genetic distance ranged from 0.4746 to 0.6025. The results of this study indicate moderate genetic variation and polymorphism across the loci. Appropriate breeding strategies should be designed for deer in captivity.     

Characterisation of 12 microsatellite loci in the Vietnamese commercial clam <Emphasis Type="Italic">Lutraria rhynchaena</Emphasis> Jonas 1844 (Heterodonta: Bivalvia: Mactridae) through next-generation sequencing

The marine clam Lutraria rhynchaena is gaining popularity as an aquaculture species in Asia. Lutraria populations are present in the wild throughout Vietnam and several stocks have been established and translocated for breeding and aquaculture grow-out purposes. In this study, we demonstrate the feasibility of utilising Illumina next-generation sequencing technology to streamline the identification and genotyping of microsatellite loci from this clam species. Based on an initial partial genome scan, 48 microsatellite markers with similar melting temperatures were identified and characterised. The 12 most suitable polymorphic loci were then genotyped using 51 individuals from a population in Quang Ninh Province, North Vietnam. Genetic variation was low (mean number of alleles per locus = 2.6; mean expected heterozygosity = 0.41). Two loci showed significant deviation from Hardy–Weinberg equilibrium (HWE) and the presence of null alleles, but there was no evidence of linkage disequilibrium among loci. Three additional populations were screened (n = 7–36) to test the geographic utility of the 12 loci, which revealed 100 % successful genotyping in two populations from central Vietnam (Nha Trang). However, a second population from north Vietnam (Co To) could not be successfully genotyped and morphological evidence and mitochondrial variation suggests that this population represents a cryptic species of Lutraria. Comparisons of the Qang Ninh and Nha Trang populations, excluding the 2 loci out of HWE, revealed statistically significant allelic variation at 4 loci. We reported the first microsatellite loci set for the marine clam Lutraria rhynchaena and demonstrated its potential in differentiating clam populations. Additionally, a cryptic species population of Lutraria rhynchaena was identified during initial loci development, underscoring the overlooked diversity of marine clam species in Vietnam and the need to genetically characterise population representatives prior to microsatellite development. The rapid identification and validation of microsatellite loci using next-generation sequencing technology warrant its integration into future microsatellite loci development for key aquaculture species in Vietnam and more generally, aquaculture countries in the South East Asia region.     

Fifteen new polymorphic microsatellite loci for the meadow brown butterfly,Maniola jurtina

We characterized fifteen microsatellite markers for the butterfly Maniola jurtina. For the six studied populations (96 samples) the total number of alleles per locus ranged from 3 to 55 and mean overall expected heterozygosity across all loci was 0.74. In spite of a high frequency of null alleles detected in part of the loci, a recurrent phenomenon in Lepidopteron, the estimation of pairwise F_ST seems rather insensitive to the presence of these null alleles as shown by the high correlation between F_ST calculated after correction for the presence of null alleles and non-corrected F_ST, indicating that the loci may be usable in population genetics, more specifically for the study of populations genetics structure.     

Effect of balancing selection on spatial genetic structure within populations: theoretical investigations on the self-incompatibility locus and empirical studies in Arabidopsis halleri

The effect of selection on patterns of genetic structure within and between populations may be studied by contrasting observed patterns at the genes targeted by selection with those of unlinked neutral marker loci. Local directional selection on target genes will produce stronger population genetic structure than at neutral loci, whereas the reverse is expected for balancing selection. However, theoretical predictions on the intensity of this signal under precise models of balancing selection are still lacking. Using negative frequency-dependent selection acting on self-incompatibility systems in plants as a model of balancing selection, we investigated the effect of such selection on patterns of spatial genetic structure within a continuous population. Using numerical simulations, we tested the effect of the type of self-incompatibility system, the number of alleles at the self-incompatibility locus and the dominance interactions among them, the extent of gene dispersal, and the immigration rate on spatial genetic structure at the selected locus and at unlinked neutral loci. We confirm that frequency-dependent selection is expected to reduce the extent of spatial genetic structure as compared to neutral loci, particularly in situations with low number of alleles at the self-incompatibility locus, high frequency of codominant interactions among alleles, restricted gene dispersal and restricted immigration from outside populations. Hence the signature of selection on spatial genetic structure is expected to vary across species and populations, and we show that empirical data from the literature as well as data reported here on three natural populations of the herb Arabidopsis halleri confirm these theoretical results.     

Genetic Polymorphism of Microsatellite DNA in Two Populations of Northern Sheatfish (Silurus soldatovi)

In this article, population variations and genetic structures of two populations of northern sheatfish (Silurus soldatovi) were analyzed using 24 microsatellite loci enriched from southern catfish (S. meriaionalis Chen) by magnetic beads. Gene frequency (P), observed heterozygosity (Ho), expected heterozygosity (He), polymorphism information contents (PIC), and number of effective alleles (Ne) were determined. One population was wild, ripe individuals collected from Heilongjiang River (HNS); the other was cultured fry collected from Songhuajiang River (SNS). The Hardy-Weinberg equilibrium (HWE) was tested by the genetic departure index (d). The coefficient of gene differentiation G_ST and Φ_ST by AMOVA (Analysis of Molecular Variety) was imputed using Arlequin software in this study. In addition, a phylogenetic tree was constructed by UPGMA method based on the pairwise Nei's standard distances using PHYLIP. A total of 1 357 fragments with sizes ranging between 102 bp and 385 bp were acquired by PCR amplifications. The average number of alleles of the two populations was 8.875. Results indicated that these microsatellite loci were highly polymorphic and could be used as genetic markers. The mean values of the parameters P, Ho, He, PIC, and Ne were 0.165, 0.435, 0.758, 0.742, and 5.019 for HNS and 0.147, 0.299, 0.847, 0.764, and 5.944 for SNS, respectively. Although there were differences, there were no significant differentiations except for the locus HLJcf37. These populations to a certain extent deviated from HWE, such as excessive and deficient heterozygote numbers. The value of G_ST was 0.078 and above 98% of the variation were differences among individuals within the population, so the variation between populations was insignificant. Cluster analysis also showed that the relationships among individuals were very close. In conclusion, the microsatellite markers that were developed through this study are useful for genetic analysis and the genetic culture that was proposed in this study has no significant impact on S. soldatovi.     

Characterization of 20 microsatellite loci by multiplex PCR in swimming crab, Portunus trituberculatus

The swimming crab, Portunus trituberculatus (Crustacea: Decapoda: Brachyura), is one of the most important fishery resources in East Asia. The aim of this work is to develop microsatellite multiplex PCR systems for 20 microsatellites which were isolated from a microsatellite-enriched genomic library and to explore the genetic diversity of P. trituberculatus captured in the coastal regions of the Yellow Sea in Korea. We prepared four multiplex PCR systems (two hexplex and tetraplex systems each) which could be used for the genotyping of 20 markers. Most markers in this system are highly informative with mean polymorphism information contents and observed heterozygosity of 0.913 and 0.930. The total number of alleles observed in this study was 724, and the mean allelic number per locus was 36.2. Significant deviations from the Hardy–Weinberg equilibrium (HWE) were observed at six loci in the Wonsando Island population, but all the loci deviated from the HWE in the Jeonjangpo population. The genetic differentiation between the two populations and the phylogenetic tree constructed by the unweighted pair group method with the arithmetic average suggested that there is no significant genetic difference between Wonsando Island and Jeonjangpo populations. The prepared multiplex PCR systems and population genetic data will be helpful to study phylogeographic analysis as well as to prepare strategies of stock management and aquaculture of P. trituberculatus.     

Natural variation in the genes responsible for maturity loci E1, E2, E3 and E4 in soybean

Background and Aims

The timing of flowering has a direct impact on successful seed production in plants. Flowering of soybean (Glycine max) is controlled by several E loci, and previous studies identified the genes responsible for the flowering loci E1, E2, E3 and E4. However, natural variation in these genes has not been fully elucidated. The aims of this study were the identification of new alleles, establishment of allele diagnoses, examination of allelic combinations for adaptability, and analysis of the integrated effect of these loci on flowering.

Methods

The sequences of these genes and their flanking regions were determined for 39 accessions by primer walking. Systematic discrimination among alleles was performed using DNA markers. Genotypes at the E1–E4 loci were determined for 63 accessions covering several ecological types using DNA markers and sequencing, and flowering times of these accessions at three sowing times were recorded.

Key Results

A new allele with an insertion of a long interspersed nuclear element (LINE) at the promoter of the E1 locus (e1-re) was identified. Insertion and deletion of 36 bases in the eighth intron (E2-in and E2-dl) were observed at the E2 locus. Systematic discrimination among the alleles at the E1–E3 loci was achieved using PCR-based markers. Allelic combinations at the E1–E4 loci were found to be associated with ecological types, and about 62–66 % of variation of flowering time could be attributed to these loci.

Conclusions

The study advances understanding of the combined roles of the E1–E4 loci in flowering and geographic adaptation, and suggests the existence of unidentified genes for flowering in soybean,     

Genetic basis of natural variation in body pigmentation in Drosophila melanogaster

Body pigmentation in insects and other organisms is typically variable within and between species and is often associated with fitness. Regulatory variants with large effects at bab1, t and e affect variation in abdominal pigmentation in several populations of Drosophila melanogaster. Recently, we performed a genome wide association (GWA) analysis of variation in abdominal pigmentation using the inbred, sequenced lines of the Drosophila Genetic Reference Panel (DGRP). We confirmed the large effects of regulatory variants in bab1, t and e; identified 81 additional candidate genes; and validated 17 candidate genes (out of 28 tested) using RNAi knockdown of gene expression and mutant alleles. However, these analyses are imperfect proxies for the effects of segregating variants. Here, we describe the results of an extreme quantitative trait locus (xQTL) GWA analysis of female body pigmentation in an outbred population derived from light and dark DGRP lines. We replicated the effects on pigmentation of 28 genes implicated by the DGRP GWA study, including bab1, t and e and 7 genes previously validated by RNAi and/or mutant analyses. We also identified many additional loci. The genetic architecture of Drosophila pigmentation is complex, with a few major genes and many other loci with smaller effects.     

Characterization of 11 novel polymorphic microsatellite loci in the threatened Korean loach,Iksookimia koreensis,isolated using a next-generation sequencing method

To investigate the genetic diversity and population structure of Iksookimia koreensis, we characterized 11 novel polymorphic microsatellite markers developed using next-generation sequencing. The number of alleles per locus ranged from 4 to 10 (average = 6.26). Observed and expected heterozygosities ranged from 0.333 to 0.866 and from 0.375 to 0.866, respectively. No loci showed significant deviation from Hardy–Weinberg equilibrium (HWE). These loci were also used successfully to study the genetic diversity of a closely related species, Iksookimia longicorpa. Four of the 11 loci amplified in the two species showed different allele frequency and distribution, indicating deep genetic divergence between I. koreensis and I. longicorpa. The newly developed microsatellite markers reported here will provide a useful tool for examining gene flow, population genetic structure, and genetic diversity of these species.