Influence of Storage Time on the DNA Integrity and Viability of Spermatozoa of the Spider Crab Maja brachydactyla.

Natural populations of the spider crab Maja brachydactyla constitute a fishery resource of great economic importance in many countries. As in the rest of eubrachyurans, the females of this species have ventral-type seminal receptacles where they store sperm from copulations. Sperm can be stored in these structures for months and even years before egg fertilisation, with the consequent degradation of the sperm cells during the time. In this work, we analyse the viability and the possible genetic damage in sperm accumulated in the seminal receptacles of M. brachydactyla females as a function of the storage time (from 0 to 14 months) using the comet assay technique. On one hand, we developed an algorithm for comet image analysis that improves the comet segmentation compared with the free software Open comet v1.3.1 (97% vs. 76% of detection). In addition, our software allows the manual modification of the contours wrongly delimited via the automatic tool. On the other hand, our data show a sharp decline in sperm viability and DNA integrity in the first four months of storage, which could lead to a decrease in the fecundity rate and/or viability of the embryos or larvae from the second and third clutches of the annual cycle if the repair capacity in these gametic cells is low.

Chromosome evolution in Iberolacerta, a genus that deviates from the standard karyotype formula of Lacertidae.

This paper describes the preparation of flow-sorted chromosome paints from the Iberian Rock lizard Iberolacerta monticola, exemplifying their subsequent use in cross-species comparisons of chromosome painting. We carried out comparative analyses of chromosome evolution in the congeneric species I. galani and I. bonnali, as well as in two other species of Lacertini (Lacerta schreiberi and Timon lepidus) whose sex chromosomes were also studied through comparative genomic hybridization. Most species of Lacertini possess a diplod number of 2n = 38, with 36 acrocentric macrochromosomes and 2 microchromosomes. However, the nine species included in the genus Iberolacerta do not possess microchromosomes. Furthermore, very conspicuous differences from the standard Lacertini karyotype were observed in the three Pyrenean species of this genus, which included several biarmed metacentrics and a Z1Z2W multiple sex-chromosome system. With the possible exception of L. schreiberi, all the species of the family Lacertidae described to date appear to share homologous Z chromosomes, which date back to the last common ancestor of the whole group. We provide conclusive evidence that L. schreiberi should no longer be considered an exception to this rule, and demonstrate that the loss of microchromosomes in Iberolacerta was produced by their fusion to a middle-sized chromosome. Furthermore, we show that the multiple sex-chromosome system of the Pyrenean species of Iberolacerta originated from the fusion of the ancestral W chromosome with one of the shortest autosomes, and provide additional evidence of the fast evolution of DNA sequences linked to the W chromosome in Lacertini.

Mitogenome phylogenetics in the genus Palaemon (Crustacea: Decapoda) sheds light on species crypticism in the rockpool shrimp P. elegans.

The genus Palaemon comprises worldwide marine and freshwater shrimps and prawns, and some of them are ecologically or commercially important species. Palaemon is not currently a monophyletic group, so phylogenetics and systematics are constantly changing. Species crypticism has been pointed out in several Palaemon species, being the clearest evidence in the European rockpool shrimp P. elegans. Here we sequenced and described seven European Palaemon mitochondrial genomes. The mitochondrial protein-coding genes were used, along with those of three other Palaemon species, to perform mitogenome phylogenetic analyses to clarify the evolutionary relationships within the genus, and particularly to shed light on the cryptic species found within P. elegans. The Messinian Salinity Crisis (5.3-5.9 Ma, late Miocene) was proposed to be the origin of this cryptic species and it was used as aged constraint for calibration analysis. We provide the largest and the first time-calibrated mitogenome phylogeny of the genus Palaemon and mitogenome substitution rate was estimated (1.59% per million years) in Decapoda for the first time. Our results highlighted the need for future systematics changes in Palaemon and crypticism in P. elegans was confirmed. Mitochondrial genome and cox1 (1.41%) substitution rate estimates matched those published elsewhere, arguing that the Messinian Salinity Crisis was a plausible event driving the split between P. elegans and its cryptic species. Molecular dating suggested that Pleistocene glaciations were likely involved in the differentiation between the Atlantic and Mediterranean populations of P. elegans. On the contrary, the divergence between the Atlantic and Mediterranean populations of the common littoral shrimp P. serratus was greater and dated to be much older (4.5-12.3 Ma, Plio-Miocene), so we considered that they could represent two separated species. Therefore, species crypticism in the genus Palaemon seems to be a common phenomenon.

