Giant lungfish genome elucidates the conquest of land by vertebrates.

Lungfishes belong to lobe-fined fish (Sarcopterygii) that, in the Devonian period, 'conquered' the land and ultimately gave rise to all land vertebrates, including humans1-3. Here we determine the chromosome-quality genome of the Australian lungfish (Neoceratodus forsteri), which is known to have the largest genome of any animal. The vast size of this genome, which is about 14× larger than that of humans, is attributable mostly to huge intergenic regions and introns with high repeat content (around 90%), the components of which resemble those of tetrapods (comprising mainly long interspersed nuclear elements) more than they do those of ray-finned fish. The lungfish genome continues to expand independently (its transposable elements are still active), through mechanisms different to those of the enormous genomes of salamanders. The 17 fully assembled lungfish macrochromosomes maintain synteny to other vertebrate chromosomes, and all microchromosomes maintain conserved ancient homology with the ancestral vertebrate karyotype. Our phylogenomic analyses confirm previous reports that lungfish occupy a key evolutionary position as the closest living relatives to tetrapods4,5, underscoring the importance of lungfish for understanding innovations associated with terrestrialization. Lungfish preadaptations to living on land include the gain of limb-like expression in developmental genes such as hoxc13 and sall1 in their lobed fins. Increased rates of evolution and the duplication of genes associated with obligate air-breathing, such as lung surfactants and the expansion of odorant receptor gene families (which encode proteins involved in detecting airborne odours), contribute to the tetrapod-like biology of lungfishes. These findings advance our understanding of this major transition during vertebrate evolution.

Crossroads in the evolution of plant specialized metabolism.

The monophyletic group of embryophytes (land plants) stands out among photosynthetic eukaryotes: they are the sole constituents of the macroscopic flora on land. In their entirety, embryophytes account for the majority of the biomass on land and constitute an astounding biodiversity. What allowed for the massive radiation of this particular lineage? One of the defining features of all land plants is the production of an array of specialized metabolites. The compounds that the specialized metabolic pathways of embryophytes produce have diverse functions, ranging from superabundant structural polymers and compounds that ward off abiotic and biotic challenges, to signaling molecules whose abundance is measured at the nanomolar scale. These specialized metabolites govern the growth, development, and physiology of land plants-including their response to the environment. Hence, specialized metabolites define the biology of land plants as we know it. And they were likely a foundation for their success. It is thus intriguing to find that the closest algal relatives of land plants, freshwater organisms from the grade of streptophyte algae, possess homologs for key enzymes of specialized metabolic pathways known from land plants. Indeed, some studies suggest that signature metabolites emerging from these pathways can be found in streptophyte algae. Here we synthesize the current understanding of which routes of the specialized metabolism of embryophytes can be traced to a time before plants had conquered land.

Widespread occurrence of covalent lysine-cysteine redox switches in proteins.

We recently reported the discovery of a lysine-cysteine redox switch in proteins with a covalent nitrogen-oxygen-sulfur (NOS) bridge. Here, a systematic survey of the whole protein structure database discloses that NOS bridges are ubiquitous redox switches in proteins of all domains of life and are found in diverse structural motifs and chemical variants. In several instances, lysines are observed in simultaneous linkage with two cysteines, forming a sulfur-oxygen-nitrogen-oxygen-sulfur (SONOS) bridge with a trivalent nitrogen, which constitutes an unusual native branching cross-link. In many proteins, the NOS switch contains a functionally essential lysine with direct roles in enzyme catalysis or binding of substrates, DNA or effectors, linking lysine chemistry and redox biology as a regulatory principle. NOS/SONOS switches are frequently found in proteins from human and plant pathogens, including severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), and also in many human proteins with established roles in gene expression, redox signaling and homeostasis in physiological and pathophysiological conditions.

Phylogeny and evolution of streptophyte algae.

The Streptophyta emerged about a billion years ago. Nowadays, this branch of the green lineage is most famous for one of its clades, the land plants (Embryophyta). While Embryophyta make up the major share of species numbers in Streptophyta, there is a diversity of likely more than 5000 species of streptophyte algae that form a paraphyletic grade next to land plants. Here, we focus on the deep divergences that gave rise to the diversity of streptophytes-and thus, particularly on the streptophyte algae. Phylogenomic efforts have not only clarified the position of streptophyte algae to land plants but recent efforts have also begun to unravel the relationships and major radiations within streptophyte algal diversity. We illustrate how new phylogenomic perspectives have changed our view on the evolutionary emergence of key traits such as intricate signaling networks that are intertwined with multicellular growth and the chemodiverse hotbed from which they emerged. These traits are key for the biology of land plants-but were bequeathed from their algal progenitors.

