Mitochondrial Sequence Data Research Articles

Comparative Analysis of Mitogenomes in Leafhopper Tribe Deltocephalini (Hemiptera: Cicadellidae: Deltocephalinae): Structural Conservatism and Phylogeny

ABSTRACTPrevious studies on the gene order and composition of leafhopper mitochondrial genomes have revealed a high level of conservation in overall genome structure. Some members of Deltocephalinae, the largest subfamily, exhibit tRNA gene rearrangements; however, few mitogenomes have been sequenced in this group and the degree of structural variation within tribes remains unclear. In this study, we sequenced the complete mitogenomes of 14 species belonging to four genera of tribe Deltocephalini from China and compared them with the two previously reported mitogenomes for this tribe. The studied mitogenomes showed a high degree of similarity to most other leafhopper mitogenomes in overall structure, mostly varying in the total length (14,961–15,416 bp) and number of non‐coding A + T‐rich regions. Gene size, order, arrangement, base composition, codon usage, and secondary structure of tRNAs in the newly sequenced mitogenomes were highly conserved in Deltocephalini, and variations in start/stop codon usage and tRNA secondary structure mostly matched those of other leafhoppers. Phylogenetic analysis of different combinations of protein‐coding and ribosomal genes using maximum likelihood and Bayesian inference under different models using either amino acid or nucleotide sequences were generally consistent and agreed with the previous nuclear and partial mitochondrial gene sequence data, indicating that complete mitochondrial genomes are phylogenetically informative at different levels of divergence within Deltocephalini and among different leafhoppers species. In addition to Deltocephalini, Deltocephalinae included members of Athysanini and Opsiini formed monophyletic groups. Maximum likelihood and Bayesian inference analyses consistently grouped Graminella nigrinota with Paralimnini, rendering Deltocephalini polyphyletic. The topology consistently divided Deltocephalini into two major branches, with Alobaldia tobae and Polyamia penistenuis forming a well‐supported sister group to the remaining species of the tribe.

Assessment of genomic diversity within and between two cryptic shiners, the West Texas shiner (Notropis megalops) and the Texas shiner (Notropis amabilis).