Genetic diversity and population structure of the rockpool shrimp Palaemon elegans based on microsatellites: evidence for a cryptic species and differentiation across the Atlantic-Mediterranean transition.

The rockpool shrimp Palaemon elegans is an ecologically important crustacean species within the European coastline fauna. In the present study, genetic diversity and population structure and connectivity were assessed by examining 21 polymorphic microsatellite loci at 13 sampling sites located along the natural distribution range. All localities showed similar levels of genetic variability. Significant deficits of heterozygosity were recorded, most likely due to the presence of null alleles. Genetic structure analyses revealed two clearly genetically distinct groups within P. elegans but without following any geographical or oceanographic basis. Thus, our results provided nuclear evidence for the existence of a Mediterranean cryptic species within P. elegans, highlighting the need to revise its taxonomic status. Regarding P. elegans sensu stricto, population structuring was reported across the Atlantic-Mediterranean transition area, where the Almería-Orán Front restricts the gene flow between the Atlantic and the Mediterranean population. Moreover, while population connectivity was suggested between all Mediterranean localities, some substructure was found within the Atlantic group. Canary Islands exhibited a weak but significant genetic differentiation from all Atlantic mainland localities, consistent with the isolation-by-distance pattern detected throughout the Atlantic population. Overall, all these findings provided new insights into the population biology of P. elegans complex.

High incidence of heteroplasmy in the mtDNA of a natural population of the spider crab Maja brachydactyla.

Mitochondria are mostly inherited by maternal via, that is, only mitochondria from eggs are retained in the embryos. However, this general assumption of uniparentally transmitted, homoplasmic and non-recombining mitochondrial genomes is becoming more and more controversial. The presence of different sequences of mtDNA within a cell or individual, known as heteroplasmy, is increasingly reported in several taxon of animals, such as molluscs, arthropods and vertebrates. In this work, a considerable frequency of heteroplasmy were detected in the COI and 16S genes of the spider crab Maja brachydactyla, possibly associated to hybridisation with the congeneric species Maja squinado. This finding is a fact to keep in mind before addressing molecular analyses based on mitochondrial markers, since the assumption of maternal inheritance could lead to erroneous results. As M. brachydactyla is a commercial species, heteroplasmy is an important aspect to take into account for the fisheries management of this resource, since effective population size could be overestimated.

De novo gonad transcriptome analysis of the common littoral shrimp Palaemon serratus: novel insights into sex-related genes.

BACKGROUND: The common littoral shrimp Palaemon serratus is an economically important decapod resource in some European communities. Aquaculture practices prevent the genetic deterioration of wild stocks caused by overfishing and at the same time enhance the production. The biotechnological manipulation of sex-related genes has the proved potential to improve the aquaculture production but the scarcity of genomic data about P. serratus hinders these applications. RNA-Seq analysis has been performed on ovary and testis samples to generate a reference gonadal transcriptome. Differential expression analyses were conducted between three ovary and three testis samples sequenced by Illumina HiSeq 4000 PE100 to reveal sex-related genes with sex-biased or sex-specific expression patterns. RESULTS: A total of 224.5 and 281.1 million paired-end reads were produced from ovary and testis samples, respectively. De novo assembly of ovary and testis trimmed reads yielded a transcriptome with 39,186 transcripts. The 29.57% of the transcriptome retrieved at least one annotation and 11,087 differentially expressed genes (DEGs) were detected between ovary and testis replicates. Six thousand two hundred seven genes were up-regulated in ovaries meanwhile 4880 genes were up-regulated in testes. Candidate genes to be involved in sexual development and gonadal development processes were retrieved from the transcriptome. These sex-related genes were discussed taking into account whether they were up-regulated in ovary, up-regulated in testis or not differentially expressed between gonads and in the framework of previous findings in other crustacean species. CONCLUSIONS: This is the first transcriptome analysis of P. serratus gonads using RNA-Seq technology. Interesting findings about sex-related genes from an evolutionary perspective (such as Dmrt1) and for putative future aquaculture applications (Iag or vitellogenesis genes) are reported here. We provide a valuable dataset that will facilitate further research into the reproductive biology of this shrimp.

