Selection on different genes with equivalent functions: the convergence story told by Hox genes along the evolution of aquatic mammalian lineages

HOXA9 versus HOXB9; particular focus on their controversial role in tumor pathogenesis

The Homeobox (HOX) gene family is essential to regulating cellular processes because it maintains the exact coordination required for tissue homeostasis, cellular differentiation, and embryonic development. The most distinctive feature of this class of genes is the presence of the highly conserved DNA region known as the homeobox, which is essential for controlling their regulatory activities. Important players in the intricate process of genetic regulation are the HOX genes. Many diseases, especially in the area of cancer, are linked to their aberrant functioning. Due to their distinctive functions in biomedical research-particularly in the complex process of tumor advancement-HOXA9 and HOXB9 have drawn particular attention. HOXA9 and HOXB9 are more significant than what is usually connected with HOX genes since they have roles in the intricate field of cancer and beyond embryonic processes. The framework for a focused study of the different effects of HOXA9 and HOXB9 in the context of tumor biology is established in this study.

Repatterning of mammalian backbone regionalization in cetaceans

Cetacean reinvasion of the aquatic realm is an iconic ecological transition that led to drastic modifications of the mammalian body plan, especially in the axial skeleton. Relative to the vertebral column of other mammals that is subdivided into numerous anatomical regions, regional boundaries of the cetacean backbone appear obscured. Whether the traditional mammalian regions are present in cetaceans but hard to detect due to anatomical homogenization or if regions have been entirely repatterned remains unresolved. Here we combine a segmented linear regression approach with spectral clustering to quantitatively investigate the number, position, and homology of vertebral regions across 62 species from all major cetacean clades. We propose the Nested Regions hypothesis under which the cetacean backbone is composed of six homologous modules subdivided into six to nine post-cervical regions, with the degree of regionalization dependent on vertebral count and ecology. Compared to terrestrial mammals, the cetacean backbone is less regionalized in the precaudal segment but more regionalized in the caudal segment, indicating repatterning of the vertebral column associated with the transition from limb-powered to axial-driven locomotion.

Natural selection and convergent evolution of the HOX gene family in Carnivora

HOX genes play a central role in the development and regulation of limb patterns. For mammals in the order Carnivora, limbs have evolved in different forms, and there are interesting cases of phenotypic convergence, such as the pseudothumb of the giant and red pandas, and the flippers or specialized limbs of the pinnipeds and sea otter. However, the molecular bases of limb development remain largely unclear. Here, we studied the molecular evolution of the HOX9 ~ 13 genes of 14 representative species in Carnivora and explored the molecular evolution of other HOX genes. We found that only one limb development gene, HOXC10, underwent convergent evolution between giant and red pandas and was thus an important candidate gene related to the development of pseudothumbs. No signals of amino acid convergence and natural selection were found in HOX9 ~ 13 genes between pinnipeds and sea otter, but there was evidence of positive selection and rapid evolution in four pinniped species. Overall, few HOX genes evolve via natural selection or convergent evolution, and these could be important candidate genes for further functional validation. Our findings provide insights into potential molecular mechanisms of the development of specialized pseudothumbs and flippers (or specialized limbs).

Juxtapose: a gene-embedding approach for comparing co-expression networks

BACKGROUND

Gene co-expression networks (GCNs) are not easily comparable due to their complex structure. In this paper, we propose a tool, Juxtapose, together with similarity measures that can be utilized for comparative transcriptomics between a set of organisms. While we focus on its application to comparing co-expression networks across species in evolutionary studies, Juxtapose is also generalizable to co-expression network comparisons across tissues or conditions within the same species.

METHODS

A word embedding strategy commonly used in natural language processing was utilized in order to generate gene embeddings based on walks made throughout the GCNs. Juxtapose was evaluated based on its ability to embed the nodes of synthetic structures in the networks consistently while also generating biologically informative results. Evaluation of the techniques proposed in this research utilized RNA-seq datasets from GTEx, a multi-species experiment of prefrontal cortex samples from the Gene Expression Omnibus, as well as synthesized datasets. Biological evaluation was performed using gene set enrichment analysis and known gene relationships in literature.

