Growth Performance and Transcriptomic Response of Warm-Acclimated Hybrid Abalone Haliotis rufescens (♀) × H. corrugata (♂)

Synthesis of transcriptomic studies reveals a core response to heat stress in abalone (genus Haliotis)

BACKGROUND

As climate change causes marine heat waves to become more intense and frequent, marine species increasingly suffer from heat stress. This stress can result in reduced growth, disrupted breeding cycles, vulnerability to diseases and pathogens, and increased mortality rates. Abalone (genus Haliotis) are an ecologically significant group of marine gastropods and are among the most highly valued seafood products. However, heat stress events have had devastating impacts on both farmed and wild populations. Members of this genus are among the most susceptible marine species to climate change impacts, with over 40% of all abalone species listed as threatened with extinction. This has motivated researchers to explore the genetics linked to heat stress in abalone. A substantial portion of publicly available studies has employed transcriptomic approaches to investigate abalone genetic response to heat stress. However, to date, no meta-analysis has been conducted to determine the common response to heat stress (i.e. the core response) across the genus. This study uses a standardized bioinformatic pipeline to reanalyze and compare publicly available RNA-seq datasets from different heat stress studies on abalone.

RESULTS

Nine publicly available RNA-seq datasets from nine different heat-stress studies on abalone from seven different abalone species and three hybrids were included in the meta-analysis. We identified a core set of 74 differentially expressed genes (DEGs) in response to heat stress in at least seven out of nine studies. This core set of DEGs mainly included genes associated with alternative splicing, heat shock proteins (HSPs), Ubiquitin-Proteasome System (UPS), and other protein folding and protein processing pathways.

CONCLUSIONS

The detection of a consistent set of genes that respond to heat stress across various studies, despite differences in experimental design (e.g. stress intensity, species studied-geographical distribution, preferred temperature range, etc.), strengthens our proposal that these genes are key elements of the heat stress response in abalone. The identification of the core response to heat stress in abalone lays an important foundation for future research. Ultimately, this study will aid conservation efforts and aquaculture through the identification of resilient populations, genetic-based breeding programs, possible manipulations such as early exposure to stress, gene editing and the use of immunostimulants to enhance thermal tolerance.

Transcriptomic analysis reveals the molecular mechanisms of heterosis in low-temperature tolerance in the hybrids of Argopecten scallops

The F1 hybrid of Argopecten irradians irradians (♀) × A. purpuratus (♂) exhibits significant heterosis in growth performance and mid-parent heterosis in low-temperature tolerance. This study presents a comparative transcriptomic analysis of A. irradians irradians (Ai), A. purpuratus (Ap), and the hybrid A. irradians irradians♀ × A. purpuratus♂ (Aip) to explore the mechanisms underlying low-temperature tolerance heterosis in Aip. A total of 33,376 differentially expressed genes (DEGs) were identified between F1 hybrids and purebreds under cold stress. In Aip, 80.32 % of DEGs exhibited non-additive expression patterns, with over-dominant expression observed in 30.65 % of these genes. Pairwise comparisons among the transcriptomes of Ai, Ap, and Aip revealed 14,959 alternative splicing events, affecting 8169 genes. KEGG pathway analysis indicated substantial enrichment of overlapping genes from common DEGs and non-additively expressed genes (NAGs) in apoptosis, longevity regulation, ABC transporters, and spliceosome pathways. Furthermore, analysis of DEGs, DAGs (Differentially Alternative Splicing genes), and NAGs identified 6 genes undergoing alternative splicing. These pathways and genes may be crucial in Aip's response to low-temperature stress and offer insights for advancing scallop cross-breeding strategies.

Diet-driven transcriptional changes in weaning red abalone (Haliotis rufescens) and its hybrid (H. rufescens [♀] x H. fulgens [♂])

