Signatures of Relaxed Selection in the CYP8B1 Gene of Birds and Mammals

CYP8B1 Catalyzes 12alpha-Hydroxylation of C27 Bile Acid: In Vitro Conversion of Dihydroxycoprostanic Acid into Trihydroxycoprostanic Acid

Sterol 12α-hydroxylase (CYP8B1) is the unique P450 enzyme with sterol 12-oxidation activity, playing an exclusive role in 12α-hydroxylating intermediates along the bile acid (BA) synthesis pathway. Despite the long history of BA metabolism studies, it is unclear whether CYP8B1 catalyzes 12α-hydroxylation of C27 BAs, the key intermediates shuttling between mitochondria and peroxisomes. This work provides robust in vitro evidence that both microsomal and recombinant CYP8B1 enzymes catalyze the 12α-hydroxylation of dihydroxycoprostanic acid (DHCA) into trihydroxycoprostanic acid (THCA). On the one hand, DHCA 12α-hydroxylation reactivity is conservatively detected in liver microsomes of both human and preclinical animals. The reactivity of human tissue fractions conforms well with the selectivity of CYP8B1 mRNA expression, while the contribution of P450 enzymes other than CYP8B1 is excluded by reaction phenotyping in commercial recombinant enzymes. On the other hand, we prepared functional recombinant human CYP8B1 proteins according to a recently published protocol. Titration of the purified CYP8B1 proteins with either C4 (7α-hydroxy-4-cholesten-3-one) or DHCA yields expected blue shifts of the heme Soret peak (type I binding). The recombinant CYP8B1 proteins efficiently catalyze 12α-hydroxylation of both DHCA and C4, with substrate concentration occupying half of the binding sites of 3.0 and 1.9 μM and kcat of 3.2 and 2.6 minutes-1, respectively. In summary, the confirmed role of CYP8B1 in 12α-hydroxylation of C27 BAs has furnished the forgotten passageway in the BA synthesis pathway. The present finding might have opened a new window to consider the biology of CYP8B1 in glucolipid metabolism and to evaluate CYP8B1 inhibition as a therapeutic approach of crucial interest for metabolic diseases. SIGNIFICANCE STATEMENT: The academic community has spent approximately 90 years interpreting the synthesis of bile acids. However, the 12α-hydroxylation of intermediates catalyzed by CYP8B1 is not completely mapped on the classic pathway, particularly for the C27 bile acids, the pivotal intermediates shuttling between mitochondria and peroxisomes. This work discloses the forgotten 12α-hydroxylation pathway from dihydroxycoprostanic acid into trihydroxycoprostanic acid. The present finding may facilitate evaluating CYP8B1 inhibition as a therapeutic approach of crucial interest for metabolic diseases.

Adaptive strategies of aquatic mammals: Exploring the role of the HIF pathway and hypoxia tolerance

Aquatic mammals (marine and freshwater species) share significant and similar adaptations, enabling them to tolerate hypoxia during regular breath-hold diving. Despite the established importance of HIF1A, a master regulator in the molecular mechanism of hypoxia response, and other associated genes, their role in the evolutionary adaptation of aquatic mammals is not fully understood. In this study, we investigated this topic by employing a candidate gene approach to analyze 11 critical genes involved in the HIF1A signaling pathway in aquatic mammals. Our gene analyses included evaluating positive and negative selection, relaxation or constriction of selection, and molecular convergence compared to other terrestrial mammals, including subterranean mammals. Evidence of selection suggested a significant role of negative selection, as well as relaxation of the selective regime in cetaceans for most of these genes. We found that the glutamine 68 variant in the HIF3α protein is unique to cetaceans and initial evaluations indicated a destabilizing effect on protein structure. However, further analyses are necessary to evaluate its functional impact and adaptive relevance in this taxon.

Pseudogenization of the Hair-Related Genes PADI3 and S100A3 in Cetaceans and Hippopotamus amphibius

