Mind the gap! The mitochondrial control region and its power as a phylogenetic marker in echinoids

Tracking Rare Single Donor and Recipient Immune and Leukemia Cells after Allogeneic Hematopoietic Cell Transplantation Using Mitochondrial DNA Mutations

Combined tracking of clonal evolution and chimeric cell phenotypes could enable detection of the key cellular populations associated with response following therapy, including after allogeneic hematopoietic stem cell transplantation (HSCT). We demonstrate that mitochondrial DNA (mtDNA) mutations coevolve with somatic nuclear DNA mutations at relapse post-HSCT and provide a sensitive means to monitor these cellular populations. Furthermore, detection of mtDNA mutations via single-cell assay for transposase-accessible chromatin with select antigen profiling by sequencing (ASAP-seq) simultaneously determines not only donor and recipient cells but also their phenotype at frequencies of 0.1% to 1%. Finally, integration of mtDNA mutations, surface markers, and chromatin accessibility profiles enables the phenotypic resolution of leukemic populations from normal immune cells, thereby providing fresh insights into residual donor-derived engraftment and short-term clonal evolution following therapy for post-transplant leukemia relapse. As throughput evolves, we envision future development of single-cell sequencing-based post-transplant monitoring as a powerful approach for guiding clinical decision-making. Significance: mtDNA mutations enable single-cell tracking of leukemic clonal evolution and donor-recipient origin following allogeneic HSCT. This provides unprecedented insight into chimeric cellular phenotypes of early immune reconstitution, incipient relapse, and quality of donor engraftment with immediate translational potential for future clinical post-transplant monitoring and decision-making.

Structures and genetic information of control region in mitogenomes of Odonata

Mitogenome data of Odonata is accumulating and widely used in phylogenetic analysis. However, noncoding regions, especially control region, were usually omitted from the phylogenetic reconstruction. In an effort to uncover the phylogenetic insights offered by the control region, we have amassed 65 Odonata mitogenomes and conducted an examination of their control regions. Our analysis discovered that species belonging to Anisoptera and Anisozygoptera exhibited a stem-loop structure, which was formed by a conserved polyC-polyG stretch located near the rrns gene (encoding 12S rRNA). Conversely, the polyC-polyG region was not a conserved fragment in Zygoptera. The length and number of repetitions within the control region were identified as the primary determinants of its overall length. Further, sibling species within Odonata, particularly those in the genus Euphaea, displayed similar patterns of repetition in their control region. Collectively, our research delineates the structural variations within the control region of Odonata and suggests the potential utility of this region in elucidating phylogenetic relationships among closely related species.

Mitogenomic Insights into the Evolution, Divergence Time, and Ancestral Ranges of Coturnix Quails

The Old-World quails, Coturnix coturnix (common quail) and Coturnix japonica (Japanese quail), are morphologically similar yet occupy distinct geographic ranges. This study aimed to elucidate their evolutionary trajectory and ancestral distribution patterns through a thorough analysis of their mitochondrial genomes. Mitogenomic analysis revealed high structural conservation, identical translational mechanisms, and similar evolutionary pressures in both species. Selection analysis revealed significant evidence of positive selection across the Coturnix lineage for the nad4 gene tree owing to environmental changes and acclimatization requirements during its evolutionary history. Divergence time estimations imply that diversification among Coturnix species occurred in the mid-Miocene (13.89 Ma), and their current distributions were primarily shaped by dispersal rather than global vicariance events. Phylogenetic analysis indicates a close relationship between C. coturnix and C. japonica, with divergence estimated at 2.25 Ma during the Pleistocene epoch. Ancestral range reconstructions indicate that the ancestors of the Coturnix clade were distributed over the Oriental region. C. coturnix subsequently dispersed to Eurasia and Africa, and C. japonica to eastern Asia. We hypothesize that the current geographic distributions of C. coturnix and C. japonica result from their unique dispersal strategies, developed to evade interspecific territoriality and influenced by the Tibetan Plateau's geographic constraints. This study advances our understanding of the biogeographic and evolutionary processes leading to the diversification of C. coturnix and C. japonica, laying important groundwork for further research on this genus.

Genetic Diversity and Population Structure of Brown Croaker (Miichthys miiuy) in Korea and China Inferred from mtDNA Control Region

Brown croaker (Miichthys miiuy), a species of fish with significant commercial value, is found in the coastal seas of Korea, China, and Japan. The genetic diversity and population structure of a representative sample of brown croaker specimens were assessed based on the control region of their mitochondrial DNA (mtDNA). Samples from a total of 115 individuals were collected from three separate locations, one in China (Lianyungang) and two in Korea (Mokpo and Gyeongnyeolbiyeoldo Island). Analysis of the 436-base-pair mtDNA control region revealed that the haplotype diversity ranged from 0.973 ± 0.025 to 0.988 ± 0.008, while the nucleotide diversity ranged from 0.012 ± 0.006 to 0.017 ± 0.009. The level of genetic diversity, star-shaped haplotype network, significant Fu's Fs test, and analysis of the mismatch distribution all suggested that this species has experienced population expansion. Fixation index analysis indicated that the population collected at the site in China differed significantly from the two populations obtained in Korea. The findings of this study extend the general understanding of the population structure of M. miiuy and can be used to develop strategies for effective resource management.

Ancient DNA re-opens the question of the phylogenetic position of the Sardinian pika Prolagus sardus (Wagner, 1829), an extinct lagomorph

Palaeogenomics is contributing to refine our understanding of many major evolutionary events at an unprecedented resolution, with relevant impacts in several fields, including phylogenetics of extinct species. Few extant and extinct animal species from Mediterranean regions have been characterised at the DNA level thus far. The Sardinian pika, Prolagus sardus (Wagner, 1829), was an iconic lagomorph species that populated Sardinia and Corsica and became extinct during the Holocene. There is a certain scientific debate on the phylogenetic assignment of the extinct genus Prolagus to the family Ochotonidae (one of the only two extant families of the order Lagomorpha) or to a separated family Prolagidae, or to the subfamily Prolaginae within the family Ochotonidae. In this study, we successfully reconstructed a portion of the mitogenome of a Sardinian pika dated to the Neolithic period and recovered from the Cabaddaris cave, an archaeological site in Sardinia. Our calibrated phylogeny may support the hypothesis that the genus Prolagus is an independent sister group to the family Ochotonidae that diverged from the Ochotona genus lineage about 30 million years ago. These results may contribute to refine the phylogenetic interpretation of the morphological peculiarities of the Prolagus genus already described by palaeontological studies.

First Report of Rickettsia conorii in Hyalomma kumari Ticks

