An improved and general streamlined phylogenetic protocol applied to the fatty acid desaturase family

Phylogenomics of the Ecdysteroid Kinase-like (EcKL) Gene Family in Insects Highlights Roles in Both Steroid Hormone Metabolism and Detoxification

The evolutionary dynamics of large gene families can offer important insights into the functions of their individual members. While the ecdysteroid kinase-like (EcKL) gene family has previously been linked to the metabolism of both steroid molting hormones and xenobiotic toxins, the functions of nearly all EcKL genes are unknown, and there is little information on their evolution across all insects. Here, we perform comprehensive phylogenetic analyses on a manually annotated set of EcKL genes from 140 insect genomes, revealing the gene family is comprised of at least 13 subfamilies that differ in retention and stability. Our results show the only two genes known to encode ecdysteroid kinases belong to different subfamilies and therefore ecdysteroid metabolism functions must be spread throughout the EcKL family. We provide comparative phylogenomic evidence that EcKLs are involved in detoxification across insects, with positive associations between family size and dietary chemical complexity, and we also find similar evidence for the cytochrome P450 and glutathione S-transferase gene families. Unexpectedly, we find that the size of the clade containing a known ecdysteroid kinase is positively associated with host plant taxonomic diversity in Lepidoptera, possibly suggesting multiple functional shifts between hormone and xenobiotic metabolism. Our evolutionary analyses provide hypotheses of function and a robust framework for future experimental studies of the EcKL gene family. They also open promising new avenues for exploring the genomic basis of dietary adaptation in insects, including the classically studied coevolution of butterflies with their host plants.

A novel Tetrahymena thermophila sterol C-22 desaturase belongs to the Fatty Acid Hydroxylase/Desaturase superfamily

Sterols in eukaryotic cells play important roles in modulating membrane fluidity and in cell signaling and trafficking. During evolution, a combination of gene losses and acquisitions gave rise to an extraordinary diversity of sterols in different organisms. The sterol C-22 desaturase identified in plants and fungi as a cytochrome P-450 monooxygenase evolved from the first eukaryotic cytochrome P450 and was lost in many lineages. Although the ciliate Tetrahymena thermophila desaturates sterols at the C-22 position, no cytochrome P-450 orthologs are present in the genome. Here, we aim to identify the genes responsible for the desaturation as well as their probable origin. We used gene knockout and yeast heterologous expression approaches to identify two putative genes, retrieved from a previous transcriptomic analysis, as sterol C-22 desaturases. Furthermore, we demonstrate using bioinformatics and evolutionary analyses that both genes encode a novel type of sterol C-22 desaturase that belongs to the large fatty acid hydroxylase/desaturase superfamily and the genes originated by genetic duplication prior to functional diversification. These results stress the widespread existence of nonhomologous isofunctional enzymes among different lineages of the tree of life as well as the suitability for the use of T. thermophila as a valuable model to investigate the evolutionary process of large enzyme families.

A non-canonical Δ9-desaturase synthesizing palmitoleic acid identified in the thraustochytrid Aurantiochytrium sp. T66

Abstract

Thraustochytrids are oleaginous marine eukaryotic microbes currently used to produce the essential omega-3 fatty acid docosahexaenoic acid (DHA, C22:6 n-3). To improve the production of this essential fatty acid by strain engineering, it is important to deeply understand how thraustochytrids synthesize fatty acids. While DHA is synthesized by a dedicated enzyme complex, other fatty acids are probably synthesized by the fatty acid synthase, followed by desaturases and elongases. Which unsaturated fatty acids are produced differs between different thraustochytrid genera and species; for example, Aurantiochytrium sp. T66, but not Aurantiochytrium limacinum SR21, synthesizes palmitoleic acid (C16:1 n-7) and vaccenic acid (C18:1 n-7). How strain T66 can produce these fatty acids has not been known, because BLAST analyses suggest that strain T66 does not encode any Δ9-desaturase-like enzyme. However, it does encode one Δ12-desaturase-like enzyme. In this study, the latter enzyme was expressed in A. limacinum SR21, and both C16:1 n-7 and C18:1 n-7 could be detected in the transgenic cells. Our results show that this desaturase, annotated T66Des9, is a Δ9-desaturase accepting C16:0 as a substrate. Phylogenetic studies indicate that the corresponding gene probably has evolved from a Δ12-desaturase-encoding gene. This possibility has not been reported earlier and is important to consider when one tries to deduce the potential a given organism has for producing unsaturated fatty acids based on its genome sequence alone.

Key points

• In thraustochytrids, automatic gene annotation does not always explain the fatty acids produced.

• T66Des9 is shown to synthesize palmitoleic acid (C16:1 n-7).

• T66des9 has probably evolved from Δ12-desaturase-encoding genes.

Supplementary Information

The online version contains supplementary material available at 10.1007/s00253-021-11425-5.

A new phylogenetic protocol: dealing with model misspecification and confirmation bias in molecular phylogenetics

Molecular phylogenetics plays a key role in comparative genomics and has increasingly significant impacts on science, industry, government, public health and society. In this paper, we posit that the current phylogenetic protocol is missing two critical steps, and that their absence allows model misspecification and confirmation bias to unduly influence phylogenetic estimates. Based on the potential offered by well-established but under-used procedures, such as assessment of phylogenetic assumptions and tests of goodness of fit, we introduce a new phylogenetic protocol that will reduce confirmation bias and increase the accuracy of phylogenetic estimates.

Gene identification and functional characterization of a Δ12 fatty acid desaturase in Tetrahymena thermophila and its influence in homeoviscous adaptation to low temperature

Homeoviscous adaptation in poikilotherms is based in the regulation of the level of desaturation of fatty acids, variation in phospholipids head groups and sterol content in the membrane lipids, in order to maintain the membrane fluidity in response to changes in environmental temperature. Increased proportion of unsaturated fatty acids is thought to be the main response to low-temperature acclimation, which is mostly achieved by fatty acid desaturases. Genome analysis of the ciliate Tetrahymena thermophila and a gene knockout approach has allowed us to identify one Δ12 FAD and to study its activity in the original host and in a yeast heterologous expression system. The "PUFA index" -relative content of polyunsaturated fatty acids compared to the sum of saturated and monounsaturated fatty acid content- was ~57% lower at 15 °C and 35 °C in the Δ12 FAD gene knockout strain (KOΔ12) compared to WT strain. We characterized the role of T. thermophila Δ12 FAD on homeoviscous adaptation and analyzed its involvement in cellular growth, cold stress response, and membrane fluidity, as well as its expression pattern during temperature shifts. Although these alterations allowed normal growth in the KOΔ12 strain at 30 °C or higher temperatures, growth was impaired at temperatures of 20 °C or lower, where homeoviscous adaptation is impaired. These results stress the importance of Δ12 FAD in the regulation of cold adaptation processes, as well as the suitability of T. thermophila as a valuable model to investigate the regulation of membrane lipids and evolutionary conservation and divergence of the underlying mechanisms.