Improved annotation of protein-coding genes boundaries in metazoan mitochondrial genomes

Comparative mitochondrial genomics of endemic Mexican vesper yellow bats genus Rhogeessa (Chiroptera: Vespertilionidae) and insights into internal relationships in the family Vespertilionidae

In the species-rich family Vespertilionidae, vesper yellow bats in the genus Rhogeessa include eleven species, three of them endemic to Mexico. These insectivorous bats provide important ecosystem services, including pest control. Even though some aspects of their biology are well- known, only a few genomic resources are available for these species, which limits our understanding of their biology. In this study, we assembled and annotated the mitochondrial genome of four species: R. aenea, R. genowaysi, R. mira, and R. parvula. We generated a phylomitogenomic hypothesis based on translated protein-coding genes for a total of 52 species in the family Vespertilionidae and examined the phylogenetic position of the genus Rhogeessa and species within the family. The AT-rich mitogenomes of R. aenea, R. genowaysi, R. mira, and R. parvula are 16,763, 16,781, 16,807, and 16,794 pb in length, respectively. Each studied mitogenome encodes 13 Protein Coding Genes (PCGs), 22 transfer RNA genes, and 2 rRNA genes, and contains a putative control region (CR). All tRNAs exhibit a 'cloverleaf' secondary structure, except tRNA-Serine-1 that lacked the DHU arm in all studied mitogenomes. Selective pressure analyses indicated that all protein-coding genes are exposed to purifying selection. The phylomitogenomic analysis supported the monophyletic status of the family Vespertilionidae, confirmed the placement of Rhogeessa within the tribe Antrozoini, and clarified phylogenetic relationships within and among subfamilies and tribes in this family. Our results indicate that phylomitogenomics are useful to explore the evolutionary history of vesper bats. The assembly and comprehensive analysis of mitochondrial genomes offer the potential to generate molecular references and resources beneficial for genetic analyses aimed at understanding the ecology and evolution of these remarkable bats.

The mitochondrial genome of Homoneura picta (Diptera: Lauxaniidae) and its phylogenetic analysis

Homoneura picta belongs to the Homoneurinae subfamily of Lauxaniidae, and it is widely distributed and common in China. This study reports the newly sequenced mitochondrial genome of H. picta. The sequence is 15,469 bp long and contains 37 genes (13 protein-coding, 22 tRNA, and 2 rRNA genes) and a control region. The overall base composition is 38.4% for A, 37.7% for T, 14.1% for C, and 9.8% for G, with a bias toward A + T (76.1%). Phylogenetic analysis show that Homoneura is a sister genus of Cestrotus. We have successfully sequenced the mitochondrial genome of H. picta, which can be useful in investigating the phylogenetic status of Homoneurinae. Our results provide data for further studies of phylogeny in Diptera.

Genome assembly and variant analysis of two Saccharomyces cerevisiae strains isolated from stingless bee pollen

Products from stingless bees are rich reservoirs of microbial diversity, including yeasts with fermentative potential. Previously, two Saccharomyces cerevisiae strains, JP14 and IP9, were isolated from Jataí (Tetragonisca angustula) and Iraí (Nannotrigona testaceicornis) bees, respectively, aiming at mead production. Both strains presented great osmotic and sulfite tolerance, and ethanol production, although they have a high free amino nitrogen demand. Herein, their genomes were sequenced, assembled, and annotated, and the variants were compared to the S. cerevisiae S288c reference strain. The final assembly of IP9 and JP14 presented high N50 and BUSCO scores, and more than 6430 protein-coding genes. Additionally, nQuire predicted the ploidy of IP9 as diploid, but the results were not enough to determine the ploidy of JP14. The mitochondrial genomes of IP9 and JP14 presented the same gene content as S288c but the genes were rearranged and fragmented in different patterns. Meanwhile, the genes with mutations of high impact (e.g., indels, gain of stop codon) for both yeasts were enriched for transmembrane transport, electron transfer, oxidoreductase, heme binding, fructose, mannose, and glucose transport, activities related to the respiratory chain and sugar metabolism. The IP9 strain presented copy number gains in genes related to sugar transport and cell morphogenesis; in JP14, genes were enriched for disaccharide metabolism and transport, response to reactive oxygen species, and polyamine transport. On the other hand, IP9 presented copy number losses related to disaccharide, thiamine, and aldehyde metabolism, while JP14 presented depletions related to disaccharide, oligosaccharide, asparagine, and aspartate metabolism. Notably, both strains presented a killer toxin gene, annotated from the assembling of unmapped reads, representing a potential mechanism for the control of other microorganisms population in the environment. Therefore, the annotated genomes of JP14 and IP9 presented a high selective pressure for sugar and nitrogen metabolism and stress response, consistent with their isolation source and fermentative properties.

New mitogenomes of Runcinidae and Facelinidae: two understudied heterobranch families (Mollusca: Gastropoda)

Here, we present the mitochondrial sequences of two sea slugs (Heterobranchia): Runcina aurata and Facelina auriculata, the latter being the type species of the family. The mitochondrial genomes are 14,282 and 14,171bp in length, respectively, with a complete set of 13 PCGs, 2 rRNAs, and 22 tRNAs. None of the mitogenomes show gene reorganization, keeping the standard mitogenomic structure of Heterobranchia. Nucleotide composition differs significantly between them, with R. aurata showing the most AT-rich mitogenome (25.7% GC content) reported to date in Heterobranchia, and F. auriculata showing a rich GC content (35%) compared with other heterobranch mitochondrial genomes.

Characterization of the mitochondrial genomes of the Mexican endemic bats Corynorhinus mexicanus and Corynorhinus leonpaniaguae (Chiroptera: Vespertilionidae)

BACKGROUND

The genus Corynorhinus is composed of four recognized species: C. rafinesquii, C. townsendii, C. mexicanus, and C. leonpaniaguae, the latter two being endemic to Mexico. According to the IUCN, C. mexicanus is considered "Near Threatened", as its populations are dwindling and habitats are affected by anthropogenic disturbance. Corynorhinus leonpaniaguae has not been assigned to an IUCN Red List risk category due to its recent description.

METHODS AND RESULTS

In this study, the mitochondrial genomes of C. mexicanus and C. leonpaniaguae were assembled and characterized in detail. The mitochondrial genomes (mtDNA) of C. mexicanus and C. leonpaniaguae have lengths of 16,470 and 16,581 bp respectively, with a predominant nucleotide usage of adenine (31.670% and 31.729%, respectively) and thymine (26.15% and 26.18%, respectively). The mtDNA of C. mexicanus and C. leonpaniaguae is composed of 37 coding and non-coding elements: 22 transfer RNAs (tRNA), 13 protein-coding genes (PCGs), two ribosomal RNAs and a non-coding region, the control region, which has a length of 933 bp and 1,149 bp, respectively. All tRNAs exhibited a cloverleaf secondary structure, with the exception of trn-Ser1 which showed a deletion of the dihydrouridine arm in the two species. All PCGs are subjected to purifying selection, with atp8 being the gene showing the highest Ka/Ks value.

CONCLUSIONS

These are the first whole mitogenomic resources developed for C. mexicanus and C. leonpaniaguae and enhance our knowledge of the ecology of these species and aid in their conservation.

High heteroplasmy is associated with low mitochondrial copy number and selection against non-synonymous mutations in the snail Cepaea nemoralis

Molluscan mitochondrial genomes are unusual because they show wide variation in size, radical genome rearrangements and frequently show high variation (> 10%) within species. As progress in understanding this variation has been limited, we used whole genome sequencing of a six-generation matriline of the terrestrial snail Cepaea nemoralis, as well as whole genome sequences from wild-collected C. nemoralis, the sister species C. hortensis, and multiple other snail species to explore the origins of mitochondrial DNA (mtDNA) variation. The main finding is that a high rate of SNP heteroplasmy in somatic tissue was negatively correlated with mtDNA copy number in both Cepaea species. In individuals with under ten mtDNA copies per nuclear genome, more than 10% of all positions were heteroplasmic, with evidence for transmission of this heteroplasmy through the germline. Further analyses showed evidence for purifying selection acting on non-synonymous mutations, even at low frequency of the rare allele, especially in cytochrome oxidase subunit 1 and cytochrome b. The mtDNA of some individuals of Cepaea nemoralis contained a length heteroplasmy, including up to 12 direct repeat copies of tRNA-Val, with 24 copies in another snail, Candidula rugosiuscula, and repeats of tRNA-Thr in C. hortensis. These repeats likely arise due to error prone replication but are not correlated with mitochondrial copy number in C. nemoralis. Overall, the findings provide key insights into mechanisms of replication, mutation and evolution in molluscan mtDNA, and so will inform wider studies on the biology and evolution of mtDNA across animal phyla.

The nuclear and mitochondrial genome assemblies of Tetragonisca angustula (Apidae: Meliponini), a tiny yet remarkable pollinator in the Neotropics

BACKGROUND