As a vector of wide range of pathogenic agents, ticks pose health threats to wild and domestic animals, and humans. Information is unavailable about the prevalence and spatial survey of Hyalomma kumari ticks and associated Rickettsia spp. in Pakistan. Concerning this knowledge gap, the present study aimed to molecularly detect Rickettsia species associated with H. kumari infesting small ruminants in Khyber Pakhtunkhwa (KP), Pakistan. A total of 409 H. kumari ticks were collected from 163/295 infested hosts with an infestation rate of 55.25%. A total of 204 females, 158 males, and 47 nymphs were collected. Goats were heavily infested by 224 ticks having an infestation rate of 58.33% (98/168), whereas sheep were infested by 185 ticks having a lesser infestation rate of 51.18% (65/127). Genomic DNA extracted from ticks was used for the amplification of tick (cox I, 16S rRNA, ITS-2) species and Rickettsia (gltA, ompA, and ompB) partial genes. Eighty-three ticks were subjected to PCR, and 8/83 (9.6%) were found positive for rickettsial agents. The cox I and 16S rRNA sequences of H. kumari showed 98.90-99.74% identity with H. kumari sequences reported from Pakistan, and phylogenetically clustered to the corresponding species reported from Pakistan and India. The obtained rickettsial gltA, ompA, and ompB sequences showed 100% identity with Rickettsia sp. of the Rickettsia conorii reported from Pakistan. In the phylogenetic trees, rickettsial sequences clustered with uncharacterized Rickettsia sp. from Pakistan and R. conorii from Israel, Russia, South Africa, and India. The present molecular based detection of H. kumari-associated R. conorii will facilitate effective surveillance in the region.

Comparative phylogeography reveals widespread cryptic diversity driven by ecology in Panamanian birds

Widespread species often harbor unrecognized genetic diversity, and investigating the factors associated with such cryptic variation can help us better understand the forces driving diversification. Here, we identify potential cryptic species based on a comprehensive dataset of COI mitochondrial DNA barcodes from 2,333 individual Panamanian birds across 429 species, representing 391 (59%) of the 659 resident landbird species of the country, as well as opportunistically sampled waterbirds. We complement this dataset with additional publicly available mitochondrial loci, such as ND2 and cytochrome b, obtained from whole mitochondrial genomes from 20 taxa. Using barcode identification numbers (BINs), we find putative cryptic species in 19% of landbird species, highlighting hidden diversity in the relatively well-described avifauna of Panama. Whereas some of these mitochondrial divergence events corresponded with recognized geographic features that likely isolated populations, such as the Cordillera Central highlands, the majority (74%) of lowland splits were between eastern and western populations. The timing of these splits are not temporally coincident across taxa, suggesting that historical events, such as the formation of the Isthmus of Panama and Pleistocene climatic cycles, were not the primary drivers of cryptic diversification. Rather, we observed that forest species, understory species, insectivores, and strongly territorial species-all traits associated with lower dispersal ability-were all more likely to have multiple BINs in Panama, suggesting strong ecological associations with cryptic divergence. Additionally, hand-wing index, a proxy for dispersal capability, was significantly lower in species with multiple BINs, indicating that dispersal ability plays an important role in generating diversity in Neotropical birds. Together, these results underscore the need for evolutionary studies of tropical bird communities to consider ecological factors along with geographic explanations, and that even in areas with well-known avifauna, avian diversity may be substantially underestimated.

LAY SUMMARY

- What factors are common among bird species with cryptic diversity in Panama? What role do geography, ecology, phylogeographic history, and other factors play in generating bird diversity?- 19% of widely-sampled bird species form two or more distinct DNA barcode clades, suggesting widespread unrecognized diversity.- Traits associated with reduced dispersal ability, such as use of forest understory, high territoriality, low hand-wing index, and insectivory, were more common in taxa with cryptic diversity. Filogeografía comparada revela amplia diversidad críptica causada por la ecología en las aves de Panamá.

RESUMEN

Especies extendidas frecuentemente tiene diversidad genética no reconocida, y investigando los factores asociados con esta variación críptica puede ayudarnos a entender las fuerzas que impulsan la diversificación. Aquí, identificamos especies crípticas potenciales basadas en un conjunto de datos de códigos de barras de ADN mitocondrial de 2,333 individuos de aves de Panama en 429 especies, representando 391 (59%) de las 659 especies de aves terrestres residentes del país, además de algunas aves acuáticas muestreada de manera oportunista. Adicionalmente, complementamos estos datos con secuencias mitocondriales disponibles públicamente de otros loci, tal como ND2 o citocroma b, obtenidos de los genomas mitocondriales completos de 20 taxones. Utilizando los números de identificación de código de barras (en ingles: BINs), un sistema taxonómico numérico que proporcina una estimación imparcial de la diversidad potencial a nivel de especie, encontramos especies crípticas putativas en 19% de las especies de aves terrestres, lo que destaca la diversidad oculta en la avifauna bien descrita de Panamá. Aunque algunos de estos eventos de divergencia conciden con características geográficas que probablemente aislaron las poblaciones, la mayoría (74%) de la divergencia en las tierras bajas se encuentra entre las poblaciones orientales y occidentales. El tiempo de esta divergencia no coincidió entre los taxones, sugiriendo que eventos históricos tales como la formación del Istmo de Panamá y los ciclos climáticos del pleistoceno, no fueron los principales impulsores de la especiación. En cambio, observamos asociaciones fuertes entre las características ecológicas y la divergencia mitocondriale: las especies del bosque, sotobosque, con una dieta insectívora, y con territorialidad fuerte mostraton múltiple BINs probables. Adicionalmente, el índice mano-ala, que está asociado a la capacidad de dispersión, fue significativamente menor en las especies con BINs multiples, sugiriendo que la capacidad de dispersión tiene un rol importamente en la generación de la diversidad de las aves neotropicales. Estos resultos demonstran la necesidad de que estudios evolutivos de las comunidades de aves tropicales consideren los factores ecológicos en conjunto con las explicaciones geográficos. Palabras clave: biodiversidad tropical, biogeografía, códigos de barras, dispersión, especies crípticas.

Comparative phylogeography reveals widespread cryptic diversity driven by ecology in Panamanian birds