Sequence similarity networks bear out hierarchical relationships of green cytochrome P450.

Land plants have diversified enzyme families. One of the most prominent is the cytochrome P450 (CYP or CYP450) family. With over 443,000 CYP proteins sequenced across the tree of life, CYPs are ubiquitous in archaea, bacteria, and eukaryotes. Here, we focused on land plants and algae to study the role of CYP diversification. CYPs, acting as monooxygenases, catalyze hydroxylation reactions crucial for specialized plant metabolic pathways, including detoxification and phytohormone production; the CYPome consists of one enormous superfamily that is divided into clans and families. Their evolutionary history speaks of high substrate promiscuity; radiation and functional diversification have yielded numerous CYP families. To understand the evolutionary relationships within the CYPs, we employed sequence similarity network analyses. We recovered distinct clusters representing different CYP families, reflecting their diversified sequences that we link to the prediction of functionalities. Hierarchical clustering and phylogenetic analysis further elucidated relationships between CYP clans, uncovering their shared deep evolutionary history. We explored the distribution and diversification of CYP subfamilies across plant and algal lineages, uncovering novel candidates and providing insights into the evolution of these enzyme families. This identified unexpected relationships between CYP families, such as the link between CYP82 and CYP74, shedding light on their roles in plant defense signaling pathways. Our approach provides a methodology that brings insights into the emergence of new functions within the CYP450 family, contributing to the evolutionary history of plants and algae. These insights can be further validated and implemented via experimental setups under various external conditions.

A seed-like proteome in oil-rich tubers.

There are numerous examples of plant organs or developmental stages that are desiccation-tolerant and can withstand extended periods of severe water loss. One prime example are seeds and pollen of many spermatophytes. However, in some plants, also vegetative organs can be desiccation-tolerant. One example are the tubers of yellow nutsedge (Cyperus esculentus), which also store large amounts of lipids similar to seeds. Interestingly, the closest known relative, purple nutsedge (Cyperus rotundus), generates tubers that do not accumulate oil and are not desiccation-tolerant. We generated nanoLC-MS/MS-based proteomes of yellow nutsedge in five replicates of four stages of tuber development and compared them to the proteomes of roots and leaves, yielding 2257 distinct protein groups. Our data reveal a striking upregulation of hallmark proteins of seeds in the tubers. A deeper comparison to the tuber proteome of the close relative purple nutsedge (C. rotundus) and a previously published proteome of Arabidopsis seeds and seedlings indicates that indeed a seed-like proteome was found in yellow but not purple nutsedge. This was further supported by an analysis of the proteome of a lipid droplet-enriched fraction of yellow nutsedge, which also displayed seed-like characteristics. One reason for the differences between the two nutsedge species might be the expression of certain transcription factors homologous to ABSCISIC ACID INSENSITIVE3, WRINKLED1, and LEAFY COTYLEDON1 that drive gene expression in Arabidopsis seed embryos.

A Phylogenomic Backbone for Gastropod Molluscs.

Gastropods have survived several mass extinctions during their evolutionary history resulting in extraordinary diversity in morphology, ecology, and developmental modes, which complicate the reconstruction of a robust phylogeny. Currently, gastropods are divided into six subclasses: Caenogastropoda, Heterobranchia, Neomphaliones, Neritimorpha, Patellogastropoda, and Vetigastropoda. Phylogenetic relationships among these taxa historically lack consensus, despite numerous efforts using morphological and molecular information. We generated sequence data for transcriptomes derived from 12 taxa belonging to clades with little or no prior representation in previous studies in order to infer the deeper cladogenetic events within Gastropoda and, for the first time, infer the position of the deep-sea Neomphaliones using a phylogenomic approach. We explored the impact of missing data, homoplasy, and compositional heterogeneity on the inferred phylogenetic hypotheses. We recovered a highly supported backbone for gastropod relationships that is congruent with morphological and mitogenomic evidence, in which Patellogastropoda, true limpets, are the sister lineage to all other gastropods (Orthogastropoda) which are divided into two main clades 1) Vetigastropoda $s.l.$ (including Pleurotomariida $+$ Neomphaliones) and 2) Neritimorpha $+$ (Caenogastropoda $+$ Heterobranchia). As such, our results support the recognition of five subclasses (or infraclasses) in Gastropoda: Patellogastropoda, Vetigastropoda, Neritimorpha, Caenogastropoda, and Heterobranchia. [Compositional heterogeneity; fast-evolving; long-branch attraction; missing data; Mollusca; phylogenetics; systematic error.].

