Identification of candidate AFLP markers for shell color of the Pacific oyster (Crassostrea gigas) under artificial selection

The shell color of the Pacific oyster (Crassostrea gigas) is a desirable trait, but only a few genetic studies on shell color have been documented. Through successive selective breeding, four shell color variants of white (W), gold (G), black (B) and purple (P) C. gigas have been developed. The amplified fragment length polymorphism (AFLP) technique was used to scan the genomes of the four variants with different shell colors and one wild population (C) to identify candidate markers for shell polymorphism. Fifteen AFLP primer combinations were used, 1079 loci were scored as polymorphic loci, and the percentage of polymorphic bands was 95.5%. In the gold, white, black, purple and wild populations, the percentages of polymorphic loci were estimated to be 90.5% (G), 90.0% (W), 91.1% (B), 95.3% (P) and 93.2% (C); the expected heterozygosity values were 0.3115 (G), 0.3044 (W), 0.3102 (B), 0.3285 (P) and 0.3105 (C). The white shell variant was observed to have slightly lower genetic diversity than others, with a F_ST value of 0.1483. These results indicated that the four different shell color variants had high genetic diversity and that the genetic differentiation of populations mostly results from genetic diversity of individuals within populations. Furthermore, 11 outlier loci were considered candidate markers for shell color. This work provides new insights on relationships among color variants of C. gigas.     

Genetic Mapping and QTL Analysis of Growth-Related Traits in Pinctada fucata Using Restriction-Site Associated DNA Sequencing

The pearl oyster, Pinctada fucata (P. fucata), is one of the marine bivalves that is predominantly cultured for pearl production. To obtain more genetic information for breeding purposes, we constructed a high-density linkage map of P. fucata and identified quantitative trait loci (QTL) for growth-related traits. One F1 family, which included the two parents, 48 largest progeny and 50 smallest progeny, was sampled to construct a linkage map using restriction site-associated DNA sequencing (RAD-Seq). With low coverage data, 1956.53 million clean reads and 86,342 candidate RAD loci were generated. A total of 1373 segregating SNPs were used to construct a sex-average linkage map. This spanned 1091.81 centimorgans (cM), with 14 linkage groups and an average marker interval of 1.41 cM. The genetic linkage map coverage, Coa, was 97.24%. Thirty-nine QTL-peak loci, for seven growth-related traits, were identified using the single-marker analysis, nonparametric mapping Kruskal-Wallis (KW) test. Parameters included three for shell height, six for shell length, five for shell width, four for hinge length, 11 for total weight, eight for soft tissue weight and two for shell weight. The QTL peak loci for shell height, shell length and shell weight were all located in linkage group 6. The genotype frequencies of most QTL peak loci showed significant differences between the large subpopulation and the small subpopulation (P<0.05). These results highlight the effectiveness of RAD-Seq as a tool for generation of QTL-targeted and genome-wide marker data in the non-model animal, P. fucata, and its possible utility in marker-assisted selection (MAS).     

Comparison of Two Pearl Sacs Formed in the Same Recipient Oyster with Different Genetic Background Involved in Yellow Pigmentation in <Emphasis Type="Italic">Pinctada fucata</Emphasis>

Color is one of the most important factors determining the commercial value of pearls. Pinctada fucata is a well-known pearl oyster producing high-quality Akoya pearls. Phenotypic variation in amount of yellow pigmentation produces white and yellowish pearls. It has been reported that polymorphism of yellow pigmentation of Akoya pearls is genetically regulated, but the responsible gene(s) has remained unknown. Here, we prepared pearl sac pairs formed in the same recipient oyster but coming from donor oysters that differ in their color. These two pearl sacs produced pearls with different yellowness even in the same recipient oyster. Yellow tone of produced pearls was consistent with shell nacre color of donor oysters from which mantle grafts were prepared, indicating that donor oysters strongly contribute to the yellow coloration of Akoya pearls. We also conducted comparative RNA-seq analysis and retrieved several candidate genes involved in the pearl coloration. Whole gene expression patterns of pair sacs were not grouped by pearl color they produced, but grouped by recipient oysters in which they were grown, suggesting that the number of genes involved in the yellow coloration is quite small, and that recipient oyster affects gene expression of the majority of genes in the pearl sac.     