Widespread species often harbor unrecognized genetic diversity, and investigating the factors associated with such cryptic variation can help us better understand the forces driving diversification. Here, we identify potential cryptic species based on a comprehensive dataset of COI mitochondrial DNA barcodes from 2,333 individual Panamanian birds across 429 species, representing 391 (59%) of the 659 resident landbird species of the country, as well as opportunistically sampled waterbirds. We complement this dataset with additional publicly available mitochondrial loci, such as ND2 and cytochrome b, obtained from whole mitochondrial genomes from 20 taxa. Using barcode identification numbers (BINs), we find putative cryptic species in 19% of landbird species, highlighting hidden diversity in the relatively well-described avifauna of Panama. Whereas some of these mitochondrial divergence events corresponded with recognized geographic features that likely isolated populations, such as the Cordillera Central highlands, the majority (74%) of lowland splits were between eastern and western populations. The timing of these splits are not temporally coincident across taxa, suggesting that historical events, such as the formation of the Isthmus of Panama and Pleistocene climatic cycles, were not the primary drivers of cryptic diversification. Rather, we observed that forest species, understory species, insectivores, and strongly territorial species-all traits associated with lower dispersal ability-were all more likely to have multiple BINs in Panama, suggesting strong ecological associations with cryptic divergence. Additionally, hand-wing index, a proxy for dispersal capability, was significantly lower in species with multiple BINs, indicating that dispersal ability plays an important role in generating diversity in Neotropical birds. Together, these results underscore the need for evolutionary studies of tropical bird communities to consider ecological factors along with geographic explanations, and that even in areas with well-known avifauna, avian diversity may be substantially underestimated.

Phylogenetic Diversity of the Red Swamp Crayfish Procambarus clarkii and Its Dispersal Pattern in Northern and Central Italy

The red swamp crayfish Procambarus clarkii is one of the most threatening freshwater species in the world. The aim of this study is to provide a better understanding of the phylogeography and the invasion routes of P. clarkii populations in the Italian Peninsula through the analysis of mitochondrial phylogeny. Mitochondrial control region and cytochrome c oxidase subunit I (COI) sequences of 153 samples collected from six Italian basins were analyzed and compared to worldwide data. Except for the lakes Bolsena and Posta Fibreno, a high genetic variability was found in the other basins. The mitochondrial DNA pattern of P. clarkii from the lakes Candia and Massaciuccoli confirmed the hypothesis of double introduction events. Another entry point could be represented by Lake Trasimeno, which shows haplotypes originating from Louisiana and not shared with other Italian basins. Moreover, unique lineages were also found in the Stella River, thus enhancing the hypothesis that multiple introductions of P. clarkii occurred in northern and Central Italy and strengthening the idea that knowledge about the dispersion routes of this alien species can be useful to predict its invasiveness and elaborate control strategies to preserve biodiversity.

Low genetic diversity in a widespread whistling alien: A comparison of Eleutherodactylus johnstonei Barbour, 1914 (Eleutherodactylidae) and congeners in native and introduced ranges

There is no clear empirical evidence to support the general assumption that genetic diversity favours successful invasions. Many invading species disperse and establish successfully despite low genetic diversity, a phenomenon known as the genetic paradox of biological invasion. Model systems that allow comparison of genetic patterns between exotic and native source populations are still scarce. This is particularly true for amphibians. Here we compare genetic patterns of the widely introduced Johnstone’s Whistling Frog, Eleutherodactylus johnstonei, with its successful alien congener E. antillensis and the single island endemic E. portoricensis. Genetic diversity and population differentiation in native and introduced populations of the three taxa were inferred from mitochondrial D-loop sequences (235 bp). Our results reveal that exotic populations of the two alien taxa, E. johnstonei and E. antillensis, are not only genetically impoverished due to founder effects, but that, moreover, their native range source-populations exhibit low genetic diversity and inter-population differentiation in the first place. Populations of the endemic E. portoricensis, on the other hand, are genetically more diverse and show marked inter-population differentiation. These observed genetic patterns are consistent with geological processes and invasion histories. We argue that the establishment success of the alien taxa in our model system is better explained by ecological factors and anthropogenic drivers than by genetic diversity. As these factors provide more parsimonious explanations, they should be given priority in management decisions. However, molecular studies with higher resolution are needed to fully test possible genetic and epigenetic components that could promote the invasion process.

In silico identification of multiple conserved motifs within the control region of Culicidae mitogenomes

Mosquitoes are important vectors for human and animal diseases. Genetic markers, like the mitochondrial COI gene, can facilitate the taxonomic classification of disease vectors, vector-borne disease surveillance, and prevention. Within the control region (CR) of the mitochondrial genome, there exists a highly variable and poorly studied non-coding AT-rich area that contains the origin of replication. Although the CR hypervariable region has been used for species differentiation of some animals, few studies have investigated the mosquito CR. In this study, we analyze the mosquito mitogenome CR sequences from 125 species and 17 genera. We discovered four conserved motifs located 80 to 230 bp upstream of the 12S rRNA gene. Two of these motifs were found within all 392 Anopheles (An.) CR sequences while the other two motifs were identified in all 37 Culex (Cx.) CR sequences. However, only 3 of the 304 non-Culicidae Dipteran mitogenome CR sequences contained these motifs. Interestingly, the short motif found in all 37 Culex sequences had poly-A and poly-T stretch of similar length that is predicted to form a stable hairpin. We show that supervised learning using the frequency chaos game representation of the CR can be used to differentiate mosquito genera from their dipteran relatives.

The Quality of Sequence Data Affects Biodiversity and Conservation Perspectives in the Neotropical Damselfly Megaloprepus caerulatus

Ideally, the footprint of the evolutionary history of a species is drawn from integrative studies including quantitative and qualitative taxonomy, biogeography, ecology, and molecular genetics. In today’s research, species delimitations and identification of conservation units is often accompanied by a set of—at minimum—two sequence markers appropriate for the systematic level under investigation. Two such studies re-evaluated the species status in the world’s largest Odonata, the Neotropical damselfly Megaloprepus caerulatus. The species status of the genus Megaloprepus has long been debated. Despite applying a highly similar set of sequence markers, the two studies reached different conclusions concerning species status and population genetic relationships. In this study, we took the unique opportunity to compare the two datasets and analyzed the reasons for those incongruences. The two DNA sequence markers used (16S rDNA and CO1) were re-aligned using a strict conservative approach and the analyses used in both studies were repeated. Going step by step back to the first line of data handling, we show that a high number of unresolved characters in the sequence alignments as well as internal gaps are responsible for the different outcomes in terms of species delimitations and population genetic relationships. Overall, this study shows that high quality raw sequence data are an indispensable requirement, not only in odonate research.

Strategy of micro-environmental adaptation to cold seep among different brittle stars’ colonization