The presence of cryptic species can hinder effective conservation planning and implementation, as has been the case for speciose groups of freshwater fishes that are difficult to differentiate due to conserved morphologies. The West Texas shiner Notropis megalops and the Texas shiner Notropis amabilis are a cryptic pair of leuciscids (minnows) that co-occur in spring-fed tributaries of the Rio Grande in Texas and Mexico. Both N. megalops and N. amabilis are listed as Species of Greatest Conservation Need by the Texas Parks and Wildlife Department. Notropis amabilis is widespread and listed as apparently secure by Texas Parks and Wildlife Department whereas N. megalops has a very limited distribution and has not been ranked by Texas Parks and Wildlife Department because of data deficiency. Morphological differences between these species have been described; however, proper identification in situ remains problematic. Furthermore, given their range of overlap there is potential for hybridization, and limited genetic data have been collected comparing the species. Therefore, reduced representation genomic and mitochondrial sequencing data were used to reassess the distinctness of the species, screen for hybridization, and characterize their relative frequencies throughout their range of overlap. Genomic analyses recovered two distinct genetic groups corresponding to the species (F'CT = 0.89) with no evidence of admixture or introgression. The species were found to co-occur at three sampling locations, two in the Devils River and one in the Pecos, but not in equal frequencies. Overall, these results provide data and tools for further research on N. megalops needed for accurate conservation policies and management practices.

Oreichthys warjaintia, a new species of cyprinid fish from Meghalaya, India (Cypriniformes: Cyprinidae).

Oreichthys warjaintia, new species, is described from the Pyrngang stream (Surma-Meghna basin), Meghalaya, northeast India. It can be distinguished from all known species within the genus Oreichthys in having the following combination of characters: complete lateral line with 25-27 pored scales, 13 branched pectoral-fin rays, 5-5½ branched anal-fin rays, 14 rows of papillae on suborbital, ½3/1/2½ scale rows in transverse line on body, presence of a black blotch on caudal-fin base, and in the coloration of dorsal, caudal, and anal fins. Oreichthys warjaintia is distinct from all known Oreichthys species for which mitochondrial cox1 sequence data are available by a raw genetic distance ranging from 8.3% to 11.1%.

A new small-sized Calotes Cuvier (Squamata: Agamidae: Draconinae) from the Subansiri river basin, Arunachal Pradesh, India.

We describe a new species of small-sized Calotes from mid-elevation (~1270 m asl.), swidden or shifting cultivation areas in the Subansiri river basin, Upper Subansiri District, Arunachal Pradesh, India. Mitochondrial sequence data (ND2) recovers Calotes sinyik sp. nov. as the closely related sister taxon to an unnamed species from Medog, Tibet; the two falling in a broader clade including C. paulus + C. zolaiking and species of the C. emma and C. mystaceus groups. The new species is 7.8% divergent in ND2 sequence data from the unnamed lineage from Medog and deeply divergent (≥ 21.6-28.1%) from other congeners. The new species can easily be distinguished from regional congeners by its small adult body size (maximum snout to vent length of 65 mm) and heterogenous dorsal scales, and from its closest relatives C. paulus and C. zolaiking by having 54 or 55 midbody scale rows and 42 or 43 vertebral scales. It is likely that many more allied species remain to be discovered from Northeast India, which remains poorly surveyed across taxonomic groups.

Snailfishes of the genus Careproctus (Perciformes: Liparidae) with a reduced pelvic disk: three new species and new records from the western North Pacific with comments on their phenotypic diversity.

A taxonomic review of the snailfish genus Careproctus (Liparidae) with a reduced pelvic disk-the Careproctus gilberti species group-from the western North Pacific recognized three new species from the southern Sea of Okhotsk, and confirmed the validity of Careproctus gilberti, Careproctus mederi, Careproctus ostentum, and Careproctus parvidiscus. One of the newly described species, Careproctus longibarbatus sp. nov., had been previously confused with C. mederi, but was readily distinguishable on the basis of counts of vertebrae (61-66 vs. 57-63), dorsal- (56-60 vs. 52-56) and anal-fin rays (50-54 vs. 47-52), as well as length of the pectoral fin lower lobe (extending past anus vs. not reaching anus) and color of the peritoneum (dark vs. pale with black dots). Careproctus barbatulus sp. nov. and Careproctus spinulosus sp. nov. were diagnosed by pectoral fin lower lobe length, anus position, and vertebral, and dorsal- and anal-fin ray counts. The validities of the three new species were also confirmed by mitochondrial sequence data. Careproctus mederi was redescribed based on specimens from the southern Sea of Okhotsk, representing the first Japanese, as well as southernmost record of the species. Furthermore, the full description of an unidentified Careproctus species, most similar to C. gilberti, from the Pacific coast of Hokkaido was provided. Significant morphological variation within the C. gilberti species group, possibly resulting from reproductive interference among sympatric species, was also discussed.

DNA barcoding is currently unreliable for species identification in most crayfishes.

DNA barcoding is commonly used for species identification. Despite this, there has not been a comprehensive assessment of the utility of DNA barcoding in crayfishes (Decapoda: Astacidea). Here we examined the extent to which local barcoding gaps (used for species identification) and global barcoding gaps (used for species discovery) exist among crayfishes, and whether global gaps met a previously suggested 10× threshold (mean interspecific difference being 10× larger than mean intra specific difference). We examined barcoding gaps using publicly available mitochondrial COI sequence data from the National Center for Biotechnology Information's nucleotide database. We created two versions of the COI datasets used for downstream analyses: one focused on the number of unique haplotypes (N H) per species, and another that focused on total number of sequences (N S; i.e., including redundant haplotypes) per species. A total of 81 species were included, with 58 species and five genera from the family Cambaridae and 23 species from three genera from the family Parastacidae. Local barcoding gaps were present in only 30 species (20 Cambaridae and 10 Parastacidae species). We detected global barcoding gaps in only four genera (Cambarus, Cherax, Euastacus, and Tenuibranchiurus), which were all below (4.2× to 5.2×) the previously suggested 10× threshold. We propose that a ~5× threshold would be a more appropriate working hypothesis for species discovery. While the N H and N S datasets yielded largely similar results, there were some discrepant inferences. To understand why some species lacked a local barcoding gap, we performed species delimitation analyses for each genus using the N H dataset. These results suggest that current taxonomy in crayfishes may be inadequate for the majority of examined species, and that even species with local barcoding gaps present may be in need of taxonomic revisions. Currently, the utility of DNA barcoding for species identification and discovery in crayfish is quite limited, and caution should be exercised when mitochondrial-based approaches are used in place of taxonomic expertise. Assessment of the evidence for local and global barcoding gaps is important for understanding the reliability of molecular species identification and discovery, but outcomes are dependent on the current state of taxonomy. As this improves (e.g., via resolving species complexes, possibly elevating some subspecies to the species-level status, and redressing specimen misidentifications in natural history and other collections), so too will the utility of DNA barcoding.

Temporal Changes in Geographic Range and Population Size of the Santa Cruz Black Salamander

Abstract The Santa Cruz black salamander Aneides niger is a priority 3 California species of special concern with a restricted geographic range confined to the Santa Cruz Mountains. Anecdotal observations suggest that the species was relatively abundant in the early 1900s, but it has become more difficult to find in the past few decades. To better understand if the species has undergone population size fluctuations, we analyzed mitochondrial and nuclear sequence data to examine levels of genetic variation and phylogeographic structure, and test for signatures of population size change. We then reconstructed the climatic suitability for the species to 1) determine if past climate fluctuations could have influenced range size and genetic diversity, and 2) estimate the effects of future climate change on geographic range suitability as a proxy for possible future population size change. Genetic analyses detected low levels of genetic variation and a general lack of genetic structuring, suggesting a recent genetic bottleneck. While neutrality tests of individual loci were nonsignificant, skyline plot and isolation-with-migration analyses detected a relatively recent reduction in population size. Interpretation of these genetic results should consider the limited number of localities and individuals sampled for this species. Climatic suitability for Santa Cruz black salamanders was much lower during the last glacial maximum, which could be the cause of the detected historical change in population size. Future projections of climatic suitability under a high-emission scenario suggest a dramatic geographic range restriction to coastal areas. These projections highlight the need for the protection of coastal habitat patches to preserve existing coastal populations, and to maintain connectivity between coastal and inland habitats to allow the westward movement of populations and genes in response to climate change.

