PCR survey of 50 introns in animals: cross-amplification of homologous EPIC loci in eight non-bilaterian, protostome and deuterostome phyla

Exon Primed Intron Crossing (EPIC) markers provide molecular tools that are susceptible to be variable within species while remaining amplifiable by PCR using potentially universal primers. In this study we tested the possibility of obtaining PCR products from 50 EPIC markers on 23 species belonging to seven different phyla (Porifera, Cnidaria, Arthropoda, Nematoda, Mollusca, Annelida, Echinodermata) using 70 new primer pairs. A previous study had identified and tested those loci in a dozen species, including another phylum, Urochordata (Chenuil et al., 2010). Results were contrasted among species. The best results were achieved with the oyster (Mollusca) where 28 loci provided amplicons susceptible to contain an intron according to their size. This was however not the case with the other mollusk Crepidula fornicata, which seems to have undergone a reduction in intron number or intron size. In the Porifera, 13 loci appeared susceptible to contain an intron, a surprisingly high number for this phylum considering its phylogenetic distance with genomic data used to design the primers. For two cnidarian species, numerous loci (24) were obtained. Ecdysozoan phyla (arthropods and nematodes) proved less successful than others as expected considering reports of their rapid rate of genome evolution and the worst results were obtained for several arthropods. Some general patterns among phyla arose, and we discuss how the results of this EPIC survey may give new insights into genome evolution of the study species. This work confirms that this set of EPIC loci provides an easy-to-use toolbox to identify genetic markers potentially useful for population genetics, phylogeography or phylogenetic studies for a large panel of metazoan species. We then argue that obtaining diploid sequence genotypes for these loci became simple and affordable owing to Next-Generation Sequencing development. Species surveyed in this study belong to several genera (Acanthaster, Alvinocaris, Aplysina, Aurelia, Crepidula, Eunicella, Hediste, Hemimysis, Litoditis, Lophelia, Mesopodopsis, Mya, Ophiocten, Ophioderma, Ostrea, Pelagia, Platynereis, Rhizostoma, Rimicaris), two of them, belonging to the family Vesicomydae and Eunicidae, could not be determined at the genus level.

QTL for resistance to summer mortality and OsHV-1 load in the Pacific oyster (Crassostrea gigas)

Summer mortality is a phenomenon severely affecting the aquaculture production of the Pacific oyster (Crassostrea gigas). Although its causal factors are complex, resistance to mortality has been described as a highly heritable trait, and several pathogens including the virus Ostreid Herpes virus type 1 (OsHV-1) have been associated with this phenomenon. A QTL analysis for survival of summer mortality and OsHV-1 load, estimated using real-time PCR, was performed using five F(2) full-sib families resulting from a divergent selection experiment for resistance to summer mortality. A consensus linkage map was built using 29 SNPs and 51 microsatellite markers. Five significant QTL were identified and assigned to linkage groups V, VI, VII and IX. Analysis of single full-sib families revealed differential QTL segregation between families. QTL for the two-recorded traits presented very similar locations, highlighting the interest of further study of their respective genetic controls. These QTL show substantial genetic variation in resistance to summer mortality, and present new opportunities for selection for resistance to OsHV-1.

Combining two-stage testing and interval mapping strategies to detect QTL for resistance to bonamiosis in the european flat oyster Ostrea edulis

We have identified quantitative trait loci (QTL) in the flat oyster (Ostrea edulis) for resistance to Bonamia ostreae, a parasite responsible for the dramatic reduction in the aquaculture of this species. An F(2) family from a cross between a wild oyster and an individual from a family selected for resistance to bonamiosis was cultured with wild oysters injected with the parasite, leading to 20% cumulative mortality. Selective genotyping of 92 out of a total of 550 F(2) progeny (i.e., 46 heavily infected oysters that died and 46 parasite-free oysters that survived) was performed using 20 microsatellites and 34 amplification fragment length polymorphism primer pairs. Both a two-stage testing strategy and QTL interval mapping methods were used. The two-stage detection strategy had a high power with a low rate of false positives and identified nine and six probable markers linked to genes of resistance and susceptibility, respectively. Parent-specific genetic linkage maps were built for the family, spanning ten linkage groups (n = 10) with an observed genome coverage of 69-84%. Three QTL were identified by interval mapping in the first parental map and two in the second. Good concordance was observed between the results obtained after the two-stage testing strategy and QTL mapping.

Single Nucleotide polymorphisms and their relationship to codon usage bias in the Pacific oyster Crassostrea gigas

DNA sequence polymorphism and codon usage bias were investigated in a set of 41 nuclear loci in the Pacific oyster Crassostrea gigas. Our results revealed a very high level of DNA polymorphism in oysters, in the order of magnitude of the highest levels reported in animals to date. A total of 290 single nucleotide polymorphisms (SNPs) were detected, 76 of which being localised in exons and 214 in non-coding regions. Average density of SNPs was estimated to be one SNP every 60 bp in coding regions and one every 40 bp in non-coding regions. Non-synonymous substitutions contributed substantially to the polymorphism observed in coding regions. The non-synonymous to silent diversity ratio was 0.16 on average, which is fairly higher to the ratio reported in other invertebrate species recognised to display large population sizes. Therefore, purifying selection does not appear to be as strong as it could have been expected for a species with a large effective population size. The level of non-synonymous diversity varied greatly from one gene to another, in accordance with varying selective constraints. We examined codon usage bias and its relationship with DNA polymorphism. The table of optimal codons was deduced from the analysis of an EST dataset, using EST counts as a rough assessment of gene expression. As recently observed in some other taxa, we found a strong and significant negative relationship between codon bias and non-synonymous diversity suggesting correlated selective constraints on synonymous and non-synonymous substitutions. Codon bias as measured by the frequency of optimal codons for expression might therefore provide a useful indicator of the level of constraint upon proteins in the oyster genome.

