Correcting Coalescent Analyses for Panel-Based SNP Ascertainment

Single-nucleotide polymorphism (SNP) data are routinely obtained by sequencing a region of interest in a small panel, constructing a chip with probes specific to sites found to vary in the panel, and using the chip to assay subsequent samples. The size of the chip is often reduced by removing low-frequency alleles from the set of SNPs. Using coalescent estimation of the scaled population size parameter, Θ, as a test case, we demonstrate the loss of information inherent in this procedure and develop corrections for coalescent analysis of SNPs obtained via a panel. We show that more accurate Θ-estimates can be recovered if the panel size is known, but at considerable computational cost as the panel individuals must be explicitly modeled in the analysis. We extend this technique to apply to the case where rare alleles have been omitted from the SNP panel. We find that when appropriate corrections for panel ascertainment and rare-allele omission are used, the biases introduced by ascertainment are largely correctable, but recovered estimates are less accurate than would be obtained with fully sequenced data. This method is then applied to recombinant multiple population data to investigate the effects of recombination and migration on the estimate of Θ.     

Inference of population structure and patterns of gene flow in canine heartworm (Dirofilaria immitis)

Understanding the genetic variation within a parasitic species is crucial to implementing successful control programs and preventing the dispersal of drug resistance alleles. We examined the population genetics and structure of canine heartworm (Dirofilaria immitis) by developing a panel of 11 polymorphic microsatellite loci for this abundant parasite. In total, 192 individual nematodes were opportunistically sampled from 9 geographic regions in the United States and Mexico and genotyped. Population genetic analyses indicate the presence of 4 genetic clusters. The canine heartworm samples used in this study were characterized by low heterozygosity, with eastern and central North America experiencing high levels of reciprocal gene flow. Geographic barriers impede the movement of vectors and infected hosts west of the Rocky Mountains and south of the Central Mexican Plateau. This, combined with corridors of contiguous habitat, could influence the spread of drug resistance alleles.     

Imputation of genotypes from low- to high-density genotyping platforms and implications for genomic selection

The objective of this study was to quantify the accuracy achievable from imputing genotypes from a commercially available low-density marker panel (2730 single nucleotide polymorphisms (SNPs) following edits) to a commercially available higher density marker panel (51 602 SNPs following edits) in Holstein-Friesian cattle using Beagle, a freely available software package. A population of 764 Holstein-Friesian animals born since 2006 were used as the test group to quantify the accuracy of imputation, all of which had genotypes for the high-density panel; only SNPs on the low-density panel were retained with the remaining SNPs to be imputed. The reference population for imputation consisted of 4732 animals born before 2006 also with genotypes on the higher density marker panel. The concordance between the actual and imputed genotypes in the test group of animals did not vary across chromosomes and was on average 95%; the concordance between actual and imputed alleles was, on average, 97% across all SNPs. Genomic predictions were undertaken across a range of production and functional traits for the 764 test group animals using either their real or imputed genotypes. Little or no mean difference in the genomic predictions was evident when comparing direct genomic values (DGVs) using real or imputed genotypes. The average correlation between the DGVs estimated using the real or imputed genotypes for the 15 traits included in the Irish total merit index was 0.97 (range of 0.92 to 0.99), indicating good concordance between proofs from real or imputed genotypes. Results show that a commercially available high-density marker panel can be imputed from a commercially available lower density marker panel, which will also have a lower cost, thereby facilitating a reduction in the cost of genomic selection. Increased available numbers of genotyped and phenotyped animals also has implications for increasing the accuracy of genomic prediction in the entire population and thus genetic gain using genomic selection.     

Multiplexed systems of microsatellite markers for genetic analysis of mahogany,Swietenia macrophylla King (Meliaceae), a threatened neotropical timber species

Mahogany (Swietenia macrophylla King [Meliaceae]) is the most valuable hardwood species in the neotropics. Its conservation status has been the subject of increasing concern due to overexploitation and habitat destruction. In this work we report the development and characterization of 10 highly variable microsatellite loci for S. macrophylla. Twenty-nine percent of the 126 sequenced mahogany clones yielded useful microsatellite loci. Three high-throughput genotyping systems were developed based on polymerase chain reaction (PCR) multiplexing of these mahogany loci. We identified a total of 158 alleles in 121 adult individuals of S. macrophylla, with an average of 15.8 alleles (range 11-25) per locus. All loci showed Mendelian inheritance in open-pollinated half-sib families. The mean expected heterozygosity was 0.84 and the mean observed heterozygosity was 0.73. The combined probability of identity-the probability that two individuals selected at random from a population would have identical genotypes--was 7.0 x 10(-15), and combined probability of paternity exclusion was 0.999998 overall loci. These microsatellite loci permit precise estimates of parameters such as gene flow, mating system, and paternity, thus providing important insights into the population genetics and conservation of S. macrophylla.     

