Single-cell mtDNA heteroplasmy in colorectal cancer.

Human mitochondria can be genetically distinct within the same individual, a phenomenon known as heteroplasmy. In cancer, this phenomenon seems exacerbated, and most mitochondrial mutations seem to be heteroplasmic. How this genetic variation is arranged within and among normal and tumor cells is not well understood. To address this question, here we sequenced single-cell mitochondrial genomes from multiple normal and tumoral locations in four colorectal cancer patients. Our results suggest that single cells, both normal and tumoral, can carry various mitochondrial haplotypes. Remarkably, this intra-cell heteroplasmy can arise before tumor development and be maintained afterward in specific tumoral cell subpopulations. At least in the colorectal patients studied here, the somatic mutations in the single-cells do not seem to have a prominent role in tumorigenesis.

PseudoChecker: an integrated online platform for gene inactivation inference.

The rapid expansion of high-quality genome assemblies, exemplified by ongoing initiatives such as the Genome-10K and i5k, demands novel automated methods to approach comparative genomics. Of these, the study of inactivating mutations in the coding region of genes, or pseudogenization, as a source of evolutionary novelty is mostly overlooked. Thus, to address such evolutionary/genomic events, a systematic, accurate and computationally automated approach is required. Here, we present PseudoChecker, the first integrated online platform for gene inactivation inference. Unlike the few existing methods, our comparative genomics-based approach displays full automation, a built-in graphical user interface and a novel index, PseudoIndex, for an empirical evaluation of the gene coding status. As a multi-platform online service, PseudoChecker simplifies access and usability, allowing a fast identification of disruptive mutations. An analysis of 30 genes previously reported to be eroded in mammals, and 30 viable genes from the same lineages, demonstrated that PseudoChecker was able to correctly infer 97% of loss events and 95% of functional genes, confirming its reliability. PseudoChecker is freely available, without login required, at http://pseudochecker.ciimar.up.pt.

PhyloToL: A Taxon/Gene-Rich Phylogenomic Pipeline to Explore Genome Evolution of Diverse Eukaryotes.

Estimating multiple sequence alignments (MSAs) and inferring phylogenies are essential for many aspects of comparative biology. Yet, many bioinformatics tools for such analyses have focused on specific clades, with greatest attention paid to plants, animals, and fungi. The rapid increase in high-throughput sequencing (HTS) data from diverse lineages now provides opportunities to estimate evolutionary relationships and gene family evolution across the eukaryotic tree of life. At the same time, these types of data are known to be error-prone (e.g., substitutions, contamination). To address these opportunities and challenges, we have refined a phylogenomic pipeline, now named PhyloToL, to allow easy incorporation of data from HTS studies, to automate production of both MSAs and gene trees, and to identify and remove contaminants. PhyloToL is designed for phylogenomic analyses of diverse lineages across the tree of life (i.e., at scales of >100 My). We demonstrate the power of PhyloToL by assessing stop codon usage in Ciliophora, identifying contamination in a taxon- and gene-rich database and exploring the evolutionary history of chromosomes in the kinetoplastid parasite Trypanosoma brucei, the causative agent of African sleeping sickness. Benchmarking PhyloToL's homology assessment against that of OrthoMCL and a published paper on superfamilies of bacterial and eukaryotic organellar outer membrane pore-forming proteins demonstrates the power of our approach for determining gene family membership and inferring gene trees. PhyloToL is highly flexible and allows users to easily explore HTS data, test hypotheses about phylogeny and gene family evolution and combine outputs with third-party tools (e.g., PhyloChromoMap, iGTP).

Complete Inactivation of Sebum-Producing Genes Parallels the Loss of Sebaceous Glands in Cetacea.