Low genetic differentiation among widely separated populations of the pearl oyster Pinctada fucata as revealed by AFLP

Genetic variation within and among five populations of the pearl oyster Pinctada fucata, from China (Daya Bay, Sanya Bay and Beibu Bay), Japan (Mie Prefecture) and Australia (Port Stephens) was studied using AFLP. Three primer pairs generated 184 loci among which 91.8-97.3% is polymorphic. An overall genetic diversity of 0.38 among populations and an average of 0.37 within populations (ranging from 0.35 in Japanese population to 0.39 in Beibu Bay population) were observed. Genetic differentiation among the five populations is low but significant as indicated by pairwise G_ST (0.0079-0.0404). AMOVA further shows that differentiation is significant among the five populations but is not significant at a broader geographical scale, among the three groups of Chinese, Japanese and Australian populations or among the two groups of Australian and north Pacific populations. The low level of genetic differentiation indicated that P. fucata populations in the west Pacific are genetically linked. Among the five populations, the Australian one is more differentiated from the others, based on both pairwise AMOVA and G_ST analyses, and is genetically isolated by distance as indicated by Mantel test. However, genetic differences among the three Chinese populations are not correlated with the geographic distances, suggesting that Hainan Island and Leizhou Peninsula may act as barriers blocking gene flow.     

Molecular identification of AFLP fragments associated with elytra color variation of the Asian ladybird beetle,Harmonia axyridis (Coleoptera: Coccinellidae)

The Asian ladybird beetle, Harmonia axyridis shows polymorphism in elytra color patterns. However, it is uncertain whether these color patterns are regulated by genetic factors. This investigation used amplified fragment length polymorphism (AFLP) analysis to determine any genetic causes of the variability of color patterns. Using four individuals of each group, AFLP analysis produced 37 polymorphic bands. Among several polymorphic bands, six AFLP markers were associated with elytra color patterns after further analysis using six additional individuals of each group. These polymorphic sites were sequenced but did not match DNA sequence data deposited in GenBank. Based on the color-associated AFLP markers, SCAR primers were designed for PCR amplification of genomic DNA. These primers (SCAR 12 and SCAR 44) were used to analyze color-associated loci and/or alleles of H. axyridis DNA. SCAR 12 primers designed from a Spectabilis type-specific fragment (AFLP 12) amplified a specific band of 530 bp in four Spectabilis individuals, but not in the insects with other color patterns.     

A Novel Acidic Matrix Protein,PfN44, Stabilizes Magnesium Calcite to Inhibit the Crystallization of Aragonite

Magnesium is widely used to control calcium carbonate deposition in the shell of pearl oysters. Matrix proteins in the shell are responsible for nucleation and growth of calcium carbonate crystals. However, there is no direct evidence supporting a connection between matrix proteins and magnesium. Here, we identified a novel acidic matrix protein named PfN44 that affected aragonite formation in the shell of the pearl oyster Pinctada fucata. Using immunogold labeling assays, we found PfN44 in both the nacreous and prismatic layers. In shell repair, PfN44 was repressed, whereas other matrix proteins were up-regulated. Disturbing the function of PfN44 by RNAi led to the deposition of porous nacreous tablets with overgrowth of crystals in the nacreous layer. By in vitro circular dichroism spectra and fluorescence quenching, we found that PfN44 bound to both calcium and magnesium with a stronger affinity for magnesium. During in vitro calcium carbonate crystallization and calcification of amorphous calcium carbonate, PfN44 regulated the magnesium content of crystalline carbonate polymorphs and stabilized magnesium calcite to inhibit aragonite deposition. Taken together, our results suggested that by stabilizing magnesium calcite to inhibit aragonite deposition, PfN44 participated in P. fucata shell formation. These observations extend our understanding of the connections between matrix proteins and magnesium.     

Identification and Mapping of Amplified Fragment Length Polymorphism Markers Linked to Shell Color in Bay Scallop, Argopecten irradians irradians (Lamarck, 1819)