Development of novel tetra- and trinucleotide microsatellite markers for giant grouper <Emphasis Type="Italic">Epinephelus lanceolatus</Emphasis> using 454 pyrosequencing

Giant grouper (Epinephelus lanceolatus) is a commercially important species, but its wild population has recently been classified as vulnerable. This species has significant potential for use in aquaculture, though a greater understanding of population genetics is necessary for selective breeding programs to minimize kinship for genetically healthy individuals. High-throughput pyrosequencing of genomic DNA was used to identify and characterize novel tetra- and trinucleotide microsatellite markers in giant grouper from Sabah, Malaysia. In total, of 62,763 sequences containing simple sequence repeats (SSRs) were obtained, and 78 SSR loci were selected to possibly contain tetra- and trinucleotide repeats. Of these loci, 16 had tetra- and 8 had trinucleotide repeats, all of which exhibited polymorphisms within easily genotyped regions. A total of 143 alleles were identified with an average of 5.94 alleles per locus, with mean observed and expected heterozygosities of 0.648 and 0.620, respectively. Among of them, 15 microsatellite markers were identified without null alleles and with Hardy–Weinberg equilibrium. These alleles showed a combined non-exclusion probability of 0.01138. The probability of individual identification (P_ID) value combined with in descending order 12 microsatellite markers was 0.00008, which strongly suggests that the use of the microsatellite markers developed in this study in various combinations would result in a high resolution method for parentage analysis and individual identification. These markers could be used to establish a broodstock management program for giant grouper and to provide a foundation for genetic studies such as population structure, parentage analysis, and kinship selection.     

Isolation and characterization of 13 microsatellite loci for Neochamaelea pulverulenta (Cneoraceae)

We report 13 polymorphic nuclear microsatellite loci for Neochamaelea pulverulenta (Cneoraceae) using an enriched-library approach. Although this plant species is tetraploid, expected patterns for tetrasomic segregation were completely absent, and all loci analysed showed a diploid pattern of inheritance. We detected a total of 102 alleles in 57 individuals genotyped (mean number of alleles per locus was 7.85). The values of observed and expected heterozygosities ranged from 0.193 to 0.737 and 0.425 to 0.812 respectively. Disomic segregation, levels of polymorphism and the exclusionary power of the developed markers render them readily applicable for parentage assignment of dispersed seeds, and for analyses of spatial genetic structure and population connectivity.     

A large panel of novel microsatellite markers for the bank vole (Myodes glareolus)

We describe a set of 66 highly polymorphic microsatellite loci isolated from the bank vole, Myodes (Clethrionomys) glareolus. These microsatellites were characterized for a long-term study on periodically fluctuating density of the bank vole population in Central Finland. We detected six to 38 alleles per locus in the population sampled at two different density phases, and the levels of observed and expected heterozygosities varied between 0.17 and 1.00, and between 0.72 and 0.95, respectively. This microsatellite panel serves as an informative tool for population and molecular genetic studies.     

BARCSoySNP23: a panel of 23 selected SNPs for soybean cultivar identification

This report describes a set of 23 informative SNPs (BARCSoySNP23) distributed on 19 of the 20 soybean linkage groups that can be used for soybean cultivar identification. Selection of the SNPs to include in this set was made based upon the information provided by each SNP for distinguishing a diverse set of soybean genotypes as well as the linkage map position of each SNP. The genotypes included the ancestors of North American cultivars, modern North American cultivars and a group of Korean cultivars. The procedure used to identify this subset of highly informative SNP markers resulted in a significant increase in the power of identification versus any other randomly selected set of equal number. This conclusion was supported by a simulation which indicated that the 23-SNP panel can uniquely distinguish 2,200 soybean cultivars, whereas sets of randomly selected 23-SNP panels allowed the unique identification of only about 50 cultivars. The 23-SNP panel can efficiently distinguish each of the genotypes within four maturity group sets of additional cultivars/lines that have identical classical pigmentation and morphological traits. Comparatively, the 13 trinucleotide SSR set published earlier (BARCSoySSR13) has more power on a per locus basis because of the multi-allelic nature of SSRs. However, the assay of bi-allelic SNP loci can be multi-plexed using non-gel based techniques allowing for rapid determination of the SNP alleles present in soybean genotypes, thereby compensating for their relatively low information content. Both BARCSoySNP23 and BARCSoySSR13 were highly congruent relative to identifying genotypes and for estimating population genetic differences.     