The field of bee genomics has considerably advanced in recent years, however, the most diverse group of honey producers on the planet, the stingless bees, are still largely neglected. In fact, only eleven of the ~ 600 described stingless bee species have been sequenced, and only three using a long-read (LR) sequencing technology. Here, we sequenced the nuclear and mitochondrial genomes of the most common, widespread and broadly reared stingless bee in Brazil and other neotropical countries-Tetragonisca angustula (popularly known in Brazil as jataí).

RESULTS

A total of 48.01 Gb of DNA data were generated, including 2.31 Gb of Pacific Bioscience HiFi reads and 45.70 Gb of Illumina short reads (SRs). Our preferred assembly comprised 683 contigs encompassing 284.49 Mb, 62.84 Mb of which (22.09%) corresponded to 445,793 repetitive elements. N50, L50 and complete BUSCOs reached 1.02 Mb, 91 contigs and 97.1%, respectively. We predicted that the genome of T. angustula comprises 17,459 protein-coding genes and 4,108 non-coding RNAs. The mitogenome consisted of 17,410 bp, and all 37 genes were found to be on the positive strand, an unusual feature among bees. A phylogenomic analysis of 26 hymenopteran species revealed that six odorant receptor orthogroups of T. angustula were found to be experiencing rapid evolution, four of them undergoing significant contractions.

CONCLUSIONS

Here, we provided the first nuclear and mitochondrial genome assemblies for the ecologically and economically important T. angustula, the fourth stingless bee species to be sequenced with LR technology thus far. We demonstrated that even relatively small amounts of LR data in combination with sufficient SR data can yield high-quality genome assemblies for bees.

The complete mitochondrial genome of Potomida acarnanica (Kobelt, 1879)

Freshwater mussels (Bivalvia, Unionida) play essential roles in the well-functioning of ecosystems, even providing essential services to humans. However, these bivalves face numerous threats (e.g. habitat loss and fragmentation, pollution, introduction of invasive species, and climate change) which have already led to the extinction of many populations. This underscores the need to fully characterize the biology of these species, particularly those, such as Potomida acarnanica, that are still poorly studied. This study presents the first mitogenome of P. acarnanica (Kobelt, 1879), an endemic species of Greece with a distribution limited to only two river basins. The mitochondrial genome of a P. acarnanica specimen, collected at Pamisos River (Peloponnese, Greece), was sequenced by Illumina high-throughput sequencing. This mitogenome (16,101 bp) is characterized by 13 protein-coding genes, 22 transfer RNA and 2 ribosomal RNA genes. The size of this mitogenome is within the range of another Potomida mitogenome already published for the species Potomida littoralis. In the phylogenetic inference, P. acarnanica was recovered as monophyletic with P. littoralis mitogenome in the Lamprotulini tribe, as expected. This genomic resource will assist in genetically characterizing the species, potentially benefiting future evolutionary studies and conservation efforts.

Japanese Planocerid Flatworms: Difference in Composition of Tetrodotoxin and Its Analogs and the Effects of Ingestion by Toxin-Bearing Fishes in the Ryukyu Islands, Japan

Tetrodotoxin (TTX), known as pufferfish toxin, is a potent neurotoxin blocking sodium channels in muscle and nerve tissues. TTX has been detected in various taxa other than pufferfish, including marine polyclad flatworms, suggesting that pufferfish toxin accumulates in fish bodies via food webs. The composition of TTX and its analogs in the flatworm Planocera multitentaculata was identical to those in wild grass puffer Takifugu alboplumbeus. Previously, Planocera sp. from Okinawa Island, Japan, were reported to possess high level of TTX, but no information was available on TTX analogs in this species. Here we identified TTX and analogs in the planocerid flatworm using high-resolution liquid chromatography-mass spectrometry, and compared the composition of TTX and analogs with those of another toxic and non-toxic planocerid species. We show that the composition of TTX and several analogs, such as 5,6,11-trideoxyTTX, dideoxyTTXs, deoxyTTXs, and 11-norTTX-6(S)-ol, of Planocera sp. was identical to those of toxic species, but not to its non-toxic counterpart. The difference in the toxin composition was reflected in the phylogenetic relationship based on the mitochondrial genome sequence. A toxification experiment using predatory fish and egg plates of P. multitentaculata demonstrated that the composition of TTX and analogs in wild T. alboplumbeus juveniles was reproduced in artificially toxified pufferfish. Additionally, feeding on the flatworm egg plates enhanced the signal intensities of all TTX compounds in Chelonodon patoca and that of deoxyTTXs in Yongeichthys criniger.

Mitonuclear interactions impact aerobic metabolism in hybrids and may explain mitonuclear discordance in young, naturally hybridizing bird lineages

Understanding genetic incompatibilities and genetic introgression between incipient species are major goals in evolutionary biology. Mitochondrial genes evolve rapidly and exist in dense gene networks with coevolved nuclear genes, suggesting that mitochondrial respiration may be particularly susceptible to disruption in hybrid organisms. Mitonuclear interactions have been demonstrated to contribute to hybrid dysfunction between deeply divergent taxa crossed in the laboratory, but there are few empirical examples of mitonuclear interactions between younger lineages that naturally hybridize. Here, we use controlled hybrid crosses and high-resolution respirometry to provide the first experimental evidence in a bird that inter-lineage mitonuclear interactions impact mitochondrial aerobic metabolism. Specifically, respiration capacity of the two mitodiscordant backcrosses (with mismatched mitonuclear combinations) differs from one another, although they do not differ significantly from the parental groups or mitoconcordant backcrosses as we would expect of mitonuclear disruptions. In the wild hybrid zone between these subspecies, the mitochondrial cline centre is shifted west of the nuclear cline centre, which is consistent with the direction of our experimental results. Our results therefore demonstrate asymmetric mitonuclear interactions that impact the capacity of cellular mitochondrial respiration and may help to explain the geographic discordance between mitochondrial and nuclear genomes observed in the wild.

Characteristics and Comparative Analysis of Mitochondrial Genomes of the Aphid Genus Hyalopterus Koch (Hemiptera: Aphididae: Aphidinae)

Using Illumina sequencing technology, we generated complete mitochondrial genomes (mitogenomes) of three constituent species of the aphid genus Hyalopterus Koch, Hyalopterus amygdali (Blanchard), Hyalopterus arundiniformis Ghulamullah, and Hyalopterus pruni (Geoffroy). The sizes of the Hyalopterus mitogenomes range from 15,306 to 15,410 bp, primarily due to variations in the length of non-coding regions. The Hyalopterus mitogenomes consist of 37 coding genes arranged in the order of the ancestral insect mitogenome, a control region, and a repeat region between trnE and trnF. According to the COI-based analysis, one previously reported mitogenome of H. pruni should be assigned to H. arundiniformis. The gene order, nucleotide composition, and codon usage in the Hyalopterus mitogenomes are highly conserved and similar to those of other species of Aphidinae. The tandem repeat units differ in nucleotide composition, length, and copy number across three Hyalopterus species. Within the widespread Eurasian species H. arundiniformis, variation in repeat units among different geographic populations is observed, indicating that the repeat region may provide valuable insights for studying the intraspecific diversification of aphids. Phylogenetic analyses based on 28 complete mitogenomes of Aphidinae supported the monophyly of Aphidinae, Aphidini, Macrosiphini, and two subtribes of Aphidini. Hyalopterus was monophyletic. H. amygdali and H. pruni formed a sister group, while H. arundiniformis was placed basally. Characterization of the mitogenomes of Hyalopterus provides valuable resources for further comparative studies and for advancing our understanding of the aphid mitogenome architecture.

Mitochondrial metagenomics reveal the independent colonization of the world's coasts by intertidal oribatid mites (Acari, Oribatida, Ameronothroidea)

Oribatid mites are an ancient group that already roamed terrestrial ecosystems in the early and middle Devonian. The superfamily of Ameronothroidea, a supposedly monophyletic lineage, represents the only group of oribatid mites that has successfully invaded the marine coastal environment. By using mitogenome data and nucleic ribosomal RNA genes (18S, 5.8S, 28S), we show that Ameronothroidea are a paraphyletic assemblage and that the land-to-sea transition happened three times independently. Common ancestors of the tropical Fortuyniidae and Selenoribatidae were the first to colonize the coasts and molecular calibration of our phylogeny dates this event to a period in the Triassic and Jurassic era (225-146 mya), whereas present-day distribution indicates that this event might have happened early in this period during the Triassic, when the supercontinent Pangaea still existed. The cold temperate northern hemispheric Ameronothridae colonized the marine littoral later in the late Jurassic-Early Cretaceous and had an ancient distribution on Laurasian coasts. The third and final land-to-sea transition happened in the same geological period, but approx. 30 my later when ancestors of Podacaridae invaded coastal marine environments of the Gondwanan landmasses.

Different Evolutionary Trends of Galloanseres: Mitogenomics Analysis

