Combined Transcriptomics and Metabolomics Identify Regulatory Mechanisms of Porcine Vertebral Chondrocyte Development In Vitro

Porcine body length is closely related to meat production, growth, and reproductive performance, thus playing a key role in the profitability of the pork industry. Cartilage development is critical to longitudinal elongation of individual vertebrae. This study isolated primary porcine vertebral chondrocytes (PVCs) to clarify the complex mechanisms of elongation. We used transcriptome and target energy metabolome technologies to confirm crucial genes and metabolites in primary PVCs at different differentiation stages (0, 4, 8, and 12 days). Pairwise comparisons of the four stages identified 4566 differentially expressed genes (DEGs). Time-series gene cluster and functional analyses of these DEGs revealed four clusters related to metabolic processes, cartilage development, vascular development, and cell cycle regulation. We constructed a transcriptional regulatory network determining chondrocyte maturation. The network indicated that significantly enriched transcription factor (TF) families, including zf-C2H2, homeobox, TF_bZIP, and RHD, are important in cell cycle and differentiation processes. Further, dynamic network biomarker (DNB) analysis revealed that day 4 was the tipping point for chondrocyte development, consistent with morphological and metabolic changes. We found 24 DNB DEGs, including the TFs NFATC2 and SP7. Targeted energy metabolome analysis showed that most metabolites were elevated throughout chondrocyte development; notably, 16 differentially regulated metabolites (DRMs) were increased at three time points after cell differentiation. In conclusion, integrated metabolome and transcriptome analyses highlighted the importance of amino acid biosynthesis in chondrocyte development, with coordinated regulation of DEGs and DRMs promoting PVC differentiation via glucose oxidation. These findings reveal the regulatory mechanisms underlying PVC development and provide an important theoretical reference for improving pork production.

A chromosome-level reference genome assembly of the Reeve's moray eel (Gymnothorax reevesii)

Due to potentially hostile behaviors and elusive habitats, moray eels (Muraenidae) as one group of apex predators in coral reefs all across the globe have not been well investigated. Here, we constructed a chromosome-level genome assembly for the representative Reeve's moray eel (Gymnothorax reevesii). This haplotype genome assembly is 2.17 Gb in length, and 97.87% of the sequences are anchored into 21 chromosomes. It contains 56.34% repetitive sequences and 23,812 protein-coding genes, of which 96.77% are functionally annotated. This sequenced marine species in Anguilliformes makes a good complement to the genetic resource of eel genomes. It not only provides a genetic resource for in-depth studies of the Reeve's moray eel, but also enables deep-going genomic comparisons among various eels.

Integrated Analysis of Transcriptome Expression Profiles Reveals miRNA-326-NKX3.2-Regulated Porcine Chondrocyte Differentiation

The porcine body length trait is an essential factor affecting meat production and reproductive performance. It is evident that the development/lengthening of individual vertebrae is one of the main reasons for increases in body length; however, the underlying molecular mechanism remains unclear. In this study, RNA-seq analysis was used to profile the transcriptome (lncRNA, mRNA, and miRNA) of the thoracic intervertebral cartilage (TIC) at two time points (1 and 4 months) during vertebral column development in Yorkshire (Y) and Wuzhishan pigs (W). There were four groups: 1- (Y1) and 4-month-old (Y4) Yorkshire pigs and 1- (W1) and 4-month-old (W4) Wuzhishan pigs. In total, 161, 275, 86, and 126 differentially expressed (DE) lncRNAs, 1478, 2643, 404, and 750 DE genes (DEGs), and 74,51, 34, and 23 DE miRNAs (DE miRNAs) were identified in the Y4 vs. Y1, W4 vs. W1, Y4 vs. W4, and Y1 vs. W1 comparisons, respectively. Functional analysis of these DE transcripts (DETs) demonstrated that they had participated in various biological processes, such as cellular component organization or biogenesis, the developmental process, the metabolic process, bone development, and cartilage development. The crucial bone development-related candidate genes NK3 Homeobox 2 (NKX3.2), Wnt ligand secretion mediator (WLS), gremlin 1 (GREM1), fibroblast growth factor receptor 3 (FGFR3), hematopoietically expressed homeobox (HHEX), (collagen type XI alpha 1 chain (COL11A1), and Wnt Family Member 16 (WNT16)) were further identified by functional analysis. Moreover, lncRNA, miRNA, and gene interaction networks were constructed; a total of 55 lncRNAs, 6 miRNAs, and 7 genes formed lncRNA-gene, miRNA-gene, and lncRNA-miRNA-gene pairs, respectively. The aim was to demonstrate that coding and non-coding genes may co-regulate porcine spine development through interaction networks. NKX3.2 was identified as being specifically expressed in cartilage tissues, and it delayed chondrocyte differentiation. miRNA-326 regulated chondrocyte differentiation by targeting NKX3.2. The present study provides the first non-coding RNA and gene expression profiles in the porcine TIC, constructs the lncRNA-miRNA-gene interaction networks, and confirms the function of NKX3.2 in vertebral column development. These findings contribute to the understanding of the potential molecular mechanisms regulating pig vertebral column development. They expand our knowledge about the differences in body length between different pig species and provide a foundation for future studies.

