Multi-omics Investigation of Freeze Tolerance in the Amur Sleeper, an Aquatic Ectothermic Vertebrate.

Freeze tolerance, the ability of an organism to survive internal ice formation, is a striking survival strategy employed by some ectotherms living in cold environments. However, the genetic bases of this remarkable adaptation are largely unknown. The Amur sleeper (Perccottus glenii), the only known freeze-tolerant fish species, can overwinter with its entire body frozen in ice. Here, we sequenced the chromosome-level genome of the Amur sleeper and performed comparative genomic, transcriptomic, and metabolomic analyses to investigate its strategies for surviving freezing. Evolutionary analysis suggested that the Amur sleeper diverged from its closest non-cold-hardy relative about 15.07 million years ago and has experienced a high rate of protein evolution. Transcriptomic and metabolomic data identified a coordinated and tissue-specific regulation of genes and metabolites involved in hypometabolism, cellular stress response, and cryoprotectant accumulation involved in freezing and thawing. Several genes show evidence of accelerated protein sequence evolution or family size expansion were found as adaptive responses to freezing-induced stresses. Specifically, genetic changes associated with cytoskeleton stability, cryoprotectant synthesis, transmembrane transport, and neuroprotective adaptations were identified as potentially key innovations that aid in freezing survival. Our work provides valuable resources and opportunities to unveil the molecular adaptations supporting freeze tolerance in ectothermic vertebrates.

Pseudo-chromosome-length genome assembly for a deep-sea eel Ilyophis brunneus sheds light on the deep-sea adaptation.

High hydrostatic pressure, low temperature, and scarce food supply are the major factors that limit the survival of vertebrates in extreme deep-sea environments. Here, we constructed a high-quality genome of the deep-sea Muddy arrowtooth eel (MAE, Ilyophis brunneus, captured below a depth of 3,500 m) by using Illumina, PacBio, and Hi-C sequencing. We compare it against those of shallow-water eel and other outgroups to explore the genetic basis that underlies the adaptive evolution to deep-sea biomes. The MAE genome was estimated to be 1.47 Gb and assembled into 14 pseudo-chromosomes. Phylogenetic analyses indicated that MAE diverged from its closely related shallow-sea species, European eel, ∼111.9 Mya and experienced a rapid evolution. The genome evolutionary analyses primarily revealed the following: (i) under high hydrostatic pressure, the positively selected gene TUBGCP3 and the expanded family MLC1 may improve the cytoskeleton stability; ACOX1 may enhance the fluidity of cell membrane and maintain transport activity; the expansion of ABCC12 gene family may enhance the integrity of DNA; (ii) positively selected HARS likely maintain the transcription ability at low temperatures; and (iii) energy metabolism under a food-limited environment may be increased by expanded and positively selected genes in AMPK and mTOR signaling pathways.

Extensive sequence divergence between the reference genomes of two zebrafish strains, Tuebingen and AB.

The zebrafish (Danio rerio) is one of the most widely used model organisms for studying vertebrate gene function and human disease, given the 70% conserved protein-coding genes between zebrafish and human. Two of the most common laboratory zebrafish strains are Tuebingen and AB. Despite the fact that the zebrafish reference genome is derived from the Tuebingen strain, the AB strain is still widely used although a high-quality genome comparable to Tuebingen is lacking. Here, we report a 1.40-Gb representative de novo genome assembly of the AB strain (DrAB1), with contig N50 length of 21 Mb, by integrating Illumina short-read sequencing, Nanopore long-read sequencing and HiC-based chromatin mapping. Compared with the published zebrafish Zv11 reference genome (GRCz11), this genome assembly shows considerable improvements in both contiguity and completeness. In addition, substantial structural differences and extensive sequence divergence of unprecedented resolution have been uncovered, especially with respect to 9,029,929 single nucleotide polymorphisms, 2,376,812 InDels, 32,623 insertions, 22,089 deletions and 220 inversions, which constitute ~2.6% of the DrAB1 genome. Many of these variants may have potential functional effects on phenotype, which should be considered in further experimental designs. Consequently, our study provides additional genomic resources and a high-resolution structural variation map based on whole-genome alignment for the zebrafish community, which could also be an indispensable reference genome from a model species in future research on fish phylogenetic genomics, comparative genomics and adaptive evolution.