The two existing clades of Galloanseres, orders Galliformes (landfowl) and Anseriformes (waterfowl), exhibit dramatically different evolutionary trends. Mitochondria serve as primary sites for energy production in organisms, and numerous studies have revealed their role in biological evolution and ecological adaptation. We assembled the complete mitogenome sequences of two species of the genus Aythya within Anseriformes: Aythya baeri and Aythya marila. A phylogenetic tree was constructed for 142 species within Galloanseres, and their divergence times were inferred. The divergence between Galliformes and Anseriformes occurred ~79.62 million years ago (Mya), followed by rapid evolution and diversification after the Middle Miocene (~13.82 Mya). The analysis of selective pressure indicated that the mitochondrial protein-coding genes (PCGs) of Galloanseres species have predominantly undergone purifying selection. The free-ratio model revealed that the evolutionary rates of COX1 and COX3 were lower than those of the other PCGs, whereas ND2 and ND6 had faster evolutionary rates. The CmC model also indicated that most PCGs in Anseriformes exhibited stronger selective constraints. Our study suggests that the distinct evolutionary trends and energy requirements of Galliformes and Anseriformes drive different evolutionary patterns in the mitogenome.

An atlas of spider development at single-cell resolution provides new insights into arthropod embryogenesis

Spiders are a diverse order of chelicerates that diverged from other arthropods over 500 million years ago. Research on spider embryogenesis, particularly studies using the common house spider Parasteatoda tepidariorum, has made important contributions to understanding the evolution of animal development, including axis formation, segmentation, and patterning. However, we lack knowledge about the cells that build spider embryos, their gene expression profiles and fate. Single-cell transcriptomic analyses have been revolutionary in describing these complex landscapes of cellular genetics in a range of animals. Therefore, we carried out single-cell RNA sequencing of P. tepidariorum embryos at stages 7, 8 and 9, which encompass the establishment and patterning of the body plan, and initial differentiation of many tissues and organs. We identified 20 cell clusters, from 18.5 k cells, which were marked by many developmental toolkit genes, as well as a plethora of genes not previously investigated. We found differences in the cell cycle transcriptional signatures, suggestive of different proliferation dynamics, which related to distinctions between endodermal and some mesodermal clusters, compared with ectodermal clusters. We identified many Hox genes as markers of cell clusters, and Hox gene ohnologs were often present in different clusters. This provided additional evidence of sub- and/or neo-functionalisation of these important developmental genes after the whole genome duplication in an arachnopulmonate ancestor (spiders, scorpions, and related orders). We also examined the spatial expression of marker genes for each cluster to generate a comprehensive cell atlas of these embryonic stages. This revealed new insights into the cellular basis and genetic regulation of head patterning, hematopoiesis, limb development, gut development, and posterior segmentation. This atlas will serve as a platform for future analysis of spider cell specification and fate, and studying the evolution of these processes among animals at cellular resolution.

From Caves to the Savannah, the Mitogenome History of Modern Lions (Panthera leo) and Their Ancestors

Lions (Panthera leo) play a crucial ecological role in shaping and maintaining fragile ecosystems within Africa. Conservation efforts should focus on genetic variability within wild populations when considering reintroduction attempts. We studied two groups of lions from two conservation sites located in Zambia and Zimbabwe to determine their genetic make-up, information that is usually unknown to the sites. In this study, we analysed 17 specimens for cytb and seven microsatellite markers to ascertain family relationships and genetic diversity previously obtained by observational studies. We then produced a standardised haplogroup phylogeny using all available entire mitogenomes, as well as calculating a revised molecular clock. The modern lion lineage diverged ~151 kya and was divided into two subspecies, both containing three distinct haplogroups. We confirm that Panthera leo persica is not a subspecies, but rather a haplogroup of the northern P.l. leo that exited Africa at least ~31 kya. The progenitor to all lions existed ~1.2 Mya, possibly in SE Africa, and later exited Africa and split into the two cave lion lineages ~175 kya. Species demography is correlated to major climactic events. We now have a detailed phylogeny of lion evolution and an idea of their conservation status given the threat of climate change.

The contribution of an X chromosome QTL to non-Mendelian inheritance and unequal chromosomal segregation in Auanema freiburgense

Auanema freiburgense is a nematode with males, females, and selfing hermaphrodites. When XO males mate with XX females, they typically produce a low proportion of XO offspring because they eliminate nullo-X spermatids. This process ensures that most sperm carry an X chromosome, increasing the likelihood of X chromosome transmission compared to random segregation. This occurs because of an unequal distribution of essential cellular organelles during sperm formation, likely dependent on the X chromosome. Some sperm components are selectively segregated into the X chromosome's daughter cell, while others are discarded with the nullo-X daughter cell. Intriguingly, the interbreeding of 2 A. freiburgense strains results in hybrid males capable of producing viable nullo-X sperm. Consequently, when these hybrid males mate with females, they yield a high percentage of male offspring. To uncover the genetic basis of nullo-spermatid elimination and X chromosome drive, we generated a genome assembly for A. freiburgense and genotyped the intercrossed lines. This analysis identified a quantitative trait locus spanning several X chromosome genes linked to the non-Mendelian inheritance patterns observed in A. freiburgense. This finding provides valuable clues to the underlying factors involved in asymmetric organelle partitioning during male meiotic division and thus non-Mendelian transmission of the X chromosome and sex ratios.

Genome Evolution and Introgression in the New Zealand mud Snails Potamopyrgus estuarinus and Potamopyrgus kaitunuparaoa

We have sequenced, assembled, and analyzed the nuclear and mitochondrial genomes and transcriptomes of Potamopyrgus estuarinus and Potamopyrgus kaitunuparaoa, two prosobranch snail species native to New Zealand that together span the continuum from estuary to freshwater. These two species are the closest known relatives of the freshwater species Potamopyrgus antipodarum-a model for studying the evolution of sex, host-parasite coevolution, and biological invasiveness-and thus provide key evolutionary context for understanding its unusual biology. The P. estuarinus and P. kaitunuparaoa genomes are very similar in size and overall gene content. Comparative analyses of genome content indicate that these two species harbor a near-identical set of genes involved in meiosis and sperm functions, including seven genes with meiosis-specific functions. These results are consistent with obligate sexual reproduction in these two species and provide a framework for future analyses of P. antipodarum-a species comprising both obligately sexual and obligately asexual lineages, each separately derived from a sexual ancestor. Genome-wide multigene phylogenetic analyses indicate that P. kaitunuparaoa is likely the closest relative to P. antipodarum. We nevertheless show that there has been considerable introgression between P. estuarinus and P. kaitunuparaoa. That introgression does not extend to the mitochondrial genome, which appears to serve as a barrier to hybridization between P. estuarinus and P. kaitunuparaoa. Nuclear-encoded genes whose products function in joint mitochondrial-nuclear enzyme complexes exhibit similar patterns of nonintrogression, indicating that incompatibilities between the mitochondrial and the nuclear genome may have prevented more extensive gene flow between these two species.

High-Quality, Chromosome-Level Reference Genomes of the Viviparous Caribbean Skinks Spondylurus nitidus and S. culebrae

New World mabuyine skinks are a diverse radiation of morphologically cryptic lizards with unique reproductive biologies. Recent studies examining population-level data (morphological, ecological, and genomic) have uncovered novel biodiversity and phenotypes, including the description of dozens of new species and insights into the evolution of their highly complex placental structures. Beyond the potential for this diverse group to serve as a model for the evolution of viviparity in lizards, much of the taxonomic diversity is concentrated in regions experiencing increasing environmental instability from climate and anthropogenic change. Consequently, a better understanding of genome structure and diversity will be an important tool in the adaptive management and conservation of this group. Skinks endemic to Caribbean islands are particularly vulnerable to global change with several species already considered likely extinct and several remaining species either endangered or threatened. Combining PacBio long-read sequencing, Hi-C, and RNAseq data, here we present the first genomic resources for this group by describing new chromosome-level reference genomes for the Puerto Rican Skink Spondylurus nitidus and the Culebra Skink S. culebrae. Results indicate two high quality genomes, both ∼1.4 Gb, assembled nearly telomere to telomere with complete mitochondrion assembly and annotation.

Mitochondrial genome fragmentation is correlated with increased rates of molecular evolution

While mitochondrial genome content and organization is quite diverse across all Eukaryotes, most bilaterian animal mitochondrial genomes (mitogenomes) exhibit highly conserved gene content and organisation, with genes typically encoded on a single circular chromosome. However, many species of parasitic lice (Insecta: Phthiraptera) are among the notable exceptions, having mitogenomes fragmented into multiple circular chromosomes. To better understand the process of mitogenome fragmentation, we conducted a large-scale genomic study of a major group of lice, Amblycera, with extensive taxon sampling. Analyses of the evolution of mitogenome structure across a phylogenomic tree of 90 samples from 53 genera revealed evidence for multiple independent origins of mitogenome fragmentation, some inferred to have occurred less than five million years ago. We leveraged these many independent origins of fragmentation to compare the rates of DNA substitution and gene rearrangement, specifically contrasting branches with fragmented and non-fragmented mitogenomes. We found that lineages with fragmented mitochondrial genomes had significantly higher rates of mitochondrial sequence evolution. In addition, lineages with fragmented mitochondrial genomes were more likely to have mitogenome gene rearrangements than those with single-chromosome mitochondrial genomes. By combining phylogenomics and mitochondrial genomics we provide a detailed portrait of mitogenome evolution across this group of insects with a remarkably unstable mitogenome structure, identifying processes of molecular evolution that are correlated with mitogenome fragmentation.