Hair-related genes in mammals play important roles in the development and maintenance of hair and other keratinous structures in mammals. The peptidyl arginine deiminase 3 (PADI3) gene encodes an enzyme that catalyzes the conversion of arginine residues to citrulline. The S100 calcium binding protein A3 (S100A3) gene encodes a protein that is highly expressed in the hair cuticle and contains arginine residues that are converted to citrullines by PADI enzymes. In this study, we investigated the pseudogenization events of PADI3 and S100A3 in cetaceans and Hippopotamus amphibius. We found that PADI3 underwent three independent pseudogenization events during cetacean evolution, in baleen whales, toothed cetaceans other than Physeter catodon, and P. catodon. Notably, the entire PADI3 gene is absent in the baleen whales. Pseudogenization of S100A3 occurred independently in cetaceans and H. amphibius. Interestingly, we found that in cetaceans S100A3 underwent pseudogenization before PADI3, suggesting that differential selection pressures were acting on the two genes. Our findings provide valuable insights into the molecular evolution of these genes in cetaceans and hippopotamuses, highlighting their importance for understanding the evolution of hair-related genes.

Transcriptome and Gut Microbiota Profiling Analysis of ANIT-Induced Cholestasis and the Effects of Da-Huang-Xiao-Shi Decoction Intervention

Cholestasis is characterized by bile acid (BA) circulation disorders, which is usually related to damage of hepatocyte barrier function. Currently, patients with cholestasis face several obstacles in seeking diagnosis and therapy. Da-Huang-Xiao-Shi decoction (DHXSD) is an ancient classic formula that has been used clinically for cholestasis treatment. Nevertheless, the underlying biological activities and therapeutic mechanisms remain unclear. In this study, an alpha-naphthylisothiocyanate (ANIT)-induced cholestasis rat model was established to examine the anticholestatic effects of DHXSD using histopathological and molecular analyses. Transcriptomic analysis combined with 16S rRNA gene sequencing analysis was systematically applied to study the mechanism of action of DHXSD. Simultaneously, the effect of DHXSD on gut microbiota, short-chain fatty acids (SCFAs), and intestinal barrier function were evaluated based on the ANIT-induced cholestasis model in rats. The results showed that DHXSD effectively attenuated ANIT-induced cholestasis by reducing liver function indicators (alanine transaminase [ALT], P < 0.05; alkaline phosphatase [ALP], P < 0.05; total bile acid [TBA], P < 0.01; γ-glutamyl transpeptidase [GGT], P < 0.001) and levels of hepatotoxicity-related enzymes (P < 0.05), thus improving the recovery of histopathological injuries, and regulating levels of inflammatory cytokines (P < 0.05). In addition, 16S rRNA gene sequencing analysis combined with intestinal barrier function analysis revealed that the DHXSD significantly ameliorated ANIT-induced gut microbiota dysbiosis. Significantly altered genes in the model and treatment groups were screened using transcriptomic analysis. Sixty-eight genes and four microbial genera were simultaneously altered with opposing trends in variation after ANIT and DHXSD treatments. We built a framework for predicting targets and host-microbe interaction mechanisms, as well as identifying alternative treatment for cholestasis, which should be validated further for clinical application. In conclusion, DHXSD appears to be a promising agent for protection against liver injury. IMPORTANCE Cholestasis is a serious manifestation of liver diseases resulting in liver injury, fibrosis, and liver failure with limited therapies. To date, only ursodeoxycholic acid (UDCA) has been approved by the U.S. Food and Drug Administration for the treatment of cholestasis. However, approximately one-third of patients with cholestasis are unresponsive to UDCA. Therefore, it is urgent to search for appropriate therapeutic agents for restoring stoppage status of the bile components to treat cholestasis. In this study, we investigated how the microbiome and transcriptome data sets correlated with each other to clarify the role of microbiome alterations in host metabolism. In combination, this research offers potential molecular biomarkers that should be validated for more accurate diagnosis of cholestasis and the clinical utilisation of gut microbiota as a target for treatment.