A review of Nearctic Lathrobium (Coleoptera, Staphylinidae), with revision and descriptions of new flightless species from the mountains of the southeastern U.S.

Species of the genus Lathrobium Gravenhorst (Coleoptera: Staphylinidae: Paederinae) from North America north of Mexico are reviewed and 41 species are recognized. Morphology and mitochondrial COI sequence data were used to guide species designations in three flightless lineages endemic to the southern Appalachian Mountains, a biologically diverse region known for cryptic diversity. Using a combination of phylogeny, algorithm-based species delimitation analyses, and genitalic morphology, five new cryptic species are described and possible biogeographic scenarios for their speciation hypothesized: L.balsamense Haberski & Caterino, sp. nov., L.camplyacra Haberski & Caterino, sp. nov., L.islae Haberski & Caterino, sp. nov., L.lividum Haberski & Caterino, sp. nov., L.smokiense Haberski & Caterino, sp. nov. Five additional species are described: L.absconditum Haberski & Caterino, sp. nov., L.hardeni Haberski & Caterino, sp. nov., L.lapidum Haberski & Caterino, sp. nov., L.solum Haberski & Caterino, sp. nov., and L.thompsonorum Haberski & Caterino, sp. nov. Two species are transferred from Lathrobium to Pseudolathra Casey: Pseudolathraparcum (LeConte, 1880), comb. nov. and Pseudolathratexana (Casey, 1905), comb. nov. Twenty-six names are reduced to synonymy. Lectotypes are designated for 47 species. Larvae are described where known, and characters of possible diagnostic value are summarized. Species diagnoses, distributions, illustrations of male and female genitalia, and a key to Lathrobium species known from the Nearctic region (including several introduced species) are provided.

Mitogenomic Analysis and Phylogenetic Implications for the Deltocephaline Tribe Chiasmini (Hemiptera: Cicadellidae: Deltocephalinae).