Sympatric or micro-allopatric speciation in a glacial lake? Genomic islands support neither.

Apparent cases of sympatric speciation may actually be due to micro-allopatric or micro-parapatric speciation. One way to distinguish between these models is to examine the existence and nature of genomic islands of divergence, wherein divergent DNA segments are interspersed with low-divergence segments. Such islands should be rare or absent under micro-allopatric speciation but common in cases of speciation with gene flow. Sympatric divergence of endemic fishes is known from isolated saline, crater, postglacial, and ancient lakes. Two morphologically distinct cyprinid fishes, Gymnocypris eckloni scoliostomus (GS) and G. eckloni eckloni (GE), in a small glacial lake on the Qinghai-Tibet Plateau, Lake Sunmcuo, match the biogeographic criteria of sympatric speciation. In this study, we examined genome-wide variation in 46 individuals from these two groups. The divergence time between the GS and GE lineages was estimated to be 20-60 Kya. We identified 54 large genomic islands (≥100 kb) of speciation, which accounted for 89.4% of the total length of all genomic islands. These islands harboured divergent genes related to olfactory receptors and olfaction signals that may play important roles in food selection and assortative mating in fishes. Although the genomic islands clearly indicated speciation with gene flow and rejected micro-allopatric speciation, they were too large to support the hypothesis of sympatric speciation. Theoretical and recent empirical studies suggested that continual gene flow in sympatry should give rise to many small genomic islands (as small as a few kilobases in size). Thus, the observed pattern is consistent with the extensive evidence on parapatric speciation, in which adjacent habitats facilitate divergent selection but also permit gene flow during speciation. We suggest that many, if not most, of the reported cases of sympatric speciation are likely to be micro-parapatric speciation.

The new chimeric chiron genes evolved essential roles in zebrafish embryonic development by regulating NAD+ levels.

The origination of new genes is important for generating genetic novelties for adaptive evolution and biological diversity. However, their potential roles in embryonic development, evolutionary processes into ancient networks, and contributions to adaptive evolution remain poorly investigated. Here, we identified a novel chimeric gene family, the chiron family, and explored its genetic basis and functional evolution underlying the adaptive evolution of Danioninae fishes. The ancestral chiron gene originated through retroposition of nampt in Danioninae 48-54 million years ago (Mya) and expanded into five duplicates (chiron1-5) in zebrafish 1-4 Mya. The chiron genes (chirons) likely originated in embryonic development and gradually extended their expression in the testis. Functional experiments showed that chirons were essential for zebrafish embryo development. By integrating into the NAD+ synthesis pathway, chirons could directly catalyze the NAD+ rate-limiting reaction and probably impact two energy metabolism genes (nmnat1 and naprt) to be under positive selection in Danioninae fishes. Together, these results mainly demonstrated that the origin of new chimeric chiron genes may be involved in adaptive evolution by integrating and impacting the NAD+ biosynthetic pathway. This coevolution may contribute to the physiological adaptation of Danioninae fishes to widespread and varied biomes in Southeast Asian.

Transcriptome Analysis of In Vitro Fertilization and Parthenogenesis Activation during Early Embryonic Development in Pigs.