Diffusing fluid from methane seepage in cold seep field creates zones with physicochemical gradients and divergent ecosystems like the mussel beds and clam beds. Three species of brittle stars (Ophiuroidea) were discovered in the Haima cold seep fields, of which Ophiophthalmus serratus and Histampica haimaensis were found on top of or within mussel beds and clam beds, whereas Amphiura sp. was only collected from muds in the clam bed assemblage. Here, we evaluated the genetic signatures of micro-environmental adaptation of brittle stars to cold seep through the comparison of mitogenomes. This study provided two complete mitogenome sequences of O. serratus and Amphiura sp. and compared with those of H. haimaensis and other non-seep species. We found that the split events of the seep and non-seep species were as ancient as the Cretaceous period (∼148–98 Mya). O. serratus and H. haimaensis display rapid residue mutation and mitogenome rearrangements compared to their shallow or deep-sea relatives, in contrast, Amphiura sp. only show medium, regardless of nucleotide mutation rate or mitogenome rearrangement, which may correlate with their adaptation to one or two micro-ecosystems. Furthermore, we identified 10 positively selected residues in ND4 in the Amphiura sp. lineage, suggesting important roles of the dehydrogenase complex in Amphiura sp. adaptive to the cold seep environment. Our results shed light on the different evolutionary strategies during colonization in different micro-environments.

Mitogenome-wise codon usage pattern from comparative analysis of the first mitogenome of Blepharipa sp. (Muga uzifly) with other Oestroid flies

Uziflies (Family: Tachinidae) are dipteran endoparasites of sericigenous insects which cause major economic loss in the silk industry globally. Here, we are presenting the first full mitogenome of Blepharipa sp. (Acc: KY644698, 15,080 bp, A + T = 78.41%), a dipteran parasitoid of Muga silkworm (Antheraea assamensis) found in the Indian states of Assam and Meghalaya. This study has confirmed that Blepharipa sp. mitogenome gene content and arrangement is similar to other Tachinidae and Sarcophagidae flies of Oestroidea superfamily, typical of ancestral Diptera. Although, Calliphoridae and Oestridae flies have undergone tRNA translocation and insertion, forming unique intergenic spacers (IGS) and overlapping regions (OL) and a few of them (IGS, OL) have been conserved across Oestroidea flies. The Tachinidae mitogenomes exhibit more AT content and AT biased codons in their protein-coding genes (PCGs) than the Oestroidea counterpart. About 92.07% of all (3722) codons in PCGs of this new species have A/T in their 3rd codon position. The high proportion of AT and repeats in the control region (CR) affects sequence coverage, resulting in a short CR (Blepharipa sp.: 168 bp) and a smaller tachinid mitogenome. Our research unveils those genes with a high AT content had a reduced effective number of codons, leading to high codon usage bias. The neutrality test shows that natural selection has a stronger influence on codon usage bias than directed mutational pressure. This study also reveals that longer PCGs (e.g., nad5, cox1) have a higher codon usage bias than shorter PCGs (e.g., atp8, nad4l). The divergence rates increase nonlinearly as AT content at the 3rd codon position increases and higher rate of synonymous divergence than nonsynonymous divergence causes strong purifying selection. The phylogenetic analysis explains that Blepharipa sp. is well suited in the family of insectivorous tachinid maggots. It's possible that biased codon usage in the Tachinidae family reduces the effective number of codons, and purifying selection retains the core functions in their mitogenome, which could help with efficient metabolism in their endo-parasitic life style and survival strategy.

The Carniolan Honeybee from Slovenia—A Complete and Annotated Mitochondrial Genome with Comparisons to Closely Related Apis mellifera Subspecies

Simple Summary

The western honeybee, Apis mellifera, is a globally distributed bee species with many recognised subspecies, one of which is Apis mellifera carnica, the Carniolan honeybee. Apis m. carnica is native to southern Central Europe and parts of the Balkans, with the locus classicus in Slovenia. It is also widely popular with beekeepers in parts of Central and Northern Europe and other parts of the world, including the USA, Canada, and even New Zealand. In Slovenia, A. m. carnica is protected, with measures to conserve the subspecies’ autochthonous domestic population in place. Such efforts depend heavily upon genomic and phylogenetic information. In this study, we sequenced and annotated the mitochondrial genome of a specimen from Slovenia and compared the obtained data with a previously published sample of the A. m. carnica from Austria and the closely related Italian honeybee A. m. ligustica. We found several features unique to the new mitochondrial genome. We also phylogenetically analyzed the relationship between our sequence and the selected available A. mellifera mitochondrial sequences. The acquired position of the sequenced A. m. carnica from Slovenia on the phylogenetic tree brings new evidence for close relationships among C and O lineages and reflects their recent historical matrilinear ancestry.

Abstract

The complete mitochondrial genome of the Carniolan honeybee (Apis mellifera carnica) from Slovenia, a homeland of this subspecies, was acquired in two contigs from WGS data and annotated. The newly obtained mitochondrial genome is a circular closed loop of 16,447 bp. It comprises 37 genes (13 protein coding genes, 22 tRNA genes, and 2 rRNA genes) and an AT-rich control region. The order of the tRNA genes resembles the order characteristic of A. mellifera. The mitogenomic sequence of A. m. carnica from Slovenia contains 44 uniquely coded sites in comparison to the closely related subspecies A. m. ligustica and to A. m. carnica from Austria. Furthermore, 24 differences were recognised in comparison between A. m. carnica and A. m. ligustica subspecies. Among them, there are three SNPs that affect translation in the nd2, nd4, and cox2 genes, respectively. The phylogenetic placement of A. m. carnica from Slovenia within C lineage deviates from the expected position and changes the perspective on relationship between C and O lineages. The results of this study represent a valuable addition to the information available in the phylogenomic studies of A. mellifera—a pollinator species of worldwide importance. Such genomic information is essential for this local subspecies’ conservation and preservation as well as its breeding and selection.

Concerted and Independent Evolution of Control Regions 1 and 2 of Water Monitor Lizards (Varanus salvator macromaculatus) and Different Phylogenetic Informative Markers

Simple Summary

The evolutionary patterns and phylogenetic utility of duplicate control regions (CRs) in 72 individuals of Varanus salvator macromaculatus and other varanids have been observed. Divergence of the two CRs from each individual revealed a pattern of independent evolution in CRs of varanid lineage. This study is a first step towards developing new phylogenetic evolutionary models of the varanid lineage, with accurate evolutionary inferences to provide basic insights into the biology of mitogenomes.

Abstract

Duplicate control regions (CRs) have been observed in the mitochondrial genomes (mitogenomes) of most varanids. Duplicate CRs have evolved in either concerted or independent evolution in vertebrates, but whether an evolutionary pattern exists in varanids remains unknown. Therefore, we conducted this study to analyze the evolutionary patterns and phylogenetic utilities of duplicate CRs in 72 individuals of Varanus salvator macromaculatus and other varanids. Sequence analyses and phylogenetic relationships revealed that divergence between orthologous copies from different individuals was lower than in paralogous copies from the same individual, suggesting an independent evolution of the two CRs. Distinct trees and recombination testing derived from CR1 and CR2 suggested that recombination events occurred between CRs during the evolutionary process. A comparison of substitution saturation showed the potential of CR2 as a phylogenetic marker. By contrast, duplicate CRs of the four examined varanids had similar sequences within species, suggesting typical characteristics of concerted evolution. The results provide a better understanding of the molecular evolutionary processes related to the mitogenomes of the varanid lineage.