AFLP-based genetic linkage maps of the blue mussel (Mytilus edulis)

We report the construction of the first genetic linkage map in the blue mussel, Mytilus edulis. AFLP markers were used in 86 full-sib progeny from a controlled pair mating, applying a double pseudo-test cross strategy. Thirty-six primer pairs generated 2354 peaks, of which 791 (33.6%) were polymorphic in the mapping family. Among those, 341 segregated through the female parent, 296 through the male parent (type 1:1) and 154 through both parents (type 3:1). Chi-square goodness-of-fit tests revealed that 71% and 73% of type 1:1 and 3:1 markers respectively segregated according to Mendelian inheritance. Sex-specific linkage maps were built with mapmaker 3.0 software. The female framework map consisted of 121 markers ordered into 14 linkage groups, spanning 862.8 cM, with an average marker spacing of 8.0 cM. The male framework map consisted of 116 markers ordered into 14 linkage groups, spanning 825.2 cM, with an average marker spacing of 8.09 cM. Genome coverage was estimated to be 76.7% and 75.9% for the female and male framework maps respectively, rising to 85.8% (female) and 86.2% (male) when associated markers were included. Twelve probable homologous linkage group pairs were identified and a consensus map was built for nine of these homologous pairs based on multiple and parallel linkages of 3:1 markers, spanning 816 cM, with joinmap 4.0 software.

Phylogeography of the green turtle, Chelonia mydas, in the Southwest Indian Ocean

Patterns of mitochondrial DNA (mtDNA) variation were used to analyse the population genetic structure of southwestern Indian Ocean green turtle (Chelonia mydas) populations. Analysis of sequence variation over 396 bp of the mtDNA control region revealed seven haplotypes among 288 individuals from 10 nesting sites in the Southwest Indian Ocean. This is the first time that Atlantic Ocean haplotypes have been recorded among any Indo-Pacific nesting populations. Previous studies indicated that the Cape of Good Hope was a major biogeographical barrier between the Atlantic and Indian Oceans because evidence for gene flow in the last 1.5 million years has yet to emerge. This study, by sampling localities adjacent to this barrier, demonstrates that recent gene flow has occurred from the Atlantic Ocean into the Indian Ocean via the Cape of Good Hope. We also found compelling genetic evidence that green turtles nesting at the rookeries of the South Mozambique Channel (SMC) and those nesting in the North Mozambique Channel (NMC) belong to separate genetic stocks. Furthermore, the SMC could be subdivided in two different genetic stocks, one in Europa and the other one in Juan de Nova. We suggest that this particular genetic pattern along the Mozambique Channel is attributable to a recent colonization from the Atlantic Ocean and is maintained by oceanic conditions in the northern and southern Mozambique Channel that influence early stages in the green turtle life cycle.

Reduced female gene flow in the European flat oyster Ostrea edulis

The geographical structure of 15 natural populations of the flat oyster (Ostrea edulis L.) was assessed by single-strand conformation polymorphism (SSCP) of a 313-base-pair (bp) fragment of the mitochondrial 12S-rRNA gene. Fourteen haplotypes were observed, with one being dominant in the Mediterranean samples and another one in the Atlantic populations. The geographically extreme populations sampled in Norway and the Black Sea appeared differentiated by exhibiting the dominance of a third group of haplotypes. The results were compared to available microsatellite data at five loci. The Atlantic/Mediterranean differentiation pattern was qualitatively the same with both types of markers, confirming an isolation-by-distance pattern. The average mitochondrial haplotypic diversity displayed a high among populations variance, reflecting small effective population size in some locations. Additionally, a 10-fold quantitative difference was observed in Fst between the mitochondrial and the nuclear genomes, which could be due to an unbalanced sex ratio or sex-biased differential reproductive success between males and females (or both).

Trans-atlantic distribution of a mangrove oyster species revealed by 16S mtDNA and karyological analyses

Three species of mangrove oysters, Crassostrea rhizophorae, C. brasiliana, and C. gasar, have been described along the Atlantic shores of South America and Africa. Because the distribution of these molluscs is of great biological and commercial interest, their taxonomy and distribution deserve further clarification. Therefore, 15 populations were sampled from both continents. Their 16S mitochondrial polymorphism was studied by sequencing and PCR-RFLP analysis. Two haplotypes were identified. Haplotype a was the only one observed in Africa, but it was also observed in South America together with haplotype b. Because C. gasar is the only mangrove oyster identified on the west coast of Africa, haplotype a was attributed to this species, which has thus been shown to occur in South America. Haplotype b is attributed to C. rhizophorae. The karyotypes of specimens of C. gasar, from Africa and from South America, were very similar, and both species were observed at the same location in Brazil. The occurrence of C. gasar in South America adds a third species-in addition to C. rhizophorae and C. brasiliana-to the list of species present along these coasts. The predominant surface circulation patterns in this part of the Atlantic Ocean favor the hypothesis that C. gasar was transported from Africa to America. Finally, a phylogenetic tree built with seven 16S sequences from Crassostrea and Saccostrea species showed that C. gasar is intermediate between the American Crassostrea species (C. virginica and C. rhizophorae) and the Asian species (C. gigas and C. ariakensis).