Two Lineages of Papillomaviruses Identified from Caracals (Caracal caracal) in South Africa

Papillomaviruses (PV) infect epithelial cells and can cause hyperplastic or neoplastic lesions. In felids, most described PVs are from domestic cats (Felis catus; n = 7 types), with one type identified in each of the five wild felid species studied to date (Panthera uncia, Puma concolor, Leopardus wiedii, Panthera leo persica and Lynx rufus). PVs from domestic cats are highly diverse and are currently classified into three genera (Lambdapapillomavirus, Dyothetapapillomavirus, and Taupapillomavirus), whereas those from wild felids, although diverse, are all classified into the Lambdapapillomavirus genus. In this study, we used a metagenomic approach to identify ten novel PV genomes from rectal swabs of five deceased caracals (Caracal caracal) living in the greater Cape Town area, South Africa. These are the first PVs to be described from caracals, and represent six new PV types, i.e., Caracal caracal papillomavirus (CcarPV) 1-6. These CcarPV fall into two phylogenetically distinct genera: Lambdapapillomavirus, and Treisetapapillomavirus. Two or more PV types were identified in a single individual for three of the five caracals, and four caracals shared at least one of the same PV types with another caracal. This study broadens our understanding of wild felid PVs and provides evidence that there may be several wild felid PV lineages.

Complete Mitochondrial Genome for Lucilia cuprina dorsalis (Diptera: Calliphoridae) from the Northern Territory, Australia

The Australian sheep blowfly, Lucilia cuprina dorsalis, is a major sheep ectoparasite causing subcutaneous myiasis (flystrike), which can lead to reduced livestock productivity and, in severe instances, death of the affected animals. It is also a primary colonizer of carrion, an efficient pollinator, and used in maggot debridement therapy and forensic investigations. In this study, we report the complete mitochondrial (mt) genome of L. c. dorsalis from the Northern Territory (NT), Australia, where sheep are prohibited animals, unlike the rest of Australia. The mt genome is 15,943 bp in length, comprising 13 protein-coding genes (PCGs), two ribosomal RNAs (rRNAs), 22 transfer RNAs (tRNAs), and a non-coding control region. The gene order of the current mt genome is consistent with the previously published L. cuprina mt genomes. Nucleotide composition revealed an AT bias, accounting for 77.5% of total mt genome nucleotides. Phylogenetic analyses of 56 species/taxa of dipterans indicated that L. c. dorsalis and L. sericata are the closest among all sibling species of the genus Lucilia, which helps to explain species evolution within the family Luciliinae. This study provides the first complete mt genome sequence for L. c. dorsalis derived from the NT, Australia to facilitate species identification and the examination of the evolutionary history of these blowflies.

Microbiome and mitogenomics of the chigger mite Pentidionis agamae: potential role as an Orientia vector and associations with divergent clades of Wolbachia and Borrelia

BACKGROUND

Trombiculid mites are globally distributed, highly diverse arachnids that largely lack molecular resources such as whole mitogenomes for the elucidation of taxonomic relationships. Trombiculid larvae (chiggers) parasitise vertebrates and can transmit bacteria (Orientia spp.) responsible for scrub typhus, a zoonotic febrile illness. Orientia tsutsugamushi causes most cases of scrub typhus and is endemic to the Asia-Pacific Region, where it is transmitted by Leptotrombidium spp. chiggers. However, in Dubai, Candidatus Orientia chuto was isolated from a case of scrub typhus and is also known to circulate among rodents in Saudi Arabia and Kenya, although its vectors remain poorly defined. In addition to Orientia, chiggers are often infected with other potential pathogens or arthropod-specific endosymbionts, but their significance for trombiculid biology and public health is unclear.

RESULTS

Ten chigger species were collected from rodents in southwestern Saudi Arabia. Chiggers were pooled according to species and screened for Orientia DNA by PCR. Two species (Microtrombicula muhaylensis and Pentidionis agamae) produced positive results for the htrA gene, although Ca. Orientia chuto DNA was confirmed by Sanger sequencing only in P. agamae. Metagenomic sequencing of three pools of P. agamae provided evidence for two other bacterial associates: a spirochaete and a Wolbachia symbiont. Phylogenetic analysis of 16S rRNA and multi-locus sequence typing genes placed the spirochaete in a clade of micromammal-associated Borrelia spp. that are widely-distributed globally with no known vector. For the Wolbachia symbiont, a genome assembly was obtained that allowed phylogenetic localisation in a novel, divergent clade. Cytochrome c oxidase I (COI) barcodes for Saudi Arabian chiggers enabled comparisons with global chigger diversity, revealing several cases of discordance with classical taxonomy. Complete mitogenome assemblies were obtained for the three P. agamae pools and almost 50 SNPs were identified, despite a common geographic origin.

CONCLUSIONS

P. agamae was identified as a potential vector of Ca. Orientia chuto on the Arabian Peninsula. The detection of an unusual Borrelia sp. and a divergent Wolbachia symbiont in P. agamae indicated links with chigger microbiomes in other parts of the world, while COI barcoding and mitogenomic analyses greatly extended our understanding of inter- and intraspecific relationships in trombiculid mites.

The complete mitochondrial genome of Phymateus saxosus (Coquerel, 1861) (Orthoptera: Pyrgomorphidae) and phylogenetic analysis

Phymateus saxosus is a member of the family Pyrgomorphidae, Orthoptera. In this study, the complete mitochondrial genome (mitogenome) of P. saxosus was determined and analyzed. Assembled mitogenome sequence of P. saxosus is 15,672 bp in size, containing 37 genes and a control region. The gene orientation and arrangement of P. saxosus are identical to other species in the Pyrgomorphoidea family. The overall nucleotide composition is as follows: A (43.6%) > T (30.2%) > C (16.1%) > G (10.1%). Phylogenetic analysis suggested that P. saxosus forms sister groups with P. morbillosus, and the monophyly of Pyrgomorphidae is supported. In general, this study provided valuable genetic information for P. saxosus and explored the phylogenetic relationships in the family Pyrgomorphidae.

The complete mitochondrial genome of Meller's mongoose (Rhynchogale melleri)

Meller's mongoose (Rhynchogale melleri) is a member of the family Herpestidae (Mammalia: Carnivora) and the sole species in the genus Rhynchogale. It is primarily found in savannas and open woodlands of eastern sub-Saharan Africa. Here, we report the first complete mitochondrial genome for a female Meller's mongoose collected in Tanzania, generated using a genome-skimming approach. The mitogenome had a final length of 16,644 bp and a total of 37 annotated genes. Phylogenetic analysis validated the placement of this species in the herpestid subfamily Herpestinae. Ultimately, the outcomes of this research offer a genetic foundation for future studies of Meller's mongoose.

The genome of the rayed Mediterranean limpet Patella caerulea (Linnaeus, 1758)

Patella caerulea (Linnaeus, 1758) is a mollusc limpet species of the class Gastropoda. Endemic to the Mediterranean Sea, it is considered a keystone species due to its primary role in structuring and regulating the ecological balance of tidal and subtidal habitats. It is currently being used as a bioindicator to assess the environmental quality of coastal marine waters and as a model species to understand adaptation to ocean acidification. Here, we provide a high-quality reference genome assembly and annotation for P. caerulea. We generated ∼30 Gb of Pacific Biosciences high-fidelity data from a single individual and provide a final 749.8 Mb assembly containing 62 contigs, including the mitochondrial genome (14,938 bp). With an N50 of 48.8 Mb and 98% of the assembly contained in the 18 largest contigs, this assembly is near chromosome-scale. Benchmarking Universal Single-Copy Orthologs scores were high (Mollusca, 87.8% complete; Metazoa, 97.2% complete) and similar to metrics observed for other chromosome-level Patella genomes, highlighting a possible bias in the Mollusca database for Patellids. We generated transcriptomic Illumina data from a second individual collected at the same locality and used it together with protein evidence to annotate the genome. A total of 23,938 protein-coding gene models were found. By comparing this annotation with other published Patella annotations, we found that the distribution and median values of exon and gene lengths was comparable with other Patella species despite different annotation approaches. The present high-quality P. caerulea reference genome, available on GenBank (BioProject: PRJNA1045377; assembly: GCA_036850965.1), is an important resource for future ecological and evolutionary studies.

Boring systematics: A genome skimmed phylogeny of ctenostome bryozoans and their endolithic family Penetrantiidae with the description of one new species