Body shape transformations by alternate anatomical adaptive peak shifts in blenniiform fishes

Extreme body elongation has occurred repeatedly in the evolutionary history of ray-finned fishes. Lengthening of the anterior-posterior body axis relative to depth and width can involve changes in the cranial skeleton and vertebral column, but to what extent is anatomical evolution determined by selective factors and intrinsic constraints that are shared broadly among closely related lineages? In this study, we fit adaptive (Ornstein-Uhlenbeck) evolutionary models to body shape and its anatomical determinants and identified two instances of extreme elongation by divergent anatomical peak shifts in the Blenniiformes, a radiation of small-bodied substrate-associated marine teleost fishes. Species in the genus Xiphasia (hairtail blennies) evolved toward a peak defined by a highly elongated caudal vertebral region but ancestral cranial and precaudal vertebral morphology. In contrast, a clade that includes the genera Chaenopsis and Lucayablennius (pike and arrow blennies) evolved toward a peak with a long slender skull but ancestral axial skeletal anatomy. Neither set of anatomical peak shifts aligns closely with the major axis of anatomical diversification in other blenniiform fishes. These results provide little evidence that ancestral constraints have affected body shape transformation, and instead suggest that extreme elongation arose with distinct shifts in selective factors and development.

Repeated ecological and life cycle transitions make salamanders an ideal model for evolution and development

Observations on the ontogeny and diversity of salamanders provided some of the earliest evidence that shifts in developmental trajectories have made a substantial contribution to the evolution of animal forms. Since the dawn of evo-devo there have been major advances in understanding developmental mechanisms, phylogenetic relationships, evolutionary models, and an appreciation for the impact of ecology on patterns of development (eco-evo-devo). Molecular phylogenetic analyses have converged on strong support for the majority of branches in the Salamander Tree of Life, which includes 764 described species. Ancestral reconstructions reveal repeated transitions between life cycle modes and ecologies. The salamander fossil record is scant, but key Mesozoic species support the antiquity of life cycle transitions in some families. Colonization of diverse habitats has promoted phenotypic diversification and sometimes convergence when similar environments have been independently invaded. However, unrelated lineages may follow different developmental pathways to arrive at convergent phenotypes. This article summarizes ecological and endocrine-based causes of life cycle transitions in salamanders, as well as consequences to body size, genome size, and skeletal structure. Salamanders offer a rich source of comparisons for understanding how the evolution of developmental patterns has led to phenotypic diversification following shifts to new adaptive zones.

Structural and environmental constraints on reduction of paired appendages among vertebrates

Burrowing habits or complex environments have generally been considered as potential drivers acting on reduction and loss of the appendicular skeleton among vertebrates. Herein, we suggest that this might be the case for lissamphibians and squamates, but that fin loss in fishes is usually prevented by important structural constraints, because pectoral fins are commonly used to control rolling and pitching. We provide an overview of the distribution of paired appendage reduction across vertebrates while examining the ecological affinities of finless and limbless clades. We analysed the correlation between lifestyle and fin or limb loss using the discrete comparative analysis. The resulting Bayesian factors indicate strong evidence of correlation between: (1) pectoral-fin loss and coexistence of anguilliform elongation and burrowing habits or complex habitat in teleost fishes; and (2) limb loss and a burrowing or grass-swimming lifestyle in squamate reptiles and lissamphibians. These correlations suggest that a complex environment or a fossorial habit is a driving force leading to appendage loss. The only style of locomotion that is functional even in the absence of paired appendages is the undulatory one, which is typical of all elongated reptiles and lissamphibians, but certainly less common in teleost fishes.

Frequent nonrandom shifts in the temporal sequence of developmental landmark events during teleost evolutionary diversification