The abalone (Haliotis sp) aquaculture industry is on the rise fueled by its high commercial demand. The use of abalone interspecific hybrids is gaining attention due to their improved growth and tolerance to environmental challenges. However, hybrids may respond differently to dietary inputs compared to their parental species, which could be optimized to achieve maximum growth potential. Accordingly, we previously reported that the postlarvae of a hybrid cross obtained by mating red and green abalone (H. rufescens [♀] x H. fulgens [♂]) presented higher growth and survival rates than red abalone postlarvae when fed with Macrocystis pyrifera (M) and a mixture of Ulva onhoi and Navicula incerta (UN). Here, to gain a deeper insight into the physiological and metabolic responses of both crosses, we report the results of a transcriptome sequencing analysis to unveil the diet-driven physiological changes within each cross. Our findings indicate that the diet induces modulation in the transcriptional machinery related to macromolecules breakdown, especially, the mixed diet UN induced higher expression of lysosomal enzymes, especially in enzymes involved in glycosaminoglycan catabolism, such as arylsulfatases, glucosidases, and alpha-fucosidases, which might indicate a higher ability to break down complex carbohydrates from algae. These transcriptional responses might explain the growth difference between pure and hybrid as well as for abalone fed with M and UN.

Comparative analysis of the growth differences between hybrid Ningdu Yellow chickens and their parentals

Although the local high quality chicken breed in China has excellent flavor, its growth rate is inferior to that of foreign breeds. To improve the growth rate of local chicken breeds, it is crucial to study the mechanism of chicken muscle development. Herein, Ningdu Yellow chicken was used as the research object, and a new hybrid breed (W) was obtained by crossing the G, H and D lines, which combined the excellent physiological characteristics of its parents. Combined analysis of Ningdu Yellow chickens and their parents was carried out. Chickens from 105-day-old lines (W, G, H) were selected, and breast meat and serum were extracted for transcriptome sequencing and metabolome determination to study their growth differences. The live weight, carcass weight, half-eviscerated weight, eviscerated weight, and breast muscle weight of W were significantly higher than those of G and H. Differential expression analysis identified 1700 differentially expressed genes (DEG), and gene ontology and kyoto encyclopedia of genes and genomes (KEGG) analyses identified 33 and 1 pathways related to growth and development and steroid biosynthesis, respectively. Next, pairwise analysis identified 57 KEGG pathways, among which the MAPK signaling, steroid hormone biosynthesis, tight junction, and PPAR signaling pathways were involved in growth and development. Cluster analysis found that genes highly expressed in the W group were associated with regulation of the actin cytoskeleton, riboflavin metabolism, steroid biosynthesis, and glycerophospholipid metabolism. The top 2 clusters obtained by protein-protein interaction analysis were important for the growth and development of chickens. Finally, the metabolomic analysis found key differentially accumulated metabolites (DAM) that might be account for the growth differences. Further integrated analysis identified key DEGs and DAMs that might be responsible for the observed growth differences. This study identified genes governing growth traits in Ningdu Yellow chickens, laying a theoretical foundation for the development of chicken breeding, the utilization of hybrid supporting lines, and promotion of the Chinese chicken industry.

From heterosis to outbreeding depression: genotype-by-environment interaction shifts hybrid fitness in opposite directions

In F1 hybrids, phenotypic values are expected to be near the parental means under additive effects or close to one parent under dominance. However, F1 traits can fall outside the parental range, and outbreeding depression occurs when inferior fitness is observed in hybrids. Another possible outcome is heterosis, a phenomenon that interspecific hybrids or intraspecific crossbred F1s exhibit improved fitness compared to both parental species or strains. As an application of heterosis, hybrids between channel catfish females and blue catfish males are superior in feed conversion efficiency, carcass yield, and harvestability. Over 20 years of hybrid catfish production in experimental settings and farming practices generated abundant phenotypic data, making it an ideal system to investigate heterosis. In this study, we characterized fitness in terms of growth and survival longitudinally, revealing environment-dependent heterosis. In ponds, hybrids outgrow both parents due to an extra rapid growth phase of 2-4 months in year 2. This bimodal growth pattern is unique to F1 hybrids in pond culture environments only. In sharp contrast, the same genetic types cultured in tanks display outbreeding depression, where hybrids perform poorly, while channel catfish demonstrate superiority in growth throughout development. Our findings represent the first example, known to the authors, of opposite fitness shifts in response to environmental changes in interspecific vertebrate hybrids, suggesting a broader fitness landscape for F1 hybrids. Future genomic studies based on this experiment will help understand genome-environment interaction in shaping the F1 progeny fitness in the scenario of environment-dependent heterosis and outbreeding depression.