Molecular Regulation of Bile Acid Homeostasis

Bile acids have been known for decades to aid in the digestion and absorption of dietary fats and fat-soluble vitamins in the intestine. The development of gene knockout mice models and transgenic humanized mouse models have helped us understand other functions of bile acids, such as their role in modulating fat, glucose, and energy metabolism, and in the molecular regulation of the synthesis, transport, and homeostasis of bile acids. The G-protein coupled receptor TGR5 regulates the bile acid induced alterations of intermediary metabolism, whereas the nuclear receptor FXR regulates bile acid synthesis and homeostasis. However, this review indicates that unidentified factors in addition to FXR must exist to aid in the regulation of bile acid synthesis and homeostasis. SIGNIFICANCE STATEMENT: This review captures the present understanding of bile acid synthesis, the role of bile acid transporters in the enterohepatic circulation of bile acids, the role of the nuclear receptor FXR on the regulation of bile acid synthesis and bile acid transporters, and the importance of bile acids in activating GPCR signaling via TGR5 to modify intermediary metabolism. This information is useful for developing drugs for the treatment of various hepatic and intestinal diseases, as well as the metabolic syndrome.

Toxicological Implications of Biological Heterogeneity

From a micro to macro scale of biological organization, macromolecular diversity and biological heterogeneity are fundamental properties of biological systems. Heterogeneity may result from genetic, epigenetic, and non-genetic characteristics (e.g., tissue microenvironment). Macromolecular diversity and biological heterogeneity are tolerated as long as the sustenance and propagation of life are not disrupted. They also provide the raw materials for microevolutionary changes that may help organisms adapt to new selection pressures arising from the environment. Sequence evolution, functional divergence, and positive selection of gene and promoter dosage play a major role in the evolution of life's diversity including complex metabolic networks, which is ultimately reflected in changes in the allele frequency over time. Robustness in evolvable biological systems is conferred by functional redundancy that is often created by macromolecular diversity and biological heterogeneity. The ability to investigate biological macromolecules at an increasingly finer level has uncovered a wealth of information in this regard. Therefore, the dynamics of biological complexity should be taken into consideration in biomedical research.

Recurrent erosion of COA1/MITRAC15 exemplifies conditional gene dispensability in oxidative phosphorylation

Skeletal muscle fibers rely upon either oxidative phosphorylation or the glycolytic pathway with much less reliance on oxidative phosphorylation to achieve muscular contractions that power mechanical movements. Species with energy-intensive adaptive traits that require sudden bursts of energy have a greater dependency on glycolytic fibers. Glycolytic fibers have decreased reliance on OXPHOS and lower mitochondrial content compared to oxidative fibers. Hence, we hypothesized that gene loss might have occurred within the OXPHOS pathway in lineages that largely depend on glycolytic fibers. The protein encoded by the COA1/MITRAC15 gene with conserved orthologs found in budding yeast to humans promotes mitochondrial translation. We show that gene disrupting mutations have accumulated within the COA1 gene in the cheetah, several species of galliform birds, and rodents. The genomic region containing COA1 is a well-established evolutionary breakpoint region in mammals. Careful inspection of genome assemblies of closely related species of rodents and marsupials suggests two independent COA1 gene loss events co-occurring with chromosomal rearrangements. Besides recurrent gene loss events, we document changes in COA1 exon structure in primates and felids. The detailed evolutionary history presented in this study reveals the intricate link between skeletal muscle fiber composition and the occasional dispensability of the chaperone-like role of the COA1 gene.

Coelacanth SERINC2 Inhibits HIV-1 Infectivity and Is Counteracted by Envelope Glycoprotein from Foamy Virus