Characterization of the complete mitochondrial genome of Haemadipsa tianmushana Song 1977 (Hirudiniformes, Haemadipsidae) and its phylogenetic analysis

The complete mitochondrial genome of Haemadipsa tianmushana Song 1977 from China has been determined and reported for the first time in this study. It was 14,625 bp in length and consisted of 22 tRNA genes, 2 rRNA genes, 13 protein-coding genes (PCGs), and 3 control regions. The nucleotide base content of the complete mitogenome for this species was 35.1% A, 10.5% C, 11.6% G, and 42.8% T. The tRNA genes were ranged from 57 bp (Ser^TCT) to 66 bp (Gln^TTG) in length. The phylogenetic analyses indicated that Hirudinea is a mono-phyletic clade. And it includes Whitmania acranulata, Whitmania pigra, Whitmania laevis, Zeylanicobdella arugamensis, Ozobranchus jantseanus and Placobdella lamothei. In Hirudiniformes, H. tianmushana and three species of Haemopidae were obviously clustered into two independent branches. This result is consistent with a taxonomy that they all belong to the same suborder. This study adds to the genetic resources currently available for the species.

The investigation of intraspecific characteristics and comparative analyses of the complete mitochondrial genome of Stegobium paniceum (Linnaeus, 1758) (Coleoptera: Ptinidae) assembled from public NGS raw reads of the black truffle, Tuber melanosporum

Due to the rapid development of NGS technologies, a huge amount of NGS raw reads have been accumulated in public repositories, such as the Short Read Archive of NCBI. We successfully rescued the complete mitochondrial genome of Stegobium paniceum, a drug store beetle, from public NGS raw reads of truffle generated from the whole genome project. The circular mitogenome of S. paniceum is 15,474 bp long including 13 protein-coding genes, two ribosomal RNA genes, 22 transfer RNAs, and a single large non-coding region of 803 bp. Intraspecific transfer RNAs structure and sequence variations were investigated and simple sequence repeats identified from three S. paniceum mitochondrial genomes were compared showing their diversities as fundamental data to utilize them in various aspects including developing efficient molecular markers in the family, Ptinidae. Phylogenetic analysis of 23 Bostrichoidea mitochondrial genomes presented better species identification based on phylogenetic analyses and the optimal options for constructing phylogenetic trees based on Bostrichoidea mitochondrial genomes. Our results present not only utilization of public NGS raw read sequences but also intraspecific features of S. paniceum mitochondrial genomes and comparative analysis of Bostrichoidea mitochondrial genomes in various aspects.

Complete mitochondrial genome sequencing of Oxycarenus laetus (Hemiptera: Lygaeidae) from two geographically distinct regions of India

Oxycarenus laetus is a seed-sap sucking pest affecting a variety of crops, including cotton plants. Rising incidence and pesticide resistance by O. laetus have been reported from India and neighbouring countries. In this study, O. laetus samples were collected from Bhatinda and Coimbatore (India). Pure mtDNA was isolated and sequenced using Illumina MiSeq. Both the samples were found to be identical species (99.9%), and the complete genome was circular (15,672 bp), consisting of 13 PCGs, 2 rRNA, 23 tRNA genes, and a 962 bp control region. The mitogenome is 74.1% AT-rich, 0.11 AT, and - 0.19 GC skewed. All the genes had ATN as the start codon except cox1 (TTG), and an additional trnT was predicted. Nearly all tRNAs folded into the clover-leaf structure, except trnS1 and trnV. The intergenic space between trnH and nad4, considered as a synapomorphy of Lygaeoidea, was displaced. Two 5 bp motifs AATGA and ACCTA, two tandem repeats, and a few microsatellite sequences, were also found. The phylogenetic tree was constructed using 36 mitogenomes from 7 super-families of Hemiptera by employing rigorous bootstrapping and ML. Ours is the first study to sequence the complete mitogenome of O. laetus or any Oxycarenus species. The findings from this study would further help in the evolutionary studies of Lygaeidae.

Non-coding Regions of Mitochondrial DNA and the cox1 Gene Reveal Genetic Variability Among Local Belarusian Populations of the Causative Agent of Cercarial Dermatitis, Bird Schistosome Trichobilharzia szidati (Digenea: Schistosomatidae)

INTRODUCTION

The cercariae of avian blood flukes Trichobilharzia szidati (Digenea, Schistosomatidae) are known to cause cercarial allergic dermatitis ("swimmer's itch") in humans. Global epidemics can have significant impacts on local tourism-related economies in recreational areas. Little is known about the genetic polymorphism of the parasite population, or about the variability of the non-coding regions of mitochondrial DNA (mtDNA) and the possibility of using this as a genetic marker.

MATERIALS AND METHODS

The T. szidati cercariae were collected over 7 years from 33 naturally infected Lymnaea stagnalis snails from five sites at two neighboring lakes in Belarus. We investigated the variability of the short (SNR) and long (LNR) non-coding regions of mt DNA and the genetic diversity within the 1125-bp sequences of the gene for subunit 1 of cytochrome c oxidase (cox1).

RESULTS

In the SNR sequences, we found only length variability caused by changes in the number of bases in the mononucleotide tracts T6-T8. LNR demonstrates high variability in nucleotide sequence length (182-260 bp) depending on the presence of two long deletions of 59 and 78 nucleotides. Both mitochondrial loci (LNR and cox1) are characterized by high haplotype diversity (H = 0.922 and H = 1.0, respectively); the nucleotide diversity is significantly higher for LNR (π = 1.926 ± 0.443) compared to cox1 (π = 0.704 ± 0.059). Phylogenetic reconstructions based on the variability of each of the loci (LNR and cox1) and their concatenated sequences revealed their shallow structure and the absence of a correlation between the distribution of single-nucleotide polymorphisms and the geographic origin of parasites from two Belarusian lakes. We identified at last four weakly sublineages in the phylogenetic pattern of T. szidati. The carriers of each deletion have specific patterns for each of the two loci and form their own phylogeographic sublineages. An association between two fixed LNR substitutions and a fixed non-synonymous substitution in cox1 was found in four representatives of one lineage that had a short deletion in the LNR.

CONCLUSIONS

This study clarified the phylogeographic structure of the Belarusian population of T. szidati. Our data provide the basis for the use two mt markers in large-scale population studies of the parasite, as well as for studying the molecular evolution of coding and non-coding mtDNA in trematodes.