Genomes are dynamic biological units, with processes of gene duplication and loss triggering evolutionary novelty. The mammalian skin provides a remarkable case study on the occurrence of adaptive morphological innovations. Skin sebaceous glands (SGs), for instance, emerged in the ancestor of mammals serving pivotal roles, such as lubrication, waterproofing, immunity, and thermoregulation, through the secretion of sebum, a complex mixture of various neutral lipids such as triacylglycerol, free fatty acids, wax esters, cholesterol, and squalene. Remarkably, SGs are absent in a few mammalian lineages, including the iconic Cetacea. We investigated the evolution of the key molecular components responsible for skin sebum production: Dgat2l6, Awat1, Awat2, Elovl3, Mogat3, and Fabp9. We show that all analyzed genes have been rendered nonfunctional in Cetacea species (toothed and baleen whales). Transcriptomic analysis, including a novel skin transcriptome from blue whale, supports gene inactivation. The conserved mutational pattern found in most analyzed genes, indicates that pseudogenization events took place prior to the diversification of modern Cetacea lineages. Genome and skin transcriptome analysis of the common hippopotamus highlighted the convergent loss of a subset of sebum-producing genes, notably Awat1 and Mogat3. Partial loss profiles were also detected in non-Cetacea aquatic mammals, such as the Florida manatee, and in terrestrial mammals displaying specialized skin phenotypes such as the African elephant, white rhinoceros and pig. Our findings reveal a unique landscape of "gene vestiges" in the Cetacea sebum-producing compartment, with limited gene loss observed in other mammalian lineages: suggestive of specific adaptations or specializations of skin lipids.

Mesozoic mitogenome rearrangements and freshwater mussel (Bivalvia: Unionoidea) macroevolution.

Using a new fossil-calibrated mitogenome-based approach, we identified macroevolutionary shifts in mitochondrial gene order among the freshwater mussels (Unionoidea). We show that the early Mesozoic divergence of the two Unionoidea clades, Margaritiferidae and Unionidae, was accompanied by a synchronous split in the gene arrangement in the female mitogenome (i.e., gene orders MF1 and UF1). Our results suggest that this macroevolutionary jump was completed within a relatively short time interval (95% HPD 201-226 Ma) that coincided with the Triassic-Jurassic mass extinction. Both gene orders have persisted within these clades for ~200 Ma. The monophyly of the so-called "problematic" Gonideinae taxa was supported by all the inferred phylogenies in this study using, for the first time, the M- and F-type mitogenomes either singly or combined. Within Gonideinae, two additional splits in the gene order (UF1 to UF2, UF2 to UF3) occurred in the Mesozoic and have persisted for ~150 and ~100 Ma, respectively. Finally, the mitogenomic results suggest ancient connections between freshwater basins of East Asia and Europe near the Cretaceous-Paleogene boundary, probably via a continuous paleo-river system or along the Tethys coastal line, which are well supported by at least three independent but almost synchronous divergence events.

Convergent inactivation of the skin-specific C-C motif chemokine ligand 27 in mammalian evolution.

The appearance of mammalian-specific skin features was a key evolutionary event contributing for the elaboration of physiological processes such as thermoregulation, adequate hydration, locomotion, and inflammation. Skin inflammatory and autoimmune processes engage a population of skin-infiltrating T cells expressing a specific C-C chemokine receptor (CCR10) which interacts with an epidermal CC chemokine, the skin-specific C-C motif chemokine ligand 27 (CCL27). CCL27 is selectively produced in the skin by keratinocytes, particularly upon inflammation, mediating the adhesion and homing of skin-infiltrating T cells. Here, we examined the evolution and coding condition of Ccl27 in 112 placental mammalian species. Our findings reveal that a number of open reading frame inactivation events such as insertions, deletions, and start and stop codon mutations independently occurred in Cetacea, Pholidota, Sirenia, Chiroptera, and Rodentia, totalizing 18 species. The diverse habitat settings and lifestyles of Ccl27-eroded lineages probably implied distinct evolutionary triggers rendering this gene unessential. For example, in Cetacea, the rapid renewal of skin layers minimizes the need for an elaborate inflammatory mechanism, mirrored by the absence of epidermal scabs. Our findings suggest that the convergent and independent loss of Ccl27 in mammalian evolution concurred with unique adaptive roads for skin physiology.

Cetacea are natural knockouts for IL20.

The Cetacea infraorder comprises a very unique group within the mammalian lineage. While sharing common ancestors with terrestrial mammals, their exclusive dependence on aquatic environments makes them attractive models to explore the landscape of molecular shifts in radical habitat transitions. Among their diverse anatomical and physiological solutions, we find detectable genetic remodeling of the immune system. In agreement, here we show that the gene sequence of interleukin-20 (IL20) displays unambiguous signs of inactivation with several disruptive mutations, including stop codons, insertions, and a conserved trans-species mutation abolishing a canonical splice site, in nine analyzed cetacean genomes. Considering the suggested role of IL20 in skin immunity processes, including inflammation, epithelization, and remodeling, we propose that gene inactivation follows specific adaptations of cetacean skin to the aquatic environment, in frame with the less-is-more hypothesis.