SERINC5 restricts nef-defective HIV-1 by affecting early steps of the virus life cycle. Distantly related retroviruses with a wide host range encode virulent factors in response to challenge by SERINC5. However, the evolutionary origins of this antiretroviral activity, its prevalence among the paralogs, and its ability to target retroviruses remain understudied. In agreement with previous studies, we found that four human SERINC paralogs inhibit nef-defective HIV-1, with SERINC2 being an exception. Here, we demonstrate that this lack of activity in human SERINC2 is associated with its post-whole-genome duplication (post-WGD) divergence, as evidenced by the ability of pre-WGD orthologs from Saccharomyces cerevisiae and flies and a post-WGD-proximate SERINC2 from coelacanths to inhibit the virus. Intriguingly, Nef is unable to counter coelacanth SERINC2, indicating that such activity was directed toward other retroviruses found in coelacanths (like foamy viruses). However, foamy virus-derived vectors are intrinsically resistant to the action of SERINC2, and we show that the foamy virus envelope confers this resistance by affecting its steady-state levels. Our study highlights an ancient origin of antiretroviral activity in SERINCs and a hitherto-unknown interaction with a foamy virus.

IMPORTANCESERINC5 constitutes a critical barrier to the propagation of retroviruses, as highlighted by parallel emergence of anti-SERINC5 activities among distant retroviral lineages. Therefore, understanding the origin and evolution of these host factors will provide key information about virus-host relationships that can be exploited for future drug development. Here, we show that SERINC5-mediated nef-defective HIV-1 infection inhibition is evolutionarily conserved. SERINC2 from coelacanth restricts HIV-1, and it was functionally adapted to target foamy viruses. Our findings provide insights into the evolutionary origin of antiretroviral activity in the SERINC gene family and uncover the role of SERINCs in shaping the long-term conflicts between retroviruses and their hosts.

Evidence for the loss of plasminogen receptor KT gene in chicken

The loss of conserved genes has the potential to alter phenotypes drastically. Screening of vertebrate genomes for lineage-specific gene loss events has identified numerous natural knockouts associated with specific phenotypes. We provide evidence for the loss of a multi-exonic plasminogen receptor KT (PLGRKT) protein-encoding gene located on the Z chromosome in chicken. Exons 1 and 2 are entirely missing; remnants of exon 3 and a mostly intact exon 4 are identified in an assembly gap-free region in chicken with conserved synteny across species and verified using transcriptome and genome sequencing. PLGRKT gene disrupting changes are present in representative species from all five galliform families. In contrast to this, the presence of an intact transcriptionally active PLGRKT gene in species such as mallard, swan goose, and Anolis lizard suggests that gene loss occurred in the galliform lineage sometime between 68 and 80 Mya. The presence of galliform specific chicken repeat 1 (CR1) insertion at the erstwhile exon 2 of PLGRKT gene suggests repeat insertion-mediated loss. However, at least nine other independent PLGRKT coding frame disrupting changes in other bird species are supported by genome sequencing and indicate a role for relaxed purifying selection before CR1 insertion. The recurrent loss of a conserved gene with a role in the regulation of macrophage migration, efferocytosis, and blood coagulation is intriguing. Hence, we propose potential candidate genes that might be compensating the function of PLGRKT based on the presence of a C-terminal lysine residue, transmembrane domains, and gene expression patterns.

Convergent Cortistatin losses parallel modifications in circadian rhythmicity and energy homeostasis in Cetacea and other mammalian lineages

The ancestors of Cetacea underwent profound morpho-physiological alterations. By displaying an exclusive aquatic existence, cetaceans evolved unique patterns of locomotor activity, vigilant behaviour, thermoregulation and circadian rhythmicity. Deciphering the molecular landscape governing many of these adaptations is key to understand the evolution of phenotypes. Here, we investigate Cortistatin (CORT), a neuropeptide displaying an important role in mammalian biorhythm regulation. This neuropeptide is a known neuroendocrine factor, stimulating slow-wave sleep, but also involved in the regulation of energy metabolism and hypomotility inducement. We assessed the functional status of CORT in 359 mammalian genomes (25 orders), including 30 species of Cetacea. Our findings indicate that cetaceans and other mammals with atypical biorhythms, thermal constraints and/or energy metabolism, have accumulated deleterious mutations in CORT. In light of the pleiotropic action of this neuropeptide, we suggest that this inactivation contributed to a plethora of phenotypic adjustments to accommodate adaptive solutions to specific ecological niches.

