Genome and Phylogenetic Analysis of Genes Involved in the Immune System of Solea senegalensis - Potential Applications in Aquaculture

Genomic Characterization of hox Genes in Senegalese Sole (Solea senegalensis, Kaup 1858): Clues to Evolutionary Path in Pleuronectiformes

The Senegalese sole (Solea senegalensis, Kaup 1858), a marine flatfish, belongs to the Pleuronectiformes order. It is a commercially important species for fisheries and aquaculture. However, in aquaculture, several production bottlenecks have still to be resolved, including skeletal deformities and high mortality during the larval and juvenile phase. The study aims to characterize the hox gene clusters in S. senegalensis to understand better the developmental and metamorphosis process in this species. Using a BAC library, the clones that contain hox genes were isolated, sequenced by NGS and used as BAC-FISH probes. Subsequently the hox clusters were studied by sequence analysis, comparative genomics, and cytogenetic and phylogenetic analysis. Cytogenetic analysis demonstrated the localization of four BAC clones on chromosome pairs 4, 12, 13, and 16 of the Senegalese sole cytogenomic map. Comparative and phylogenetic analysis showed a highly conserved organization in each cluster and different phylogenetic clustering in each hox cluster. Analysis of structural and repetitive sequences revealed accumulations of polymorphisms mediated by repetitive elements in the hoxba cluster, mainly retroelements. Therefore, a possible loss of the hoxb7a gene can be established in the Pleuronectiformes lineage. This work allows the organization and regulation of hox clusters to be understood, and is a good base for further studies of expression patterns.

Characterization of TRIM16, a member of the fish-specific finTRIM family, in olive flounder Paralichthys olivaceus

Tripartite motif-containing (TRIM) proteins are conserved throughout the metazoan kingdom, and the TRIM subset finTRIM is highly diversified in fish. We isolated TRIM16 cDNA, a member of the finTRIM family, from the olive flounder Paralichthys olivaceus (PoTRIM16). PoTRIM16 contained a 1,725-bp coding sequence encoding a 574-amino acid polypeptide, which in turn contained a really interesting new gene (RING) finger domain, B-box-type zinc finger (B-BOX), nuclease SbcCD subunit C (SbcC), structural maintenance of chromosome (SMC prok B), and stonustoxin (SNTX) subunit alpha (SPRY-PRY-SNTX). Multiple alignment of related sequences revealed that PoTRIM16 showed 86.63-97.40% identity with fish orthologues, and a phylogenetic tree was constructed of vertebrates. PoTRIM16 mRNA was detected in all tissues examined; levels were highest in the eye and ovary. PoTRIM16 mRNA expression was investigated during early development. Under VHSV infection, PoTRIM16 mRNA was downregulated in the liver of P. olivaceus. This is the first study to characterize fish-specific finTRIM in P. olivaceus, which may play a role in the immune response against virus infection.

Integration of Maps Enables a Cytogenomics Analysis of the Complete Karyotype in Solea senegalensis

The Pleuronectiformes order, which includes several commercially-important species, has undergone extensive chromosome evolution. One of these species is Solea senegalensis, a flatfish with 2n = 42 chromosomes. In this study, a cytogenomics approach and integration with previous maps was applied to characterize the karyotype of the species. Synteny analysis of S. senegalensis was carried out using two flatfish as a reference: Cynoglossus semilaevis and Scophthalmus maximus. Most S. senegalensis chromosomes (or chromosome arms for metacentrics and submetacentrics) showed a one-to-one macrosyntenic pattern with the other two species. In addition, we studied how repetitive sequences could have played a role in the evolution of S. senegalensis bi-armed (3, and 5–9) and acrocentric (11, 12 and 16) chromosomes, which showed the highest rearrangements compared with the reference species. A higher abundance of TEs (Transposable Elements) and other repeated elements was observed adjacent to telomeric regions on chromosomes 3, 7, 9 and 16. However, on chromosome 11, a greater abundance of DNA transposons was detected in interstitial BACs. This chromosome is syntenic with several chromosomes of the other two flatfish species, suggesting rearrangements during its evolution. A similar situation was also found on chromosome 16 (for microsatellites and low complexity sequences), but not for TEs (retroelements and DNA transposons). These differences in the distribution and abundance of repetitive elements in chromosomes that have undergone remodeling processes during the course of evolution also suggest a possible role for simple repeat sequences in rearranged regions.

Toll-Like Receptor Evolution: Does Temperature Matter?

