Genotyping the clonal structure of a gorgonian coral, Junceella juncea (Anthozoa: Octocorallia), using microsatellite loci

The identification of different clones is fundamental to the study of population structure among organisms with mixed reproductive modes such as cnidarians. However, due to the low genetic variation of coral mtDNA and contamination by zooxanthellate DNA, very few molecular markers are available for studying the clonal structure of cnidarians. Herein we used four polymorphic loci of microsatellite DNA isolated from a zooxanthellae-free octocoral, Junceella juncea, to study its clonal structure in seven populations collected from three localities in Taiwan. In total, 40 multilocus genotypes were found among 152 colonies, and the number of genotypes (clones) identified in the seven populations ranged from 2 to 16. Each of the 40 multilocus genotypes was restricted to a single population, even where adjacent populations were only 100 m distant. The ratio of observed to expected genotypic diversity (G_o:G_e) ranged from 0.217 to 0.650, and G_o showed a significant departure from G_e (p<0.05) at each site indicating that asexual fragmentation may play a major role in the maintenance of established populations. Mean relatedness (R) values showed that genotypes within reefs were more closely related than those between regions. The results indicate that microsatellites are useful for discerning the clonal structures among and within populations at different spatial scales. Electronic supplement: Unique multilocus genotypes (clones) revealed by 4 polymorphic loci for Junceella juncea colonies collected from Xiashuijui (Reefs A, B, C, Transplant, Transect), Nanwan and Shicheng     

Genotypic diversity,molecular markers and spatial distribution of genets in clonal plants,a literature survey

We present a literature survey of studies using molecular markers to investigate genet diversity and structure in clonal plants. The data from 40 studies comprised 45 species of which only two were studied by DNA methods, the rest by isozyme markers. Less than one third of the studies provided information about the spatial distribution of individual genets within populations, and only 12.5% of the studies used mapping of all ramets within plots or part of the population in combination with identification of multilocus genotypes. We also present two case studies. InGlechoma hederacea morphological criteria were used to select clones. Multi-samples of ramets within these “clones” turned out to be variable using isozymes indicating that these “clones” in most cases consisted of several genets. One individual multilocus genotype covered tens of square metres. InHylocomium splendens samples from 10×10 cm plots usually consisted of a mixture of multilocus genotypes, but occasionally the whole plot consisted of one genet.     

Evolution and virulence of serogroup 6 pneumococci on a global scale

To study the evolution and virulence of pneumococcal populations, we used multilocus sequence typing to identify the major clones among 212 carriage and invasive isolates expressing capsular serogroup 6 from 39 countries. The global population consisted of 8 major complexes and 6 minor complexes of related clones and 32 clones of diverse origin. Surprisingly, serotype 6A clones evolved by mutation nearly as often as by recombination, whereas serotype 6B clones evolved almost exclusively by recombination (P = 0.0029). This is the first report of population genetic differences among serotypes of this species. The largest clonal complex was associated with invasive disease (P = 0.019) and included a common ancestor for five previously identified drug-resistant clones. The putative ancestors of the major clonal complexes were represented by a greater proportion of carriage isolates than were their descendents (P = 0.001), and the ancestors tended to be less virulent than their descendents in a mouse model of infection. These data suggested that virulent serogroup 6 clones have evolved multiple times from less-virulent ancestral clones.     

Population structure and evolutionary dynamics of pathogenic bacteria

Evidence concerning the significance of recombination within natural bacterial populations has historically come from two main sources: multilocus enzyme electrophoresis (MLEE) and nucleotide sequence data. Here we discuss evidence from a third method, multilocus sequence typing (MLST), which is a development of MLEE based on nucleotide sequencing that combines the advantages of both approaches. MLST has confirmed both the existence of clones and the high rates of recombination for several bacterial pathogens. The data are consistent with "epidemic" population structures, where clones are superimposed upon a backdrop of frequent recombination, thus, in the short term, resisting the homogenising effect of recombination. The nature of the selective advantage of clones, however, and how this advantage relates to virulence are unclear. The current evidence also has broader implications concerning bacterial species definition, the management of antibiotic-resistant bacteria and the assessment of the dangers of releasing genetically modified organisms into the environment.     