Ctenostomes are a group of gymnolaemate bryozoans with an uncalcified chitinous body wall having few external, skeletal characters. Hence, species identification is challenging and their systematics remain poorly understood, even more so when they exhibit an endolithic (boring) lifestyle. Currently, there are four Recent families of endolithic bryozoans that live inside mineralized substrates like mollusk shells. In particular, Penetrantiidae Silén, 1946 has received considerable attention and its systematic affinity to either cheilostomes or ctenostomes has been debated. Species delimitation of penetrantiids remains difficult, owing to a high degree of colonial and zooidal plasticity. Consequently, an additional molecular approach is essential to unravel the systematics of penetrantiids, their phylogenetic placement and their species diversity. We therefore sequenced the mitochondrial (mt) genomes and two nuclear markers of 27 ctenostome species including nine penetrantiids. Our phylogeny supports the Penetrantiidae as a monophyletic group placed as sister taxon to the remaining ctenostomes alongside paludicellids, arachnidioids and terebriporids. The boring family Terebriporidae d'Orbigny, 1847 were previously considered to be among vesicularioids, but our results suggest an arachnidioid affinity instead. Ctenostome paraphyly is supported by our data, as the cheilostomes nest within them. A Multiporata clade is also well supported, including the former victorelloid genus Sundanella. Altogether, this study provides new insights into ctenostome systematics, assists with species delimitation and contributes to our understanding of the bryozoan tree of life.

The first complete mitochondrial genome data of the Afghan pika Ochotona rufescens (Lagomorpha, Ochotonidae), near the type locality

The Afghan pika Ochotona rufescens (Gray, 1842) is widely distributed across the mountains of Afghanistan, Iran, Pakistan, and southwestern Turkmenistan, most often at elevations between 2,000 and 3,000 m. Here we present, for the first time, the complete mitochondrial genomes of two specimens of the nominotypical subspecies Ochotona rufescens rufescens, de novo assembled from Illumina short reads of fragmented probe-enriched DNA. The lengths of the circular mitogenomes are 16,408 bp and 16,407 bp, respectively. Both mitogenomes contain 13 protein-coding genes (PCGs), two ribosomal RNAs (16S rRNA and 12S rRNA), 22 transfer RNA genes, and a control region. The gene NAD6 and the tRNA (Gln), tRNA (Ala), tRNA (Asn), tRNA (Cys), tRNA (Tyr), tRNA (Ser), tRNA (Glu), and tRNA (Pro) are encoded on the light strand while the rest are encoded on the heavy strand. The overall nucleotide composition was ∼30% for A, 25% for T, 15% for G, and 29% for C. The mitogenome data are available in the GenBank under the accession numbers ON859136 and ON859137.

Molecular Identification of Acrossocheilus jishouensis (Teleostei: Cyprinidae) and Its Complete Mitochondrial Genome

Molecular identification, such as DNA barcoding, is a useful tool that is widely applied in distinguishing species. To identify the cyprinid Acrossocheilus jishouensis, which was previously known to be restricted to only its type locality, we conducted molecular identification of this species based on 23 samples in five localities. Molecular identification based on the mitochondrial COI gene sequence showed that the morphologically similar samples from the five populations were all A. jishouensis, as the mean genetic distances between populations were very small (0.1-1.6%); thus, the distribution of this species was substantially expanded. The whole mitochondrial genome of one sample was also assembled, which was 16,594 bp in length and consisted of 13 protein-coding genes (PCGs), two rRNA genes, 22 tRNA genes, and one control region. All PCGs began with ATG except the COI gene, which started with GTG; seven PCGs used the complete stop codon TAA, while four terminated in T(AA) and two ended with TAG. The overall base composition reflected a higher proportion of A+T than G+C and a positive AT-skew and negative GC-skew pattern except for the opposite in ND6. Phylogenetic relationships inferred using BI and ML methods revealed that both Acrossocheilus and Onychostoma were nonmonophyletic, which indicated that the traditional diagnoses between these two genera need to be assessed further. The results of this study not only expanded the known distribution ranges of A. jishouensis, but also provided a valuable data resource for future molecular and evolutionary studies of Acrossocheilus and other cyprinids in Barbinae.

The updated genome of the Hungarian population of Aedes koreicus

Vector-borne diseases pose a potential risk to human and animal welfare, and understanding their spread requires genomic resources. The mosquito Aedes koreicus is an emerging vector that has been introduced into Europe more than 15 years ago but only a low quality, fragmented genome was available. In this study, we carried out additional sequencing and assembled and characterized the genome of the species to provide a background for understanding its evolution and biology. The updated genome was 1.1 Gbp long and consisted of 6099 contigs with an N50 value of 329,610 bp and a BUSCO score of 84%. We identified 22,580 genes that could be functionally annotated and paid particular attention to the identification of potential insecticide resistance genes. The assessment of the orthology of the genes indicates a high turnover at the terminal branches of the species tree of mosquitoes with complete genomes, which could contribute to the adaptation and evolutionary success of the species. These results could form the basis for numerous downstream analyzes to develop targets for the control of mosquito populations.

Complete Mitochondrial Genome of the Eggplant Fruit and Shoot Borer, Leucinodes orbonalis Guenée (Lepidoptera: Crambidae), and Comparison with Other Pyraloid Moths

The eggplant fruit and shoot borer (EFSB) (Leucinodes orbonalis Guenée) is a devastating lepidopteran pest of eggplant (Solanum melongena L.) in the Philippines. Management of an insect pest like the EFSB requires an understanding of its biology, evolution, and adaptations. Genomic resources provide a starting point for understanding EFSB biology, as the resources can be used for phylogenetics and population structure studies. To date, genomic resources are scarce for EFSB; thus, this study generated its complete mitochondrial genome (mitogenome). The circular mitogenome is 15,244 bp-long. It contains 37 genes, namely 13 protein-coding, 22 tRNA, and 2 rRNA genes, and has conserved noncoding regions, motifs, and gene syntenies characteristic of lepidopteran mitogenomes. Some protein-coding genes start and end with non-canonical codons. The tRNA genes exhibit a conserved cloverleaf structure, with the exception in trnS1. Partitioned phylogenetic analysis using 72 pyraloids generated highly supported maximum likelihood and Bayesian inference trees revealing expected basal splits between Crambidae and Pyralidae, and Spilomelinae and Pyraustinae. Spilomelinae was recovered to be paraphyletic, with the EFSB robustly placed before the split of Spilomelinae and Pyraustinae. Overall, the EFSB mitogenome resource will be useful for delineations within Spilomelinae and population structure analysis.

First mitogenomic characterization of Macromotettixoides (Orthoptera, Tetrigidae), with the descriptions of two new species

Classification of species is commonly based on morphological, molecular, and distribution depending on the taxa. Macromotettixoides Zheng, Wei & Jiang, 2005 and Formosatettix Tinkham, 1937 are both wingless types of Tetrigidae with extremely similar morphological characteristics, and in the current taxonomic system they are placed in two different subfamilies, Metrodorinae and Tetriginae, respectively. It is difficult to clearly identify the species of these two genera by morphological characteristics, and molecular data is often needed to assist identification. Here, the complete mitogenomes of two new species were sequenced and assembled, with that of Macromotettixoidesorthomargina. Molecular data of species of Formosatettix were used to test the monophyly of Macromotettixoides and to re-assess the generic characters, and also to test whether Macromotettixoides belongs to the Asian Metrodorinae or Tetriginae. Furthermore, mitochondrial characteristics were analyzed and the phylogeny of the Tetrigidae reconstructed based on mitochondrial protein-coding genes (PCGs). The results indicated that the two new species were clustered with Macromotettixoides rather than Formosatettix, and the anterior margin of the fastigium and pronotum of the two new species usually had the humeral angle different from that of Formosatettix. Therefore, after integrating morphological and molecular data, the two new species were placed in the genus Macromotettixoides, M.maoershanensissp. nov. and M.brachycornasp. nov. Finally, a phylogenetic reconstruction supported Macromotettixoides being assigned to Tetriginae rather than Metrodorinae, in contrast to the previous classification of this genus.

The complete mitochondrial genome of Rhynchocinetes brucei Okuno 1994 (Decapoda: Rhynchocinetidae)

We report the complete mitochondrial genome of Rhynchocinetes brucei Okuno 1994. The mitogenome was found to contain 16158 bp with 13 protein-coding genes (PCGs), 22 tRNA genes (tRNAs), 2 rRNA genes (rRNAs), and 1 putative control region. Phylogenetic analysis indicated that R. brucei was sister to Rhynchocinetes durbanensis (PP= 1), of the same family Rhynchocinetidae. These results are helpful for research on the phylogenetic and evolutionary studies of this group.

Complete mitogenome of Lepidopygopsis typus, an evolutionarily-distinct, endangered cyprinid fish from the Western Ghats Biodiversity Hotspot: Phylogenetic relationships and implications for conservation