Toll-like receptors (TLRs) recognize conserved pathogen-associated molecular patterns (PAMPs) and are an ancient and well-conserved group of pattern recognition receptors (PRRs). The isolation of the Antarctic continent and its unique teleost fish and microbiota prompted the present investigation into Tlr evolution. Gene homologues of tlr members in teleosts from temperate regions were present in the genome of Antarctic Nototheniidae and the non-Antarctic sister lineage Bovichtidae. Overall, in Nototheniidae apart from D. mawsoni, no major tlr gene family expansion or contraction occurred. Instead, lineage and species-specific changes in the ectodomain and LRR of Tlrs occurred, particularly in the Tlr11 superfamily that is well represented in fish. Positive selective pressure and associated sequence modifications in the TLR ectodomain and within the leucine-rich repeats (LRR), important for pathogen recognition, occurred in Tlr5, Tlr8, Tlr13, Tlr21, Tlr22, and Tlr23 presumably associated with the unique Antarctic microbiota. Exposure to lipopolysaccharide (Escherichia coli O111:B4) Gram negative bacteria did not modify tlr gene expression in N. rossii head–kidney or anterior intestine, although increased water temperature (+4°C) had a significant effect.

Genome-wide in silico identification of phospholipase D (PLD) gene family from Corchorus capsularis and Corchorus olitorius: reveals their responses to plant stress

Background

Plant grows in nature facing various types of abiotic stresses for their normal growth and development. During abiotic stress, plants evolve different types of mechanisms to survive in a hostile environment. Phospholipase D (PLD) plays important role in the regulation of diverse cellular processes including stress responses in plants. Member of PLD genes are well studied in different model plants; however, their functions in the jute are not clear yet.

Result

In the present study, a total of 12 and 11 PLD genes were identified in the genome of C. capsularis and C. olitorius, respectively. The presence of the two conserved HKD motifs in PLD genes except for CoPLDδ-2 in jute suggests their strong lipase activity. Twenty different motifs were found in the identified PLD genes, and PLD-β1, PLD-γ1, and all members of PLD-δ1 of both jute species contained the highest number of motifs. Phylogenetic analysis showed the close evolutionary relationship among the five groups of jute PLD proteins along with the PLD proteins from Arabidopsis. Tissue-specific expression pattern of PLDα1-2, PLD-α2, PLDβ1, PLDγ1, and PLDδ1 of two jute species suggested their involvement in plant growth and development. However, the expression pattern of PLDα1-2, PLDα1-3, PLD-α4, PLDδ1, and PLDδ3 indicated their association during waterlogging stress. In addition, PLD-α2, PLDβ1, and PLDδ2 seemed to be involved in drought stress as well as salinity stress.

Conclusion

This genome-wide identification of jute PLD genes from C. capsularis and C. olitorius will help to further functional characterization of the PLD genes for developing stress-tolerant jute variety.

Cytogenomics Unveil Possible Transposable Elements Driving Rearrangements in Chromosomes 2 and 4 of Solea senegalensis

Cytogenomics, the integration of cytogenetic and genomic data, has been used here to reconstruct the evolution of chromosomes 2 and 4 of Solea senegalensis. S. senegalensis is a flat fish with a karyotype comprising 2n = 42 chromosomes: 6 metacentric + 4 submetacentric + 8 subtelocentric + 24 telocentric. The Fluorescence in situ Hybridization with Bacterial Artificial Chromosomes (FISH-BAC) technique was applied to locate BACs in these chromosomes (11 and 10 BACs in chromosomes 2 and 4, respectively) and to generate integrated maps. Synteny analysis, taking eight reference fish species (Cynoglossus semilaevis, Scophthalmus maximus, Sparus aurata, Gasterosteus aculeatus, Xiphophorus maculatus, Oryzias latipes, Danio rerio, and Lepisosteus oculatus) for comparison, showed that the BACs of these two chromosomes of S. senegalensis were mainly distributed in two principal chromosomes in the reference species. Transposable Elements (TE) analysis showed significant differences between the two chromosomes, in terms of number of loci per Mb and coverage, and the class of TE (I or II) present. Analysis of TE divergence in chromosomes 2 and 4 compared to their syntenic regions in four reference fish species (C. semilaevis, S. maximus, O. latipes, and D. rerio) revealed differences in their age of activity compared with those species but less notable differences between the two chromosomes. Differences were also observed in peaks of divergence and coverage of TE families for all reference species even in those close to S. senegalensis, like S. maximus and C. semilaevis. Considered together, chromosomes 2 and 4 have evolved by Robertsonian fusions, pericentric inversions, and other chromosomal rearrangements mediated by TEs.