Amplified fragment length polymorphisms (AFLP) were used to study the inheritance of shell color in Argopecten irradians. Two scallops, one with orange and the other with white shells, were used as parents to produce four F₁ families by selfing and outcrossing. Eighty-eight progeny, 37 orange and 51 white, were randomly selected from one of the families for segregation and mapping analysis with AFLP and microsatellite markers. Twenty-five AFLP primer pairs were screened, yielding 1138 fragments, among which 148 (13.0%) were polymorphic in two parents and segregated in progeny. Six AFLP markers showed significant (P < 0.05) association with shell color. All six loci were mapped to one linkage group. One of the markers, F1f335, is completely linked to the gene for orange shell, which we designated as Orange1, without any recombination in the progeny we sampled. The marker was amplified in the orange parent and all orange progeny, but absent in the white parent and all the white progeny. The close linkage between F1f335 and Orange1 was validated using bulk segregation analysis in two natural populations, and all our data indicate that F1f335 is specific for the shell color gene, Orange1. The genomic mapping of a shell color gene in bay scallop improves our understanding of shell color inheritance and may contribute to the breeding of molluscs with desired shell colors.     

Challenges and pitfalls in the characterization of anonymous outlier AFLP markers in non-model species: lessons from an ocellated lizard genome scan

In the last few years, dozens of studies have documented the detection of loci influenced by selection from genome scans in a wide range of non-model species. Many of those studies used amplified fragment length polymorphism (AFLP) markers, which became popular for being easily applicable to any organism. However, because they are anonymous markers, AFLPs impose many challenges for their isolation and identification. Most recent AFLP genome scans used capillary electrophoresis (CE), which adds even more obstacles to the isolation of bands with a specific size for sequencing. These caveats might explain the extremely low number of studies that moved from the detection of outlier AFLP markers to their actual isolation and characterization. We document our efforts to characterize a set of outlier AFLP markers from a previous genome scan with CE in ocellated lizards (Lacerta lepida). Seven outliers were successfully isolated, cloned and sequenced. Their sequences are noncoding and show internal indels or polymorphic repetitive elements (microsatellites). Three outliers were converted into codominant markers by using specific internal primers to sequence and screen population variability from undigested DNA. Amplification in closely related lizard species was also achieved, revealing remarkable interspecific conservation in outlier loci sequences. We stress the importance of following up AFLP genome scans to validate selection signatures of outlier loci, but also report the main challenges and pitfalls that may be faced during the process.     

Differential community development of fouling species on the pearl oysters Pinctada fucata,Pteria penguin and Pteria chinensis (Bivalvia,Pteriidae)

Abstract

A field experiment documented the development of fouling communities on two shell regions, the lip and hinge, of the pearl oyster species Pinctada fucata, Pteria penguin and Pteria chinensis. Fouling communities on the three species were not distinct throughout the experiment. However, when each species was analysed separately, fouling communities on the lip and hinge of P. penguin and P. chinensis were significantly different during the whole sampling period and after 12 weeks, respectively, whereas no significant differences could be detected for P. fucata. There was no significant difference in total fouling cover between shell regions of P. fucata and P. chinensis after 16 weeks; however, the hinge of P. penguin was significantly more fouled than the lip. The most common fouling species (the hydroid Obelia bidentata, the bryozoan Parasmittina parsevalii, the bivalve Saccostrea glomerata and the ascidian Didemnum sp.) showed species-specific fouling patterns with differential fouling between shell regions for each species. The role of the periostracum in determining the community development of fouling species was investigated by measuring the presence and structure of the periostracum at the lip and hinge of the three pearl oyster species. The periostracum was mainly present at the lip of the pearl oysters, while the periostracum at the hinge was absent and the underlying prismatic layer eroded. The periostracum of P. fucata lacked regular features, whereas the periostracum of P. penguin and P. chinensis consisted of a regular strand-like structure with mean amplitudes of 0.84 μm and 0.65 μm, respectively. Although the nature and distribution of fouling species on the pearl oysters was related to the presence of the periostracum, the periostracum does not offer a fouling-resistant surface for these pearl oyster species.     

Development and characterization of six new microsatellite markers for the silver‐ or gold‐lipped pearl oyster,Pinctada maxima (Pteriidae)