The grassland leafhopper tribe Chiasmini (Cicadellidae: Deltocephalinae) presently comprises 324 described species worldwide, with the highest species diversity occurring in the Nearctic region but a greater diversity of genera occurring in the Old World. In China, this tribe comprises 39 described species in 11 genera, but the fauna remains understudied. The complete mitogenomes of three species of this tribe have been sequenced previously. In order to better understand the phylogenetic position of Chiasmini within the subfamily Deltocephalinae and to investigate relationships among Chiasmini genera and species, we sequenced and analyzed the complete mitogenomes of 13 species belonging to seven genera from China. Comparison of the newly sequenced mitogenomes reveals a closed circular double-stranded structure containing 37 genes with a total length of 14,805 to 16,269 bp and a variable number of non-coding A + T-rich regions. The gene size, gene order, gene arrangement, base composition, codon usage, and secondary structure of tRNAs of the newly sequenced mitogenomes of these 13 species are highly conserved in Chiasmini. The ATN codon is commonly used as the start codon in protein-coding genes (PCGs), except for ND5 in Doratura sp. and ATP6 in Nephotettix nigropictus, which use the rare GTG start codon. Most protein-coding genes have TAA or TAG as the stop codon, but some genes have an incomplete T stop codon. Except for the tRNA for serine (trnS1(AGN)), the secondary structure of the other 21 tRNAs is a typical cloverleaf structure. In addition to the primary type of G-U mismatch, five other types of tRNA mismatches were observed: A-A, A-C, A-G, U-C, and U-U. Chiasmini mitochondrial genomes exhibit gene overlaps with three relatively stable regions: the overlapping sequence between trnW and trnC is AAGTCTTA, the overlapping sequence between ATP8 and ATP6 is generally ATGATTA, and the overlapping sequence between ND4 and ND4L is generally TTATCAT. The largest non-coding region is the control region, which exhibits significant length and compositional variation among species. Some Chiasmini have tandem repeat structures within their control regions. Unlike some other deltocephaline leafhoppers, the sequenced Chiasmini lack mitochondrial gene rearrangements. Phylogenetic analyses of different combinations of protein-coding and ribosomal genes using maximum likelihood and Bayesian methods under different models, using either amino acid or nucleotide sequences, are generally consistent and also agree with results of prior analyses of nuclear and partial mitochondrial gene sequence data, indicating that complete mitochondrial genomes are phylogenetically informative at different levels of divergence within Chiasmini and among leafhoppers in general. Apart from Athysanini and Opsiini, most of the deltocephaline tribes are recovered as monophyletic. The results of ML and BI analyses show that Chiasmini is a monophyletic group with seven monophyletic genera arranged as follows: ((Zahniserius + (Gurawa + (Doratura + Aconurella))) + (Leofa + (Exitianus + Nephotettix))).

First molecular description of Neorhadinorhynchus nudus (Acanthocephala: Cavisomidae) from fish in the pacific coast of Vietnam, with notes on biogeography.

Neorhadinorhynchus nudus (Harada, 1938) Yamaguti, 1939 (Cavisomidae) was morphologically described from the frigate tuna Auxis thazard (Lacépède) (Scombridae) in Nha Trang, Pacific south Vietnam. Females of N. nudus were fully described for the first time in the Pacific. Its original inadequate description as Rhadinorhynchus nudus (Harada, 1938) was corrected in material from Fiji Island, the Red Sea and Pacific Vietnam and errors in the text and line drawings of Harada were repeated in subsequent major publications where it underwent considerable nomenclature changes. New descriptive and biogeographical notes are included. We also provided here the molecular characterization of the nuclear gene (18S) and the mitochondrial cytochromecoxidase subunit 1 (cox1) sequence data ofN. nudus. Furthermore, to elucidate the phylogenetic relationship of N. nudus within the family Cavisomidae and with other isolates were performed incorporating nuclear (18S) and mitochondrial (cox1) sequence data using maximum likelihood (ML) and Bayesian inference (BI). The phylogenetic results showed thatN. nudus has a relationship withother isolates of the same species and the median-joining network showed the pattern of haplotypes that reflected the structure of the populations.

Morphological Ontogeny, Ecology, and Biogeography of Fuscozetes fuscipes (Acari, Oribatida, Ceratozetidae).