Microsatellite markers polymorphism in the breeding nutria (<Emphasis Type="Italic">Myocastor coypus</Emphasis>) population in Poland

The aim of the research was to establish a microsatellite panel to determine the genetic diversity within the breeding nutria population in Poland. In the study, 92 animals representing six color forms were used. Ten fluorescently labeled microsatellite markers were investigated by multicolored capillary electrophoresis. All the microsatellites were polymorphic. The average heterozygosity observed among the population was 41%. The mean number of alleles per locus was 9.2. The average heterozygosity observed in the whole population was lower than expected. This implies that the nutria population deviates from the Hardy–Weinberg equilibrium. Low M values (from 0.078 to 0.545) of the Garza–Williamson index reveal a reduction of genetic variation in the investigated population and suggest that the breeding nutria population is remnant.     

Characterization of 14 microsatellite loci in a tropical palm, Attalea phalerata (Arecaceae)

? Premise of the study: We developed microsatellite primers for the widely distributed tropical palm Attalea phalerata for studies on the dispersal and spatial genetic structure of palm populations. ? Methods and Results: Fourteen di-, tri-, and tetra-nucleotide microsatellite primer pairs were identified. The number of alleles in the population tested ranged between 3 and 25, with a mean of 12.1. Ten microsatellite loci exhibited no significant deviations from Hardy-Weinberg Equilibrium or presence of null alleles, and their combined probability of exclusion was 0.998. ? Conclusions: These microsatellite loci will be useful in parentage analysis and population genetics studies of Attalea phalerata.     

Comparison of genotype imputation strategies using a combined reference panel for chicken population

Using whole-genome sequence (WGS) data are supposed to be optimal for genome-wide association studies and genomic predictions. However, sequencing thousands of individuals of interest is expensive. Imputation from single nucleotide polymorphisms panels to WGS data is an attractive approach to obtain highly reliable WGS data at low cost. Here, we conducted a genotype imputation study with a combined reference panel in yellow-feather dwarf broiler population. The combined reference panel was assembled by sequencing 24 key individuals of a yellow-feather dwarf broiler population (internal reference panel) and WGS data from 311 chickens in public databases (external reference panel). Three scenarios were investigated to determine how different factors affect the accuracy of imputation from 600 K array data to WGS data, including: genotype imputation with internal, external and combined reference panels; the number of internal reference individuals in the combined reference panel; and different reference sizes and selection strategies of an external reference panel. Results showed that imputation accuracy from 600 K to WGS data were 0.834±0.012, 0.920±0.007 and 0.982±0.003 for the internal, external and combined reference panels, respectively. Increasing the reference size from 50 to 250 improved the accuracy of genotype imputation from 0.848 to 0.974 for the combined reference panel and from 0.647 to 0.917 for the external reference panel. The selection strategies for the external reference panel had no impact on the accuracy of imputation using the combined reference panel. However, if only an external reference panel with reference size >50 was used, the selection strategy of minimizing the average distance to the closest leaf had the greatest imputation accuracy compared with other methods. Generally, using a combined reference panel provided greater imputation accuracy, especially for low-frequency variants. In conclusion, the optimal imputation strategy with a combined reference panel should comprehensively consider genetic diversity of the study population, availability and properties of external reference panels, sequencing and computing costs, and frequency of imputed variants. This work sheds light on how to design and execute genotype imputation with a combined external reference panel in a livestock population.     

Polymorphic heterologous microsatellite loci for population genetics studies of the white-faced ibis Plegadis chihi (Vieillot, 1817) (Pelecaniformes, Threskiornithidae)

We screened 44 heterologous microsatellites isolated in species of the families Threskiornithidae, Ciconiidae and Ardeidae for their use in a migratory waterbird, the white-faced ibis Plegadis chihi (Vieillot, 1817) (Threskiornithidae). Of the screened loci, 57% amplified successfully and 24% were polymorphic. In two breeding colonies from southern Brazil (N = 131) we detected 32 alleles (2-10 alleles/locus). Average He over all loci and colonies was 0.55, and the combined probability of excluding false parents, 98%. There was no departure from HWE in any loci or population. Eru6 and Eru4 loci were in non-random association in the Alvorada colony, and NnNF5 and Eru5 in both populations. AMOVA analysis indicated that most of the genetic diversity was contained within populations. Structure analysis suggested a single population, and F(ST) value showed weak genetic structuring (F(ST) = 0.009, p = 0.05). The two populations are apparently connected through gene-flow. The panel of six microsatellites optimized here was sufficiently informative for characterizing the genetic diversity and structure in these natural populations of the white-faced ibis. The information generated could be useful in future studies of genetic diversity, relatedness and the mating system in Plegadis chihi and related species.     