Parthenogenesis activation (PA), as an important artificial breeding method, can stably preserve the dominant genotype of a species. However, the delayed development of PA embryos is still overly severe and largely leads to pre-implantation failure in pigs. The mechanisms underlying the deficiencies of PA embryos have not been completely understood. For further understanding of the molecular mechanism behind PA embryo failure, we performed transcriptome analysis among pig oocytes (meiosis II, MII) and early embryos at three developmental stages (zygote, morula, and blastocyst) in vitro fertilization (IVF) and PA group. Totally, 11,110 differentially expressed genes (DEGs), 4694 differentially expressed lincRNAs (DELs) were identified, and most DEGs enriched the regulation of apoptotic processes. Through cis- and trans-manner functional prediction, we found that hub lincRNAs were mostly involved in abnormal parthenogenesis embryonic development. In addition, twenty DE imprinted genes showed that some paternally imprinted genes in IVF displayed higher expression than that in PA. Notably, we identified that three DELs of imprinted genes (MEST, PLAGL1, and DIRAS3) were up regulated in IVF, and there was no significant change in PA group. Disordered expression of key genes for embryonic development might play key roles in abnormal parthenogenesis embryonic development. Our study indicates that embryos derived from different production techniques have varied in vitro development to the blastocyst stage, and they also affect the transcription level of corresponding genes, such as imprinted genes. This work will help future research on these genes and molecular-assisted breeding for pig parthenotes.

Characterization and Functional Analysis of Polyadenylation Sites in Fast and Slow Muscles.

Many increasing documents have proved that alternative polyadenylation (APA) events with different polyadenylation sites (PAS) contribute to posttranscriptional regulation. However, little is known about the detailed molecular features of PASs and its role in porcine fast and slow skeletal muscles through microRNAs (miRNAs) and RNA binding proteins (RBPs). In this study, we combined single-molecule real-time sequencing and Illumina RNA-seq datasets to comprehensively analyze polyadenylation in pigs. We identified a total of 10,334 PASs, of which 8734 were characterized by reference genome annotation. 32.86% of PAS-associated genes were determined to have more than one PAS. Further analysis demonstrated that tissue-specific PASs between fast and slow muscles were enriched in skeletal muscle development pathways. In addition, we obtained 1407 target genes regulated by APA events through potential binding 69 miRNAs and 28 RBPs in variable 3' UTR regions and some are involved in myofiber transformation. Furthermore, the de novo motif search confirmed that the most common usage of canonical motif AAUAAA and three types of PASs may be related to the strength of motifs. In summary, our results provide a useful annotation of PASs for pig transcriptome and suggest that APA may serve as a role in fast and slow muscle development under the regulation of miRNAs and RBPs.

MEG3 Promotes Differentiation of Porcine Satellite Cells by Sponging miR-423-5p to Relieve Inhibiting Effect on SRF.

Although thousands of long noncoding RNAs (lncRNAs) have been identified in porcine growth and development, the regulation mechanisms of functional lncRNAs have not been well explored. In this study, using 5'- and 3'-rapid amplification of cDNA ends (RACE) assays, we obtained two different variants of lncRNA maternally expressed gene 3 (MEG3), namely, MEG3 v1 and MEG3 v2, that were both highly expressed in porcine skeletal muscle and in the early stage of the differentiation of porcine satellite cells. Moreover, we identified the core transcript MEG3 v2. Functional analyses showed that MEG3 overexpression could effectively arrest myoblasts in the G1 phase, inhibit DNA replication, and promote myoblast differentiation, whereas MEG3 knockdown resulted in the opposite effects. Interestingly, the expression of serum response factor (SRF), a crucial transcription factor for myogenesis process, remarkably increased and decreased in mRNA and protein levels with the respective overexpression and knockdown of MEG3. Dual luciferase reporter assay showed that MEG3 could attenuate the decrease of luciferase activity of SRF induced by miR-423-5p in a dose-dependent manner. MEG3 overexpression could relieve the inhibitory effect on SRF and myoblast differentiation induced by miR-423-5p. In addition, results of RNA immunoprecipitation analysis suggested that MEG3 could act as a ceRNA for miR-423-5p. Our findings initially established a novel connection among MEG3, miR-423-5p, and SRF in porcine satellite cell differentiation. This novel role of MEG3 may shed new light on understanding of molecular regulation of lncRNA in porcine myogenesis.

RAI14 in the blood feather regulates chicken pigmentation.