Single-cell genomics reveals new rozellid lineages and supports their sister relationship to Microsporidia.

The phylum Rozellomycota has been proposed for a group of early-branching holomycotan lineages representing obligate parasites and hyperparasites of zoosporic fungi, oomycotes or phytoplankton. Given their predominantly intracellular lifestyle, rozellids are typically known from environmental ribosomal DNA data, except for the well-studied Rozella species. To date, the phylogenetic relationship between rozellids and microsporidians (Microsporidia) is not fully understood and most reliable hypotheses are based on phylogenomic analyses that incorporate the only publicly available rozellid genome of Rozella allomycis. Here, we provide genomic data of three new rozellid lineages obtained by single-cell sequencing from environmental samples and show with a phylogenomic approach that rozellids form a monophyletic group that is sister to microsporidians, corroborating the previously proposed phylum Rozellomycota. Whereas no mitochondrial genes coding for the respiratory Complex I could be found, we discovered a gene coding for a nucleotide phosphate transporter in one of the three draft genomes. The scattered absence of Complex I genes and scattered presence of nucleotide transporter genes across diverse microsporidian and rozellid lineages suggest that these adaptations to a parasitic lifestyle, which reduce the parasite's capability to synthesize ATP but enables it to steal ATP from its host, evolved independently in microsporidians and rozellids.

Divergent responses in desiccation experiments in two ecophysiologically different Zygnematophyceae.

Water scarcity can be considered a major stressor on land, with desiccation being its most extreme form. Land plants have found two different solutions to this challenge: avoidance and tolerance. The closest algal relatives to land plants, the Zygnematophyceae, use the latter, and how this is realized is of great interest for our understanding of the conquest of land. Here, we worked with two representatives of the Zygnematophyceae, Zygnema circumcarinatum SAG 698-1b and Mesotaenium endlicherianum SAG 12.97, who differ in habitats and drought resilience. We challenged both algal species with severe desiccation in a laboratory setup until photosynthesis ceased, followed by a recovery period. We assessed their morphological, photophysiological, and transcriptomic responses. Our data pinpoint global differential gene expression patterns that speak of conserved responses, from calcium-mediated signaling to the adjustment of plastid biology, cell envelopes, and amino acid pathways, between Zygnematophyceae and land plants despite their strong ecophysiological divergence. The main difference between the two species appears to rest in a readjustment of the photobiology of Zygnema, while Mesotaenium experiences stress beyond a tipping point.

The evolution of the phenylpropanoid pathway entailed pronounced radiations and divergences of enzyme families.

Land plants constantly respond to fluctuations in their environment. Part of their response is the production of a diverse repertoire of specialized metabolites. One of the foremost sources for metabolites relevant to environmental responses is the phenylpropanoid pathway, which was long thought to be a land-plant-specific adaptation shaped by selective forces in the terrestrial habitat. Recent data have, however, revealed that streptophyte algae, the algal relatives of land plants, have candidates for the genetic toolkit for phenylpropanoid biosynthesis and produce phenylpropanoid-derived metabolites. Using phylogenetic and sequence analyses, we here show that the enzyme families that orchestrate pivotal steps in phenylpropanoid biosynthesis have independently undergone pronounced radiations and divergence in multiple lineages of major groups of land plants; sister to many of these radiated gene families are streptophyte algal candidates for these enzymes. These radiations suggest a high evolutionary versatility in the enzyme families involved in the phenylpropanoid-derived metabolism across embryophytes. We suggest that this versatility likely translates into functional divergence, and may explain the key to one of the defining traits of embryophytes: a rich specialized metabolism.

Extensive introgression at late stages of species formation: Insights from grasshopper hybrid zones.