Subtyping for TNFa microsatellite sequence variation

The microsatellite locus TNFa is frequently used as an additional genetic marker in studies of the major histocompatibility complex (MHC). Novel sequence variations at the TNFa locus have been described, and which may have implications for genetic analyses. In this study, we set up a nested polymerase chain reaction-sequence-specific primer (PCR-SSP) approach to type for these TNFa sequence variations. First, sequencing analysis of workshop B lymphoblastoid cell lines (n=13) showed the presence of three sequence variations upstream of the dinucleotide repeat at TNFa. Using nested PCR-SSP, we were able to detect these variations in a larger B lymphoblastoid cell line panel (n=34). Furthermore, we were able to show that TNFa alleles a7 and a10 are present in two distinct conformations leading to "splitting" of TNFa alleles exhibiting identical fragment lengths. To establish the frequency of the TNFa alleles and their variants, we performed microsatellite typing of a large panel of random individuals from the Dutch population (n=272). Subsequent nested PCR-SSP typing showed the presence of three previously described sequence variations in the Dutch population. Furthermore, the presence of a fourth subtype was established. The described variations of allele TNFa7 and TNFa10 are present in the random population with significant frequencies. Haplotyping analysis between HLA-DR, TNFa, and HLA-B showed that allele TNFa7.2 is present in an extended DR7-TNFa7.2-B13 haplotype. In this way, we were able to show that the additional sequence variations behave like distinct TNFa alleles.     

广西毛南族DXS7133、DXS8378、DXS6789和DXS7423基因座遗传多态性

采用PCR-STR及基因分型技术,对广西毛南族167名(女57,男110)健康无关个体4个X-STR基因座(DXS7133、DXS8378、DXS6789和DXS7423)的遗传多态性进行研究。结果显示4个X-STR基因座分别检出4、5、9、3个等位基因和5、9、18、5种基因型,4个X-STR基因座女性的基因型频率分布均符合Hardy-Weinberg平衡定律(P>0.05)。群体遗传多态性指标为:多态信息含量(PIC)0.9611、男性个体识别率(DPmale)0.9771、女性个体识别率(DPfemale)0.9980、父-母-女三联体非父排除率(MECtrio)0.9611、父-女二联体非父排除率(MECduo)0.8821,显示上述4个X-STR基因座均具有较高多态性,在法医学个人识别、亲权鉴定及群体遗传学研究中有重要应用价值,同时也为人类群体遗传学、法医学等研究提供了广西毛南族群体X-STR基因座的基础数据,丰富了中华民族基因数据库。     

Isolation and characterization of microsatellite markers for the Korean rockfish, Sebastes schlegeli

The Korean rockfish (Sebastes schlegeli) is an important commercial fish that is widely used in aquaculture. We isolated and characterized 18 polymorphic microsatellite loci from the Korean rockfish using a (GT)(13)-enriched genomic library. Polymorphism was assessed in 48 individuals from a single population collected from the northern coastal waters of the Yellow Sea. The observed and expected heterozygosities ranged from 0.0244 to 0.7660 (mean 0.4194) and 0.0244 to 0.8758 (mean 0.5002), respectively. Polymorphism at these loci indicated from two to 15 alleles (mean 5.7); 14 of 18 loci conformed to Hardy-Weinberg equilibrium. These markers should be useful for management and conservation studies of this species.     

Novel mutations and DNA-based screening in non-Jewish carriers of Tay-Sachs disease.

We have evaluated the feasibility of using PCR-based mutation screening for non-Jewish enzyme-defined carriers identified through Tay-Sachs disease-prevention programs. Although Tay-Sachs mutations are rare in the general population, non-Jewish individuals may be screened as spouses of Jewish carriers or as relatives of probands. In order to define a panel of alleles that might account for the majority of mutations in non-Jewish carriers, we investigated 26 independent alleles from 20 obligate carriers and 3 affected individuals. Eighteen alleles were represented by 12 previously identified mutations, 7 that were newly identified, and 1 that remains unidentified. We then investigated 46 enzyme-defined carrier alleles: 19 were pseudodeficiency alleles, and five mutations accounted for 15 other alleles. An eighth new mutation was detected among enzyme-defined carriers. Eleven alleles remain unidentified, despite the testing for 23 alleles. Some may represent false positives for the enzyme test. Our results indicate that predominant mutations, other than the two pseudodeficiency alleles (739C-->T and 745C-->T) and one disease allele (IVS9+1G-->A), do not occur in the general population. This suggests that it is not possible to define a collection of mutations that could identify an overwhelming majority of the alleles in non-Jews who may require Tay-Sachs carrier screening. We conclude that determination of carrier status by DNA analysis alone is inefficient because of the large proportion of rare alleles. Notwithstanding the possibility of false positives inherent to enzyme screening, this method remains an essential component of carrier screening in non-Jews. DNA screening can be best used as an adjunct to enzyme testing to exclude known HEXA pseudodeficiency alleles, the IVS9+1G-->A disease allele, and other mutations relevant to the subject's genetic heritage.     