Metabolic and transcriptional responses demonstrating enhanced thermal tolerance in domesticated abalone

Global warming threatens aquatic systems and organisms. Many studies have focused on the vulnerability and stress responses of aquaculture organisms to future thermal conditions. However, it may be of more practical significance to reveal their acclimation potential and mechanisms. In this study, the physiological, metabolic, and transcriptional responses to long-term temperature acclimation of northern and southern populations of Pacific abalone Haliotis discus hannai, a commercially important gastropod sensitive to environmental changes, were compared. This study conducted two common-garden experiments, including a thermostatic experiment in the lab and an aquaculture experiment on the farm. The abalone population cultured in warmer southern waters was tolerant of ongoing high temperatures, whereas the abalone population originally cultured in cooler northern waters exhibited vulnerability to high temperatures but could enhance its thermal tolerance through the process of natural selection in warmer southern waters. This difference was linked to divergence in the metabolic and transcriptional processes of the two populations. The tolerant population exhibited a greater capacity for carbohydrate and amino acid metabolism regulation and energy redistribution to cope with heat stress. This capacity may have been selected for, and accumulated, over many generations because the tolerant population originated from the intolerant population over two decades ago. This work provides insight into the vulnerability and acclimation potential of abalone to heat stress and discloses the molecular and metabolic traits underlying this phenomenon. Future research on the ability of abalone and other commercial shellfish species to acclimate to global warming should take this potential into account.

Transcriptome Analysis of Reciprocal Hybrids Between Crassostrea gigas and C. angulata Reveals the Potential Mechanisms Underlying Thermo-Resistant Heterosis

Heterosis, also known as hybrid vigor, is widely used in aquaculture, but the molecular causes for this phenomenon remain obscure. Here, we conducted a transcriptome analysis to unveil the gene expression patterns and molecular bases underlying thermo-resistant heterosis in Crassostrea gigas ♀ × Crassostrea angulata ♂ (GA) and C. angulata ♀ × C. gigas ♂ (AG). About 505 million clean reads were obtained, and 38,210 genes were identified, of which 3779 genes were differentially expressed between the reciprocal hybrids and purebreds. The global gene expression levels were toward the C. gigas genome in the reciprocal hybrids. In GA and AG, 95.69% and 92.00% of the differentially expressed genes (DEGs) exhibited a non-additive expression pattern, respectively. We observed all gene expression modes, including additive, partial dominance, high and low dominance, and under- and over-dominance. Of these, 77.52% and 50.00% of the DEGs exhibited under- or over-dominance in GA and AG, respectively. The over-dominance DEGs common to reciprocal hybrids were significantly enriched in protein folding, protein refolding, and intrinsic apoptotic signaling pathway, while the under-dominance DEGs were significantly enriched in cell cycle. As possible candidate genes for thermo-resistant heterosis, GRP78, major egg antigen, BAG, Hsp70, and Hsp27 were over-dominantly expressed, while MCM6 and ANAPC4 were under-dominantly expressed. This study extends our understanding of the thermo-resistant heterosis in oysters.

The expanded inhibitor of apoptosis gene family in oysters possesses novel domain architectures and may play diverse roles in apoptosis following immune challenge

Background

Apoptosis plays important roles in a variety of functions, including immunity and response to environmental stress. The Inhibitor of Apoptosis (IAP) gene family of apoptosis regulators is expanded in molluscs, including eastern, Crassostrea virginica, and Pacific, Crassostrea gigas, oysters. The functional importance of IAP expansion in apoptosis and immunity in oysters remains unknown.

Results