The systematic status of Fuscozetes Sellnick, 1928, is not clear in the literature. Therefore, the morphological ontogeny of F. fuscipes (C.L. Koch, 1844), the type species of this genus, was investigated and compared with its congeners in this study, and a new diagnosis of Fuscozetes is given. The juveniles of F. fuscipes are light brown, with a brown prodorsum, sclerites, epimeres, and legs. In all juveniles, a humeral organ and a humeral macrosclerite are present. The gastronotum of the larva has 12 pairs of setae (h3 is present), whereas the nymphs have 15 pairs. In the larva, the gastronotal shield is weakly developed, and most gastronotal setae are short except for a slightly longer h2. Most of the gastronotal setae are inserted on the microsclerites except for h3, and several other macrosclerites and many microsclerites are present on the hysterosoma. In the nymphs, the gastronotal shield is well developed, with 10 pairs of setae (d-, l-, and h-series, and p1), and setae p2 and p3 are located on a large posteroventral macrosclerite. In all the instars, femora I and II are oval in cross-section, without a large ventral carina. Mitochondrial COI sequence data revealed a deep split between the Nearctic and Palearctic populations of F. fuscipes, and a less, but significant, divergence within each continent. These strong geographical barriers were contrasted with multiple cases of shared haplotypes over long distances in the Palearctic, indicating high migration rates in modern times.

MSeqDR Quick-Mitome (QM): Combining Phenotype-Guided Variant Interpretation and Machine Learning Classifiers to Aid Primary Mitochondrial Disease Genetic Diagnosis.

The international Mitochondrial Disease Sequence Data Resource Consortium (MSeqDR) Quick-Mitome (QM) is a web-based platform enabling automated variant interpretation of whole-exome sequencing (WES) datasets for the genetic diagnosis of primary mitochondrial diseases (PMD). Designed specifically to address the unique dual genome nature of PMD etiologies, QM includes features for both nuclear and mitochondrial DNA (mtDNA) genome analysis. QM requires VCF variant lists, HPO ID clinical phenotypes, and pedigree files for multiple-sample VCF inputs. QM maps phenotypes to HPO terms before analysis. QM analysis requires 2 to 20 min for 100,000 variants on an 8-vCPU AWS server using Exomiser's "PASS_ONLY" mode for nuclear variants. QM ranks variants based on allele frequency, phenotype-gene association, functional impact, and inheritance mode. Variants are further annotated with multiple data sources such as OMIM, ClinVar, dbNSFP, gnoMAD, MITOMAP, and MSeqDR. In addition to standard Exomiser results, QM generates an Analysis Report and QM Integrated Report with add-on mtDNA-specific analyses, including haplogroup prediction with Phy-Mer, heteroplasmy calculation, and mvTool annotations. We developed the Mitochondrial Disease Variant (MDV) classifier using XGBoost to predict variant pathogenicity for PMD. The MDV classifier was trained on >120 features and performance benchmarking showed that it correctly classified >98% of nuclear gene variants as being pathogenic or benign, and predicted PMD-causing variants with 94% precision. The MSeqDR QM server is an open-access resource for phenotype-driven dual-genome analyses for PMD diagnosis by the global mitochondrial disease community. It is publicly available for non-commercial, non-clinical research use at https://mseqdr.org/quickmitome.php. © 2024 Wiley Periodicals LLC. Basic Protocol 1: Standardizing clinical phenotypes into human phenotype ontology (HPO) terms as the phenotype input for Quick-Mitome (QM) Basic Protocol 2: Prepare the pedigree input for multiple-sample VCF Basic Protocol 3: Quick-Mitome (QM) analysis Basic Protocol 4: Reviewing and understanding the QM Integrated Report and Analysis Report.

Genetic diversity, variation, and structure of two populations of bigfin reef squid (Sepioteuthis lessoniana d’Orbigny) in Con Dao and Phu Quoc islands, Vietnam