The process of species formation is characterized by the accumulation of multiple reproductive barriers. The evolution of hybrid male sterility, or Haldane's rule, typically characterizes later stages of species formation, when reproductive isolation is strongest. Yet, understanding how quickly reproductive barriers evolve and their consequences for maintaining genetic boundaries between emerging species remains a challenging task because it requires studying taxa that hybridize in nature. Here, we address these questions using the meadow grasshopper Pseudochorthippus parallelus, where populations that show multiple reproductive barriers, including hybrid male sterility, hybridize in two natural hybrid zones. Using mitochondrial data, we infer that such populations diverged some 100,000 years ago, at the beginning of the last glacial cycle in Europe. Nuclear data show that contractions at multiple glacial refugia, and post-glacial expansions have facilitated genetic differentiation between lineages that today interact in hybrid zones. We find extensive introgression throughout the sampled species range, irrespective of the current strength of reproductive isolation. Populations exhibiting hybrid male sterility in two hybrid zones show repeatable patterns of genomic differentiation, consistent with shared genomic constraints affecting ancestral divergence or with the role of those regions in reproductive isolation. Together, our results suggest that reproductive barriers that characterize late stages of species formation can evolve relatively quickly, particularly when associated with strong demographic changes. Moreover, we show that such barriers persist in the face of extensive gene flow, allowing future studies to identify associated genomic regions.

Phylogenomic Insights into the Origin of Primary Plastids.

The origin of plastids was a major evolutionary event that paved the way for an astonishing diversification of photosynthetic eukaryotes. Plastids originated by endosymbiosis between a heterotrophic eukaryotic host and cyanobacteria, presumably in a common ancestor of the primary photosynthetic eukaryotes (Archaeplastida). A single origin of primary plastids is well supported by plastid evidence but not by nuclear phylogenomic analyses, which have consistently failed to recover the monophyly of Archaeplastida hosts. Importantly, plastid monophyly and nonmonophyletic hosts could be explained under scenarios of independent or serial eukaryote-to-eukaryote endosymbioses. Here, we assessed the strength of the signal for the monophyly of Archaeplastida hosts in four available phylogenomic data sets. The effect of phylogenetic methodology, data quality, alignment trimming strategy, gene and taxon sampling, and the presence of outlier genes were investigated. Our analyses revealed a lack of support for host monophyly in the shorter individual data sets. However, when analyzed together under rigorous data curation and complex mixture models, the combined nuclear data sets supported the monophyly of primary photosynthetic eukaryotes (Archaeplastida) and recovered a putative association with plastid-lacking Picozoa. This study represents an important step toward better understanding deep eukaryotic evolution and the origin of plastids. [Archaeplastida; Bayesian; chloroplast; maximum likelihood; mixture model; ortholog; outlier loci; paralog; protist.].

Functional genomics of abiotic environmental adaptation in lacertid lizards and other vertebrates.

Understanding the genomic basis of adaptation to different abiotic environments is important in the context of climate change and resulting short-term environmental fluctuations. Using functional and comparative genomics approaches, we here investigated whether signatures of genomic adaptation to a set of environmental parameters are concentrated in specific subsets of genes and functions in lacertid lizards and other vertebrates. We first identify 200 genes with signatures of positive diversifying selection from transcriptomes of 24 species of lacertid lizards and demonstrate their involvement in physiological and morphological adaptations to climate. To understand how functionally similar these genes are to previously predicted candidate functions for climate adaptation and to compare them with other vertebrate species, we then performed a meta-analysis of 1,100 genes under selection obtained from -omics studies in vertebrate species adapted to different abiotic factors. We found that the vertebrate gene set formed a tightly connected interactome, which was to 23% enriched in previously predicted functions of adaptation to climate, and to a large part (18%) involved in organismal stress response. We found a much higher degree of identical genes being repeatedly selected among different animal groups (43.6%), and of functional similarity and post-translational modifications than expected by chance, and no clear functional division between genes used for ectotherm and endotherm physiological strategies. In total, 171 out of 200 genes of Lacertidae were part of this network. These results highlight an important role of a comparatively small set of genes and their functions in environmental adaptation and narrow the set of candidate pathways and markers to be used in future research on adaptation and stress response related to climate change.

Unexpected cryptic species among streptophyte algae most distant to land plants.

Streptophytes are one of the major groups of the green lineage (Chloroplastida or Viridiplantae). During one billion years of evolution, streptophytes have radiated into an astounding diversity of uni- and multicellular green algae as well as land plants. Most divergent from land plants is a clade formed by Mesostigmatophyceae, Spirotaenia spp. and Chlorokybophyceae. All three lineages are species-poor and the Chlorokybophyceae consist of a single described species, Chlorokybus atmophyticus. In this study, we used phylogenomic analyses to shed light into the diversity within Chlorokybus using a sampling of isolates across its known distribution. We uncovered a consistent deep genetic structure within the Chlorokybus isolates, which prompted us to formally extend the Chlorokybophyceae by describing four new species. Gene expression differences among Chlorokybus species suggest certain constitutive variability that might influence their response to environmental factors. Failure to account for this diversity can hamper comparative genomic studies aiming to understand the evolution of stress response across streptophytes. Our data highlight that future studies on the evolution of plant form and function can tap into an unknown diversity at key deep branches of the streptophytes.