Phylogenetic analysis of IAP genes in 10 molluscs identified lineage specific gene expansion in bivalve species. Greater IAP gene family expansion was observed in C. virginica than C. gigas (69 vs. 40), resulting mainly from tandem duplications. Functional domain analysis of oyster IAP proteins revealed 3 novel Baculoviral IAP Repeat (BIR) domain types and 14 domain architecture types across gene clusters, 4 of which are not present in model organisms. Phylogenetic analysis of bivalve IAPs suggests a complex history of domain loss and gain. Most IAP genes in oysters (76% of C. virginica and 82% of C. gigas), representing all domain architecture types, were expressed in response to immune challenge (Ostreid Herpesvirus OsHV-1, bacterial probionts Phaeobacter inhibens and Bacillus pumilus, several Vibrio spp., pathogenic Aliiroseovarius crassostreae, and protozoan parasite Perkinsus marinus). Patterns of IAP and apoptosis-related differential gene expression differed between the two oyster species, where C. virginica, in general, differentially expressed a unique set of IAP genes in each challenge, while C. gigas differentially expressed an overlapping set of IAP genes across challenges. Apoptosis gene expression patterns clustered mainly by resistance/susceptibility of the oyster host to immune challenge. Weighted Gene Correlation Network Analysis (WGCNA) revealed unique combinations of transcripts for 1 to 12 IAP domain architecture types, including novel types, were significantly co-expressed in response to immune challenge with transcripts in apoptosis-related pathways.

Conclusions

Unprecedented diversity characterized by novel BIR domains and protein domain architectures was observed in oyster IAPs. Complex patterns of gene expression of novel and conserved IAPs in response to a variety of ecologically-relevant immune challenges, combined with evidence of direct co-expression of IAP genes with apoptosis-related transcripts, suggests IAP expansion facilitates complex and nuanced regulation of apoptosis and other immune responses in oysters.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-08233-6.

Integration of mRNA and miRNA Profiling Reveals Heterosis in Oreochromis niloticus × O. aureus Hybrid Tilapia

Heterosis is a widespread biological phenomenon in fishes, in which hybrids have superior traits to parents. However, the underlying molecular basis for heterosis remains uncertain. Heterosis in growth and survival rates is apparent in hybrid tilapia (Oreochromis niloticus ♀ × O. aureus ♂). Comparisons of growth and hematological biochemical characteristics and mRNA and miRNA transcriptional analyses were performed in hybrid and parents tilapia stocks to investigate the underlying molecular basis for heterosis. Growth characteristics and hematological glucose and cholesterol parameters were significantly improved in hybrids. Of 3097 differentially expressed genes (DEGs) and 120 differentially expressed miRNAs (DEMs) identified among three stocks (O. niloticus, O. aureus, and hybrids), 1598 DEGs and 62 DEMs were non-additively expressed in hybrids. Both expression level dominance and overdominance patterns occurred for DEGs and DEMs, indicating that dominance and overdominance models are widespread in the transcriptional and post-transcriptional regulation of genes involved in growth, metabolism, immunity, and antioxidant capacity in hybrid tilapia. Moreover, potential negative regulation networks between DEMs and predicted target DEGs revealed that most DEGs from miRNA-mRNA pairs are up-regulated. Dominance and overdominance models in levels of transcriptome and miRNAome facilitate the integration of advantageous parental alleles into hybrids, contributing to heterosis of growth and improved survival. The present study provides new insights into molecular heterosis in hybrid tilapia, advancing our understanding of the complex mechanisms involved in this phenomenon in aquatic animals.

Environment-Dependent Heterosis and Transgressive Gene Expression in Reciprocal Hybrids between the Channel Catfish Ictalurus punctatus and the Blue Catfish Ictalurus furcatus

The hybrid between female channel catfish (Ictalurus punctatus) and male blue catfish (Ictalurus furcatus) is superior in feed conversion, disease resistance, carcass yield, and harvestability compared to both parental species. However, heterosis and heterobeltiosis only occur in pond culture, and channel catfish grow much faster than the other genetic types in small culture units. This environment-dependent heterosis is intriguing, but the underlying genetic mechanisms are not well understood. In this study, phenotypic characterization and transcriptomic analyses were performed in the channel catfish, blue catfish, and their reciprocal F₁s reared in tanks. The results showed that the channel catfish is superior in growth-related morphometrics, presumably due to significantly lower innate immune function, as investigated by reduced lysozyme activity and alternative complement activity. RNA-seq analysis revealed that genes involved in fatty acid metabolism/transport are significantly upregulated in channel catfish compared to blue catfish and hybrids, which also contributes to the growth phenotype. Interestingly, hybrids have a 40-80% elevation in blood glucose than the parental species, which can be explained by a phenomenon called transgressive expression (overexpression/underexpression in F₁s than the parental species). A total of 1140 transgressive genes were identified in F₁ hybrids, indicating that 8.5% of the transcriptome displayed transgressive expression. Transgressive genes upregulated in F₁s are enriched for glycan degradation function, directly related to the increase in blood glucose level. This study is the first to explore molecular mechanisms of environment-dependent heterosis/heterobeltiosis in a vertebrate species and sheds light on the regulation and evolution of heterosis vs. hybrid incompatibility.