BackgroundBigfin squid is one of the economically important seafood resources in Vietnam’s fisheries and the waters around Con Dao and Phu Quoc islands are two major fishing grounds where this species has been actively exploited. The start codon targeted polymorphism (SCoT) and CAAT box–derived polymorphism (CBDP) techniques were used to generate DNA fingerprinting data to analyze the genetic diversity, variation, and structure of the two populations in the waters surrounding Phu Quoc and Con Dao islands together with mitochondrial cytochrome C oxidase subunit I (COI) gene sequence data. ResultsCon Dao population possessed a higher diversity [expected heterozygosity (He) = 0.2254, Shannon index (I) = 0.3459, percentage of polymorphic bands (PPB) = 80.14%, nucleotide diversity (π) = 0.0336, haplotype diversity (h) = 0.910 with 16 haplotypes] than Phu Quoc population (He = 0.1854, I = 0.2873, PPB = 70.38%, π = 0.0246, h = 0.838 with 14 haplotypes). The genetic diversity at species level in the investigated region was at level of He = 0.2169, I = 0.3399, PPB = 86.41, π = 0.0289, and h = 0.892 with 24 haplotypes. Based on DNA fingerprinting data, the pairwise genetic similarity coefficients among individuals of the Con Dao population were lower (average of 0.7977) than the Phu Quoc population (average of 0.8316). Based on mitochondrial COI data, the pairwise genetic distances among individuals of the Con Dao population were higher (average of 0.0361) than the Phu Quoc population (average of 0.0263). Gene differentiation (GST) between two investigated populations was 0.0316 and 0.0310 leading to the genetic distance was 0.0573 and 0.0213 and the gene flow between them was Nm = 8.2209 and 11.4700 migrants per generation among populations based on DNA fingerprinting and based on COI gene sequence data, respectively. Genetic variation within individuals of both populations (WP) played the key role in total genetic variation at species level in surveyed region. ConclusionsFor the bigfin reef squid species in the surveyed region, the Con Dao population had the higher genetic diversity than the Phu Quoc population, between them existed a low to moderate genetic differentiation and a genetic exchange via gene flow. The DNA fingerprinting data better revealed the genetic differentiation between the two surveyed populations while the mitochondrial COI gene sequence data could show the phylogenetic relationship among the surveyed individuals and the other from the sea regions in Southeast Asia. Based on the results obtained, fisheries management strategies are suggested toward the conservation and sustainable exploitation of this species.

Adulis and the transshipment of baboons during classical antiquity.

Adulis, located on the Red Sea coast in present-day Eritrea, was a bustling trading centre between the first and seventh centuries CE. Several classical geographers--Agatharchides of Cnidus, Pliny the Elder, Strabo-noted the value of Adulis to Greco--Roman Egypt, particularly as an emporium for living animals, including baboons (Papio spp.). Though fragmentary, these accounts predict the Adulite origins of mummified baboons in Ptolemaic catacombs, while inviting questions on the geoprovenance of older (Late Period) baboons recovered from Gabbanat el-Qurud ('Valley of the Monkeys'), Egypt. Dated to ca. 800-540 BCE, these animals could extend the antiquity of Egyptian-Adulite trade by as much as five centuries. Previously, Dominy et al. (2020) used stable istope analysis to show that two New Kingdom specimens of P. hamadryas originate from the Horn of Africa. Here, we report the complete mitochondrial genomes from a mummified baboon from Gabbanat el-Qurud and 14 museum specimens with known provenance together with published georeferenced mitochondrial sequence data. Phylogenetic assignment connects the mummified baboon to modern populations of Papio hamadryas in Eritrea, Ethiopia, and eastern Sudan. This result, assuming geographical stability of phylogenetic clades, corroborates Greco-Roman historiographies by pointing toward present-day Eritrea, and by extension Adulis, as a source of baboons for Late Period Egyptians. It also establishes geographic continuity with baboons from the fabled Land of Punt (Dominy et al., 2020), giving weight to speculation that Punt and Adulis were essentially the same trading centres separated by a thousand years of history.

Phylogenetic review of the comb-tooth blenny genus Hypleurochilus in the northwest Atlantic and Gulf of Mexico

As some of the smallest vertebrates, yet largest producers of consumed reef biomass, cryptobenthic reef fishes serve a disproportionate role in reef ecosystems and are one of the most poorly understood groups of fish. The blenny genera Hypleurochilus and Parablennius are currently considered paraphyletic and the interrelationships of Parablennius have been the focus of recent phylogenetic studies. However, the interrelationships of Hypleurochilus remain understudied. This genus is transatlantically distributed and comprises 11 species with a convoluted taxonomic history. In this study, relationships for ten Hypleurochilus species are resolved using multi-locus nuclear and mtDNA sequence data, morphological data, and mined COI barcode data. Mitochondrial and nuclear sequence data from 61 individuals collected from the western Atlantic and northern Gulf of Mexico (N. GoM) delimit seven species into a temperate clade, a tropical clade, and a third distinct lineage. This lineage, herein referred to as H. cf. aequipinnis, may represent a species of Hypleurochilus whose range has expanded into the N. GoM. Inclusion of publicly available COI sequence for an additional three species provides further phylogenetic resolution. H. bananensis forms a new eastern Atlantic clade with H. cf. aequipinnis, providing further evidence for a western Atlantic range expansion. Single marker COI delimitation was unable to elucidate the relationships between H. springeri/H. pseudoaequipinnis and between H. multifilis/H. caudovittatus due to incomplete lineage sorting. Mitochondrial data are also unable to accurately resolve the placement of H. bermudensis. However, a comprehensive approach using multi-locus phylogenetic and species delimitation methods was able to resolve these relationships. While mining publicly available sequence data allowed for the inclusion of an increased number of species in the analysis and a more comprehensive phylogeny, it was not without drawbacks, as a handful of sequences are potentially mis-identified. Overall, we find that the recent divergence of some species within this genus and potential introgression events confound the results of single locus delimitation methods, yet a combination of single and multi-locus analyses has allowed for insights into the biogeography of this genus and uncovered a potential transatlantic range expansion.