Phylogenomics and evolutionary history of Oreobates (Anura: Craugastoridae) Neotropical frogs along elevational gradients.

Mountain ranges offer opportunities for understanding how species evolved and diversified across different environmental conditions. Neotropical frogs of the genus Oreobates (Anura: Craugastoridae) are adapted to highland and lowland habitats along the Andes, but many aspects of their evolution remain unknown. We studied their evolutionary history using ~18,000 exons enriched by targeted sequence-capture. Since capture success was very variable across samples, we evaluated to what degree differing data filtering produced robust inferences. The inferred evolutionary framework evidenced phylogenetic discordances among lowland species that can be explained by taxonomic misidentification or admixture of ancestral lineages. Highland species showed smaller effective populations than lowland frogs, probably due to greater habitat fragmentation in montane environments. Stronger genetic drift likely decreased the power of purifying selection and led to an increased proportion of nonsynonymous mutations in highland populations that could play an important role in their adaptation. Overall, our work sheds light on the evolutionary history and diversification of this group of Neotropical frogs along elevational gradients in the Andes as well as on their patterns of intraspecific diversity.

Phylotranscriptomic evidence for pervasive ancient hybridization among Old World salamanders.

Hybridization can leave genealogical signatures in an organism's genome, originating from the parental lineages and persisting over time. This potentially confounds phylogenetic inference methods that aim to represent evolution as a strictly bifurcating tree. We apply a phylotranscriptomic approach to study the evolutionary history of, and test for inter-lineage introgression in the Salamandridae, a Holarctic salamanders group of interest in studies of toxicity and aposematism, courtship behavior, and molecular evolution. Although the relationships between the 21 currently recognized salamandrid genera have been the subject of numerous molecular phylogenetic studies, some branches have remained controversial and sometimes affected by discordances between mitochondrial vs. nuclear trees. To resolve the phylogeny of this family, and understand the source of mito-nuclear discordance, we generated new transcriptomic (RNAseq) data for 20 salamandrids and used these along with published data, including 28 mitochondrial genomes, to obtain a comprehensive nuclear and mitochondrial perspective on salamandrid evolution. Our final phylotranscriptomic data set included 5455 gene alignments for 40 species representing 17 of the 21 salamandrid genera. Using concatenation and species-tree phylogenetic methods, we find (1) Salamandrina sister to the clade of the "True Salamanders" (consisting of Chioglossa, Mertensiella, Lyciasalamandra, and Salamandra), (2) Ichthyosaura sister to the Near Eastern genera Neurergus and Ommatotriton, (3) Triturus sister to Lissotriton, and (4) Cynops paraphyletic with respect to Paramesotriton and Pachytriton. Combining introgression tests and phylogenetic networks, we find evidence for introgression among taxa within the clades of "Modern Asian Newts" and "Modern European Newts". However, we could not unambiguously identify the number, position, and direction of introgressive events. Combining evidence from nuclear gene analysis with the observed mito-nuclear phylogenetic discordances, we hypothesize a scenario with hybridization and mitochondrial capture among ancestral lineages of (1) Lissotriton into Ichthyosaura and (2) Triturus into Calotriton, plus introgression of nuclear genes from Triturus into Lissotriton. Furthermore, both mitochondrial capture and nuclear introgression may have occurred among lineages assigned to Cynops. More comprehensive genomic data will, in the future, allow testing this against alternative scenarios involving hybridization with other, extinct lineages of newts.

A mitogenomic phylogeny of chitons (Mollusca: Polyplacophora).

BACKGROUND: Polyplacophora, or chitons, have long fascinated malacologists for their distinct and rather conserved morphology and lifestyle compared to other mollusk classes. However, key aspects of their phylogeny and evolution remain unclear due to the few morphological, molecular, or combined phylogenetic analyses, particularly those addressing the relationships among the major chiton lineages. RESULTS: Here, we present a mitogenomic phylogeny of chitons based on 13 newly sequenced mitochondrial genomes along with eight available ones and RNAseq-derived mitochondrial sequences from four additional species. Reconstructed phylogenies largely agreed with the latest advances in chiton systematics and integrative taxonomy but we identified some conflicts that call for taxonomic revisions. Despite an overall conserved gene order in chiton mitogenomes, we described three new rearrangements that might have taxonomic utility and reconstructed the most likely scenario of gene order change in this group. Our phylogeny was time-calibrated using various fossils and relaxed molecular clocks, and the robustness of these analyses was assessed with several sensitivity analyses. The inferred ages largely agreed with previous molecular clock estimates and the fossil record, but we also noted that the ambiguities inherent to the chiton fossil record might confound molecular clock analyses. CONCLUSIONS: In light of the reconstructed time-calibrated framework, we discuss the evolution of key morphological features and call for a continued effort towards clarifying the phylogeny and evolution of chitons.