Characterization of a C-Type Lectin Domain-Containing Protein with Antibacterial Activity from Pacific Abalone (Haliotis discus hannai)

Genes that influence the growth of Pacific abalone (Haliotis discus hannai) may improve the productivity of the aquaculture industry. Previous research demonstrated that the differential expression of a gene encoding a C-type lectin domain-containing protein (CTLD) was associated with a faster growth in Pacific abalone. We analyzed this gene and identified an open reading frame that consisted of 145 amino acids. The sequence showed a significant homology to other genes that encode CTLDs in the genus Haliotis. Expression profiling analysis at different developmental stages and from various tissues showed that the gene was first expressed at approximately 50 days after fertilization (shell length of 2.47 ± 0.13 mm). In adult Pacific abalone, the gene was strongly expressed in the epipodium, gill, and mantle. Recombinant Pacific abalone CTLD purified from Escherichia coli exhibited antimicrobial activity against several Gram-positive bacteria (Bacillus subtilis, Streptococcus iniae, and Lactococcus garvieae) and Gram-negative bacteria (Vibrio alginolyticus and Vibrio harveyi). We also performed bacterial agglutination assays in the presence of Ca²⁺, as well as bacterial binding assays in the presence of the detergent dodecyl maltoside. Incubation with E. coli and B. subtilis cells suggested that the CTLD stimulated Ca²⁺-dependent bacterial agglutination. Our results suggest that this novel Pacific abalone CTLD is important for the pathogen recognition in the gastropod host defense mechanism.

Transcriptome analysis reveals the molecular mechanisms of heterosis on thermal resistance in hybrid abalone

BACKGROUND

Heterosis has been exploited for decades in different animals and crops due to it resulting in dramatic increases in yield and adaptability. Hybridization is a classical breeding method that can effectively improve the genetic characteristics of organisms through heterosis. Abalone has become an increasingly economically important aquaculture resource with high commercial value. However, due to changing climate, abalone is now facing serious threats of high temperature in summer. Interspecific hybrid abalone (Haliotis gigantea ♀ × H. discus hannai ♂, SD) has been cultured at large scale in southern China and has been shown high survival rates under heat stress in summer. Therefore, SD has become a good model material for heterosis research, but the molecular basis of heterosis remains elusive.

RESULTS

Heterosis in thermal tolerance of SD was verified through Arrhenius break temperatures (ABT) of cardiac performance in this study. Then RNA-Sequencing was conducted to obtain gene expression patterns and alternative splicing events at control temperature (20 °C) and heat stress temperature (30 °C). A total of 356 (317 genes), 476 (435genes), and 876 (726 genes) significantly diverged alternative splicing events were identified in H. discus hannai (DD), H. gigantea (SS), and SD in response to heat stress, respectively. In the heat stress groups, 93.37% (20,512 of 21,969) of the expressed genes showed non-additive expression patterns, and over-dominance expression patterns of genes account for the highest proportion (40.15%). KEGG pathway enrichment analysis showed that the overlapping genes among common DEGs and NAGs were significantly enriched in protein processing in the endoplasmic reticulum, mitophagy, and NF-κB signaling pathway. In addition, we found that among these overlap genes, 39 genes had undergone alternative splicing events in SD. These pathways and genes may play an important role in the thermal resistance of hybrid abalone.

CONCLUSION

More alternative splicing events and non-additive expressed genes were detected in hybrid under heat stress and this may contribute to its thermal heterosis. These results might provide clues as to how hybrid abalone has a better physiological regulation ability than its parents under heat stress, to increase our understanding of heterosis in abalone.