Effects of sample size,number of markers,and allelic richness on the detection of spatial genetic pattern

The influence of study design on the ability to detect the effects of landscape pattern on gene flow is one of the most pressing methodological gaps in landscape genetic research. To investigate the effect of study design on landscape genetics inference, we used a spatially‐explicit, individual‐based program to simulate gene flow in a spatially continuous population inhabiting a landscape with gradual spatial changes in resistance to movement. We simulated a wide range of combinations of number of loci, number of alleles per locus and number of individuals sampled from the population. We assessed how these three aspects of study design influenced the statistical power to successfully identify the generating process among competing hypotheses of isolation‐by‐distance, isolation‐by‐barrier, and isolation‐by‐landscape resistance using a causal modelling approach with partial Mantel tests. We modelled the statistical power to identify the generating process as a response surface for equilibrium and non‐equilibrium conditions after introduction of isolation‐by‐landscape resistance. All three variables (loci, alleles and sampled individuals) affect the power of causal modelling, but to different degrees. Stronger partial Mantel r correlations between landscape distances and genetic distances were found when more loci were used and when loci were more variable, which makes comparisons of effect size between studies difficult. Number of individuals did not affect the accuracy through mean equilibrium partial Mantel r, but larger samples decreased the uncertainty (increasing the precision) of equilibrium partial Mantel r estimates. We conclude that amplifying more (and more variable) loci is likely to increase the power of landscape genetic inferences more than increasing number of individuals.     

Isolation and characterization of novel microsatellite loci for the short‐beaked common dolphin (Delphinus delphis) and cross‐amplification in other cetacean species

Seven tetranucleotide and three dinucleotide microsatellite loci were isolated from the short‐beaked common dolphin (Delphinus delphis). Number of alleles and observed heterozygosity ranged from three to twelve (mean 8.4) and from 0.087 to 0.978 (mean 0.741), respectively. All loci were tested in four other odontocete species and showed polymorphisms in most locus–species combinations. This panel of nuclear markers will be useful for the investigation of population ecology and social interactions in common dolphins and other cetacean species.     

Birth weight and neonatal survival of harbour seal pups are positively correlated with genetic variation measured by microsatellites.

We examined the relations between fitness-related traits of wild harbour seal (Phoca vitulina) pups with microsatellite heterozygosity, and with a measure of genomic diversity based on the mean squared distance between microsatellite alleles within an individual, mean d2. Birth weight was positively influenced by maternal age, pup sex, and either mean d2 or individual heterozygosity in separate multiple regression models. The association of birth weight with mean d2 was stronger than that with heterozygosity, however. The factors maternal age, pup sex, and mean d2 combined to account for 36.8% of the variation in birth weight, with mean d2 accounting for the greatest explanatory power (52.3% of the variance explained). Pups which survived until weaning had significantly higher mean d2 than pups which died, independent of birth weight. These effects are consistent with heterosis resulting from recent population mixing, and/or inbreeding depression in this population. Mean d2 thus provides (i) a better measure of individual genetic variability than heterozygosity for microsatellite data; and (ii) a convenient tool for assessing the effects of inbreeding and outbreeding in natural populations.     

Microsatellite markers for the tambaqui (Colossoma macropomum, Serrasalmidae, Characiformes), an economically important keystone species of the Amazon River floodplain

Colossoma macropomum is a keystone species of the Amazon floodplain, and is an important but severely overexploited commercial species. To provide tools for addressing ecological and conservation questions, we developed 14 highly polymorphic microsatellite markers that had between four and 21 alleles per locus in the 25 tested individuals. With the exception of comparisons involving the locus Cm1F5 that also showed heterozygosity deficiency, no pairs of loci were at linkage disequilibrium. Many of the microsatellite loci were also variable in three other serrasalmid species which span the phylogenetic depth of the Serrasalmidae.