A Comprehensive Integrated Genetic Map of the Complete Karyotype of Solea senegalensis (Kaup 1858)

Solea senegalensis aquaculture production has experienced a great increase in the last decade and, consequently, the genome knowledge of the species is gaining attention. In this sense, obtaining a high-density genome mapping of the species could offer clues to the aquaculture improvement in those aspects not resolved so far. In the present article, a review and new processed data have allowed to obtain a high-density BAC-based cytogenetic map of S. senegalensis beside the analysis of the sequences of such BAC clones to achieve integrative data. A total of 93 BAC clones were used to localize the chromosome complement of the species and 588 genes were annotated, thus almost reaching the 2.5% of the S. senegalensis genome sequences. As a result, important data about its genome organization and evolution were obtained, such as the lesser gene density of the large metacentric pair compared with the other metacentric chromosomes, which supports the theory of a sex proto-chromosome pair. In addition, chromosomes with a high number of linked genes that are conserved, even in distant species, were detected. This kind of result widens the knowledge of this species’ chromosome dynamics and evolution.

Development of whole-genome multiplex assays and construction of an integrated genetic map using SSR markers in Senegalese sole

The Senegalese sole (Solea senegalensis) is an economically important flatfish species. In this study, a genome draft was analyzed to identify microsatellite (SSR) markers for whole-genome genotyping. A subset of 224 contigs containing SSRs were preselected and validated by using a de novo female hybrid assembly. Overall, the SSR density in the genome was 886.7 markers per megabase of genomic sequences and the dinucleotide motif was the most abundant (52.4%). In silico comparison identified a set of 108 SSRs (with di-, tetra- or pentanucleotide motifs) widely distributed in the genome and suitable for primer design. A total of 106 markers were structured in thirteen multiplex PCR assays (with up to 10-plex) and the amplification conditions were optimized with a high-quality score. Main genetic diversity statistics and genotyping reliability were assessed. A subset of 40 high polymorphic markers were selected to optimize four supermultiplex PCRs (with up to 11-plex) for pedigree analysis. Theoretical exclusion probabilities and real parentage allocation tests using parent–offspring information confirmed their robustness and effectiveness for parental assignment. These new SSR markers were combined with previously published SSRs (in total 229 makers) to construct a new and improved integrated genetic map containing 21 linkage groups that matched with the expected number of chromosomes. Synteny analysis with respect to C. semilaevis provided new clues on chromosome evolution in flatfish and the formation of metacentric and submetacentric chromosomes in Senegalese sole.

Gene clusters related to metamorphosis in Solea senegalensis are highly conserved

The flatfish, Solea senegalensis has considerable scientific interest and commercial value. The metamorphosis in this species occurs between 12 and 19 days after hatching and it takes about 1 week to complete. Eleven Bacterial Artificial Chromosomes (BAC) clones containing the various candidate genes involved in the process of metamorphosis: thyroxine 5 deiodinase 3 (dio3); forkhead box protein E4 (foxe4); melatonin receptor type 1C (mel1c); calsequestrin 1b (casq1b); thyrotropin subunit beta (tshβ); thyrotropin-releasing hormone receptor 1, 2, and 3 (trhr1, trhr2, trhr3); thyroid hormone receptor α a and b (thrαa, thrαb); and thyroid hormone receptor beta (thrβ) were analyzed by multiple Fluorescence in situ Hybridization (mFISH) and Next Generation Sequencing (NGS) techniques. The mFISH technique localized the 11 BAC clones on 12 different chromosome pairs because three of them, specifically the trhr1a, trhr2 and thrβ BAC clones, showed double signals. This signal duplication indicates a duplication of the genomic region inserted within the BAC clone, which provides evidence for the Teleost-Specific Whole Genome Duplication (TS-WGD). Micro-synteny and phylogenetic analysis showed that Cynoglossus semilaevis is the nearest species to S. senegalensis and that Danio rerio is the most distant one. The tshβ BAC clone was highly conserved as the genes belonging to this BAC were located on a single chromosome in all the species studied. These genes participate in proliferation, migration and cell-death, which are key processes during metamorphosis. Overall, micro-synteny analysis showed that most candidate genes are found in conserved genomic surroundings.