Expansion and systematics redefinition of the most threatened freshwater mussel family, the Margaritiferidae.

Two Unionida (freshwater mussel) families are present in the Northern Hemisphere; the Margaritiferidae, representing the most threatened of unionid families, and the Unionidae, which include several genera of unresolved taxonomic placement. The recent reassignment of the poorly studied Lamprotula rochechouartii from the Unionidae to the Margaritiferidae motivated a new search for other potential species of margaritiferids from members of Gibbosula and Lamprotula. Based on molecular and morphological analyses conducted on newly collected specimens from Vietnam, we here assign Gibbosula crassa to the Margaritiferidae. Additionally, we reanalyzed all diagnostic characteristics of the Margaritiferidae and examined museum specimens of Lamprotula and Gibbosula. As a result, two additional species are also moved to the Margaritiferidae, i.e. Gibbosula confragosa and Gibbosula polysticta. We performed a robust five marker phylogeny with all available margaritiferid species and discuss the taxonomy within the family. The present phylogeny reveals the division of Margaritiferidae into four ancient clades with distinct morphological, biogeographical and ecological characteristics that justify the division of the Margaritiferidae into two subfamilies (Gibbosulinae and Margaritiferinae) and four genera (Gibbosula, Cumberlandia, Margaritifera, and Pseudunio). The systematics of the Margaritiferidae family is re-defined as well as their distribution, potential origin and main biogeographic patterns.

HR/MS-based lipidome analysis of rat brain modulated by tolcapone.

Lipids play key roles in the body, influencing cellular regulation, function, and signalling. Tolcapone, a potent catechol-O-methyltransferase (COMT) inhibitor described to enhance cognitive performance in healthy subjects, was previously shown to impact fatty acid ß-oxidation and oxidative phosphorylation. However, its impact on the brain lipidome remains unexplored. Hence, this study aimed to assess how tolcapone affects the lipidome of the rat pre-frontal cortex (PFC), a region of the brain highly relevant to tolcapone therapeutic effect, while evaluating its influence on operant behaviour. Tolcapone at 20 mg/kg was chronically administered to Wistar rats during a behavioural task and an untargeted liquid chromatography high-resolution mass spectrometry (LC-HR/MS) approach was employed to profile lipid species. The untargeted analysis identified 7227 features, of which only 33% underwent statistical analysis following data pre-processing. The results revealed an improved cognitive performance and a lipidome remodelling promoted by tolcapone. The lipidomic analysis showed 32 differentially expressed lipid species in tolcapone-treated animals (FC ≥ 1.2, p-value ≤ 0.1), and among these several triacylglycerols, cardiolipins and N-acylethanolamine (NAE 16:2) were found upregulated whereas fatty acids, hexosylceramides, and several phospholipids including phosphatidylcholines and phosphatidylethanolamines were downregulated. These preliminary findings shed light on tolcapone impact on lipid pathways within the brain. Although tolcapone improved cognitive performance and literature suggests the significance of lipids in cognition, this study did not conclusively establish that lipids directly drove or contributed to this outcome. Nevertheless, it underscores the importance of lipid modulation and encourages further exploration of tolcapone-associated mechanisms in the central nervous system (CNS).

Long-term effects of adjunctive eslicarbazepine acetate in adult Asian patients with refractory focal seizures: Post hoc analysis of a phase III trial.