RESULTS

We show that Juxtapose is capable of globally aligning synthesized networks as well as identifying areas that are conserved in real gene co-expression networks without reliance on external biological information. Furthermore, output from a matching algorithm that uses cosine distance between GCN embeddings is shown to be an informative measure of similarity that reflects the amount of topological similarity between networks.

CONCLUSIONS

Juxtapose can be used to align GCNs without relying on known biological similarities and enables post-hoc analyses using biological parameters, such as orthology of genes, or conserved or variable pathways.

AVAILABILITY

A development version of the software used in this paper is available at https://github.com/klovens/juxtapose.

Oxidative stress drives divergent evolution of the glutathione peroxidase (GPX) gene family in mammals

The molecular basis for adaptations to extreme environments can now be understood by interrogating the ever-increasing number of sequenced genomes. Mammals such as cetaceans, bats, and highland species can protect themselves from oxidative stress, a disruption in the balance of reactive oxygen species, which results in oxidative injury and cell damage. Here, we consider the evolution of the glutathione peroxidase (GPX) family of antioxidant enzymes by interrogating publicly available genome data from 70 mammalian species from all major clades. We identified 8 GPX subclasses ubiquitous to all mammalian groups. Mammalian GPX gene families resolved into the GPX4/7/8 and GPX1/2/3/5/6 groups and are characterized by several instances of gene duplication and loss, indicating a dynamic process of gene birth and death in mammals. Seven of the eight GPX subfamilies (all but GPX7) were under positive selection, with the residues under selection located at or close to active sites or at the dimer interface. We also reveal evidence of a correlation between ecological niches (e.g. high oxidative stress) and the divergent selection and gene copy number of GPX subclasses. Notably, a convergent expansion of GPX1 was observed in several independent lineages of mammals under oxidative stress and may be important for avoiding oxidative damage. Collectively, this study suggests that the GPX gene family has shaped the adaption of mammals to stressful environments.

Functional Studies with Primary Cells Provide a System for Genome-to-Phenome Investigations in Marine Mammals

Marine mammals exhibit some of the most dramatic physiological adaptations in their clade and offer unparalleled insights into the mechanisms driving convergent evolution on relatively short time scales. Some of these adaptations, such as extreme tolerance to hypoxia and prolonged food deprivation, are uncommon among most terrestrial mammals and challenge established metabolic principles of supply and demand balance. Non-targeted omics studies are starting to uncover the genetic foundations of such adaptations, but tools for testing functional significance in these animals are currently lacking. Cellular modeling with primary cells represents a powerful approach for elucidating the molecular etiology of physiological adaptation, a critical step in accelerating genome-to-phenome studies in organisms in which transgenesis is impossible (e.g., large-bodied, long-lived, fully aquatic, federally protected species). Gene perturbation studies in primary cells can directly evaluate whether specific mutations, gene loss, or duplication confer functional advantages such as hypoxia or stress tolerance in marine mammals. Here, we summarize how genetic and pharmacological manipulation approaches in primary cells have advanced mechanistic investigations in other non-traditional mammalian species, and highlight the need for such investigations in marine mammals. We also provide key considerations for isolating, culturing, and conducting experiments with marine mammal cells under conditions that mimic in vivo states. We propose that primary cell culture is a critical tool for conducting functional mechanistic studies (e.g., gene knockdown, over-expression, or editing) that can provide the missing link between genome- and organismal-level understanding of physiological adaptations in marine mammals.

Convergent Evolution at the Pathway Level: Predictable Regulatory Changes during Flower Color Transitions