Clonal diversity and host distribution in Bordetella bronchiseptica.

A total of 303 isolates of Bordetella bronchiseptica recovered from 11 host species were characterized by the electrophoretic mobilities of 15 metabolic enzymes, and 21 distinctive multilocus genotypes (electrophoretic types) were distinguished on the basis of allele profiles at the enzyme loci. The population structure of B. bronchiseptica is clonal, and its genetic diversity is limited in comparison with most other pathogenic bacteria, perhaps reflecting a relatively recent origin of the species. Electrophoretic types mark clones which are, in many cases, nonrandomly associated with host species. Clones differing only slightly in overall chromosomal genetic character may have pronounced differences in virulence potential. There was considerable variation among individual clones and clone families in degree of host specificity and among various species of hosts in the diversity of clones causing disease. The diversity of clones infecting dogs was an order of magnitude greater than that of clones infecting pigs. Most bordetellosis in pigs in the United States and Japan was found to be caused by strains of a single multilocus genotype.     

Zygotic associations and multilocus statistics in a nonequilibrium diploid population

The usual approach to characterizing and estimating multilocus associations in a diploid population assumes that the population is in Hardy-Weinberg equilibrium. The purpose of this study is to develop a set of summary statistics that can be used to characterize and estimate the multilocus associations in a nonequilibrium population. The concept of "zygotic associations" is first expanded to facilitate the development. The summary statistics are calculated using the distribution of a random variable, the number of heterozygous loci (K) found in diploid individuals in the population. In particular, the variance of K consists of single-locus and multilocus components with the latter being the sum of zygotic associations between pairs of loci. Simulation results show that the multilocus associations in the variance of K are detectable in a sample of moderate size (> or =30) when the sum of all pairwise zygotic associations is greater than zero and when gene frequency is intermediate. The method presented here is a generalization of the well-known development for the Hardy-Weinberg equilibrium population and thus may be of more general use in elucidating the multilocus organizations in nonequilibrium and equilibrium populations.     

Ecological genetics of a cyclical parthenogen in temporary habitats

Populations of the rotifer Brachionus plicatilis inhabiting three temporary ponds in Torreblanca Marsh (Castellón, Spain) were regularly screened for allozyme variation, sexual reproduction levels and population densities during an annual cycle. Relevant ecological parameters in the ponds were also recorded. The electrophoretic survey of the three ponds (Poza Sur, Poza Norte and Canal Central) revealed a high level of overall genetic polymorphism in four marker loci, but only 13 multilocus genotypes were found. We classified clones into three clonal groups (SS, SM, L) characterized by unique arrays of alleles in the four marker loci, and significant differences in body shape and size. Clonal group succession took place in Poza Sur; SM clones occurring in spring, SS clones from spring to fall, and L clones from fall to spring. Despite the partial overlapping of sexual periods, the absence of heterozygotes indicates that gene flow is strongly restricted between clonal groups. In the transition periods, the population was far from Hardy—Weinberg equilibrium due to hetereozygote deficiencies in all four markers. In clonal group L, the only polymorphic locus, Pgi, was in Hardy—Weinberg equilibrium. In clonal group SM, both polymorphic loci were in Hardy—Weinberg equilibrium in Poza Norte, but not in Poza Sur. Factor analysis on the limnological parameters recorded during field sampling indicates that genetically different clonal groups are also ecologically specialized. These findings suggest that the taxonomical species B. plicatilis is in fact a species complex. The complete genetic discontinuity model of rotifer population succession is supported by these data. Our results are similar to those found in cladocerans inhabiting temporary habitats.     

Clonal and spatial genetic structure within populations of a coastal plant, Carex kobomugi (Cyperaceae)