A post hoc analysis of data from Asian patients included in the study BIA-2093-304 was conducted to evaluate the long-term safety/tolerability and efficacy of adjunctive eslicarbazepine acetate (ESL) in adult Asian patients with refractory focal seizures. Part I was a randomized controlled trial, in which patients received ESL (800 or 1200 mg once daily [QD]) or placebo, assessed over a 12-week maintenance period. Patients completing Part I could enter two open-label extension periods (Part II, 1 year; Part III, ≥2 years), during which all received ESL (400-1600 mg QD). Safety/tolerability was assessed by evaluating treatment-emergent adverse events (TEAEs). Efficacy assessments included responder and seizure freedom rates. The safety population included 125, 92, and 23 Asian patients in Parts I, II, and III, respectively. Incidence of ESL-related TEAEs was 61.3%, 45.7%, and 17.4% during Parts I, II, and III, respectively. ESL-related TEAEs (most commonly, dizziness, somnolence, and headache) were consistent with ESL's known safety profile. During Part I, responder rates were higher with ESL 800 (41.7%) and 1200 mg QD (44.4%) versus placebo (32.6%), although not statistically significant. Seizure freedom rates with ESL 800 (5.5%) and 1200 mg QD (11.1%) were also higher versus placebo (0%) (p < 0.05 for ESL 1200 mg QD versus placebo). At the end of Part II, responder and seizure freedom rates were 60.3% and 14.7%, respectively. In summary, adult Asian patients with refractory focal seizures were responsive to treatment with ESL as adjunctive therapy and generally showed treatment tolerance well for up to 3 years. No new/unexpected safety findings were observed.

Population Genomics Reveals the Underlying Structure of the Small Pelagic European Sardine and Suggests Low Connectivity within Macaronesia.

The European sardine (Sardina pilchardus, Walbaum 1792) is indisputably a commercially important species. Previous studies using uneven sampling or a limited number of makers have presented sometimes conflicting evidence of the genetic structure of S. pilchardus populations. Here, we show that whole genome data from 108 individuals from 16 sampling areas across 5000 km of the species' distribution range (from the Eastern Mediterranean to the archipelago of Azores) support at least three genetic clusters. One includes individuals from Azores and Madeira, with evidence of substructure separating these two archipelagos in the Atlantic. Another cluster broadly corresponds to the center of the distribution, including the sampling sites around Iberia, separated by the Almeria-Oran front from the third cluster that includes all of the Mediterranean samples, except those from the Alboran Sea. Individuals from the Canary Islands appear to belong to the Mediterranean cluster. This suggests at least two important geographical barriers to gene flow, even though these do not seem complete, with many individuals from around Iberia and the Mediterranean showing some patterns compatible with admixture with other genetic clusters. Genomic regions corresponding to the top outliers of genetic differentiation are located in areas of low recombination indicative that genetic architecture also has a role in shaping population structure. These regions include genes related to otolith formation, a calcium carbonate structure in the inner ear previously used to distinguish S. pilchardus populations. Our results provide a baseline for further characterization of physical and genetic barriers that divide European sardine populations, and information for transnational stock management of this highly exploited species towards sustainable fisheries.

Lipid profile with eslicarbazepine acetate and carbamazepine monotherapy in adult patients with newly diagnosed focal seizures: post hoc analysis of a phase III trial and open-label extension study.

Background: Antiseizure medications can have negative effects on plasma lipid levels. Objectives: To evaluate plasma lipid changes in patients with newly diagnosed focal epilepsy treated with eslicarbazepine acetate (ESL) or controlled-release carbamazepine (CBZ-CR) monotherapy during a phase III, randomized, double-blind (DB) trial and 2 years of ESL treatment in an open-label extension (OLE). Design: Post hoc analysis of a phase III trial and OLE study. Methods: Proportions of patients with elevated levels of total cholesterol and low-density lipoprotein (LDL) cholesterol were assessed at DB baseline, OLE baseline (last visit of DB trial), and end of OLE. Results: A total of 184 patients received ESL monotherapy during the OLE: 96 received ESL monotherapy in the DB trial and 88 patients received CBZ-CR monotherapy. The proportions of patients with elevated total cholesterol and LDL cholesterol increased significantly during the DB trial in those treated with CBZ-CR monotherapy [total cholesterol, +14.9% (p < 0.001); LDL cholesterol, +11.5% (p = 0.012)] but decreased significantly after switching to ESL monotherapy in the OLE [total cholesterol, -15.3% (p = 0.008); LDL cholesterol, -11.1% (p = 0.021)]. No significant changes were observed in those treated with ESL monotherapy during the DB trial and OLE. At the end of the DB trial, between-group differences (ESL-CBZ-CR) in the proportions of patients with elevated total and LDL cholesterol were -13.6% (p = 0.037) and -12.3% (p = 0.061), respectively; at the end of the OLE, these between-group differences were -6.0% (p = 0.360) and -0.6% (p = 1.000), respectively. Conclusion: A lower proportion of patients with newly diagnosed focal epilepsy had increased levels of total and LDL cholesterol, compared to baseline, following monotherapy with ESL versus CBZ-CR; after switching from CBZ-CR to ESL, the proportions of patients with increased levels decreased significantly. Registration: ClinicalTrials.gov NCT01162460/NCT02484001; EudraCT 2009-011135-13/2015-001243-36.