Morphological transformations can be generated by evolutionary changes in the sequence of developmental events. In this study, we examined the evolutionary dynamics of the developmental sequence on a macroevolutionary scale in teleosts. Using the information from previous reports describing the development of 31 species, we extracted the developmental sequences of 19 landmark events involving the formation of phylogenetically conserved body parts; we then inferred ancestral developmental sequences by two different parsimony-based methods-event-pairing and continuous analysis. The phylogenetic comparisons of these sequences revealed event-dependent heterogeneity in the frequency of sequence changes. Most of the sequence changes occurred as exchanges of temporally neighboring events. These heterochronic changes in developmental sequences accumulated along evolutionary time, but the precise distribution of the changes over the teleostean phylogeny remains unclear due to technical limitations.

Evidence for complex life cycle constraints on salamander body form diversification

Metazoans display a tremendous diversity of developmental patterns, including complex life cycles composed of morphologically disparate stages. In this regard, the evolution of life cycle complexity promotes phenotypic diversity. However, correlations between life cycle stages can constrain the evolution of some structures and functions. Despite the potential macroevolutionary consequences, few studies have tested the impacts of life cycle evolution on broad-scale patterns of trait diversification. Here we show that larval and adult salamanders with a simple, aquatic-only (paedomorphic) life cycle had an increased rate of vertebral column and body form diversification compared to lineages with a complex, aquatic-terrestrial (biphasic) life cycle. These differences in life cycle complexity explain the variations in vertebral number and adult body form better than larval ecology. In addition, we found that lineages with a simple terrestrial-only (direct developing) life cycle also had a higher rate of adult body form evolution than biphasic lineages, but still 10-fold lower than aquatic-only lineages. Our analyses demonstrate that prominent shifts in phenotypic evolution can follow long-term transitions in life cycle complexity, which may reflect underlying stage-dependent constraints.

No support for Heincke's law in hagfish (Myxinidae): lack of an association between body size and the depth of species occurrence

This study tests for interspecific evidence of Heincke's law among hagfishes and advances the field of research on body size and depth of occurrence in fishes by including a phylogenetic correction and by examining depth in four ways: maximum depth, minimum depth, mean depth of recorded specimens and the average of maximum and minimum depths of occurrence. Results yield no evidence for Heincke's law in hagfishes, no phylogenetic signal for the depth at which species occur, but moderate to weak phylogenetic signal for body size, suggesting that phylogeny may play a role in determining body size in this group.

Prolongation of somitogenesis in two anguilliform species, the Japanese eel Anguilla japonica and pike eel Muraenesox cinereus, with refined descriptions of their early development

The embryonic development of the Japanese eel Anguilla japonica and pike eel Muraenesox cinereus was morphologically investigated with laboratory-reared specimens to clarify the characteristics of somitogenesis. In A. japonica, somites were first observed at 18 h post fertilization (hpf) when epiboly reached 90%. Somitogenesis progressed at a rate of 1·6 h^-1 at mean ± s.d. 22·6 ± 0·7° C and completed at 107 hpf (3 days post hatching; dph) when total number of somites (ST) reached 114, which corresponds to the species' number of vertebrae (112-119). In M. cinereus, somites were first observed at 14 hpf when epiboly completed. Somitogenesis progressed at a rate of 1·9 h^-1 at mean ± s.d. 24·4 ± 0·2° C and completed at 90 hpf (2 dph) with 149 ± 4 ST, which corresponds to the species' number of vertebrae (142-158). Both species hatched before somitogenesis was completed, at 37 hpf with 47 ST and 42 hpf with 82 ± 4 ST, respectively. The formation of other organs such as the heart, mouth and pectoral fin bud occurred during somitogenesis. Comparison with the development of zebrafish Danio rerio indicates a prolongation of somitogenesis in A. japonica and M. cinereus. Their somitogenesis rates, however, correspond well with that of D. rerio estimated at the same temperature and their developmental stages at hatching are almost equivalent to other fishes having similar yolk sizes. Therefore, the prolongation of somitogenesis in A. japonica and M. cinereus may be accounted for solely by the increased numbers of somites to be formed, not by a slow somitogenesis rate or an acceleration in organogenesis.

Evolution of the locomotory system in eels (Teleostei: Elopomorpha)

Background