The mitogenome of Lepidopygopsis typus, an evolutionarily distinct, endangered, cyprinid fish from the Western Ghats Biodiversity Hotspot, was characterized. Total length of the mitogenome was 16,729 bp, and comprised of 13 protein coding, 22 tRNA and two rRNA genes. Thirteen protein coding genes contributed a total nucleotide length of 11,408 bp, which coded for a total of 3794 amino acids and accounting for 68.2 % of the mitogenome. A maximum likelihood phylogenetic tree based on mitogenomes of members of cyprinid subfamilies Torinae, Schizothoracinae and Schizopygopsinae, recovered L. typus in a well-resolved clade of Torinae with members of the Western Ghats endemic genus Hypselobarbus as sister taxa. Selection analysis suggested that the branch of L. typus and Hypselobarbus spp. was under intensified selection with reference to other members of Torinae. There were 19 codons under diversifying selection in L. typus, which could be the result of positive selection for adapting to high-altitude, upstream tributaries of Periyar River, where the species has a restricted distribution. With respect to Hypselobarbus spp., L. typus was under relaxed selection with about 68 % of codons experiencing neutral evolution. Restricted distribution, low population size and relaxation in selection can likely trigger extinction in L. typus, and therefore urgent conservation and monitoring plans are required to secure the future of this evolutionary distinct, and globally endangered species.

Rediscovering the unusual, solitary bryozoan Monobryozoon ambulans Remane, 1936: first molecular and new morphological data clarify its phylogenetic position

BACKGROUND

One of the most peculiar groups of the mostly colonial phylum Bryozoa is the taxon Monobryozoon, whose name already implies non-colonial members of the phylum. Its peculiarity and highly unusual lifestyle as a meiobenthic clade living on sand grains has fascinated many biologists. In particular its systematic relationship to other bryozoans remains a mystery. Despite numerous searches for M. ambulans in its type locality Helgoland, a locality with a long-lasting marine station and tradition of numerous courses and workshops, it has never been reencountered until today. Here we report the first observations of this almost mythical species, Monobryozoon ambulans.

RESULTS

For the first time since 1938, we present new modern, morphological analyses of this species as well as the first ever molecular data. Our detailed morphological analysis confirms most previous descriptions, but also ascertains the presence of special ambulatory polymorphic zooids. We consider these as bud anlagen that ultimately consecutively separate from the animal rendering it pseudo-colonial. The remaining morphological data show strong ties to alcyonidioidean ctenostome bryozoans. Our morphological data is in accordance with the phylogenomic analysis, which clusters it with species of Alcyonidium as a sister group to multiporate ctenostomes. Divergence time estimation and ancestral state reconstruction recover the solitary state of M. ambulans as a derived character that probably evolved in the Late Cretaceous. In this study, we also provide the entire mitogenome of M. ambulans, which-despite the momentary lack of comparable data-provides important data of a unique and rare species for comparative aspects in the future.

CONCLUSIONS

We were able to provide first sequence data and modern morphological data for the unique bryozoan, M. ambulans, which are both supporting an alcyonidioidean relationship within ctenostome bryozoans.

De novo assembly and annotation of the Patagonian toothfish (Dissostichus eleginoides) genome

BACKGROUND

Patagonian toothfish (Dissostichus eleginoides) is an economically and ecologically important fish species in the family Nototheniidae. Juveniles occupy progressively deeper waters as they mature and grow, and adults have been caught as deep as 2500 m, living on or in just above the southern shelves and slopes around the sub-Antarctic islands of the Southern Ocean. As apex predators, they are a key part of the food web, feeding on a variety of prey, including krill, squid, and other fish. Despite its importance, genomic sequence data, which could be used for more accurate dating of the divergence between Patagonian and Antarctic toothfish, or establish whether it shares adaptations to temperature with fish living in more polar or equatorial climes, has so far been limited.

RESULTS

A high-quality D. eleginoides genome was generated using a combination of Illumina, PacBio and Omni-C sequencing technologies. To aid the genome annotation, the transcriptome derived from a variety of toothfish tissues was also generated using both short and long read sequencing methods. The final genome assembly was 797.8 Mb with a N50 scaffold length of 3.5 Mb. Approximately 31.7% of the genome consisted of repetitive elements. A total of 35,543 putative protein-coding regions were identified, of which 50% have been functionally annotated. Transcriptomics analysis showed that approximately 64% of the predicted genes (22,617 genes) were found to be expressed in the tissues sampled. Comparative genomics analysis revealed that the anti-freeze glycoprotein (AFGP) locus of D. eleginoides does not contain any AFGP proteins compared to the same locus in the Antarctic toothfish (Dissostichus mawsoni). This is in agreement with previously published results looking at hybridization signals and confirms that Patagonian toothfish do not possess AFGP coding sequences in their genome.

CONCLUSIONS

We have assembled and annotated the Patagonian toothfish genome, which will provide a valuable genetic resource for ecological and evolutionary studies on this and other closely related species.

Comparative Analysis of the Mitochondrial Genomes of Chloropidae and Their Implications for the Phylogeny of the Family

Chloropidae, commonly known as grass flies, represent the most taxonomically diverse family of Diptera Carnoidea, comprising over 3000 described species worldwide. Previous phylogenetic studies of this family have predominantly relied on morphological characters, with mitochondrial genomes being reported in a few species. This study presents 11 newly sequenced mitochondrial genomes (10 Chloropidae and 1 Milichiidae) and provides the first comprehensive comparative analysis of mitochondrial genomes for Chloropidae. Apart from 37 standard genes and the control region, three conserved intergenic sequences across Diptera Cyclorrhapha were identified in all available chloropid mitochondrial genomes. Evolutionary rates within Chloropidae exhibit significant variation across subfamilies, with Chloropinae displaying higher rates than the other three subfamilies. Phylogenetic relationships based on mitochondrial genomes were inferred using maximum likelihood and Bayesian methods. The monophyly of Chloropidae and all four subfamilies is consistently strongly supported, while subfamily relationships within Chloropidae remain poorly resolved, possibly due to rapid evolution.

Giants among Cnidaria: Large Nuclear Genomes and Rearranged Mitochondrial Genomes in Siphonophores

Siphonophores (Cnidaria: Hydrozoa) are abundant predators found throughout the ocean and are important constituents of the global zooplankton community. They range in length from a few centimeters to tens of meters. They are gelatinous, fragile, and difficult to collect, so many aspects of the biology of these roughly 200 species remain poorly understood. To survey siphonophore genome diversity, we performed Illumina sequencing of 32 species sampled broadly across the phylogeny. Sequencing depth was sufficient to estimate nuclear genome size from k-mer spectra in six specimens, ranging from 0.7 to 2.3 Gb, with heterozygosity estimates between 0.69% and 2.32%. Incremental k-mer counting indicates k-mer peaks can be absent with nearly 20× read coverage, suggesting minimum genome sizes range from 1.4 to 5.6 Gb in the 25 samples without peaks in the k-mer spectra. This work confirms most siphonophore nuclear genomes are large relative to the genomes of other cnidarians, but also identifies several with reduced size that are tractable targets for future siphonophore nuclear genome assembly projects. We also assembled complete mitochondrial genomes for 33 specimens from these new data, indicating a conserved gene order shared among nonsiphonophore hydrozoans, Cystonectae, and some Physonectae, revealing the ancestral mitochondrial gene order of siphonophores. Our results also suggest extensive rearrangement of mitochondrial genomes within other Physonectae and in Calycophorae. Though siphonophores comprise a small fraction of cnidarian species, this survey greatly expands our understanding of cnidarian genome diversity. This study further illustrates both the importance of deep phylogenetic sampling and the utility of k-mer-based genome skimming in understanding the genomic diversity of a clade.

The genome of Salmacisia buchloëana, the parasitic puppet master pulling strings of sexual phenotypic monstrosities in buffalograss

To complete its parasitic lifecycle, Salmacisia buchloëana, a biotrophic fungus, manipulates reproductive organ development, meristem determinacy, and resource allocation in its dioecious plant host, buffalograss (Bouteloua dactyloides; Poaceae). To gain insight into S. buchloëana's ability to manipulate its host, we sequenced and assembled the 20.1 Mb genome of S. buchloëana into 22 chromosome-level pseudomolecules. Phylogenetic analysis suggests that S. buchloëana is nested within the genus Tilletia and diverged from Tilletia caries and Tilletia walkeri ∼40 MYA. We find that S. buchloëana contains a novel chromosome arm with no syntenic relationship to other publicly available Tilletia genomes, and that genes on the novel arm are upregulated upon infection, suggesting that this unique chromosomal segment may have played a critical role in S. buchloëana's evolution and host specificity. Salmacisia buchloëana has one of the largest fractions of serine peptidases (1.53% of the proteome) and one of the highest GC contents (62.3%) in all classified fungi. Analysis of codon base composition indicated that GC content is controlled more by selective constraints than directional mutation, and that S. buchloëana has a unique bias for the serine codon UCG. Finally, we identify 3 inteins within the S. buchloëana genome, 2 of which are located in a gene often used in fungal taxonomy. The genomic and transcriptomic resources generated here will aid plant pathologists and breeders by providing insight into the extracellular components contributing to sex determination in dioecious grasses.