The impact of treatment with either eslicarbazepine acetate or controlled-release carbamazepine on cholesterol levels in patients with newly diagnosed focal epilepsy Patients with epilepsy have an increased risk of having cardiovascular and cerebrovascular diseases (e.g., myocardial infarction and stroke). Treatment with antiseizure medications can have a negative effect on blood cholesterol levels [such as total cholesterol and low-density lipoprotein (LDL) cholesterol], which can further increase the risk of cardiovascular and cerebrovascular diseases. We examined the impact of monotherapy treatment (i.e., treatment with only one antiseizure medication) using either eslicarbazepine acetate (ESL) or a controlled-release formulation of carbamazepine (CBZ-CR) in 184 patients with newly diagnosed focal epilepsy (ESL, 96 patients; CBZ-CR, 88 patients). Patients received monotherapy with ESL or CBZ-CR for approximately 1 year in a phase III clinical trial. After this, the patients could continue into a 2-year extension study during which they all received monotherapy with ESL. We assessed the proportions of patients with elevated levels of total cholesterol and LDL cholesterol at the beginning and end of the phase III trial, and at the end of the extension study. At the beginning of the phase III trial, the proportions of patients with elevated total cholesterol and elevated LDL cholesterol were similar between treatment groups. During the phase III trial, the proportions of patients with elevated total cholesterol and elevated LDL cholesterol increased in those treated with CBZ-CR monotherapy (total cholesterol, +14.9%; LDL cholesterol, +11.5%) but decreased after switching to ESL monotherapy in the extension study (total cholesterol, −15.3%; LDL cholesterol, −11.1%). By contrast, the proportions of patients with elevated levels of total cholesterol and LDL cholesterol remained relatively stable in those treated with ESL monotherapy during the phase III trial and extension study. These findings indicate that ESL monotherapy may be an appropriate treatment option for patients with newly diagnosed focal epilepsy who either already have, or who are at risk of developing, high levels of cholesterol, since this may reduce their likelihood of having cardiovascular and cerebrovascular diseases.

Examination of gammarid transcriptomes reveals a widespread occurrence of key metabolic genes from epibiont bdelloid rotifers in freshwater species.

Previous data revealed the unexpected presence of genes encoding for long-chain polyunsaturated fatty acid (LC-PUFA) biosynthetic enzymes in transcriptomes from freshwater gammarids but not in marine species, even though closely related species were compared. This study aimed to clarify the origin and occurrence of selected LC-PUFA biosynthesis gene markers across all published gammarid transcriptomes. Through systematic searches, we confirmed the widespread occurrence of sequences from seven elongases and desaturases involved in LC-PUFA biosynthesis, in transcriptomes from freshwater gammarids but not marine species, and clarified that such occurrence is independent from the gammarid species and geographical origin. The phylogenetic analysis established that the retrieved elongase and desaturase sequences were closely related to bdelloid rotifers, confirming that multiple transcriptomes from freshwater gammarids contain contaminating rotifers' genetic material. Using the Adineta steineri genome, we investigated the genomic location and exon-intron organization of the elongase and desaturase genes, establishing they are all genome-anchored and, importantly, identifying instances of horizontal gene transfer. Finally, we provide compelling evidence demonstrating Bdelloidea desaturases and elongases enable these organisms to perform all the reactions for de novo biosynthesis of PUFA and, from them, LC-PUFA, an advantageous trait when considering the low abundance of these essential nutrients in freshwater environments.

Mitochondrial replication's role in vertebrate mtDNA strand asymmetry.