A genome-wide association study (GWAS) was performed on a resource family consisting of white and colored chickens for identification of genes related to plumage coloration using the Fixed and random model Circulating Probability Unification (FarmCPU) package. GWAS identified three chromosomal single-nucleotide polymorphisms (SNPs), demonstrating the polygenic basis of plumage phenotypes. Herein, retinoic acid-induced protein 14 (RAI14), a developmentally regulated gene that encodes a protein containing many ankyrin repeats, was identified as a candidate gene involved in plumage color. In this study, mRNA expression profiles of chicken RAI14 were determined, indel (insertion-deletion) variants were identified, and their association was analyzed in white and colored chickens. RA114 mRNA was expressed in all tissues tested (brain, spleen, liver, heart, oviduct, kidney, lung, pituitary gland, ovary, muscle, feather bulb, and skin). A relatively high RAI14 expression in white feather bulb compared to colored feather bulb ( P < 0.01 ) indicated a potential association with plumage color. Additionally, statistical analysis revealed that a 4 bp indel genetic variation in RAI14 was associated with plumage phenotypes ( P < 0.01 ). Together, our analysis of the identification of the RAI14 gene will enable us to understand the genetic mechanisms behind chicken pigmentation.

Long Non-coding RNA H19 Regulates Porcine Satellite Cell Differentiation Through miR-140-5p/SOX4 and DBN1.

The H19 gene promotes skeletal muscle differentiation in mice, but the regulatory models and mechanisms of myogenesis regulated by H19 are largely unknown in pigs. Therefore, the regulatory modes of H19 in the differentiation of porcine skeletal muscle satellite cells (PSCs) need to be determined. We observed that H19 gene silencing could decrease the expressions of the myogenin (MYOG) gene, myogenic differentiation (MYOD), and myosin heavy chain (MYHC) in PSCs. Therefore, we constructed and sequenced 12 cDNA libraries of PSCs after knockdown of H19 at two differentiation time points to analyze the transcriptome differences. A total of 11,419 differentially expressed genes (DEGs) were identified. Among these DEGs, we found through bioinformatics analysis and protein interaction experiment that SRY-box transcription factor 4 (SOX4) and Drebrin 1 (DBN1) were the key genes in H19-regulated PSC differentiation. Functional analysis shows that SOX4 and DBN1 promote PSC differentiation. Mechanistically, H19 regulates PSC differentiation through two different pathways. On the one hand, H19 functions as a molecular sponge of miR-140-5p, which inhibits the differentiation of PSCs, thereby modulating the derepression of SOX4. On the other hand, H19 regulates PSC differentiation through directly binding with DBN1. Furthermore, MYOD binds to the promoters of H19 and DBN1. The knockdown of MYOD inhibits the expression of H19 and DBN1. We determined the function of H19 and provided a molecular model to elucidate H19's role in regulating PSC differentiation.

Identification and Functional Prediction of Long Intergenic Non-coding RNAs Related to Subcutaneous Adipose Development in Pigs.

An increasing number of studies have shown that long intergenic non-coding RNAs (lincRNAs) are a very important class of non-coding RNAs that plays a vital role in many biological processes. Adipose tissue is an important place for storing energy, but few studies on lincRNAs were related to pig subcutaneous fat development. Here, we used published RNA-seq data from subcutaneous adipose tissue of Italian Large White pigs and identified 252 putative lincRNAs, wherein 34 were unannotated. These lincRNAs had relatively shorter length, lower number of exons, and lower expression level compared with protein-coding transcripts. Gene ontology and pathway analysis indicated that the adjacent genes of lincRNAs were involved in lipid metabolism. In addition, differentially expressed lincRNAs (DELs) between low and high backfat thickness pigs were identified. Through the detection of quantitative trait locus (QTL), DELs were mainly located in QTLs related to adipose development. Based on the expression correlation of DEL genes and their differentially expressed potential target genes, we constructed a co-expression network and a potential pathway of DEL's effect on lipid metabolism. Our study identified and analyzed lincRNAs in subcutaneous adipose tissue, and results suggested that lincRNAs may be involved in the regulation of subcutaneous fat development. Our findings provided new insights into the biological function of porcine lincRNAs.