Isolation and characterization of 21 polymorphic microsatellite loci for the rockpool shrimp Palaemon elegans using Illumina MiSeq sequencing.

The rockpool shrimp Palaemon elegans is considered an important crustacean species within the European coastline fauna. This species is experiencing an ongoing geographical expansion beyond its native distribution range due to unintentional human introductions. A better knowledge of the genetic diversity, geographic structure and connectivity of its populations is necessary. In the present study, microsatellite loci were isolated using the Illumina MiSeq platform. The microsatellite-enriched library sequencing produced 3.9 million raw reads. Reads were processed and primer pairs were designed for microsatellite sequences amplification. Ninety-six microsatellite loci were preliminary screened in individuals from Atlantic and Mediterranean localities. From them, 21 loci exhibited reliable polymorphism and were thoroughly characterized in 30 individuals from a Cantabrian locality (Spain). No linkage disequilibrium between pairs of loci was detected. Number of alleles per locus ranged from 2 to 12. Observed and expected heterozygosities ranged from 0.033 to 0.833 and from 0.033 to 0.869 respectively. No significant departure from the Hardy-Weinberg equilibrium was detected in most of loci. This is the first time that microsatellite markers have been developed for P. elegans. This characterized microsatellite suite provides new suitable tools for further analyses, facilitating the understanding of population genetics both in natural and introduced populations.

Comparative cytogenetic analysis of marine Palaemon species reveals a X1X1X2X2/X1X2Y sex chromosome system in Palaemon elegans.

BACKGROUND: The maintenance of species and the promotion of speciation are closely related to chromosomal rearrangements throughout evolution. Decapoda represents the most species-rich order among crustaceans and, despite its ecological and economic importance, little is known about decapod karyology. We aim at cytogenetically characterizing two sympatric prawn species. RESULTS: Analysis of mitotic metaphases and meiotic diakinesis of the common prawn Palaemon serratus and the rockpool prawn P. elegans, revealed considerable differences between their karyotypes including chromosome numbers and sex determination systems. The cytogenetic data for P. serratus showed a diploid number of 56 and the putative absence of heteromorphic sex chromosomes. However, the diploid chromosome number in P. elegans was 90 for females and 89 for males. The karyotype of the females consisted of the three largest acrocentric pairs and 42 submetacentric and metacentric pairs, while the karyotype of the males comprised a clearly identifiable large metacentric chromosome and two acrocentric pairs as well as the smaller 42 pairs. These results highlight the presence of the X1X1X2X2/X1X2Y multiple sex chromosome system in P. elegans, which constitute the only sexual system for Decapoda reported cytogenetically using modern techniques. The origin of this sex chromosome system is discussed. We hypothesize that the chromosome evolution within the genus could involve several fusion events giving rise to a reduction on the chromosome number in P. serratus. In both species, the major ribosomal genes were located in two chromosome pairs and hybridization signals of the telomeric sequences (TTAGGG)n were visualized at the telomeres of all chromosomes. C-banding revealed that, when present, constitutive heterochromatin had a predominantly telomeric distribution and no centromeric constitutive heterochromatin was observed. CONCLUSIONS: Although more comparative cytogenetic analyses are needed to clarify our hypotheses, the findings of this work indicate that the prawns of the genus Palaemon represent a promising model among Decapoda representatives to investigate the karyotype evolution and the patterns of sex chromosome differentiation.

De novo transcriptome assembly of shrimp Palaemon serratus.