Mitogenomes are defined as compact and structurally stable over aeons. This perception results from a vertebrate-centric vision, where few types of mtDNA rearrangements are described. Here, we bring a new light to the involvement of mitochondrial replication in the strand asymmetry of the vertebrate mtDNA. Using several species of deep-sea hatchetfish (Sternoptychidae) displaying distinct mtDNA structural arrangements, we unravel the inversion of the coding direction of protein-coding genes (PCGs). This unexpected change is coupled with a strand asymmetry nucleotide composition reversal and is shown to be directly related to the strand location of the Control Region (CR). An analysis of the fourfold redundant sites of the PCGs (greater than 6000 vertebrates), revealed the rarity of this phenomenon, found in nine fish species (five deep-sea hatchetfish). Curiously, in Antarctic notothenioid fishes (Trematominae), where a single PCG inversion (the only other record in fish) is coupled with the inversion of the CR, the standard asymmetry is disrupted for the remaining PCGs but not yet reversed, suggesting a transitory state. Our results hint that a relaxation of the classic vertebrate mitochondrial structural stasis promotes disruption of the natural balance of asymmetry of the mtDNA. These findings support the long-lasting hypothesis that replication is the main molecular mechanism promoting the strand-specific compositional bias of this unique and indispensable molecule.

Phylogenomic Analyses of 2,786 Genes in 158 Lineages Support a Root of the Eukaryotic Tree of Life between Opisthokonts and All Other Lineages.

Advances in phylogenomics and high-throughput sequencing have allowed the reconstruction of deep phylogenetic relationships in the evolution of eukaryotes. Yet, the root of the eukaryotic tree of life remains elusive. The most popular hypothesis in textbooks and reviews is a root between Unikonta (Opisthokonta + Amoebozoa) and Bikonta (all other eukaryotes), which emerged from analyses of a single-gene fusion. Subsequent, highly cited studies based on concatenation of genes supported this hypothesis with some variations or proposed a root within Excavata. However, concatenation of genes does not consider phylogenetically-informative events like gene duplications and losses. A recent study using gene tree parsimony (GTP) suggested the root lies between Opisthokonta and all other eukaryotes, but only including 59 taxa and 20 genes. Here we use GTP with a duplication-loss model in a gene-rich and taxon-rich dataset (i.e., 2,786 gene families from two sets of 155 and 158 diverse eukaryotic lineages) to assess the root, and we iterate each analysis 100 times to quantify tree space uncertainty. We also contrasted our results and discarded alternative hypotheses from the literature using GTP and the likelihood-based method SpeciesRax. Our estimates suggest a root between Fungi or Opisthokonta and all other eukaryotes; but based on further analysis of genome size, we propose that the root between Opisthokonta and all other eukaryotes is the most likely.

A genome assembly of the Atlantic chub mackerel (Scomber colias): a valuable teleost fishing resource.

The Atlantic chub mackerel, Scomber colias (Gmelin, 1789), is a medium-sized pelagic fish with substantial importance in the fisheries of the Atlantic Ocean and the Mediterranean Sea. Over the past decade, this species has gained special relevance, being one of the main targets of pelagic fisheries in the NE Atlantic. Here, we sequenced and annotated the first high-quality draft genome assembly of S. colias, produced with PacBio HiFi long reads and Illumina paired-end short reads. The estimated genome size is 814 Mbp, distributed into 2,028 scaffolds and 2,093 contigs with an N50 length of 4.19 and 3.34 Mbp, respectively. We annotated 27,675 protein-coding genes and the BUSCO analyses indicated high completeness, with 97.3% of the single-copy orthologs in the Actinopterygii library profile. The present genome assembly represents a valuable resource to address the biology and management of this relevant fishery. Finally, this genome assembly ranks fourth in high-quality genome assemblies within the order Scombriformes and first in the genus Scomber.

A drastic shift in the energetic landscape of toothed whale sperm cells.

Mammalian spermatozoa are a notable example of metabolic compartmentalization.1 Energy in the form of ATP production, vital for motility, capacitation, and fertilization, is subcellularly separated in sperm cells. While glycolysis provides a local, rapid, and low-yielding input of ATP along the flagellum fibrous sheath, oxidative phosphorylation (OXPHOS), far more efficient over a longer time frame, is concentrated in the midpiece mitochondria.2 The relative weight of glycolysis and OXPHOS pathways in sperm function is variable among species and sensitive to oxygen and substrate availability.3-5 Besides partitioning energy production, sperm cell energetics display an additional singularity: the occurrence of sperm-specific gene duplicates and alternative spliced variants, with conserved function but structurally bound to the flagellar fibrous sheath.6,7 The wider selective forces driving the compartmentalization and adaptability of this energy system in mammalian species remain largely unknown, much like the impact of ecosystem resource availability (e.g., carbohydrates, fatty acids, and proteins) and dietary adaptations in reproductive physiology traits.8 Here, we investigated the Cetacea, an iconic group of fully aquatic and carnivorous marine mammals, evolutionarily related to extant terrestrial herbivores.9 In this lineage, episodes of profound trait remodeling have been accompanied by clear genomic signatures.10-14 We show that toothed whales exhibit impaired sperm glycolysis, due to gene and exon erosion, and demonstrate that dolphin spermatozoa motility depends on endogenous fatty acid ß-oxidation, but not carbohydrates. Such unique energetic rewiring substantiates the observation of large mitochondria in toothed whale spermatozoa and emphasizes the radical physiological reorganization imposed by the transition to a carbohydrate-depleted marine environment.