Living anguilliform eels represent a distinct clade of elongated teleostean fishes inhabiting a wide range of habitats. Locomotion of these fishes is highly influenced by the elongated body shape, the anatomy of the vertebral column, and the corresponding soft tissues represented by the musculotendinous system. Up to now, the evolution of axial elongation in eels has been inferred from living taxa only, whereas the reconstruction of evolutionary patterns and functional ecology in extinct eels still is scarce. Rare but excellently preserved fossil eels from the Late Cretaceous and Cenozoic were investigated here to gain a better understanding of locomotory system evolution in anguilliforms and, consequently, their habitat occupations in deep time.

Results

The number of vertebrae in correlation with the body length separates extinct and extant anguilliforms. Even if the phylogenetic signal cannot entirely be excluded, the analyses performed here reveal a continuous shortening of the vertebral column with a simultaneous increase in vertebral numbers in conjunction with short lateral tendons throughout the order. These anatomical changes contradict previous hypotheses based on extant eels solely.

Conclusions

The body curvatures of extant anguilliforms are highly flexible and can be clearly distinguished from extinct species. Anatomical changes of the vertebral column and musculotendinous system through time and between extinct and extant anguilliforms correlate with changes of the body plan and swimming performance and reveal significant shifts in habitat adaptation and thus behaviour. Evolutionary changes in the skeletal system of eels established here also imply that environmental shifts were triggered by abiotic rather than biotic factors (e.g., K/P boundary mass extinction event).

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-016-0728-7) contains supplementary material, which is available to authorized users.

Morphological variation in salamanders and their potential response to climate change

Despite the recognition that some species might quickly adapt to new conditions under climate change, demonstrating and predicting such a fundamental response is challenging. Morphological variations in response to climate may be caused by evolutionary changes or phenotypic plasticity, or both, but teasing apart these processes is difficult. Here, we built on the number of thoracic vertebrae (NTV) in ectothermic vertebrates, a known genetically based feature, to establish a link with body size and evaluate how climate change might affect the future morphological response of this group of species. First, we show that in old-world salamanders, NTV variation is strongly related to changes in body size. Secondly, using 22 salamander species as a case study, we found support for relationships between the spatial variation in selected bioclimatic variables and NTV for most of species. For 44% of species, precipitation and aridity were the predominant drivers of geographical variation of the NTV. Temperature features were dominant for 31% of species, while for 19% temperature and precipitation played a comparable role. This two-step analysis demonstrates that ectothermic vertebrates may evolve in response to climate change by modifying the number of thoracic vertebrae. These findings allow to develop scenarios for potential morphological evolution under future climate change and to identify areas and species in which the most marked evolutionary responses are expected. Resistance to climate change estimated from species distribution models was positively related to present-day species morphological response, suggesting that the ability of morphological evolution may play a role for species' persistence under climate change. The possibility that present-day capacity for local adaptation might help the resistance response to climate change can be integrated into analyses of the impact of global changes and should also be considered when planning management actions favouring species persistence.

A morphospace for reef fishes: elongation is the dominant axis of body shape evolution

Tropical reef fishes are widely regarded as being perhaps the most morphologically diverse vertebrate assemblage on earth, yet much remains to be discovered about the scope and patterns of this diversity. We created a morphospace of 2,939 species spanning 56 families of tropical Indo-Pacific reef fishes and established the primary axes of body shape variation, the phylogenetic consistency of these patterns, and whether dominant patterns of shape change can be accomplished by diverse underlying changes. Principal component analysis showed a major axis of shape variation that contrasts deep-bodied species with slender, elongate forms. Furthermore, using custom methods to compare the elongation vector (axis that maximizes elongation deformation) and the main vector of shape variation (first principal component) for each family in the morphospace, we showed that two thirds of the families diversify along an axis of body elongation. Finally, a comparative analysis using a principal coordinate analysis based on the angles among first principal component vectors of each family shape showed that families accomplish changes in elongation with a wide range of underlying modifications. Some groups such as Pomacentridae and Lethrinidae undergo decreases in body depth with proportional increases in all body regions, while other families show disproportionate changes in the length of the head (e.g., Labridae), the trunk or caudal region in all combinations (e.g., Pempheridae and Pinguipedidae). In conclusion, we found that evolutionary changes in body shape along an axis of elongation dominates diversification in reef fishes. Changes in shape on this axis are thought to have immediate implications for swimming performance, defense from gape limited predators, suction feeding performance and access to some highly specialized habitats. The morphological modifications that underlie changes in elongation are highly diverse, suggesting a role for a range of developmental processes and functional consequences.