The shrimp Palaemon serratus is a coastal decapod crustacean with a high commercial value. It is harvested for human consumption. In this study, we used Illumina sequencing technology (HiSeq 2000) to sequence, assemble and annotate the transcriptome of P. serratus. RNA was isolated from muscle of adults individuals and, from a pool of larvae. A total number of 4 cDNA libraries were constructed, using the TruSeq RNA Sample Preparation Kit v2. The raw data in this study was deposited in NCBI SRA database with study accession number of SRP090769. The obtained data were subjected to de novo transcriptome assembly using Trinity software, and coding regions were predicted by TransDecoder. We used Blastp and Sma3s to annotate the identified proteins. The transcriptome data could provide some insight into the understanding of genes involved in the larval development and metamorphosis. SPECIFICATIONS: [Table: see text].

Evolutionary dynamics of two satellite DNA families in rock lizards of the genus Iberolacerta (Squamata, Lacertidae): different histories but common traits.

Satellite DNAs compose a large portion of all higher eukaryotic genomes. The turnover of these highly repetitive sequences is an important element in genome organization and evolution. However, information about the structure and dynamics of reptilian satellite DNA is still scarce. Two satellite DNA families, HindIII and TaqI, have been previously characterized in four species of the genus Iberolacerta. These families showed different chromosomal locations, abundances, and evolutionary rates. Here, we extend the study of both satellite DNAs (satDNAs) to the remaining Iberolacerta species, with the aim to investigate the patterns of variability and factors influencing the evolution of these repetitive sequences. Our results revealed disparate patterns but also common traits in the evolutionary histories of these satellite families: (i) each satellite DNA is made up of a library of monomer variants or subfamilies shared by related species; (ii) species-specific profiles of satellite repeats are shaped by expansions and/or contractions of different variants from the library; (iii) different turnover rates, even among closely related species, result in great differences in overall sequence homogeneity and in concerted or non-concerted evolution patterns, which may not reflect the phylogenetic relationships among taxa. Contrasting turnover rates are possibly related to genomic constraints such as karyotype architecture and the interspersed organization of diverging repeat variants in satellite arrays. Moreover, rapid changes in copy number, especially in the centromeric HindIII satDNA, may have been associated with chromosomal rearrangements and even contributed to speciation within Iberolacerta.

Isolation and characterization of two satellite DNAs in some Iberian rock lizards (Squamata, Lacertidae).

Satellite DNAs represent a large portion of all high eukaryotic genomes. They consist of numerous very similar repeated sequences, tandemly arranged in large clusters up to 100 million base pairs in length, usually located in the heterochromatic parts of chromosomes. The biological significance of satDNAs is still under discussion, but most of their proposed functions are related to heterochromatin and/or centromere formation and function. Because information about the structure of reptilian satDNA is far from exhaustive, we present a molecular and cytogenetic characterization of two satDNA families in four lacertid species. Two families of tandemly repeated DNAs, namely TaqI and HindIII satDNAs, have been cloned and sequenced from four species belonging to the genus Iberolacerta. These satDNAs are characterized by a monomer length of 171-188 and 170-172 bp, and by an AT content of 60.5% and 58.1%, respectively. FISH experiments with TaqI satDNA probe produced bright signals in pericentromeric regions of a subset of chromosomes whereas all the centromeres were marked by HindIII probe. The results obtained in this study suggest that chromosome location and abundance of satDNAs influence the evolution of these elements, with centromeric families evolving tenfold faster than interstitial/pericentromeric ones. Such different rates render different satellites useful for phylogenetic investigation at different taxonomic ranks.

Molecular organization and phylogenetic analysis of 5S rDNA in crustaceans of the genus Pollicipes reveal birth-and-death evolution and strong purifying selection.