The predictability of evolution, or whether lineages repeatedly follow the same evolutionary trajectories during phenotypic convergence remains an open question of evolutionary biology. In this study, we investigate evolutionary convergence at the biochemical pathway level and test the predictability of evolution using floral anthocyanin pigmentation, a trait with a well-understood genetic and regulatory basis. We reconstructed the evolution of floral anthocyanin content across 28 species of the Andean clade Iochrominae (Solanaceae) and investigated how shifts in pigmentation are related to changes in expression of seven key anthocyanin pathway genes. We used phylogenetic multivariate analysis of gene expression to test for phenotypic and developmental convergence at a macroevolutionary scale. Our results show that the four independent losses of the ancestral pigment delphinidin involved convergent losses of expression of the three late pathway genes (F3'5'h, Dfr, and Ans). Transitions between pigment types affecting floral hue (e.g., blue to red) involve changes to the expression of branching genes F3'h and F3'5'h, while the expression levels of early steps of the pathway are strongly conserved in all species. These patterns support the idea that the macroevolution of floral pigmentation follows predictable evolutionary trajectories to reach convergent phenotype space, repeatedly involving regulatory changes. This is likely driven by constraints at the pathway level, such as pleiotropy and regulatory structure.

Accelerated Evolution of Limb-Related Gene Hoxd11 in the Common Ancestor of Cetaceans and Ruminants (Cetruminantia)

Reduced numbers of carpal and tarsal bones (wrist and ankle joints) are extensively observed in the clade of Cetacea and Ruminantia (Cetruminantia). Homebox D11 (Hoxd11) is one of the important genes required for limb development in mammals. Mutations in Hoxd11 can lead to defects in particular bones of limbs, including carpus and tarsus. To test whether evolutionary changes in Hoxd11 underlie the loss of these bones in Cetruminantia, we sequenced and analyzed Hoxd11 coding sequences and compared them with other 5' HoxA and HoxD genes in a taxonomic coverage of Cetacea, Ruminantia and other mammalian relatives. Statistical tests on the Hoxd11 sequences found an accelerated evolution in the common ancestor of cetaceans and ruminants, which coincided with the reduction of carpal and tarsal bones in this clade. Five amino acid substitutions (G222S, G227A, G229S, A240T and G261V) and one amino acid deletion (G254Del) occurred in this lineage. In contrast, other 5' HoxA and HoxD genes do not show this same evolutionary pattern, but instead display a highly conserved pattern of evolution in this lineage. Accelerated evolution of Hoxd11, but not other 5' HoxA and HoxD genes, is probably related to the reduction of the carpal and tarsal bones in Cetruminantia. Moreover, we found two amino acid substitutions (G110S and D223N) in Hoxd11 that are unique to the lineage of Cetacea, which coincided with hindlimb loss in the common ancestor of cetaceans. Our results give molecular evidence of Hoxd11 adaptive evolution in cetaceans and ruminants, which could be correlated with limb morphological adaptation.

Signatures of selection and environmental adaptation across the goat genome post-domestication

BACKGROUND

Since goat was domesticated 10,000 years ago, many factors have contributed to the differentiation of goat breeds and these are classified mainly into two types: (i) adaptation to different breeding systems and/or purposes and (ii) adaptation to different environments. As a result, approximately 600 goat breeds have developed worldwide; they differ considerably from one another in terms of phenotypic characteristics and are adapted to a wide range of climatic conditions. In this work, we analyzed the AdaptMap goat dataset, which is composed of data from more than 3000 animals collected worldwide and genotyped with the CaprineSNP50 BeadChip. These animals were partitioned into groups based on geographical area, production uses, available records on solid coat color and environmental variables including the sampling geographical coordinates, to investigate the role of natural and/or artificial selection in shaping the genome of goat breeds.

RESULTS

Several signatures of selection on different chromosomal regions were detected across the different breeds, sub-geographical clusters, phenotypic and climatic groups. These regions contain genes that are involved in important biological processes, such as milk-, meat- or fiber-related production, coat color, glucose pathway, oxidative stress response, size, and circadian clock differences. Our results confirm previous findings in other species on adaptation to extreme environments and human purposes and provide new genes that could explain some of the differences between goat breeds according to their geographical distribution and adaptation to different environments.