Fully automated annotation of mitochondrial genomes using a cluster-based approach with de Bruijn graphs.

A wide range of scientific fields, such as forensics, anthropology, medicine, and molecular evolution, benefits from the analysis of mitogenomic data. With the development of new sequencing technologies, the amount of mitochondrial sequence data to be analyzed has increased exponentially over the last few years. The accurate annotation of mitochondrial DNA is a prerequisite for any mitogenomic comparative analysis. To sustain with the growth of the available mitochondrial sequence data, highly efficient automatic computational methods are, hence, needed. Automatic annotation methods are typically based on databases that contain information about already annotated (and often pre-curated) mitogenomes of different species. However, the existing approaches have several shortcomings: 1) they do not scale well with the size of the database; 2) they do not allow for a fast (and easy) update of the database; and 3) they can only be applied to a relatively small taxonomic subset of all species. Here, we present a novel approach that does not have any of these aforementioned shortcomings, (1), (2), and (3). The reference database of mitogenomes is represented as a richly annotated de Bruijn graph. To generate gene predictions for a new user-supplied mitogenome, the method utilizes a clustering routine that uses the mapping information of the provided sequence to this graph. The method is implemented in a software package called DeGeCI (De Bruijn graph Gene Cluster Identification). For a large set of mitogenomes, for which expert-curated annotations are available, DeGeCI generates gene predictions of high conformity. In a comparative evaluation with MITOS2, a state-of-the-art annotation tool for mitochondrial genomes, DeGeCI shows better database scalability while still matching MITOS2 in terms of result quality and providing a fully automated means to update the underlying database. Moreover, unlike MITOS2, DeGeCI can be run in parallel on several processors to make use of modern multi-processor systems.

Till evolution do us part: The diversity of symbiotic associations across populations of Philaenus spittlebugs.

Symbiotic bacteria have played crucial roles in the evolution of sap-feeding insects and can strongly affect host function. However, their diversity and distribution within species are not well understood; we do not know to what extent environmental factors or associations with other species may affect microbial community profiles. We addressed this question in Philaenus spittlebugs by surveying both insect and bacterial marker gene amplicons across multiple host populations. Host mitochondrial sequence data confirmed morphology-based identification of six species and revealed two divergent clades of Philaenus spumarius. All of them hosted the primary symbiont Sulcia that was almost always accompanied by Sodalis. Interestingly, populations and individuals often differed in the presence of Sodalis sequence variants, suggestive of intra-genome 16S rRNA variant polymorphism combined with rapid genome evolution and/or recent additional infections or replacements of the co-primary symbiont. The prevalence of facultative endosymbionts, including Wolbachia, Rickettsia, and Spiroplasma, varied among populations. Notably, cytochrome I oxidase (COI) amplicon data also showed that nearly a quarter of P. spumarius were infected by parasitoid flies (Verralia aucta). One of the Wolbachia operational taxonomic units (OTUs) was exclusively present in Verralia-parasitized specimens, suggestive of parasitoids as their source and highlighting the utility of host gene amplicon sequencing in microbiome studies.

Morphological constraint obscures richness: a mitochondrial exploration of cryptic richness in Transversotrema (Trematoda: Transversotrematidae)