Insights into the nuclear and mitochondrial genome of the Lemon shark Negaprion brevirostris using low-coverage sequencing: Genome size, repetitive elements, mitochondrial genome, and phylogenetic placement

The Lemon shark Negaprion brevirostris is an important species experiencing conservation issues that is in need of genomic resources. Herein, we conducted a genome survey sequencing in N. brevirostris and determined genome size, explored repetitive elements, assembled and annotated the 45S rRNA DNA operon, and assembled and described in detail the mitochondrial genome. Lastly, the phylogenetic position of N. brevirostris in the family Carcharhinidae was examined using translated protein coding genes. The estimated haploid genome size ranged between 2.29 and 2.58 Gbp using a k-mer analysis, which is slightly below the genome size estimated for other sharks belonging to the family Carcharhinidae. Using a k-mer analysis, approx. 64-71 % of the genome of N. brevirostris was composed of repetitive elements. A relatively large proportion of the 'repeatome' could not be annotated. Taking into account only annotated repetitive elements, Class I - Long Interspersed Nuclear Element (LINE) were the most abundant repetitive elements followed by Class I - Penelope and Satellite DNA. The nuclear ribosomal operon was fully assembled. The AT-rich complete mitochondrial genome was 16,703 bp long and encoded 13 protein coding genes, 2 ribosomal RNA genes, and 22 transfer RNA genes. Negaprion brevirostris is closely related to the genera Carcharhinus, Glyphis and Lamiopsis in the family Carcharinidae. This new genomic resources will aid with the development of conservation plans for this large coastal shark.

The Complete Mitochondrial Genome of Chilo infuscatellus (Lepidoptera: Pyralidae), and Related Phylogenetic Analysis

The Chilo infuscatellus (Lepidoptera: Pyralidae) is a significant pest of sugarcane in China. The genome-level characteristics of this pest are important genetic resources for identification, phylogenetic analysis, and even management. In the present study, the complete mitogenome of C. infuscatellus was sequenced and characterized. The assembled mitochondrial genome is 15,252 bp in length and includes 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), 2 ribosomal RNA genes (rRNAs), and an A + T-rich region. Except for the CGA codon for the cox1 gene, the PCGs are initiated with ATN codons (ATG, ATT, and ATA). These PCGs are terminated with TAA or an incomplete termination codon of a single T. Except for the loss of the "DHU" arm for trnS1, the tRNA genes were folded into the typical cloverleaf structure. The A + T-rich region has a high AT content of 96.19% and contains the motifs "ATAGA" and "ATTTA", as well as a 19 bp poly-T stretch and microsatellite regions. The C. infuscatellus mitogenome exhibits a conserved gene order among lepidopteran insects, with a rearrangement of the trnM gene compared to the ancestral insect gene order. Phylogenetic analysis based on the 13 PCGs using Bayesian inference (BI) and maximum likelihood (ML) methods confirmed the monophyly of Pyralidae and Crambidae within Pyraloidea. The relationships between subfamilies in Pyralidae can be described as (Galleriinae + (Phycitinae + (Pyralinae + Epipaschiinae))). The "PS clade" and "non-PS clade" were formed within the family Crambidae. These findings provide valuable genetic resources for the identification, phylogenetic analysis, and management of sugarcane borers, contributing significantly to our understanding of the phylogeny of Pyraloidea insects and their evolution.

The complete mitochondrial genome of Aegialites californicus (Motchoulsky, 1845) (insecta: coleoptera: salpingidae)

The flightless intertidal beetle genus Aegialites (family Salpingidae) is distributed along the Northern Pacific coasts, from California to Alaska and from Northern Japan to Kamchatka. Systematics of Aegialites and its phylogenetic relationships to other members of Salpingidae are unclear, and little genetic information is available. We here present the first complete mitochondrial genome of this genus, represented by Aegialites californicus (Motchoulsky, 1845) from Sonoma County, California, U.S.A. The complete mitochondrial genome of A. californicus is 15,899 bp long and comprises 13 protein-coding (PCG), two ribosomal RNA (rRNA) and 22 transfer RNA (tRNA) genes. The phylogenetic analysis places A.californicus as sister to other members of family Salpingidae. The mitochondrial genome sequence of A. californicus will contribute to future phylogenetic and taxonomic studies of genus Aegialites, family Salpingidae and superfamily Tenebrionoidea.

Demographic dynamics and molecular evolution of the rare and endangered subsect. Gerardianae of Pinus: insights from chloroplast genomes and mitochondrial DNA markers

MAIN CONCLUSION

The divergence of subsect. Gerardianae was likely triggered by the uplift of the Qinghai-Tibetan Plateau and adjacent mountains. Pinus bungeana might have probably experienced expansion since Last Interglacial period. Historical geological and climatic oscillations have profoundly affected patterns of nucleotide variability, evolutionary history, and species divergence in numerous plants of the Northern Hemisphere. However, how long-lived conifers responded to geological and climatic fluctuations in East Asia remain poorly understood. Here, based on paternally inherited chloroplast genomes and maternally inherited mitochondrial DNA markers, we investigated the population demographic history and molecular evolution of subsect. Gerardianae (only including three species, Pinus bungeana, P. gerardiana, and P. squamata) of Pinus. A low level of nucleotide diversity was found in P. bungeana (π was 0.00016 in chloroplast DNA sequences, and 0.00304 in mitochondrial DNAs). The haplotype-based phylogenetic topology and unimodal distributions of demographic analysis suggested that P. bungeana probably originated in the southern Qinling Mountains and experienced rapid population expansion since Last Interglacial period. Phylogenetic analysis revealed that P. gerardiana and P. squamata had closer genetic relationship. The species divergence of subsect. Gerardianae occurred about 27.18 million years ago (Mya) during the middle to late Oligocene, which was significantly associated with the uplift of the Qinghai-Tibetan Plateau and adjacent mountains from the Eocene to the mid-Pliocene. The molecular evolutionary analysis showed that two chloroplast genes (psaI and ycf1) were under positive selection, the genetic lineages of P. bungeana exhibited higher transition and nonsynonymous mutations, which were involved with the strongly environmental adaptation. These findings shed light on the population evolutionary history of white pine species and provide striking insights for comprehension of their species divergence and molecular evolution.

Variability and evolution of gene order rearrangement in mitochondrial genomes of arthropods (except Hexapoda)

In the species-rich Phylum Arthropoda, the mitochondrial genome is relatively well conserved both in terms of number and order of genes. However, specific clades have a 'typical' gene order that differs from the putative arthropod ancestral arrangement. The aim of this work was to compare the rate of mitochondrial gene rearrangements at inter- and intra-taxonomic levels in the Arthropoda and to postulate the most parsimonious ancestral orders representing the four major arthropod lineages. For this purpose, we performed a comparative genomic analysis of arthropod mitochondrial genomes available in the NCBI database. Using a combination of bioinformatics methods that examined mitochondrial gene rearrangements in 464 species of arthropods from three subphyla (Chelicerata, Myriapoda, and Crustacea [except Hexapoda, previously analyzed]), we observed differences in the rate of rearrangement within major lineages. A higher rate of mitochondrial genome rearrangement was observed in Crustacea and Chelicerata compared to Myriapoda. Likewise, early branching clades exhibit less variability in mitochondrial genome order than late branching clades, within each subphylum. We identified 'hot regions' in the mitochondrial genome of each studied subphylum, and postulated the most likely ancestral gene order in each subphylum and taxonomic order. Our work provides new evidence on the evolutionary dynamics of mitochondrial genome gene order in arthropods and new mitochondrial genome architectures in different taxonomic divisions within each major lineage of arthropods.

Insights into the genome of the 'Loco' Concholepas concholepas (Gastropoda: Muricidae) from low-coverage short-read sequencing: genome size, ploidy, transposable elements, nuclear RNA gene operon, mitochondrial genome, and phylogenetic placement in the family Muricidae

BACKGROUND

The Peruvian 'chanque' or Chilean 'loco' Concholepas concholepas is an economically, ecologically, and culturally important muricid gastropod heavily exploited by artisanal fisheries in the temperate southeastern Pacific Ocean. In this study, we have profited from a set of bioinformatics tools to recover important biological information of C. concholepas from low-coverage short-read NGS datasets. Specifically, we calculated the size of the nuclear genome, ploidy, and estimated transposable elements content using an in silico k-mer approach, we discovered, annotated, and quantified those transposable elements, we assembled and annotated the 45S rDNA RNA operon and mitochondrial genome, and we confirmed the phylogenetic position of C. concholepas within the muricid subfamily Rapaninae based on translated protein coding genes.

RESULTS