CONCLUSIONS

These analyses of signatures of selection provide a comprehensive first picture of the global domestication process and adaptation of goat breeds and highlight possible genes that may have contributed to the differentiation of this species worldwide.

Deciphering the evolutionary signatures of pinnipeds using novel genome sequences: The first genomes of Phoca largha, Callorhinus ursinus, and Eumetopias jubatus

The pinnipeds, which comprise seals, sea lions, and walruses, are a remarkable group of marine animals with unique adaptations to semi-aquatic life. However, their genomes are poorly characterized. In this study, we sequenced and characterized the genomes of three pinnipeds (Phoca largha, Callorhinus ursinus, and Eumetopias jubatus), focusing on site-wise sequence changes. We detected rapidly evolving genes in pinniped lineages and substitutions unique to pinnipeds associated with amphibious sound perception. Phenotypic convergence-related sequence convergences are not common in marine mammals. For example, FASN, KCNA5, and IL17RA contain substitutions specific to pinnipeds, yet are potential candidates of phenotypic convergence (blubber, response to hypoxia, and immunity to pathogens) in all marine mammals. The outcomes of this study will provide insight into targets for future studies of convergent evolution or gene function.

Lineage-specific evolution of bitter taste receptor genes in the giant and red pandas implies dietary adaptation

Taste 2 receptors (TAS2R) mediate bitterness perception in mammals, thus are called bitter taste receptors. It is believed that these genes evolved in response to species‐specific diets. The giant panda (Ailuropoda melanoleuca) and red panda (Ailurus fulgens styani) in the order Carnivora are specialized herbivores with an almost exclusive bamboo diet (>90% bamboo). Because bamboo is full of bitter tasting compounds, we hypothesized that adaptive evolution has occurred at TAS2R genes in giant and red pandas throughout the course of their dietary shift. Here, we characterized 195 TAS2R genes in 9 Carnivora species and examined selective pressures on these genes. We found that both pandas harbor more putative functional TAS2R genes than other carnivores, and pseudogenized TAS2R genes in the giant panda are different from the red panda. The purifying selection on TAS2R1, TAS2R9 and TAS2R38 in the giant panda, and TAS2R62 in the red panda, has been strengthened throughout the course of adaptation to bamboo diet, while selective constraint on TAS2R4 and TAS2R38 in the red panda is relaxed. Remarkably, a few positively selected sites on TAS2R42 have been specifically detected in the giant panda. These results suggest an adaptive response in both pandas to a dietary shift from carnivory to herbivory, and TAS2R genes evolved independently in the 2 pandas. Our findings provide new insight into the molecular basis of mammalian sensory evolution and the process of adaptation to new ecological niches.

Evolutionary changes of Hox genes and relevant regulatory factors provide novel insights into mammalian morphological modifications

The diversity of body plans of mammals accelerates the innovation of lifestyles and the extensive adaptation to different habitats, including terrestrial, aerial and aquatic habitats. However, the genetic basis of those phenotypic modifications, which have occurred during mammalian evolution, remains poorly explored. In the present study, we synthetically surveyed the evolutionary pattern of Hox clusters that played a powerful role in the morphogenesis along the head–tail axis of animal embryos and the main regulatory factors (Mll, Bmi1 and E2f6) that control the expression of Hox genes. A deflected density of repetitive elements and lineage‐specific radical mutations of Mll have been determined in marine mammals with morphological changes, suggesting that evolutionary changes may alter Hox gene expression in these lineages, leading to the morphological modification of these lineages. Although no positive selection was detected at certain ancestor nodes of lineages, the increased ω values of Hox genes implied the relaxation of functional constraints of these genes during the mammalian evolutionary process. More importantly, 49 positively‐selected sites were identified in mammalian lineages with phenotypic modifications, indicating adaptive evolution acting on Hox genes and regulatory factors. In addition, 3 parallel amino acid substitutions in some Hox genes were examined in marine mammals, which might be responsible for their streamlined body.