Six di‐, tri‐ and tetranucleotide microsatellite loci were developed for the silver‐ or gold‐lipped pearl oyster Pinctada maxima using a linker‐ligated, magnetic bead enrichment protocol. Based on a minimum of 134 Indonesian pearl oyster samples, number of alleles and observed heterozygosity at each locus ranged from six to 17 alleles and from 0.172 to 0.813 (mean = 0.448), respectively. Mean polymorphic information content for the six loci was 0.562. These loci should be very useful in DNA parentage analyses and population differentiation of P. maxima in Australia and Indonesia.     

Assessing genetic diversity of populations of topmouth culter (Culter alburnus) in China using AFLP markers

Genetic diversity of five wild populations and a cultured population of topmouth culter (Culter alburnus) was investigated using amplified fragment length polymorphism (AFLP). A total of 373 reproducible bands amplified with seven AFLP primer combinations were obtained from 163 fish. The percentage of polymorphic loci ranged widely from 37.0% to 69.2% within distinct populations. The cultured population appeared to have a lower level of polymorphism (37.0%), gene diversity (0.121 ± 0.188) and Shannon's Information index (0.183 ± 0.263) than the wild populations. Analysis of molecular variance (AMOVA) revealed that average F_ST value overall loci was 0.2671, and the percentage of variation within population (73.29%) was larger than among populations (26.71%) (P < 0.01). The six populations were clustered into two major clades with UPGMA. The results from analysis of population pairwise gene flow indicated moderate gene flow among populations. Our study indicated that the genetic diversity of the cultured population was reduced compared with the wild populations. Geographic isolation, habitat, and artificial selection all may have played important roles in population differentiation. The information may be beneficial to future broodstock selection and defining conservation management for the different populations of topmouth culter.     

A multilocus sequence analysis of Xanthomonas campestris reveals a complex structure within crucifer-attacking pathovars of this species

Previous classification of Xanthomonas campestris has defined six pathovars (aberrans, armoraciae, barbareae, campestris, incanae, and raphani) that cause diseases on cruciferous plants. However, pathogenicity assays with a range of strains and different hosts identifies only three types of symptom: black rot, leaf spot and bacterial blight. These findings raise the question of the genetic relatedness between strains assigned to different pathovars or symptom phenotypes. Here we have addressed this issue by multilocus sequence analysis of 42 strains. The X. campestris species was polymorphic at the 8 loci analysed and had a high genetic diversity; 23 sequence types were identified of which 16 were unique. All strains that induce black rot (pathovars aberrans and campestris) were genetically close but split in two groups. Only three clonal complexes were found, all within pathovar campestris. The assignment of the genome-sequenced strain 756C to pathovar raphani suggested from disease symptoms was confirmed, although this group of strains was particularly polymorphic. Strains belonging to pathovars barbareae and incanae were closely related, but distinct from pathovar campestris. There is evidence of genetic exchanges of housekeeping genes within this species as deduced from a clear incongruence between individual gene phylogenies and from network structures from SplitsTree analysis. Overall this study showed that the high genetic diversity derived equally from recombination and point mutation accumulation. However, X. campestris remains a species with a clonal evolution driven by a differential adaptation to cruciferous hosts.     

AFLPs and Mitochondrial Haplotypes Reveal Local Adaptation to Extreme Thermal Environments in a Freshwater Gastropod

The way environmental variation shapes neutral and adaptive genetic variation in natural populations is a key issue in evolutionary biology. Genome scans allow the identification of the genetic basis of local adaptation without previous knowledge of genetic variation or traits under selection. Candidate loci for divergent adaptation are expected to show higher F_ST than neutral loci influenced solely by random genetic drift, migration and mutation. The comparison of spatial patterns of neutral markers and loci under selection may help disentangle the effects of gene flow, genetic drift and selection among populations living in contrasting environments. Using the gastropod Radix balthica as a system, we analyzed 376 AFLP markers and 25 mtDNA COI haplotypes for candidate loci and associations with local adaptation among contrasting thermal environments in Lake Mývatn, a volcanic lake in northern Iceland. We found that 2% of the analysed AFLP markers were under directional selection and 12% of the mitochondrial haplotypes correlated with differing thermal habitats. The genetic networks were concordant for AFLP markers and mitochondrial haplotypes, depicting distinct topologies at neutral and candidate loci. Neutral topologies were characterized by intense gene flow revealed by dense nets with edges connecting contrasting thermal habitats, whereas the connections at candidate loci were mostly restricted to populations within each thermal habitat and the number of edges decreased with temperature. Our results suggest microgeographic adaptation within Lake Mývatn and highlight the utility of genome scans in detecting adaptive divergence.     