BACKGROUND: The 5S ribosomal DNA (5S rDNA) is organized in tandem arrays with repeat units that consist of a transcribing region (5S) and a variable nontranscribed spacer (NTS), in higher eukaryotes. Until recently the 5S rDNA was thought to be subject to concerted evolution, however, in several taxa, sequence divergence levels between the 5S and the NTS were found higher than expected under this model. So, many studies have shown that birth-and-death processes and selection can drive the evolution of 5S rDNA. In analyses of 5S rDNA evolution is found several 5S rDNA types in the genome, with low levels of nucleotide variation in the 5S and a spacer region highly divergent. Molecular organization and nucleotide sequence of the 5S ribosomal DNA multigene family (5S rDNA) were investigated in three Pollicipes species in an evolutionary context. RESULTS: The nucleotide sequence variation revealed that several 5S rDNA variants occur in Pollicipes genomes. They are clustered in up to seven different types based on differences in their nontranscribed spacers (NTS). Five different units of 5S rDNA were characterized in P. pollicipes and two different units in P. elegans and P. polymerus. Analysis of these sequences showed that identical types were shared among species and that two pseudogenes were present. We predicted the secondary structure and characterized the upstream and downstream conserved elements. Phylogenetic analysis showed an among-species clustering pattern of 5S rDNA types. CONCLUSIONS: These results suggest that the evolution of Pollicipes 5S rDNA is driven by birth-and-death processes with strong purifying selection.

The 5S rDNA gene family in mollusks: characterization of transcriptional regulatory regions, prediction of secondary structures, and long-term evolution, with special attention to Mytilidae mussels.

Several reports on the characterization of 5S ribosomal DNA (5S rDNA) in various animal groups have been published to date, but there is a lack of studies analyzing this gene family in a much broader context. Here, we have studied 5S rDNA variation in several molluskan species, including bivalves, gastropods, and cephalopods. The degree of conservation of transcriptional regulatory regions was analyzed in these lineages, revealing a conserved TATA-like box in the upstream region. The evolution of the 120 bp coding region (5S) was also studied, suggesting the occurrence of paralogue groups in razor clams, clams, and cockles. In addition, 5S rDNA sequences from 11 species and 7 genus of Mytilidae Rafinesque, 1815 mussels were sampled and studied in detail. Four different 5S rDNA types, based on the nontranscribed spacer region were identified. The phylogenetic analyses performed within each type showed a between-species gene clustering pattern, suggesting ancestral polymorphism. Moreover, some putative pseudogenized 5S copies were also identified. Our report, together with previous studies that found high degree of intragenomic divergence in bivalve species, suggests that birth-and-death evolution may be the main force driving the evolution of 5S rDNA in these animals, even at the genus level.

Analysis of ITS1 and ITS2 sequences in Ensis razor shells: suitability as molecular markers at the population and species levels, and evolution of these ribosomal DNA spacers.

Internal transcribed spacer 1 and 2 (ITS1 and ITS2) sequences were analysed in Ensis razor shells (Mollusca: Bivalvia: Pharidae). We aimed to (1) test ITS1 and ITS2 as molecular markers at the population level in the successful alien E. directus (Conrad, 1843); (2) test these spacers at the species level in E. directus and three other Ensis species, E. siliqua (L., 1758), E. macha (Molina, 1782), and E. magnus (Schumacher, 1817); and (3) analyse the evolutionary processes that may be shaping Ensis ITS1 and ITS2 extant variation. In E. directus, despite the intragenomic divergence detected, ITS1 and ITS2 were informative in differentiating the geographic areas considered (Denmark and Canada) by means of both the insertion-deletion polymorphism and the nucleotide polymorphism. In this species, the 5.8S ribosomal gene (5.8S) showed scarce polymorphism. At the species level, maximum parsimony and maximum likelihood analyses revealed that ITS1 and ITS2 may be suitable to reconstruct Ensis phylogenetic relationships. Finally, the evolutionary models that best fit the long-term evolution of Ensis ITS1-5.8S-ITS2 are discussed. A mixed process of concerted evolution, birth-and-death evolution, and selection is chosen as an option that may reconcile the long-term evolution of Ensis ITS1-5.8S-ITS2 and 5S ribosomal DNA.

Long-term evolution of 5S ribosomal DNA seems to be driven by birth-and-death processes and selection in Ensis razor shells (Mollusca: Bivalvia).

A study of nucleotide sequence variation of 5S ribosomal DNA from six Ensis species revealed that several 5S ribosomal DNA variants, based on differences in their nontranscribed spacers (NTS), occur in Ensis genomes. The 5S rRNA gene was not very polymorphic, compared with the NTS region. The phylogenetic analyses performed showed a between-species clustering of 5S ribosomal DNA variants. Sequence divergence levels between variants were very large, revealing a lack of sequence homogenization. These results strongly suggest that the long-term evolution of Ensis 5S ribosomal DNA is driven by birth-and-death processes and selection.