Transcriptome Analysis Suggests the Roles of Long Intergenic Non-coding RNAs in the Growth Performance of Weaned Piglets.

Long intergenic non-coding RNAs (lincRNAs) have been considered to play a key regulatory role in various biological processes. An increasing number of studies have utilized transcriptome analysis to obtain lincRNAs with functions related to cancer, but lincRNAs affecting growth rates in weaned piglets are rarely described. Although lincRNAs have been systematically identified in various mouse tissues and cell lines, studies of lincRNA in pigs remain rare. Therefore, identifying and characterizing novel lincRNAs affecting the growth performance of weaned piglets is of great importance. Here, we reconstructed 101,988 lincRNA transcripts and identified 1,078 lincRNAs in two groups of longissimus dorsi muscle (LDM) and subcutaneous fat (SF) based on published RNA-seq datasets. These lincRNAs exhibit typical characteristics, such as shorter lengths and lower expression relative to protein-encoding genes. Gene ontology analysis revealed that some lincRNAs could be involved in weaned piglet related processes, such as insulin resistance and the AMPK signaling pathway. We also compared the positional relationship between differentially expressed lincRNAs (DELs) and quantitative trait loci (QTL) and found that some of DELs may play an important role in piglet growth and development. Our work details part of the lincRNAs that may affect the growth performance of weaned piglets and promotes future studies of lincRNAs for molecular-assisted development in weaned piglets.

Transcriptome Analysis Reveals the Effect of Long Intergenic Noncoding RNAs on Pig Muscle Growth and Fat Deposition.

Muscle growth and fat deposition are the two important biological processes in the development of pigs which are closely related to the pig production performance. Long intergenic noncoding RNAs (lincRNAs), with lack of coding potential and the length of at least 200nt, have been extensively studied to play important roles in many biological processes. However, the importance and molecular regulation mechanism of lincRNAs in the process of muscle growth and fat deposition in pigs are still to be further studied comprehensively. In our study, we used the data, including liver, abdominal fat, and longissimus dorsi muscle of 240 days' age of two F2 full-sib female individuals from the white Duroc and Erhualian crossbreed, to identify 581 putative lincRNAs associated with pig muscle growth and fat deposition. The 581 putative lincRNAs shared many common features with other mammalian lincRNAs, such as fewer exons, lower expression levels, and shorter transcript lengths. Cross-tissue comparisons showed that many transcripts were tissue-specific and were involved in the important biological processes in their corresponding tissues. Gene ontology and pathway analysis revealed that many potential target genes (PTGs) of putative lincRNAs were involved in pig muscle growth and fat deposition-related processes, including muscle cell proliferation, lipid metabolism, and fatty acid degradation. In Quantitative Trait Locus (QTLs) analysis, some PTGs were screened from putative lincRNAs, MRPL12 is associated with muscle growth, GCGR and SLC25A10 were associated with fat deposition, and PPP3CA, DPYD, and FGGY were related not only to muscle growth but also to fat deposition. Therefore, it implied that these lincRNAs might participate in the biological processes related to muscle growth or fat deposition through homeostatic regulation of PTGs, but the detailed molecular regulatory mechanisms still needed to be further explored. This study lays the molecular foundation for the in-depth study of the role of lincRNAs in the pig muscle growth and fat deposition and further provides the new molecular markers for understanding the complex biological mechanisms of pig muscle growth and fat deposition.

Genomic Analysis To Identify Signatures of Artificial Selection and Loci Associated with Important Economic Traits in Duroc Pigs.

Identifying genetic basis of domestication and improvement in livestock contributes to our understanding of the role of artificial selection in shaping the genome. Here we used whole-genome sequencing and the genotyping by sequencing approach to detect artificial selection signatures and identify the associated SNPs of two economic traits in Duroc pigs. A total of 38 candidate selection regions were detected by combining the fixation index and the Composite Likelihood Ratio methods. Further genome-wide association study revealed seven associated SNPs that were related with intramuscular fat content and feed conversion ratio traits, respectively. Enrichment analysis suggested that the artificial selection regions harbored genes, such as MSTN, SOD2, MC5R and CD83, which are responsible for economic traits including lean muscle mass, fertility and immunization. Overall, this study found a series of candidate genes putatively associated with the breeding improvement of Duroc pigs and the polygenic basis of adaptive evolution, which can provide important references and fundamental information for future breeding programs.