Inverted replication of vertebrate mitochondria.

After analyzing the base composition asymmetry of coding regions in vertebrate mitochondria, we identified 2 fishes, Albula glossodonta and Bathygadus antrodes, with inverted compositional patterns. Both species appear to have an unusual control region (CR), and in B. antrodes, it has switched from the light strand to the heavy strand. To our knowledge, this is the first report in vertebrates of inverted mitochondrial replication, caused by an inversion of the CR. These findings support the strand-asymmetric model of mtDNA replication and suggest that vertebrate mtDNA can tolerate globally reversed mutational pressures. In addition, we propose that nucleotide bias is not strand specific but that it depends on the location of the CR.

Systematic and phylogeographical assessment of the Acanthodactylus erythrurus group (Reptilia: Lacertidae) based on phylogenetic analyses of mitochondrial and nuclear DNA.

We have used mitochondrial 12S rRNA, 16S rRNA and nuclear beta-fibrinogen (intron 7) sequences to investigate the phylogenetic and phylogeographic relationships between Acanthodactylus erythrurus group species (except for A. boueti). The phylogenetic analyses of the Acanthodactylus genus did not cluster A. guineensis and A. savignyi with the remaining species of the group (A. blanci, A. lineomaculatus and A. erythrurus). Within the A. erythrurus group, the results of the mitochondrial DNA (mtDNA) and nuclear DNA (nDNA) showed a complex phylogeny with geographic structure, but it was not congruent with the present taxonomy. Some taxonomic units, such as A. blanci, A. lineomaculatus, A. e. atlanticus and A. e. belli did not form monophyletic genetic units. The application of a molecular clock suggested that the uplift of the Atlas Mountains in the mid-late Miocene and the reopening of the Strait of Gibraltar could be major biogeographic events responsible for the genetic differentiation in the group. Additionally, diverse micro-evolutionary patterns due to the recent contraction/expansion phases of the habitats in North Africa associated with the high dispersal capabilities of these lizards could be related to the complex phylogenetic patterns observed.

Phylogenomics suggests oxygen availability as a driving force in Thaumarchaeota evolution.

Ammonia-oxidizing archaea (AOA) of the phylum Thaumarchaeota are widespread in marine and terrestrial habitats, playing a major role in the global nitrogen cycle. However, their evolutionary history remains unexplored, which limits our understanding of their adaptation mechanisms. Here, our comprehensive phylogenomic tree of Thaumarchaeota supports three sequential events: origin of AOA from terrestrial non-AOA ancestors, colonization of the shallow ocean, and expansion to the deep ocean. Careful molecular dating suggests that these events coincided with the Great Oxygenation Event around 2300 million years ago (Mya), and oxygenation of the shallow and deep ocean around 800 and 635-560 Mya, respectively. The first transition was likely enabled by the gain of an aerobic pathway for energy production by ammonia oxidation and biosynthetic pathways for cobalamin and biotin that act as cofactors in aerobic metabolism. The first transition was also accompanied by the loss of dissimilatory nitrate and sulfate reduction, loss of oxygen-sensitive pyruvate oxidoreductase, which reduces pyruvate to acetyl-CoA, and loss of the Wood-Ljungdahl pathway for anaerobic carbon fixation. The second transition involved gain of a K+ transporter and of the biosynthetic pathway for ectoine, which may function as an osmoprotectant. The third transition was accompanied by the loss of the uvr system for repairing ultraviolet light-induced DNA lesions. We conclude that oxygen availability drove the terrestrial origin of AOA and their expansion to the photic and dark oceans, and that the stressors encountered during these events were partially overcome by gene acquisitions from Euryarchaeota and Bacteria, among other sources.