Heterogeneous nature and distribution of interruptions in dinucleotides may indicate the existence of biased substitutions underlying microsatellite evolution.

Some aspects of microsatellite evolution, such as the role of base substitutions, are far from being fully understood. To examine the significance of base substitutions underlying the evolution of microsatellites we explored the nature and the distribution of interruptions in dinucleotide repeats from the human genome. The frequencies that we inferred in the repetitive sequences were statistically different from the frequencies observed in other noncoding sequences. Additionally, we detected that the interruptions tended to be towards the ends of the microsatellites and 5'-3' asymmetry. In all the estimates nucleotides forming the same repetitive motif seem to be affected by different base substitution rates in AC and AG. This tendency itself could generate patterning and similarity in flanking sequences and reconcile these phenomena with the high mutation rate found in flanking sequences without invoking convergent evolution. Nevertheless, our data suggest that there is a regional bias in the substitution pattern of microsatellites. The accumulation of random substitutions alone cannot explain the heterogeneity and the asymmetry of interruptions found in this study or the relative frequency of different compound microsatellites in the human genome. Therefore, we cannot rule out the possibility of a mutational bias leading to convergent or parallel evolution in flanking sequences.

Genetic divergence detected by ISSR markers and characterization of microsatellite regions in Mytilus mussels.

The wide distribution of microsatellites in mussels of the Mytilus edulis complex (Mytilidae) enables the analysis of inter-simple-sequence repeat (ISSR) markers. The aim of this investigation was to assess genetic differentiation in six sampling localities distributed along the European Atlantic coast to expose the potential of these markers in genetic studies requiring the detection of low polymorphism and as a source of sequences for developing microsatellite markers. We detected low genetic structuring within each member of the Mytilus edulis complex. Nei and Li distances and AMOVA clustered the individuals studied into two groups. On the basis of these results two sampling localities coming from the M. edulis x M. galloprovincialis mosaic hybrid zone in Western Europe were assigned to one species. On the other hand, mussels of a sampling locality in the Baltic Sea were not significantly different from a pure M. edulis locality supporting an extensive introgression of M. edulis in these individuals. Finally, 148 microsatellites were found in the sequences of 51 ISSR markers, and two polymorphic microsatellite markers were developed.

Telomeric localization of the vertebrate-type hexamer repeat, (TTAGGG)n, in the wedgeshell clam Donax trunculus and other marine invertebrate genomes.

The hexamer repeat sequence (TTAGGG)(n), found at the ends of all vertebrate chromosomes, was previously identified as the main building element of one member of a HindIII satellite DNA family characterized in the genome of the bivalve mollusc Donax trunculus. It was also found in 22 perfect tandem repeats in a cloned junction region juxtaposed to the proper satellite sequence, from which the DNA tract encompassing the clustered tandem copies was excised and subcloned. Here, the chromosomal distribution of (TTAGGG)(n) sequences in the Donax was studied by the sensitivity to Bal31 exonuclease digestion, fluorescence in situ hybridization (FISH) on metaphase chromosomes and rotating-field gel electrophoresis. To verify the occurrence of the hexamer repeat in the genomes of taxonomically related molluscs and other marine invertebrates, genomic DNA from the mussel Mytilus galloprovincialis and the echinoderm Holothuria tubulosa was also analyzed. The kinetics of Bal31 hydrolysis of high molecular mass DNA from the three marine invertebrates revealed a marked decrease over time of the hybridization with the cloned (TTAGGG)(22) sequence, concomitantly with a progressive shortening of the positively reacting DNA fragments. This revealed a marked susceptibility to exonuclease consistent with terminal positioning on the respective chromosomal DNAs. In full agreement, FISH results with the (TTAGGG)(22) probe showed that the repeat appears located in telomeric regions in all chromosomes of both bivalve molluscs. The presence of (TTAGGG)(n) repeat tracts in marine invertebrate telomeres points to its wider distribution among eukaryotic organisms and suggests an ancestry older than originally presumed from its vertebrate distinctiveness.