DNA methylation changes and evolution of RNA-based duplication in Sus scrofa: based on a two-step strategy.

AIM: This study aims to couple DNA methylation changes and evolution of retrogenes. MATERIALS & METHODS: A new two-step strategy was developed to screen retrogenes. Further, reduced representation bisulfite sequencing and RNA-seq data of eight tissues were used to analyze retrogenes. RESULTS: A total of 964 retrocopies were identified and new retrocopies were available for the synthesis of glycans and lipids corresponding to pig phenotypic traits. Retrogenes were consistently hypermethylated. Hypomethylation of parental genes presented more susceptibility to retroposition. Promoter DNA methylation of retrogenes was negatively correlated with evolutionary time and played important roles in regulating retrogene tissue-specific expression pattern. CONCLUSION: A two-step procedure is effective and necessary for identifying retrogenes. DNA methylation drives origination, survival, evolution and expression of retrogenes.

Ancient duplications and functional divergence in the interferon regulatory factors of vertebrates provide insights into the evolution of vertebrate immune systems.

Interferon regulatory factors (IRFs) were first discovered as transcription factors that regulate the transcription of human interferon (IFN)-ß. Increasing evidence shows that they might be important players involved in Adaptive immune system (AIS) evolution. Although numbers of IRFs have been identified in chordates, the evolutionary history and functional diversity of this gene family during the early evolution of vertebrates have remained obscure. Using IRF HMM profile and HMMER searches, we identified 148 IRFs in 11 vertebrates and 4 protochordates. For them, we reconstructed the phylogenetic relationships, determined the synteny conservation, investigated the profile of natural selection, and analyzed the expression patterns in four "living fossil" vertebrates: lamprey, elephant shark, coelacanth and bichir. The results from phylogeny and synteny analysis imply that vertebrate IRFs evolved from three predecessors, instead of four as suggested in a previous study, as results from an ancient duplication followed by special expansions and lost during the vertebrate evolution. The profile of natural selection and expression reveals functional dynamics during the process. Together, they suggest that the 2nd whole-genome duplication (2WGD) provided raw materials for innovation in the IRF family, and that the birth of type-I IFN might be an important factor inducing the establishment of IRF-mediated immune networks. As a member involved in the AIS evolution, IRF provide insights into the process and mechanism involved in the complexity and novelties of vertebrate immune systems.

Identification and Functional Analysis of Long Intergenic Non-coding RNAs Underlying Intramuscular Fat Content in Pigs.

Intramuscular fat (IMF) content is an important trait that can affect pork quality. Previous studies have identified many genes that can regulate IMF. Long intergenic non-coding RNAs (lincRNAs) are emerging as key regulators in various biological processes. However, lincRNAs related to IMF in pig are largely unknown, and the mechanisms by which they regulate IMF are yet to be elucidated. Here we reconstructed 105,687 transcripts and identified 1,032 lincRNAs in pig longissimus dorsi muscle (LDM) of four stages with different IMF contents based on published RNA-seq. These lincRNAs show typical characteristics such as shorter length and lower expression compared with protein-coding genes. Combined with methylation data, we found that both the promoter and genebody methylation of lincRNAs can negatively regulate lincRNA expression. We found that lincRNAs exhibit high correlation with their protein-coding neighbors in expression. Co-expression network analysis resulted in eight stage-specific modules, gene ontology and pathway analysis of them suggested that some lincRNAs were involved in IMF-related processes, such as fatty acid metabolism and peroxisome proliferator-activated receptor signaling pathway. Furthermore, we identified hub lincRNAs and found six of them may play important roles in IMF development. This work detailed some lincRNAs which may affect of IMF development in pig, and facilitated future research on these lincRNAs and molecular assisted breeding for pig.