Clarification of clonal growth pattern is critical for understanding the population dynamics and reproductive system evolution of clonal plant species. The contribution of clonality to the spatial genetic structure (SGS) within populations is also an important issue. I examined the spatial distribution of genetic variability within two populations of the coastal plant Carex kobomugi using seven microsatellite loci. Genotyping of 226 and 140 ramets within 14 × 40 m and 14 × 34 m plots on two populations revealed 36 and 33 multilocus genotypes, respectively. To quantify the extent of intermingling among clones, for each genet, I calculated the dominance of ramets belonging to a particular genet within a spatial range of the genet. Furthermore, I analyzed spatial distribution of genotypes within 2 × 2 m and 1 × 2 m quadrats using second-order spatial statistics. These analyses indicated that clones are highly intermingled, suggesting a low level of spatial interaction among clones. Spatial autocorrelation analysis of kinship coefficient including all pairs of ramets showed significantly stronger SGS than analysis considering only pairs between different genets. I conclude that clonal propagation largely contributes to SGS at a fine scale.     

Evidence for a monophyletic origin of triploid clones of the Amazon molly, Poecilia formosa

Asexual reproduction in vertebrates is rare and generally considered an evolutionary dead end. Asexuality is often associated with polyploidy, and several hypotheses have been put forward to explain this relationship. So far, it remains unclear whether polyploidization in asexual organisms is a frequent or a rare event. Here we present a field study on the gynogenetic Amazon molly, Poecilia formosa. We used multilocus fingerprints and microsatellites to investigate the genetic diversity in 339 diploid and 55 triploid individuals and in 25 P. mexicana, its sexual host. Although multilocus DNA fingerprints found high clonal diversity in triploids, microsatellites revealed only two very similar clones in the triploids. Phylogenetic analysis of microsatellite data provided evidence for a monophyletic origin of the triploid clones of P. formosa. In addition, shared alleles within the triploid clones between the triploid and diploid genotypes and between asexual and sexual lineages indicate a recent origin of triploid clones in Poecilia formosa.     

Genetics of the polycross

Summary Rates and patterns of male gamete incorporation for a polycross mating design were studied for two independent years of pollination in Norway spruce, Picea abies (L) Karst. Segregation distortion in a subset of maternal clones was documented for one locus. We have proposed a model, involving the existence of a linked lethal allele, which accounts for these observations. Significant temporal and maternal clonal differences were observed in the rates at which single locus and multilocus gametes were incorporated. Striking differences in apparent fertility existed among four clones which produced unique multilocus gametes. One clone, in particular, was shown to be contributing three times as many gametes to the next generation as predicted by the hypothesis of equal clonal male contribution. These deviations from expectation were also detected in the genotypic distributions of the resultant filial generation. Ramifications of these results on family structures in the filial generation, effective size of the male population, and possible bias in inferences of genetic differences and parameter estimation are discussed.     

Ecology and genetic structure of a northern temperate Vibrio cholerae population related to toxigenic isolates

Although Vibrio cholerae is an important human pathogen, little is known about its populations in regions where the organism is endemic but where cholera disease is rare. A total of 31 independent isolates confirmed as V. cholerae were collected from water, sediment, and oysters in 2008 and 2009 from the Great Bay Estuary (GBE) in New Hampshire, a location where the organism has never been detected. Environmental analyses suggested that abundance correlates most strongly with rainfall events, as determined from data averaged over several days prior to collection. Phenotyping, genotyping, and multilocus sequence analysis (MLSA) revealed a highly diverse endemic population, with clones recurring in both years. Certain isolates were closely related to toxigenic O1 strains, yet no virulence genes were detected. Multiple statistical tests revealed evidence of recombination among strains that contributed to allelic diversity equally as mutation. This relatively isolated population discovered on the northern limit of detection for V. cholerae can serve as a model of natural population dynamics that augments predictive models for disease emergence.     

identix,a software to test for relatedness in a population using permutation methods

The computer program identix estimates relatedness in natural populations using multilocus genotypic data. Queller & Goodnight's (1989) and Lynch & Ritland's (1999) estimators of pairwise relatedness are implemented, as well as the identity index of Mathieu et al. (1990). Estimates of the confidence intervals around these pairwise values are also provided. The null hypothesis of no relatedness (multilocus genotypes are independent draws from a panmictic population) is tested using a permutation method that compares the observed distribution of the moments of pairwise relatedness coefficients to that expected in unstructured population.     