Loss of inner kinetochore genes is associated with the transition to an unconventional point centromere in budding yeast

Background

The genomic sequences of centromeres, as well as the set of proteins that recognize and interact with centromeres, are known to quickly diverge between lineages potentially contributing to post-zygotic reproductive isolation. However, the actual sequence of events and processes involved in the divergence of the kinetochore machinery is not known. The patterns of gene loss that occur during evolution concomitant with phenotypic changes have been used to understand the timing and order of molecular changes.

Methods

I screened the high-quality genomes of twenty budding yeast species for the presence of well-studied kinetochore genes. Based on the conserved gene order and complete genome assemblies, I identified gene loss events. Subsequently, I searched the intergenic regions to identify any un-annotated genes or gene remnants to obtain additional evidence of gene loss.

Results

My analysis identified the loss of four genes (NKP1, NKP2, CENPL/IML3 and CENPN/CHL4) of the inner kinetochore constitutive centromere-associated network (CCAN/also known as CTF19 complex in yeast) in both the Naumovozyma species for which genome assemblies are available. Surprisingly, this collective loss of four genes of the CCAN/CTF19 complex coincides with the emergence of unconventional centromeres in N. castellii and N. dairenensis. My study suggests a tentative link between the emergence of unconventional point centromeres and the turnover of kinetochore genes in budding yeast.

Parallel selection on ecologically relevant gene functions in the transcriptomes of highly diversifying salmonids

Background

Salmonid fishes are characterised by a very high level of variation in trophic, ecological, physiological, and life history adaptations. Some salmonid taxa show exceptional potential for fast, within-lake diversification into morphologically and ecologically distinct variants, often in parallel; these are the lake-resident charr and whitefish (several species in the genera Salvelinus and Coregonus). To identify selection on genes and gene categories associated with such predictable diversifications, we analysed 2702 orthogroups (4.82 Mbp total; average 4.77 genes/orthogroup; average 1783 bp/orthogroup). We did so in two charr and two whitefish species and compared to five other salmonid lineages, which do not evolve in such ecologically predictable ways, and one non-salmonid outgroup.

Results

All selection analyses are based on Coregonus and Salvelinus compared to non-diversifying taxa. We found more orthogroups were affected by relaxed selection than intensified selection. Of those, 122 were under significant relaxed selection, with trends of an overrepresentation of serine family amino acid metabolism and transcriptional regulation, and significant enrichment of behaviour-associated gene functions. Seventy-eight orthogroups were under significant intensified selection and were enriched for signalling process and transcriptional regulation gene ontology terms and actin filament and lipid metabolism gene sets. Ninety-two orthogroups were under diversifying/positive selection. These were enriched for signal transduction, transmembrane transport, and pyruvate metabolism gene ontology terms and often contained genes involved in transcriptional regulation and development. Several orthogroups showed signs of multiple types of selection. For example, orthogroups under relaxed and diversifying selection contained genes such as ap1m2, involved in immunity and development, and slc6a8, playing an important role in muscle and brain creatine uptake. Orthogroups under intensified and diversifying selection were also found, such as genes syn3, with a role in neural processes, and ctsk, involved in bone remodelling.

Conclusions

Our approach pinpointed relevant genomic targets by distinguishing among different kinds of selection. We found that relaxed, intensified, and diversifying selection affect orthogroups and gene functions of ecological relevance in salmonids. Because they were found consistently and robustly across charr and whitefish and not other salmonid lineages, we propose these genes have a potential role in the replicated ecological diversifications.