Historical isolation facilitates species radiation by sexual selection: Insights from Chorthippus grasshoppers.

Theoretical and empirical studies have shown that species radiations are facilitated when a trait under divergent natural selection is also involved in sexual selection. It is yet unclear how quick and effective radiations are where assortative mating is unrelated to the ecological environment and primarily results from sexual selection. We address this question using sympatric grasshopper species of the genus Chorthippus, which have evolved strong behavioural isolation while lacking noticeable ecomorphological divergence. Mitochondrial genomes suggest that the radiation is relatively recent, dating to the mid-Pleistocene, which leads to extensive incomplete lineage sorting throughout the mitochondrial and nuclear genomes. Nuclear data shows that hybrids are absent in sympatric localities but that all species have experienced gene flow, confirming that reproductive isolation is strong but remains incomplete. Demographic modelling is most consistent with a long period of geographic isolation, followed by secondary contact and extensive introgression. Such initial periods of geographic isolation might facilitate the association between male signaling and female preference, permitting the coexistence of sympatric species that are genetically, morphologically, and ecologically similar, but otherwise behave mostly as good biological species.

PREQUAL: detecting non-homologous characters in sets of unaligned homologous sequences.

Summary: Phylogenomic datasets invariably contain undetected stretches of non-homologous characters due to poor-quality sequences or erroneous gene models. The large-scale multi-gene nature of these datasets renders impractical or impossible detailed manual curation of sequences, but few tools exist that can automate this task. To address this issue, we developed a new method that takes as input a set of unaligned homologous sequences and uses an explicit probabilistic approach to identify and mask regions with non-homologous adjacent characters. These regions are defined as sharing no statistical support for homology with any other sequence in the set, which can result from e.g. sequencing errors or gene prediction errors creating frameshifts. Our methodology is implemented in the program PREQUAL, which is a fast and accurate tool for high-throughput filtering of sequences. The program is primarily aimed at amino acid sequences, although it can handle protein coding DNA sequences as well. It is fully customizable to allow fine-tuning of the filtering sensitivity. Availability and implementation: The program PREQUAL is written in C/C++ and available through a GNU GPL v3.0 at https://github.com/simonwhelan/prequal. Supplementary information: Supplementary data are available at Bioinformatics online.

New patellogastropod mitogenomes help counteracting long-branch attraction in the deep phylogeny of gastropod mollusks.

Long-branch attraction (LBA) is a well-known artifact in phylogenetic reconstruction. Sparse taxon sampling and extreme heterogeneity of evolutionary rates among lineages generate propitious situations for LBA, even defying probabilistic methods of phylogenetic inference. A clear example illustrating LBA challenges is the difficulty of reconstructing the deep gastropod phylogeny, particularly using mitochondrial (mt) genomes. Previous studies consistently obtained unorthodox phylogenetic relationships due to the LBA between the mitogenomes of patellogastropods (true limpets, represented only by Lottia digitalis), heterobranchs, and outgroup taxa. Here, we use the reconstruction of the gastropod mitogenomic phylogeny as a case exercise to test the effect of key methodological approaches proposed to counteract LBA, including the selection of slow-evolving representatives, the use of different outgroups, the application of site-heterogeneous evolutionary models, and the removal of fast-evolving sites. In this regard, we sequenced three new patellogastropod mt genomes, which displayed shorter branches than the one of Lottia as well as gene organizations more similar to that of the hypothetical gastropod ancestor. Phylogenetic analyses incorporating the mt genomes of Patella ferruginea, Patella vulgata, and Cellana radiata allowed eliminating the artificial clustering of Patellogastropoda and Heterobranchia that had prevailed in previous studies. Furthermore, the use of site-heterogeneous models with certain combinations of lineages within the outgroup allowed eliminating also the LBA between Heterobranchia and the outgroup, and recovering Apogastropoda (i.e., Caenogastropoda + Heterobranchia). Hence, for the first time, we were able to obtain a mitogenomic phylogeny of gastropods that is congruent with both morphological and nuclear datasets.