Genetic diversity of the expansive grass Brachypodium pinnatum in a changing landscape: Effect of habitat age

Perennial grasses constitute a major group of species showing a dramatic decline of biodiversity in successional plant communities. Using AFLP markers, we examined 12 populations of the expansive grass Brachypodium pinnatum differing in habitat age (30–50, ca. 100 and >300 years old) in order to determine whether clonal diversity of populations, genetic variation, and the relative importance of clonal propagation versus sexual reproduction change with grassland age. Five AFLP primer combinations gave a total of 517 bands, 79% of which were polymorphic. 314 different multilocus lineages were distinguished among the 453 samples analyzed. The number of genotypes (G) and clonal richness (R) decreased with habitat age, while the distribution of the frequency of genets changed from many clones of similar size to dominance by one or a few large clones. We consider these results to give evidence of significant role of sexual reproduction in the early phases of colonization and prevalence of clonal growth and competitive exclusion of less adapted genotypes in the later ones. However, habitat age had only marginal effect on genetic diversity, as percentage of polymorphic loci (PPL) within all the populations analyzed was similar, viz. 38.6–43.5%.     

Development and characterization of novel tetra‐ and dinucleotide microsatellite markers for the French Polynesia black‐lipped pearl oyster,Pinctada margaritifera

Ten microsatellite loci were developed for the black‐lipped pearl oyster Pinctada margaritifera with a magnetic bead enrichment protocol. These tetra‐ and dinucleotide markers were polymorphic, with 10 to 43 alleles observed in 97 individuals from two Tuamotu atoll populations. Most loci revealed significant genic differentiation between the two populations and also exhibited some degree of heterozygote deficiencies, probably due to the presence of null alleles. These loci should be very useful to describe genetic structure, genetic variability and reproductive success in the various aquaculture and wild populations of pearl oyster in French Polynesia.     

An intraspecific genetic map of velvetbean (<Emphasis Type="Italic">Mucuna</Emphasis> sp.) based on AFLP markers

Velvetbean (Mucuna sp., n=11), a self-pollinated species, is an important legume used in tropical agricultural systems in rotation with other crops for nematode management and/or soil improvement. A genetic map of velvetbean was constructed in order to identify potential molecular markers linked to important morphological and agronomic traits that would be particularly useful for developing and improving the species. Traits such as seed coat color, pod color, and pod pubescence were among the main parameters observed in a process of genetic diversity estimation. Two slightly divergent velvetbean accessions, PI364362 and Edgar Farm White, a land race from Alabama, were used to make an intraspecific F₁ hybrid. Amplified fragment length polymorphism analysis (AFLP) detected an average of six polymorphic fragments per primer pair between the two parents. As expected for dominant markers, the sum of all AFLP bands from both parents was generally observed to be present in the AFLP profiles of the F₁ progeny, indicating full penetrance and the dominant nature of AFLP markers. An F₂ population was generated by self-pollinating a single F₁ plant. Using 37 AFLP primer pairs, we detected 233 polymorphic markers of which 164 (70.4%) segregated in 3:1 Mendelian ratios, while the remaining 69 (29.6%) both segregated and were scorable. The genetic linkage map constructed from this population comprised 166 markers, including two morphological traits (pod color and pod pubescence). Twenty linkage groups were found with an average distance between markers of 34.4 cM, covering a total of 687.9 cM. The linkage groups contained from 2 to 12 loci each and the distance between two consecutive loci ranged from 0 to 21.8 cM. The newly designated morphological traits pod color (pdc) and pod pubescence (pdp) co-segregated with each other at a distance of 4.2 cM. Two DNA markers designated ACGCAG2 and ACTCTG1 were located in the same group as pdc and pdp. The AFLP linkage map provides opportunities for use in marker-assisted selection and in the detection of loci controlling morphologically important traits.Communicated by J. Dvorak