Ecological specialization of the aphid Aphis gossypii Glover on cultivated host plants

Many plant-feeding insect species considered to be polyphagous are in fact composed of genetically differentiated sympatric populations that use different hosts and between which gene flow still exists. We studied the population genetic structure of the cotton-melon aphid Aphis gossypii that is considered as one of the most polyphagous aphid species. We used eight microsatellites to analyse the genetic diversity of numerous samples of A. gossypii collected over several years at a large geographical scale on annual crops from different plant families. The number of multilocus genotypes detected was extremely low and the genotypes were found to be associated with host plants. Five host races were unambiguously identified (Cucurbitaceae, cotton, eggplant, potato and chili- or sweet pepper). These host races were dominated by asexual clones. Plant transfer experiments using several specialized clones further confirmed the existence of host-associated trade-offs. Finally, both genetic and experimental data suggested that plants of the genus Hibiscus may be used as refuge for the specialized clones. Resource abundance is discussed as a key factor involved in the process of ecological specialization in A. gossypii.     

Sex-biased dispersal in a migratory bat: a characterization using sex-specific demographic parameters

We studied the noctule bat (Nyctalus noctula), in which the mitochondrial F(ST) is about 10 times that revealed by nuclear markers, to address two questions. We first verified whether random dispersal of one sex is compatible with highly contrasted mitochondrial and nuclear population structures. Using computer simulations, we then assessed the power of multilocus population differentiation tests when the expected population structure departs only slightly from panmixia. Using an island model with sex-specific demographic parameters, we found that random male dispersal is consistent with the population structure observed in the noctule. However, other parameter combinations are also compatible with the data. We computed the minimum sex bias in dispersal (at least 69% of the dispersing individuals are males), a result that would not be available if we had used more classical population genetic models. The power of multilocus population differentiation tests was unexpectedly high, the tests being significant in almost 100% of the replicates, although the observed population structure infered from nuclear markers was extremely low (F(ST) = 0.6%).     

Population genetic analysis of Pseudomonas aeruginosa using multilocus sequence typing

To study the population genetic structure of Pseudomonas aeruginosa, we developed a multilocus sequence typing scheme. The sequences of internal fragments of seven housekeeping genes were obtained for 34 P. aeruginosa isolates from patients hospitalized in five different European cities. Twenty-six different allelic profiles were identified. The mean allelic diversity was 0.854 (range: 0.606-0.978), which was about six times greater than the results obtained with the multilocus enzyme electrophoresis method. Linkage disequilibrium was measured with the index of association. An index of 1.95+/-0.24 was calculated when all the strains were considered. This index was 1.76+/-0.27 when only one strain per sequence type was considered. Both results were different from 0, indicating linkage among loci, which means that the population structure of our set of P. aeruginosa isolates is clonal. The clonal structure of the population was also suggested by the congruence of the topology of the different trees obtained from the seven housekeeping genes. These results are in contrast to previous studies, finding a non clonal population structure. Since a small number of isolates was analyzed in this study, there might be a bias of selection which includes the possibility that they belong to widely disseminated epidemic clones. Another possibility is that recombination did not occurred homogeneously throughout the genome of P. aeruginosa, so that part of it has a clonal structure, while the remaining part of the genome is more frequently subject to recombination.     

Clonal diversity in a Rhododendron ferrugineum L. (Ericaceae) population inferred from AFLP markers

In the European Alps, Rhododendron ferrugineum can constitute dense populations with almost 100% of cover. The developmental pattern by layering and the resulting complexity of population structure make it difficult to identify distinct clones even by excavation. Therefore genotypic structure of a R. ferrugineum population, in the French Alps, was inferred from AFLP markers. In a first step, we analysed 400 samples using AFLP profiles generated by one selective primer pair. Seventeen bands out of 25 were polymorphic (68%). We identified a total of 32 multilocus genotypes. In a second step, the 32 genotypes were verified by applying two additional primer pairs to the two most distant samples from each genotype. The mean similarity (proportion of band sharing) between pairs of clones was 0.85 (range from 0.52 to 0.94). The spatial distribution of clones showed that vegetative spreading mainly occurred down a slope. Based on an annual shoot mean growth of 2.6 cm/year and the size of the widest clone, we estimated the age of the oldest individual to be at least 300 years. A single genotype can occupy a large surface and sometimes form a dense patch, suggesting that this species adopts a phalanx growth form with limited intermingling of some genets.     