pgHMA: Application of the heteroduplex mobility assay analysis in phylogenetics and population genetics

The heteroduplex mobility assay (HMA) has proven to be a robust tool for the detection of genetic variation. Here, we describe a simple and rapid application of the HMA by microfluidic capillary electrophoresis, for phylogenetics and population genetic analyses (pgHMA). We show how commonly applied techniques in phylogenetics and population genetics have equivalents with pgHMA: phylogenetic reconstruction with bootstrapping, skyline plots, and mismatch distribution analysis. We assess the performance and accuracy of pgHMA by comparing the results obtained against those obtained using standard methods of analyses applied to sequencing data. The resulting comparisons demonstrate that: (a) there is a significant linear relationship (R²  = .992) between heteroduplex mobility and genetic distance, (b) phylogenetic trees obtained by HMA and nucleotide sequences present nearly identical topologies, (c) clades with high pgHMA parametric bootstrap support also have high bootstrap support on nucleotide phylogenies, (d) skyline plots estimated from the UPGMA trees of HMA and Bayesian trees of nucleotide data reveal similar trends, especially for the median trend estimate of effective population size, and (e) optimized mismatch distributions of HMA are closely fitted to the mismatch distributions of nucleotide sequences. In summary, pgHMA is an easily-applied method for approximating phylogenetic diversity and population trends.

Characterization of Penaeus vannamei mitogenome focusing on genetic diversity

The diversity of the Penaeus vannamei mitochondrial genome has still been poorly characterized, there are no validated mitochondrial markers available for populational studies, and the heteroplasmy has not yet been investigated in this species. In this study, metagenomic reads extracted from the muscle of a single individual were used to assemble the mitochondrial genome (mtDNA). These data associated with mitochondrial genomes previously described allowed to evaluate the inter-individual variability and heteroplasmy. Comparison among 45 mtDNA control regions led to the detection of conserved and variable segments and the characterization of two hypervariable regions. The analysis of diversity revealed mostly low frequency polymorphisms, and heteroplasmy was found in practically all mitochondrial genes, with a high occurrence of indels. These results indicate that the design of mitochondrial markers for P. vannamei must be done with caution. The mapping of conserved and variable regions and the characterization of heteroplasmy presented here will contribute to increasing the efficiency of mitochondrial markers for population or individual studies.

Mitogenomes Reveal Alternative Initiation Codons and Lineage-Specific Gene Order Conservation in Echinoderms

The mitochondrial genetic code is much more varied than the standard genetic code. The invertebrate mitochondrial code, for instance, comprises six initiation codons, including five alternative start codons. However, only two initiation codons are known in the echinoderm and flatworm mitochondrial code, the canonical ATG and alternative GTG. Here, we analyzed 23 Asteroidea mitogenomes, including ten newly sequenced species and unambiguously identified at least two other start codons, ATT and ATC, both of which also initiate translation of mitochondrial genes in other invertebrates. These findings underscore the diversity of the genetic code and expand upon the suite of initiation codons among echinoderms to avoid erroneous annotations. Our analyses have also uncovered the remarkable conservation of gene order among asteroids, echinoids, and holothuroids, with only an interchange between two gene positions in asteroids over ∼500 Ma of echinoderm evolution.

Studying genetic population structure to shed light on the demographic explosion of the rare species Barbitistes vicetinus (Orthoptera, Tettigoniidae)

Insect outbreaks usually involve important ecological and economic consequences for agriculture and forestry. The short-winged bush-cricket Barbitistes vicetinus Galvagni & Fontana, 1993 is a recently described species that was considered rare until ten years ago, when unexpected population outbreaks causing severe defoliations across forests and crops were observed in north-eastern Italy. A genetic approach was used to analyse the origin of outbreak populations. The analysis of two mitochondrial regions (Cytochrome Oxidase I and II and 12S rRNA-Control Region) of 130 samples from the two disjunct ranges (Euganean and Berici Hills) showed high values of haplotype diversity and revealed a high geographical structure among populations of the two ranges. The high genetic variability observed supports the native origin of this species. In addition, results suggest that unexpected outbreaks are not a consequence of a single or few pestiferous haplotypes but rather the source of outbreaks are local populations which have experienced an increase in each area. The recent outbreaks have probably appeared independently of the genetic haplotypes whereas environmental conditions could have affected the outbreak populations. These findings contribute to a growing understanding of the status and evolutionary history of the pest that would be useful for developing and implementing biological control strategies for example by maximizing efforts to locate native natural enemies.

Complete vertebrate mitogenomes reveal widespread repeats and gene duplications

BACKGROUND

Modern sequencing technologies should make the assembly of the relatively small mitochondrial genomes an easy undertaking. However, few tools exist that address mitochondrial assembly directly.

RESULTS

As part of the Vertebrate Genomes Project (VGP) we develop mitoVGP, a fully automated pipeline for similarity-based identification of mitochondrial reads and de novo assembly of mitochondrial genomes that incorporates both long (> 10 kbp, PacBio or Nanopore) and short (100-300 bp, Illumina) reads. Our pipeline leads to successful complete mitogenome assemblies of 100 vertebrate species of the VGP. We observe that tissue type and library size selection have considerable impact on mitogenome sequencing and assembly. Comparing our assemblies to purportedly complete reference mitogenomes based on short-read sequencing, we identify errors, missing sequences, and incomplete genes in those references, particularly in repetitive regions. Our assemblies also identify novel gene region duplications. The presence of repeats and duplications in over half of the species herein assembled indicates that their occurrence is a principle of mitochondrial structure rather than an exception, shedding new light on mitochondrial genome evolution and organization.

CONCLUSIONS

Our results indicate that even in the "simple" case of vertebrate mitogenomes the completeness of many currently available reference sequences can be further improved, and caution should be exercised before claiming the complete assembly of a mitogenome, particularly from short reads alone.

Nuclear (18S-28S rRNA) and mitochondrial genome markers of Carios (Carios) vespertilionis (Argasidae) support Carios Latreille, 1796 as a lineage embedded in the Ornithodorinae: re-classification of the Carios sensu Klompen and Oliver (1993) clade into its respective subgenera