Using a k-mer approach, the haploid genome size estimated for the predicted diploid genome of C. concholepas varied between 1.83 Gbp (with kmer = 24) and 2.32 Gbp (with kmer = 36). Between half and two thirds of the nuclear genome of C. concholepas was composed of transposable elements. The most common transposable elements were classified as Long Interspersed Nuclear Elements and Short Interspersed Nuclear Elements, which were more abundant than DNA transposons, simple repeats, and Long Terminal Repeats. Less abundant repeat elements included Helitron mobile elements, 45S rRNA DNA, and Satellite DNA, among a few others.The 45S rRNA DNA operon of C. concholepas that encodes for the ssrRNA, 5.8S rRNA, and lsrRNA genes was assembled into a single contig 8,090 bp long. The assembled mitochondrial genome of C. concholepas is 15,449 bp long and encodes 13 protein coding genes, two ribosomal genes, and 22 transfer RNAs.

CONCLUSION

The information gained by this study will inform the assembly of a high quality nuclear genome for C. concholepas and will support bioprospecting and biomonitoring using environmental DNA to advance development of conservation and management plans in this overexploited marine snail.

Molecular phylogenetic relationships based on mitochondrial genomes of novel deep-sea corals (Octocorallia: Alcyonacea): Insights into slow evolution and adaptation to extreme deep-sea environments

A total of 10 specimens of Alcyonacea corals were collected at depths ranging from 905 m to 1 633 m by the manned submersible Shenhai Yongshi during two cruises in the South China Sea (SCS). Based on mitochondrial genomic characteristics, morphological examination, and sclerite scanning electron microscopy, the samples were categorized into four suborders (Calcaxonia, Holaxonia, Scleraxonia, and Stolonifera), and identified as 9 possible new cold-water coral species. Assessments of GC-skew dissimilarity, phylogenetic distance, and average nucleotide identity (ANI) revealed a slow evolutionary rate for the octocoral mitochondrial sequences. The nonsynonymous ( Ka) to synonymous ( Ks) substitution ratio ( Ka/ Ks) suggested that the 14 protein-coding genes (PCGs) were under purifying selection, likely due to specific deep-sea environmental pressures. Correlation analysis of the median Ka/ Ks values of five gene families and environmental factors indicated that the genes encoding cytochrome b (cyt b) and DNA mismatch repair protein ( mutS) may be influenced by environmental factors in the context of deep-sea species formation. This study highlights the slow evolutionary pace and adaptive mechanisms of deep-sea corals.

Complete Mitochondrial Genome and Phylogenetic Analysis of Tarsiger indicus (Aves: Passeriformes: Muscicapidae)

Tarsiger indicus (Vieillot, 1817), the White-browed Bush Robin, is a small passerine bird widely distributed in Asian countries. Here, we successfully sequenced its mitogenome using the Illumina Novaseq 6000 platform (Illumina, San Diego, CA, USA) for PE 2 × 150 bp sequencing. Combined with other published mitogenomes, we conducted the first comprehensive comparative mitogenome analysis of Muscicapidae birds and reconstructed the phylogenetic relationships between Muscicapidae and related groups. The T. indicus mitogenome was 16,723 bp in size, and it possessed the typical avian mitogenome structure and organization. Most PCGs of T. indicus were initiated strictly with the typical start codon ATG, while COX1 and ND2 were started with GTG. RSCU statistics showed that CUA, CGA, and GCC were relatively high frequency in the T. indicus mitogenome. T. cyanurus and T. indicus shared very similar mitogenomic features. All 13 PCGs of Muscicapidae mitogenomes had experienced purifying selection. Specifically, ATP8 had the highest rate of evolution (0.13296), whereas COX1 had the lowest (0.01373). The monophylies of Muscicapidae, Turdidae, and Paradoxornithidae were strongly supported. The clade of ((Muscicapidae + Turdidae) + Sturnidae) in Passeriformes was supported by both Bayesian Inference and Maximum likelihood analyses. The latest taxonomic status of many passerine birds with complex taxonomic histories were also supported. For example, Monticola gularis, T. indicus, and T. cyanurus were allocated to Turdidae in other literature; our phylogenetic topologies clearly supported their membership in Muscicapidae; Paradoxornis heudei, Suthora webbiana, S. nipalensis, and S. fulvifrons were formerly classified into Muscicapidae; we supported their membership in Paradoxornithidae; Culicicapa ceylonensis was originally classified as a member of Muscicapidae; our results are consistent with a position in Stenostiridae. Our study enriches the genetic data of T. indicus and provides new insights into the molecular phylogeny and evolution of passerine birds.

The first record of Egglestonichthys bombylios (Gobiiformes: Gobiidae) from China with its first fresh colouration information

During bottom trawl surveys carried out between 20132021, 52 specimens (33.854.0 mm SL) of Egglestones bumblebee goby Egglestonichthys bombylios were collected at a depth of 1.515 m from Dongshan Bay, Sanmen Bay, and Niushan Island, China. They represent the first records of this species from China. A full description, including fresh colouration of the species is provided as it is poorly known. The individuals collected in China agree with most morphological features of the holotype, except for the pelvic fin fraenum that was not observed or appears to be absent in most specimens. A strong relationship between E. bombylios, Larsonella pumilus, and the genus Priolepis is herein demonstrated by the mitochondrial genome sequences of E. bombylios and twenty closely related species. This study enriches the existing genetic data of the so-called Priolepis lineage and provides useful insights into the phylogenetic relationships across species of the genera Egglestonichthys, Priolepis, and Larsonella.

The complete mitochondrial genome of Slaty-backed Forktail, Enicurus schistaceus (Passeriformes: Muscicapidae) and phylogenetic analysis

The complete mitogenome of the Slaty-backed Forktail (Enicurus schistaceus) was first sequenced using next-generation sequencing. It was 17,112 bp long, with a base composition of 14.17% G, 31.77% C, 30.73% A, and 23.33% T and an AT content of 54.06%. Similar to other mitochondrial genomes within the Muscicapidae family, E. schistaceus exhibited a relatively consistent mitogenome arrangement; it consisted of 22 tRNA genes, two rRNA genes, 13 protein-coding genes, and one control region. Notably, ND6 and eight tRNA genes were encoded on the light strand. Phylogenetic analysis of the 12 Muscicapidae mitogenomes substantiated the monophyly of all genera, including E. schistaceus. Furthermore, the analysis demonstrated a close relationship between Enicurus and Myophonus.

Complete mitochondrial genome sequence of Aureoboletus raphanaceus (Boletales, Basidiomycota)

Aureoboletus raphanaceus is a member of boletoid mushroom, which is named after its distinctive radish smell. The mitochondrial genome and phylogenetic relationships with other boletes need to be investigated to gain a comprehensive understanding of it. In this study, we sequenced the mitochondrial genome of A. raphanaceus using next-generation sequencing technology and found that its mitochondrial genome is a circular DNA molecule measuring 42,157 bp. It consists of 15 core protein-coding genes, 27 transfer RNA genes, and two ribosomal RNA genes. The mitochondrial genome had a base composition of A (39.89%), C (11.06%), G (11.67%), and T (37.38%), with a GC content of 22.73%. A phylogenetic tree based on 22 mitochondrial genomes was constructed, which provided the first insights into the phylogenetic relationships of this species with related boletes.

Novel insight into lepidopteran phylogenetics from the mitochondrial genome of the apple fruit moth of the family Argyresthiidae

BACKGROUND

The order Lepidoptera has an abundance of species, including both agriculturally beneficial and detrimental insects. Molecular data has been used to investigate the phylogenetic relationships of major subdivisions in Lepidoptera, which has enhanced our understanding of the evolutionary relationships at the family and superfamily levels. However, the phylogenetic placement of many superfamilies and/or families in this order is still unknown. In this study, we determine the systematic status of the family Argyresthiidae within Lepidoptera and explore its phylogenetic affinities and implications for the evolution of the order. We describe the first mitochondrial (mt) genome from a member of Argyresthiidae, the apple fruit moth Argyresthia conjugella. The insect is an important pest on apples in Fennoscandia, as it switches hosts when the main host fails to produce crops.

RESULTS

The mt genome of A. conjugella contains 16,044 bp and encodes all 37 genes commonly found in insect mt genomes, including 13 protein-coding genes (PCGs), two ribosomal RNAs, 22 transfer RNAs, and a large control region (1101 bp). The nucleotide composition was extremely AT-rich (82%). All detected PCGs (13) began with an ATN codon and terminated with a TAA stop codon, except the start codon in cox1 is ATT. All 22 tRNAs had cloverleaf secondary structures, except trnS1, where one of the dihydrouridine (DHU) arms is missing, reflecting potential differences in gene expression. When compared to the mt genomes of 507 other Lepidoptera representing 18 superfamilies and 42 families, phylogenomic analyses found that A. conjugella had the closest relationship with the Plutellidae family (Yponomeutoidea-super family). We also detected a sister relationship between Yponomeutoidea and the superfamily Tineidae.

CONCLUSIONS

Our results underline the potential importance of mt genomes in comparative genomic analyses of Lepidoptera species and provide valuable evolutionary insight across the tree of Lepidoptera species.