POPULATION STRUCTURE OF A CLONAL GORGONIAN CORAL: THE INTERPLAY BETWEEN CLONAL REPRODUCTION AND DISTURBANCE

Clonality is a common feature of plants and benthic marine organisms. In some cases clonal propagation results in a modest increase in population density, while in other cases dense populations may be generated by the propagation of only a few clones. We analyzed the population structure of the clonal gorgonian Plexaura kuna across several reef habitats with a range of disturbance regimes in the San Blas Islands, Panama, and the Florida Keys, U.S.A. Using multilocus DNA fingerprinting to distinguish clones, we estimated that clones ranged in size from single individuals to 500 colonies. The number of genotypes identified on nine reefs ranged from three to 25. Population density and clonal structure varied markedly among reefs with G_O:G_E ranging from 0.03 to 1.00. On some reefs vegetative reproduction transformed P. kuna from a rare species to the numerically most abundant gorgonian. The effect of clonal propagation on P. kuna population structure was dependent on interactions between fragmentation and the reef environment (disturbance regime, substratum). We present a generalized model relating population structure of clonal species to disturbance and the mode of vegetative propagation. Disturbance promotes colony propagation and skews the size-frequency distribution of clones among P. kuna and many species that propagate via fragmentation. Propagation of these species is promoted by disturbance (disturbance sensitive), and they tend to have clones that are dispersed across local sites. Species that fragment and have dispersed clones, have high genotypic diversity in habitats with low levels of disturbance. Genotypic diversity then decreases at intermediate disturbance and increases again at the highest disturbance levels. Clonal species that do not rely on disturbance for vegetative propagation (disturbance insensitive) generally do not disperse and form aggregated clones. Among these taxa disturbance has a greater affect on individual survival than on propagation. Genotypic diversity is directly related to the level of disturbance until very high levels of disturbance, at which time genotypic diversity declines.     

High level of polymorphism and spatial structure in a selfing plant species,Medicago truncatula (Leguminosae), shown using RAPD markers

Using RAPD markers and one morphological marker, we studied the among- and within-population structure in a selfing annual plant species, Medicago truncatula GAERTN. About 200 individuals, sampled from four populations subdivided into three subpopulations each, were scored for 22 markers. It was found that the within-population variance component accounted for 55% of the total variance, while the among-population variance component accounted for 45%. Eighteen percent of the total variance was due to within-population structure (i.e., among subpopulations). Thus, 37% of the total variance was within subpopulations. Using a multilocus approach, it was found that no multilocus genotype was common to two populations. Two of the four studied populations were composed of few (≤6) multilocus genotypes, whereas the other two had many (≥15) multilocus genotypes. In the most polymorphic population (37 genotypes), only one genotype was found to be common to two subpopulations. Resampling experiments show that, depending on the population, three to 16 polymorphic loci were necessary and sufficient to score all multilocus genotypes in the population. When these data are compared to published results, it appears that on some occasions, the number of genotypes per population of selfing species might be larger than would be expected from the sole consideration of effective population size. The large within-subpopulation genetic variance observed in some populations could be explained by either small neighborhood sizes within subpopulations, or by outcrossing following migration through seed and/or pollen.     

spip 1.0: a program for simulating pedigrees and genetic data in age‐structured populations

We present the program spip for simulating multilocus genetic data on individuals in age‐structured populations. In addition to genetic data on sampled individuals, the pedigree connecting all individuals in the population is recorded. This allows investigation of the relationship between family structure and population parameters. We foresee that spip will be useful for evaluating multilocus estimators of pairwise relatedness and population structure, and for simulating the distribution of relatedness in populations with varying demographies. It also provides a method for simulating genetic drift in complex populations.