Argasid systematics remains controversial with widespread adherence to the Hoogstraal (1985) classification scheme, even though it does not reflect evolutionary relationships and results in paraphyly for the main genera of soft ticks (Argasidae), namely Argas and Ornithodoros. The alternative classification scheme, proposed by Klompen and Oliver (1993), has problems of its own: most notably paraphyly of the subgenus Pavlovskyella and the controversial grouping together of the subgenera Alectorobius, Antricola, Carios, Chiropterargas, Nothoaspis, Parantricola, Reticulinasus and Subparmatus into the genus Carios. Recent phylogenetic analyses of 18S/28S rRNA sequences and mitochondrial genomes agree with the scheme of Klompen and Oliver (1993), with regard to the paraphyly of Pavlovskyella, placement of Alveonasus, Ogadenus, Proknekalia and Secretargas in the Argasinae and placement of Carios and Chiropterargas in the Ornithodorinae (Mans et al., 2019). The Carios clade and its constituent subgenera remain controversial, since the phylogenetic position of its type species Carios (Carios) vespertilionis Latreille, 1796 (formerly Argas vespertilionis) has not been determined with confidence. The current study aimed to resolve Carios sensu lato Klompen and Oliver, 1993, and Carios sensu stricto Hoogstraal, 1985, by determining and analysing phylogenetic nuclear and mitochondrial markers for C. (C.) vespertilionis. Both the nuclear and mitochondrial markers support placement of Carios s.s. within the subfamily Ornithodorinae, but to the exclusion of the clade that includes the 6 other subgenera that are part of Carios s.l. Klompen and Oliver (1993), namely Alectorobius, Antricola, Nothoaspis, Parantricola, Reticulinasus and Subparmatus. These 6 subgenera form a monophyletic clade that might be placed as new subgenera within the genus Alectorobius, or elevated to genera. Given the substantial differences in biology among these subgenera, we propose that these 6 subgenera be elevated to genera. Thus, we propose to modify the classification scheme of Mans et al. (2019) so that the subfamily Argasinae now has six genera, Alveonasus, Argas (subgenera Argas and Persicargas), Navis, Ogadenus, Proknekalia and Secretargas, and the subfamily Ornithodorinae has nine genera, Alectorobius, Antricola (subgenera Antricola and Parantricola), Carios, Chiropterargas, Nothoaspis, Ornithodoros (subgenera Microargas, Ornamentum, Ornithodoros, Pavlovskyella and Theriodoros), Otobius, Reticulinasus and Subparmatus (genera indicated in bold).

Not all mitochondrial DNAs are made equal and the nucleus knows it

The oxidative phosphorylation (OXPHOS) system is the only structure in animal cells with components encoded by two genomes, maternally transmitted mitochondrial DNA (mtDNA), and biparentally transmitted nuclear DNA (nDNA). MtDNA‐encoded genes have to physically assemble with their counterparts encoded in the nucleus to build together the functional respiratory complexes. Therefore, structural and functional matching requirements between the protein subunits of these molecular complexes are rigorous. The crosstalk between nDNA and mtDNA needs to overcome some challenges, as the nuclear‐encoded factors have to be imported into the mitochondria in a correct quantity and match the high number of organelles and genomes per mitochondria that encode and synthesize their own components locally. The cell is able to sense the mito‐nuclear match through changes in the activity of the OXPHOS system, modulation of the mitochondrial biogenesis, or reactive oxygen species production. This implies that a complex signaling cascade should optimize OXPHOS performance to the cellular‐specific requirements, which will depend on cell type, environmental conditions, and life stage. Therefore, the mitochondria would function as a cellular metabolic information hub integrating critical information that would feedback the nucleus for it to respond accordingly. Here, we review the current understanding of the complex interaction between mtDNA and nDNA.

Phylogenies from mitochondrial genomes of 120 species of ticks: Insights into the evolution of the families of ticks and of the genus Amblyomma

The evolution and phylogenetic relationships of the ticks at both the family and genus levels are contested. The genus Amblyomma and its subgenera are in a state of flux; moreover, the relationships among the three tick families are controversial due to conflicting phylogenetic support for different arrangements of the three families of living ticks. With 18 newly sequenced mitochondrial (mt) genomes of ticks included, we executed the largest mt genome phylogenetic study of ticks so far. Phylogenetic trees were inferred from one sea spider mt genome, one horseshoe crab, five mite mt genomes and 146 tick mt genomes from 120 species: 153 mt genomes in total. Sixteen phylogenetic trees were inferred from 10 datasets using both maximum likelihood and Bayesian inference methods. We describe the first novel mt gene-arrangement for the metastriate Ixodidae in Amblyomma (Africaniella) transversale. Also, three unusual partial 16S rRNA gene inserts were found in the mt genome of Haemaphysalis (Alloceraea) kitaokai: we consider the possible role of past genome translocation events in the formation of these inserts. Our phylogenies revealed evidence that: (i) the genus Amblyomma is polyphyletic with respect to Amblyomma (Africaniella) transversale; (ii) the subgenus Aponomma is apparently embedded in the genus Amblyomma; (iii) Haemaphysalis (Segalia) parva and Haemaphysalis (Alloceraea) kitaokai form a clade to the exclusion of other Haemaphysalis species; and (iv) the phylogenetic position of the family Nuttalliellidae is unstable among phylogenies from different datasets.

Distant hybrids of Heliocidaris crassispina (♀) and Strongylocentrotus intermedius (♂): identification and mtDNA heteroplasmy analysis

Background

Distant hybridization between the sea urchin Heliocidaris crassispina (♀) and the sea urchin Strongylocentrotus intermedius (♂) was successfully performed under laboratory conditions. A new variety of hybrid sea urchin (HS hybrid) was obtained. However, the early-development success rates for the HS hybrids were significantly lower than those of purebred H. crassispina or S. intermedius offspring. In addition, it was difficult to distinguish the HS-hybrid adults from the pure H. crassispina adults, which might lead to confusion in subsequent breeding attempts. In this study, we attempted to develop a method to quickly and effectively identify HS hybrids, and to preliminarily investigate the molecular mechanisms underlying the poor early-development success rates in the HS hybrids.

Results

The hybrid sea urchins (HS hybrids) were identified both morphologically and molecularly. There were no significant differences in the test height to test diameter ratios between the HS hybrids and the parents. The number and arrangement of ambulacral pore pairs in the HS hybrids differed from those of the parental lines, which might serve as a useful morphological character for the identification of the HS hybrids. A primer pair that identified the HS hybrids was screened by comparing the mitochondrial genomes of the parental lines. Moreover, paternal leakage induced mitochondrial DNA heteroplasmy in the HS hybrids, which might explain the low rates of early development success in these hybrids.

Conclusions

The distant-hybrid sea urchins were accurately identified using comparative morphological and molecular genetic methods. The first evidence of mtDNA heteroplasmy after the distant hybridization of an echinoderm was also provided.

The first complete mitochondrial genome of the sand dollar Sinaechinocyamus mai (Echinoidea: Clypeasteroida)