Species of Transversotrema Witenberg, 1944 (Transversotrematidae) occupy a unique ecological niche for the Trematoda, living externally under the scales of their teleost hosts. Previous studies of the genus have been impeded partly by limited variation in ribosomal DNA sequence data between closely related species and partly by a lack of morphometrically informative characters. Here, we assess richness of the tropical Indo-west Pacific species through parallel phylogenetic and morphometric analyses, generating cytochrome c oxidase subunit 1 mitochondrial sequence data and morphometric data for hologenophore specimens from Australia, French Polynesia, Japan and Palau. These analyses demonstrate that molecular data provide the only reliable basis for species identification; host distribution, and to a lesser extent morphology, are useful for identifying just a few species of Transversotrema. We infer that a combination of morphological simplicity and infection site constraint has led to the group displaying exceptionally low morphological diversification. Phylogenetic analyses of the mitochondrial data broadly support previous systematic interpretations based on ribosomal data, but also demonstrate the presence of several morphologically and ecologically cryptic species. Ten new species are described, eight from the Great Barrier Reef, Australia (Transversotrema chrysallis n. sp., Transversotrema daphnidis n. sp., Transversotrema enceladi n. sp., Transversotrema hyperionis n. sp., Transversotrema iapeti n. sp., Transversotrema rheae n. sp., Transversotrema tethyos n. sp., and Transversotrema titanis n. sp.) and two from off Japan (Transversotrema methones n. sp. and Transversotrema panos n. sp.). There are now 26 Transversotrema species known from Australian marine fishes, making it the richest trematode genus for the fauna.

Strong Population Genetic Structure for the Endangered Micro-Trapdoor Spider Moggridgea rainbowi (Mygalomorphae, Migidae) in Unburnt Habitat after Catastrophic Bushfires

Catastrophic wildfires impacted large areas of western Kangaroo Island (KI), South Australia in 2019–2020, burning habitat for many species, including large proportions of the distributional range of the KI micro-trapdoor spider Moggridgea rainbowi, which led to it being listed as Endangered under Australia’s Environment Protection and Biodiversity Conservation Act (EPBC Act). In order to assess population genetic structure in M. rainbowi and detect diagnosable lineages and their distributional patterns across KI, 28 individuals were genotyped for 2495 loci, sampling from all known populations of the species. Population genetic and phylogenetic analyses of nuclear and mitochondrial sequence data provided strong support for three Evolutionarily Significant Units (ESUs) within M. rainbowi; two populations on eastern KI and a heavily fire-impacted western population. High levels of divergence and fixed allelic differences between 5 and 16% indicate a lack of gene flow between ESUs and long periods of isolation. Distributional patterns of these lineages match likely locations of isolation events caused by successive changes to sea level during the Quaternary (2.58 million years ago to present), which led to KI being intermittently connected to the mainland or separated into one or more islands. Our findings have strong conservation implications for M. rainbowi and highlight the importance of inclusion of population genetic structure to inform conservation strategies and to conserve lineage biodiversity at the species level and below.

Environmental and Socio-Cultural Factors Impacting the Unique Gene Pool Pattern of Mae Hong-Son Chicken.

Understanding the genetic diversity of domestic chicken breeds under the impact of socio-cultural and ecological dynamics is vital for the conservation of natural resources. Mae Hong Son chicken is a local breed of North Thai domestic chicken widely distributed in Mae Hong Son Province, Thailand; however, its genetic characterization, origin, and diversity remain poorly understood. Here, we studied the socio-cultural, environmental, and genetic aspects of the Mae Hong Son chicken breed and investigated its diversity and allelic gene pool. We genotyped 28 microsatellite markers and analyzed mitochondrial D-loop sequencing data to evaluate genetic diversity and assessed spatial habitat suitability using maximum entropy modeling. Sequence diversity analysis revealed a total of 188 genotyped alleles, with overall nucleotide diversity of 0.014 ± 0.007, indicating that the Mae Hong Son chicken population is genetically highly diverse, with 35 (M1-M35) haplotypes clustered into haplogroups A, B, E, and F, mostly in the North ecotype. Allelic gene pool patterns showed a unique DNA fingerprint of the Mae Hong Son chicken, as compared to other breeds and red junglefowl. A genetic introgression of some parts of the gene pool of red junglefowl and other indigenous breeds was identified in the Mae Hong Son chicken, supporting the hypothesis of the origin of the Mae Hong Son chicken. During domestication in the past 200-300 years after the crossing of indigenous chickens and red junglefowl, the Mae Hong Son chicken has adapted to the highland environment and played a significant socio-cultural role in the Northern Thai community. The unique genetic fingerprint of the Mae Hong Son chicken, retaining a high level of genetic variability that includes a dynamic demographic and domestication history, as well as a range of ecological factors, might reshape the adaptation of this breed under selective pressure.