Morphologic and molecular data often lead to different hypotheses of phylogenetic relationships. Such incongruence has been found in the echinoderm class Echinoidea. In particular, the phylogenetic status of the order Clypeasteroida is not well resolved. Complete mitochondrial genomes are currently available for 29 echinoid species, but no clypeasteroid had been sequenced to date. DNA extracted from a single live individual of Sinaechinocyamus mai was sequenced with 10× Genomics technology. This first complete mitochondrial genome (mitogenome) for the order Clypeasteroida is 15,756 base pairs in length. Phylogenomic analysis based on 34 ingroup taxa belonging to nine orders of the class Echinoidea show congruence between our new genetic inference and published trees based on morphologic characters, but also includes some intriguing differences that imply the need for additional investigation.

de Menezes E, Tamang P, Hamid Z, Naik A, Parsade CK, Sivaprasad S. The Detection and Partial Localisation of Heteroplasmic Mutations in the Mitochondrial Genome of Patients with Diabetic Retinopathy

Diabetic retinopathy (DR) is a common complication of diabetes and a major cause of acquired blindness in adults. Mitochondria are cellular organelles involved in energy production which contain mitochondrial DNA (mtDNA). We previously showed that levels of circulating mtDNA were dysregulated in DR patients, and there was some evidence of mtDNA damage. In the current project, our aim was to confirm the presence of, and determine the location and prevalence of, mtDNA mutation in DR. DNA isolated from peripheral blood from diabetes patients (n = 59) with and without DR was used to amplify specific mtDNA regions which were digested with surveyor nuclease S1 to determine the presence and location of heteroplasmic mtDNA mutations were present. An initial screen of the entire mtDNA genome of 6 DR patients detected a higher prevalence of mutations in amplicon P, covering nucleotides 14,443 to 1066 and spanning the control region. Further analysis of 42 subjects showed the presence of putative mutations in amplicon P in 36% (14/39) of DR subjects and in 10% (2/20) non-DR subjects. The prevalence of mutations in DR was not related to the severity of the disease. The detection of a high-prevalence of putative mtDNA mutations within a specific region of the mitochondrial genome supports the view that mtDNA damage contributes to DR. The exact location and functional impact of these mutations remains to be determined.

Aquatic Hemiptera in Southwest Cameroon: Biodiversity of Potential Reservoirs of Mycobacterium ulcerans and Multiple Wolbachia Sequence Types Revealed by Metagenomics

Buruli ulcer (BU), caused by Mycobacterium ulcerans, is a neglected tropical disease associated with freshwater habitats. A variety of limnic organisms harbor this pathogen, including aquatic bugs (Hemiptera: Heteroptera), which have been hypothesized to be epidemiologically important reservoirs. Aquatic Hemiptera exhibit high levels of diversity in the tropics, but species identification remains challenging. In this study, we collected aquatic bugs from emerging foci of BU in the Southwest Region of Cameroon, which were identified using morphological and molecular methods. The bugs were screened for mycobacterial DNA and a selection of 20 mycobacteria-positive specimens from the families Gerridae and Veliidae were subjected to next-generation sequencing. Only one individual revealed putative M. ulcerans DNA, but all specimens contained sequences from the widespread alpha-proteobacterial symbiont, Wolbachia. Phylogenetic analysis placed the Wolbachia sequences into supergroups A, B, and F. Circularized mitogenomes were obtained for seven gerrids and two veliids, the first from these families for the African continent. This study suggests that aquatic Hemiptera may have a minor role (if any) in the spread of BU in Southwest Cameroon. Our metagenomic analysis provides new insights into the incursion of Wolbachia into aquatic environments and generated valuable resources to aid molecular taxonomic studies of aquatic Hemiptera.

The AC4 Protein of a Cassava Geminivirus Is Required for Virus Infection

Geminiviruses (family Geminiviridae) are among the most devastating plant viruses worldwide, causing severe damage in crops of economic and subsistence importance. These viruses have very compact genomes and many of the encoded proteins are multifunctional. Here, we investigated the role of the East African cassava mosaic Cameroon virus (EACMCV) AC4 on virus infectivity in Nicotiana benthamiana. Results showed that plants inoculated with EACMCV containing a knockout mutation in an AC4 open reading frame displayed symptoms 2 to 3 days later than plants inoculated with wild-type virus, and these plants recovered from infection, whereas plants inoculated with the wild-type virus did not. Curiously, when an additional mutation was made in the knockout mutant, the resulting double mutant virus completely failed to cause any apparent symptoms. Interestingly, the role of AC4 on virus infectivity appeared to be dependent on an encoded N-myristoylation motif that mediates cell membrane binding. We previously showed that EACMCV containing the AC4T38I mutant produced virus progeny characterized by second-site mutations and reversion to wild-type virus. These results were confirmed in this study using additional mutations. Together, these results show involvement of EACMCV AC4 in virus infectivity; they also suggest a role for the combined action of mutation and selection, under prevailing environmental conditions, on begomovirus genetic variation and diversity.

A phylogenomic resolution of the sea urchin tree of life

BACKGROUND

Echinoidea is a clade of marine animals including sea urchins, heart urchins, sand dollars and sea biscuits. Found in benthic habitats across all latitudes, echinoids are key components of marine communities such as coral reefs and kelp forests. A little over 1000 species inhabit the oceans today, a diversity that traces its roots back at least to the Permian. Although much effort has been devoted to elucidating the echinoid tree of life using a variety of morphological data, molecular attempts have relied on only a handful of genes. Both of these approaches have had limited success at resolving the deepest nodes of the tree, and their disagreement over the positions of a number of clades remains unresolved.

RESULTS

We performed de novo sequencing and assembly of 17 transcriptomes to complement available genomic resources of sea urchins and produce the first phylogenomic analysis of the clade. Multiple methods of probabilistic inference recovered identical topologies, with virtually all nodes showing maximum support. In contrast, the coalescent-based method ASTRAL-II resolved one node differently, a result apparently driven by gene tree error induced by evolutionary rate heterogeneity. Regardless of the method employed, our phylogenetic structure deviates from the currently accepted classification of echinoids, with neither Acroechinoidea (all euechinoids except echinothurioids), nor Clypeasteroida (sand dollars and sea biscuits) being monophyletic as currently defined. We show that phylogenetic signal for novel resolutions of these lineages is strong and distributed throughout the genome, and fail to recover systematic biases as drivers of our results.

CONCLUSIONS

Our investigation substantially augments the molecular resources available for sea urchins, providing the first transcriptomes for many of its main lineages. Using this expanded genomic dataset, we resolve the position of several clades in agreement with early molecular analyses but in disagreement with morphological data. Our efforts settle multiple phylogenetic uncertainties, including the position of the enigmatic deep-sea echinothurioids and the identity of the sister clade to sand dollars. We offer a detailed assessment of evolutionary scenarios that could reconcile our findings with morphological evidence, opening up new lines of research into the development and evolutionary history of this ancient clade.