Bayesian Phylogenetic Analyses Research Articles

Rectifying a century of misidentifications: a taxonomic re-evaluation of sipo snakes (Colubridae: Chironius) on Trinidad

The island of Trinidad has several endemic reptile species that in some cases are morphologically indistinguishable, or almost so, from their mainland counterparts. In particular, many snakes from the island have not been examined thoroughly with modern scientific methods and may therefore be misidentified. At least two species of sipo snakes (Chironius carinatus, C. septentrionalis) are reported to inhabit Trinidad, both being considered species already inhabiting the mainland, though their identities are based solely on morphology. Here, we evaluate the molecular distinctiveness of these Trinidadian snakes and assess their relationships to other members of the genus. We constructed a multi-locus data set (12S rRNA, 16S rRNA, cyt b, ND4, cmos, NT3, Rag-1, Rag-2) including novel sequences and those available on GenBank to perform Bayesian and maximum likelihood phylogenetic analyses. Our phylogenetic reconstruction elucidates the identity of both Trinidad species, which have been misidentified since their discovery, and we provide a literature review of their taxonomic history. More specifically, our results suggest that the Trinidad population of Chironius carinatus is an undescribed species, that we describe herein as Chironius nigelnoriegai sp. nov., and is the sister to C. flavopictus. Our results also demonstrate that the Trinidad population previously identified as C. septentrionalis is instead the mainland species C. cochranae, a member of a clade with C. foveatus, C. laurenti, and C. multiventris. Finally, we generated a time tree and inferred that our new species separated from C. flavopictus approximately 4 million years ago in the Pliocene, a time when the island of Trinidad detached from Northern Venezuela. This work contributes to a better understanding of species diversity of Trinidad and we hope that it assists in conservation efforts towards this important endemic region. These findings support the prospect of rapid speciation on Trinidad and implies that more extensive surveys of island squamates will reveal additional cryptic diversity. http://zoobank.org/urn:lsid:zoobank.org:pub:500DB91B-62CD-4051-9831-0BEEAC64C8E2

Characterization of the Complete Mitochondrial Genome of Salassa thespis (Lepidoptera: Saturniidae) and Comparison with Other Lepidoptera Species

There are many species of Lepidoptera, but few complete mitochondrial genomes of Lepidoptera have been included in databases. Here, the complete mitochondrial genome sequence of Salassa thespis was isolated and characterized. It was 15,302 bp in length and contained 13 protein-coding genes (PCGs), two rRNA genes, 22 tRNA genes and an A + T-rich region. Among the 13 PCGs, the initiation codon of cytochrome c oxidase subunit 1 (cox1) was CGA, and the rest were ATN. The cox1 and cox2 genes had an incomplete stop codon T, while the rest terminated with TAA. Codon usage analysis showed that Phe, Ile, Leu and Asn were the most frequent amino acids, while Trp was the least. Like other Lepidopterans, some conserved motifs were found in the A + T-rich region, including a 17 bp poly-T guided by ATAGA, the AT-rich area and a poly-A element. Bayesian inference and maximum likelihood phylogenetic tree analysis based on 13 PCGs of S. thespis confirmed that it belonged to the Saturniidae family and showed the following relationship: (S. thespis + (Rhodinia fugax + Samia canningi)). The enrichment of mitochondrial DNA provides reference information for the study of the evolution and diversity of Lepidoptera insects.

HIV-1 transmission dynamics among people who inject drugs on the US/Mexico border during the COVID-19 pandemic: a prospective cohort study

We examined HIV prevalence and transmission dynamics among people who inject drugs in the U.S./Mexico border region during the COVID-19 pandemic. People who inject drugs aged ≥18 years from 3 groups were recruited: people who inject drugs who live in San Diego (SD) and engaged in cross-border drug use in Tijuana, Mexico (SD CBDUs), and people who inject drugs in SD and Tijuana (TJ) who did not engage in cross-border drug use (NCBDUs). We computed HIV prevalence at baseline and bivariate incidence-density rates (IR) at 18-month follow-up. Bayesian phylogenetic analysis was used to identify local transmission clusters, estimate their age, and effective reproductive number (Re) over time within the clusters. At baseline (n=612), 26% of participants were female, 9% engaged in sex work, and HIVprevalencewas 8% (4% SD CBDU, 4% SD NCBDU, 16% TJ NCBDU). Nine HIV seroconversions occurred over 18 months, IR: 1.357 per 100 person-years (95% CI: 0.470, 2.243); 7 in TJ NCBDU and 2 in SD CBDU. Out of 16 identified phylogenetic clusters, 9 (56%) had sequences from both the U.S. and Mexico (mixed-country). The age of three youngest mixed-country dyads (2018-2021) overlapped with the COVID-related US-Mexico border closure in 2020. One large mixed-country cluster (N=15) continued to grow during the border closure (Re=4.8, 95% Highest Posterior Density (HPD) 1.5-9.1) with 47% engaging in sex work. Amidst the COVID-19 pandemic and the border closure, cross-border HIV clusters grew. Efforts to end the HIV epidemic in the U.S. should take into account cross-border HIV-1 transmission from Tijuana. Mobile harm reduction services and coordination with municipal HIV programs to initiate anti-retroviral therapy and pre-exposure prophylaxisis are needed to reduce transmission. This research was supported by the James B. Pendleton Charitable Trust and the San Diego Center for AIDS Research.

Molecular identification and phylogenetic analysis of the mitogenome in endangered giant nuthatch Sitta magna (Passeriformes, Sittidae)

The Giant Nuthatch Sitta magna (family Sittidae) is a passerine bird, the quantification of the number of habitats and species on a global scale remains low. Most species are restricted to low elevations in southwest China, eastern Myanmar, and northern Thailand. To characterize the mitochondrial genome sequence of S. magna and its phylogenetic relationships with other members within the genus Sitta, the mitochondrial genome of S. magna was sequenced using the whole genome shotgun method. The sequencing results showed that the mitochondrial genome was 16,829 bp long and consisted of 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), 2 ribosomal RNA genes (rRNAs), and one control region (D-loop). All tRNAs were predicted to form a typical clover secondary structure. Among the 13 PCGs, only the start codon in COI was ATC, the start codon by the remaining 12 PCGs was ATG, and the stop codons were TAG, TAA, AGG, AGA, and TA. Bayesian inference and maximum likelihood phylogenetic analysis of the sequences of 17 species generated consistent well-supported phylogenies. The family Polioptilidae and the family Troglodytidae were closely related, and the family Sittidae was confined to a single branch. The genus Sitta in the family Sittidae was mainly clustered into three branches. Our findings provide new mitochondrial genomic data that could be used for phylogenetic and taxonomic studies; our results also certificate into the phylogenetic relationships within the genus Sitta ((S. himalayensi+(S. nagaensis + S. europaea))+(S. villosa + S. yunnanensis))+(S. carolinensis + S. magna).

Post-mortem detection of hemoplasmas (hemotropic Mycoplasma spp.) in South American fur seal (Arctocephalus australis) sampled in Rio Grande do Sul State, southern Brazil

Hemotropic mycoplasmas are bacteria that attaches to erythrocytes surface, which some species presents zoonotic concerns. In the suborder Pinnipedia, genera Otaria and Arctocephalus are prominent in Brazil. This study investigated the occurrence of hemoplasmas in Arctocephalus sp. and Otaria flavescens found dead along the coast of a Southern Brazilian State. DNA from 135 spleen samples were extracted and subjected to conventional PCR protocols, targeting the 16 S rRNA and 23 S rRNA gene. Three (2.22 %) Arctocephalus australis were positive in the 16 S rRNA gene, and no samples amplified in the 23 S rRNA gene. Samples from this study clustered with Zalophus californianus and Arctocephalus tropicalis mycoplasmas on a Bayesian phylogenetic analysis. Genetic diversity analysis suggested distinct genotypes, indicating A. australis as a new host for hemoplasma, and also a potential putative novel hemoplasma genotype. These findings raises future awareness for pinnipeds conservation, and adds Mycoplasma spp. to be taken into consideration when clinically evaluating rescued animals.

First Byzantinia from Afro-Arabia and the evolutionary history of extinct cricetodontine rodents investigated through Bayesian phylogenetic inference

The Cricetodontinae is a subfamily of extinct rodents from Europe, northern Africa, and Asia. Byzantinia is a cricetodontine genus that currently includes 10 species. This genus originated in Anatolia in the late Middle Miocene and became extinct at the end of the Miocene. Until now, Byzantinia was known only from Anatolia and Greece. Here, we report isolated molars of Byzantinia from a new locality close to Zahleh, Lebanon. These represent the first remains of Byzantinia from the Arabian area. All these teeth but one, which is very close to Byzantinia dardanellensis, belong to the new taxon Byzantinia rosamariae sp. nov. The evolutionary position of the two Lebanese Byzantinia species allows us to infer an age for this site of approximately 10.9–10 Ma. Bayesian tip-dating phylogenetic analyses, using morphological data for a large variety of cricetodontines, revealed that Anatolia served as the primary centre for speciation and dispersal of this rodent subfamily. This approach allowed us to infer the biogeographical history and time of divergence of the main cricetodontine lineages and sheds light on the evolutionary history of this group. The evolutionary position of the new Lebanese species is intermediate between Byzantinia sp. from Direcik I (Turkey) and B. nikosi from Biodrak (Greece), from which it diverged c. 11.85 Ma. These new Lebanese fossils provide the first physical evidence of the dispersal of these rodents into Afro-Arabia. http://zoobank.org/urn:lsid:zoobank.org:pub:2D8FD3C8-26E1-4C62-8D65-9D58666D70A3

Preliminary study on buffy coat smear and molecular detection of microfilaria in domestic chickens (Gallus gallus domesticus) raised in Southern Thailand

Background and Aim: Filarial nematode typically produces a larval stage (microfilariae) in the bloodstream of vertebrate hosts, where microfilariae reside in the blood or subcutaneous tissues. Filarial nematodes cause human diseases, such as river blindness and elephantiasis, which are widely studied. However, in avian species, they are overlooked because they are nonpathogenic. In Thailand, microfilaria can be found in wild birds and domestic chickens. Recently, an increase in the number of blood samples submitted to veterinary diagnostic laboratories may have increased the number of microfilariae. Therefore, knowledge about filarial species and reliable detection methods are important. Therefore, this study aimed to investigate the efficacy of buffy coat smear and polymerase chain reaction (PCR)-based methods for the detection of microfilaria in domestic chickens. In addition, parasites were identified using the sequence of the cytochrome c oxidase subunit 1 (COX1) gene. Materials and Methods: Giemsa-stained buffy coat smears from a previous study were reanalyzed. These available buffy coat smears were prepared from 55 domestic chickens raised as backyard free-ranging in Southern Thailand. Fifty-seven frozen genomic DNA extracted from chicken blood were used to detect the presence of the COX1 gene in Onchocercidae nematodes. The nested PCR protocol for amplification of the OnchoCOI_ R2-OnchoCOI_ R2 fragment of the COX1 gene was applied from a previous report. Sequences of COX1 were analyzed to identify Onchocercidae nematodes and if they were single or mixed infections. We constructed Bayesian phylogenetics to identify parasites and assessment of the relationship between filarial nematodes in avian species and other vertebrate hosts. Results: Buffy coat smears from 15 samples revealed microfilaria. Of these 15 samples, only eight were positive for COX1 nested-PCR amplification. The other two buffy coat-negative samples were also positive for nested-PCR. Sequencing of these 11 nested PCR-positive samples revealed that almost all of them were Onchocercidae nematodes. Bayesian phylogenetic analysis showed that chicken Onchocercidae spp. were grouped with other avian filarial nematodes. However, all chickens Onchocercidae spp. showed a double peak in the sequencing chromatogram, indicating mixed filarial infection (species or haplotypes). Therefore, no chicken Onchocercidae sequence was deposited on National Center for Biotechnology Information, GenBank. Conclusion: Giemsa-stained buffy coat smear was a reliable method for the detection of chicken microfilaria in routine veterinary diagnostic laboratories. Development of a new PCR-based method is necessary. This method may provide greater sensitivity and specificity of detection. In addition, the PCR method allowed us to access the genetic characteristics of nematodes, which helped us maximize our knowledge of nematodes. Further investigations, such as the pathogenicity of filarial nematodes in chickens and their potential vectors, are required. Keywords: Buffy coat smear, Chickens, Cytochrome c oxidase I, Microfilaria, Nested-polymerase chain reaction.

Discrepancy of flowering time between genetically close sublineages of Aegilops umbellulata Zhuk.

Aegilops umbellulata Zhuk., a wild diploid wheat-related species, has been used as a genetic resource for several important agronomic traits. However, its genetic variations have not been comprehensively studied. We sequenced RNA from 114 accessions of Ae. umbellulata to evaluate DNA polymorphisms and phenotypic variations. Bayesian clustering and phylogenetic analysis based on SNPs detected by RNA sequencing revealed two divergent lineages, UmbL1 and UmbL2. The main differences between them were in the sizes of spikes and spikelets, and culm diameter. UmbL1 is divided into two sublineages, UmbL1e and UmbL1w. These genetic differences corresponded to geographic distributions. UmbL1e, UmbL1w, and UmbL2 are found in Turkey, Iran/Iraq, and Greece, respectively. Although UmbL1e and UmbL1w were genetically similar, flowering time and other morphological traits were more distinct between these sublineages than those between the lineages. This discrepancy can be explained by the latitudinal and longitudinal differences in habitats. Specifically, latitudinal clines of flowering time were clearly observed in Ae. umbellulata, strongly correlated with solar radiation in the winter season. This observation implies that latitudinal differences are a factor in differences in the flowering times of Ae. umbellulata. Differences in flowering time could influence other morphological differences and promote genetic divergence between sublineages.

A new loach species of the genus Oreonectes (Teleostei, Cypriniformes, Nemacheilidae) from Guangxi, China.

A new loach species, Oreonectesandongensissp. nov. is described from the Guangxi Zhuang Autonomous Region, China. The new species can be differentiated from other members of the genus by combinations of characters: a developed posterior chamber of the swim bladder, 13-14 branched caudal-fin rays, 8-16 lateral-line pores, body width 12-15% of standard length (SL), interorbital width 42-47% of head length (HL), and caudal peduncle length 11-16% of SL. Bayesian inference phylogenetic analysis based on mitochondrial Cyt b provided strong support for validity of O.andongensissp. nov. (uncorrected p-distance 6.0-7.5%).

Complete genome characterization of cacao leafroll virus, a newly described cacao-infecting polerovirus.

The complete genome sequence of cacao leafroll virus (CaLRV; family Solemoviridae, genus Polerovirus) was determined by high-throughput sequencing of total RNA isolated from symptomatic cacao Theobroma cacao L. plants (n = 4). The CaLRV genome sequences ranged from 5,976 to 5,997 nucleotides (nt) in length and contained seven open reading frames (ORFs). Nucleotide and amino acid (aa) sequence comparisons showed that, among selected well-characterized poleroviruses, the CaLRV genome shared the highest nt sequence identity of 62% with that of potato leafroll virus (PLRV, NC_076505). A comparison of the predicted aa sequence of the CaLRV coat protein indicated that cotton leafroll dwarf virus (CLRDV, NC_014545) and melon aphid-borne yellows virus (MABYV, NC_010809) were the closest relatives, sharing 57% aa sequence identity. Bayesian phylogenetic analysis based on complete genome sequences showed that CaLRV grouped with well-characterized poleroviruses that cause diseases of cereal and vegetable crops. During the course of publishing this work, the nearly complete genome sequence of a member of the same polerovirus species, referred to as "cacao polerovirus" (OR605721), with which CaLRV shares 99% nt sequence identity, was reported.

Cystic echinococcosis in cattle and sheep caused by Echinococcus granulosus sensu stricto genotypes G1 and G3 in the USA

BackgroundEndemic domestic dog-ruminant cycles and human cystic echinococcosis caused by Echinococcus granulosus have been sporadically reported in the United States. However, there is a paucity of molecular data describing the genotypes and haplotypes of this important cestode in domestic ruminant hosts.MethodsNinety-four cysts from the lungs and/or livers of slaughtered beef cattle (76 samples), dairy cows (five samples) and sheep (13 samples) were collected from abattoirs in four states of the USA. Samples were genotyped at two mitochondrial loci, cox1 and nad5. Sequences were used to determine species, genotypes and haplotypes using median joining networks and Bayesian phylogenetic analyses. Cyst fertility was assessed in hematoxylin and eosin-stained sections. Additionally, previously reported autochthonous E. granulosus infections in the USA in various hosts were mapped.ResultsBased on cox1 sequences obtained from 94 cysts, 89 (94.7%) were identified as E. granulosus G1/G3, while five (5.3%) were Taenia hydatigena. Taenia hydatigena were only isolated from sheep. Based on nad5 sequences obtained from 89 hydatid cysts, 96.6% and 3.4% belonged to E. granulosus sensu stricto genotypes G1 and G3 respectively. Two haplotypes were found among E. granulosus cox1 sequences, neither of which was geographically unique. Six haplotypes were found among nad5 sequences in genotype G1, of which five were novel, while one haplotype was found in genotype G3. In the concatenated cox1-nad5 dataset, seven haplotypes were identified, of which six were geographically unique. All cysts from cattle were non-fertile. Four cysts from sheep were fertile.ConclusionsAll genotyped samples belonged to E. granulosus s.s. This is the first study to our knowledge to confirm the presence of genotypes G1 and G3 in domestic cattle and sheep intermediate hosts in the USA and provide data for future diagnostic and epidemiological studies. Sequences have been deposited in GenBank (cox1 sequences: OR398494-OR398496, nad5 sequences: OR400695-OR400702).Graphical

Terrestrial predatory leeches of the genus Orobdella (Hirudinea: Erpobdelliformes: Orobdellidae) endemic to the Southern Russian Far East: a new species of the genus from Primorsky Krai, Russia.

A quadrannulate species, Orobdella ganini sp. nov., is described from the Lazovsky Nature Reserve in Primorsky Krai, the Southern Russian Far East, Russia. Morphological features of O. ghilarovi Nakano & Prozorova, 2019 from the reserve are also provided leading to an amendment of the species diagnosis. Maximum likelihood and Bayesian inference phylogenetic analyses, which were performed using nuclear 18S rRNA, 28S rRNA, and histone H3, and mitochondrial cytochrome c oxidase subunit I, tRNACys, tRNAMet, 12S rRNA, tRNAVal, 16S rRNA, tRNALeu and NADH dehydrogenase subunit 1 markers, show that O. ganini sp. nov., O. ghilarovi and two species endemic to Hokkaido, Japan form a clade, with the new species sister to a lineage composed of the two Japanese species. A partial cytochrome c oxidase subunit I sequence obtained from a cocoon found in the Lazovsky Nature Reserve reveals that Orobdella leeches deposit cocoons somewhat similar to those deposited by terrestrial blood-sucking leeches of Haemadipsidae.

Let's get high: Cladogenesis in freshwater crabs (Decapoda: Potamonautidae: Potamonautes) supports the mountain gradient speciation hypothesis in the Cape Fold and Drakensberg Mountains, South Africa.

During the present study, the evolutionary relationship within a clade of mountain clade of freshwater crabs (Potamonautes) was examined using mtDNA sequence data for species from the Cape Fold Mountain (CFM) and Great Escarpment (Drakensberg Mountain range). We undertook phylogenetic analyses, divergence time estimation, and an ancestral area reconstruction to explore the period of cladogenesis and understand the biogeographic history in this high-altitude clade. Furthermore, we applied four species delimitation methods using ASAP, bPTP, bGMYC, and STACEY on the latter clade. Bayesian phylogenetic analyses retrieved a monophyletic freshwater crab clade comprised of two major sister clades, one comprised of the Cape Fold (clade A) and two comprised of Drakensberg Mountains (clade B) species. Divergence time estimation indicated that the two clades underwent Mio/Pliocene cladogenesis. Within the CFM clade (A), P. amathole (Amathola Mountains) was sister to P. parvispina (Cederberg and Kouebokkeveld Mountains) and the latter species were sister to P. parvicorpus (Cape Peninsula, Jonkershoek, and Helderberg Mountains) sister to P. tuerkayi (Overberg Mountains) and P. brincki (Hottentots Holland Mountains). Within the Drakensberg Mountain clade (B), we observed in situ diversification. Specimens from the southcentral Drakensberg Mountains (Dargle Forest, Injasuti, Karkloof, and Impendle) represent a new undescribed lineage Potamonautes sp. nov. 1. The second clade from the northern Drakensberg, representing P. clarus, was sister to a central Drakensberg Mountain clade that comprised P. depressus that was in turn sister to P. baziya from the Eastern Cape Province. The application of species delimitation methods generally overestimated the number of species. The biogeographic analyses indicated that the Eastern Cape Province is the most likely ancestral range area. Ecological niche modelling of representative species in clades A (Cape Fold Mountains) and B (Drakensberg Mountains) demonstrated that temperature and rainfall were the major abiotic drivers that differentiated the two clades. Our data favours the mountain gradient speciation hypothesis.

Echoes from the Cretaceous: new fossils shed light on the evolution of host detection and concealed ovipositor apparatus in the parasitoid wasp superfamily Orussoidea (Hymenoptera)

Abstract †Kryptovelona carstengroehni gen. et sp. nov. and †Orussus juttagroehnae sp. nov. are the first female members of the parasitoid wasp family Orussidae recorded from Baltic amber. We describe them, including relevant parts of the internal anatomy examined with synchrotron scanning. The fossils display a number of modifications in the antennae and foreleg correlated with the specialized host-detection mechanism, and in the ovipositor apparatus, as well as in the thorax and abdomen for accommodating the internalized ovipositor. The presence of these and other features places them as crown-group members of the Orussidae, as demonstrated by Bayesian phylogenetic analyses. The recently described stem-group orussoid fossils from Burmese amber, the probable female †Cretorussus vilhelmseni, and probable male †Burmorussus mirabilis (both placed in Burmorussidae), were also included in the dataset. By comparing the new Baltic amber taxa with †Cretorussus, it is possible to trace the progressive refinement of the echolocation mechanism through reductions in the number of antennomeres and foreleg tarsomeres. Unfortunately, †Cretorussus does not have the posterior part of the abdomen with the ovipositor preserved. Nevertheless, it is possible to infer that the putative echolocation mechanism for host detection evolved at least 100 Mya, whereas the concealed ovipositor apparatus has not been documented in fossils older than approx. 35 Mya.

Molecular evolution and phylogeographic analysis of wheat dwarf virus.

Wheat dwarf virus (WDV) has caused considerable economic loss in the global production of grain crops. Knowledge of the evolutionary biology and population history of the pathogen remain poorly understood. We performed molecular evolution and worldwide phylodynamic analyses of the virus based on the genes in the protein-coding region of the entire viral genome. Our results showed that host-driven and geography-driven adaptation are major factors that affects the evolution of WDV. Bayesian phylogenetic analysis estimates that the average WDV substitution rate was 4.240 × 10-4 substitutions/site/year (95% credibility interval, 2.828 × 10-4-5.723 × 10-4), and the evolutionary rates of genes encoding proteins with virion-sense transcripts and genes encoding proteins with complementary-sense transcripts were different. The positively selected sites were detected in only two genes encoding proteins with complementary-sense, and WDV-barley are subject to stronger purifying selection than WDV-wheat. The time since the most recent common WDV ancestor was 1746 (95% credibility interval, 1517-1893) CE. Further analyses identified that the WDV-barley population and WDV-wheat population experienced dramatic expansion-decline episodes, and the expansion time of the WDV-barley population was earlier than that of the WDV-wheat population. Our phylogeographic analysis showed that the WDV population originating in Iran was subsequently introduced to Europe, and then spread from Eastern Europe to China.

Geographic characteristics of HTLV-1 molecular subgroups and genetic substitutions in East Asia: Insights from complete genome sequencing of HTLV-1 strains isolated in Taiwan and Japan.

Although Japan is a major endemic area for human T-lymphotropic virus type 1 (HTLV-1) and the virus has been well-studied in this region, there is limited research on HTLV-1 in surrounding regions. In this study, we determined the complete genome sequences of HTLV-1 strains isolated from Taiwan and Japan and investigated the geographic characteristics of molecular subgroups and substitution mutations to understand the spread of HTLV-1 and its correlation with human migration. The complete genome sequences of 26 HTLV-1 isolates from Taiwan were determined using next-generation sequencing and were compared with those of 211 isolates from Japan in terms of subgroup and genetic mutations. In total, 15/26 (58%) isolates from Taiwan belonged to the transcontinental subgroup and 11/26 (42%) isolates belonged to the Japanese subgroup. The transcontinental subgroup was significantly more prevalent among Taiwanese isolates than Japanese isolates (58% vs 18%, P < 0.0001). The mutation rate for the complete HTLV-1 sequence was as low as 0.2%. On examining individual base substitutions, the G-to-A mutation was predominant. Bayesian phylogenetic tree analysis estimated the time to the most recent common ancestor for the transcontinental and Japanese subgroups to be 28447 years. The transcontinental subgroups from Taiwan and Japan appeared to form clusters according to their respective regions. The transcontinental subgroup of HTLV-1 is predominant in Taiwan, while the Japanese subgroup is common in Japan. The difference in subgroup distribution may be attributed to the initial spread of the transcontinental subgroup in East Asia, followed by the influx of the Japanese subgroup.

Mitogenomic Characterization and Phylogenetic Placement of African Hind, Cephalopholis taeniops: Shedding Light on the Evolution of Groupers (Serranidae: Epinephelinae).

The global exploration of evolutionary trends in groupers, based on mitogenomes, is currently underway. This research extensively investigates the structure of and variations in Cephalopholis species mitogenomes, along with their phylogenetic relationships, focusing specifically on Cephalopholis taeniops from the Eastern Atlantic Ocean. The generated mitogenome spans 16,572 base pairs and exhibits a gene order analogous to that of the ancestral teleost's, featuring 13 protein-coding genes (PCGs), two ribosomal RNA genes (rRNAs), 22 transfer RNA genes (tRNAs), and an AT-rich control region. The mitogenome of C. taeniops displays an AT bias (54.99%), aligning with related species. The majority of PCGs in the mitogenome initiate with the start codon ATG, with the exceptions being COI (GTG) and atp6 (TTG). The relative synonymous codon usage analysis revealed the maximum abundance of leucine, proline, serine, and threonine. The nonsynonymous/synonymous ratios were <1, which indicates a strong negative selection among all PCGs of the Cephalopholis species. In C. taeniops, the prevalent transfer RNAs display conventional cloverleaf secondary structures, except for tRNA-serine (GCT), which lacks a dihydrouracil (DHU) stem. A comparative examination of conserved domains and sequence blocks across various Cephalopholis species indicates noteworthy variations in length and nucleotide diversity. Maximum likelihood, neighbor-joining, and Bayesian phylogenetic analyses, employing the concatenated PCGs and a combination of PCGs + rRNAs, distinctly separate all Cephalopholis species, including C. taeniops. Overall, these findings deepen our understanding of evolutionary relationships among serranid groupers, emphasizing the significance of structural considerations in mitogenomic analyses.

First Report of Foliar Nematode (Aphelenchoides pseudobesseyi) on Soybean in Colombia

Soybean (Glycine max L.) is produced in over 70,000 ha in the Altillanura Region, eastern Colombia (Agronet 2023). From 2018 to 2020, foliar symptoms like green stem and foliar retention of soybean, which in Brazil can cause up to 100% soybean yield losses (Meyer et al. 2017), were observed in soybean fields in Colombia. During 2020, samples from symptomatic plants in reproductive stages (R1-R8) were collected from different commercial soybean fields in the Altillanura Region. Over 200 samples were processed, using an incubation method described in Coyne et al. (2014). Nematodes were recovered from photosynthetic leaf tissues and enlarged nodes/buds with population densities ranging from 13 to 132 and 36 to 936 nematodes/10g, respectively. Adult females were morphologically and molecularly characterized as Aphelenchoides pseudobesseyi (Oliveira et al. 2019; Subbotin et al. 2020). Female body length (n = 20) ranged from 653.3 to 806.3 μm (mean = 723 μm ± 52.7), stylet length from 11.0 to 12.3 μm (11.8 μm ± 0.3), body diameter from 14.8 to 17.9 μm (16.3 μm ± 1.1), post-uterine sac length from 38.7 to 51.9 μm (44.6 μm ± 5.1), vulva to anus from 145.5 to 223.2 μm (172.2 μm ± 22.4), and 26% of the vulva-anus distance. Genomic DNA was extracted (QIAGEN DNeasy® Blood &amp; Tissue kit) from a pool of nematodes. The D2A/D3B (Tenente et al. 2004) primers were used to amplify and sequence the D2/D3 expansion region of the 28S rRNA gene. PCR product (~759 bp) was purified, sequenced, deposited in GenBank (OQ930285), and compared to previously deposited sequences (e.g., KX356756, KY510840, KY510839, KY510841, KT692694, KY510842, MH187565) by means of the BLAST algorithm. Similarly, 988F and 18SR-Burs (De Jesus et al. 2016) primers were used to amplify and sequence the near full-length 18S RNA gene (SSU). PCR product was purified, sequenced, deposited in GenBank (OQ954344), and compared to previously deposited sequences (e.g., KT454962, KT943534, KT943535, KY510835, KY510836, KY510837, KY510838, MH187565). Phylogenetic Bayesian analysis (Ronquist et al. 2012) of the of the D2/D3 and 18S regions placed this nematode from Colombia in the A. pseudobesseyi clade (PP = 100). To fulfill a modified Koch’s postulates, the A. pseudobesseyi population described above was used in a greenhouse assay. In total, 120 soybean plants (cv. Flor Blanca) were infected with 200 A. pseudobesseyi (females + males)/plant. Briefly, at cotyledon stage (VC), 50 µl aliquot containing 50 A. pseudobesseyi was delivered onto each cotyledon and unifoliolate leaves (200 nematodes/plant). Sterile water was delivered to 80 plants which served as control. Plants were kept in the greenhouse at approximately 25°C and covered with clear plastic bag for 72 h to maintain over 90% relative humidity. After 15, 30, 45, and 60 days, soybean plants (n = 20) were processed, A. pseudobesseyi quantified, and the average reproduction factor (final population/initial population) was 0.1, 2.9, 14.0, and 1.8, respectively. Infected plants showed symptoms of blistering leaves with malformation (midrib vein twist), and A. pseudobesseyi was not observed in control plants. To our knowledge, this is the first report of A. pseudobesseyi parasitizing soybean buds and leaves in Colombia. Soybean is an important commodity for the Altillanura Region, and it is important to monitor the risk posed by this nematode. Furthermore, a better understanding of the nematode-host interaction and epidemiology in Colombia soybean producing regions is needed.

Transcending borders: Remarkable genetic homogeneity among Barred mudskipper (Periophthalmus argentilineatus) populations in the southwestern Indian Ocean

AbstractIn the past 24 years, the Southwestern Indian Ocean (SWIO) region has lost about 4% of its mangroves due to unsustainable extraction, land clearance for agriculture and climate change impacts. Since this loss risks fragmenting mangrove fauna, this study analysed 179 D‐loop sequences (329 base pairs) of Barred mudskipper (Periophthalmus argentilineatus) sampled from estuarine mangroves in the SWIO to test two hypotheses: (1) whether Barred mudskipper populations in the region display genetic connectivity along the fragmented mangroves and (2) whether these populations have experienced demographic declines in recent history. The populations showed low and insignificant indices of genetic differentiation (FST = −0.00068, p &gt; 0.05; ΦST = 0.005, p &gt; 0.05), and haplotypes from different localities did not cluster according to their geographical origins. Similarly, Bayesian phylogenetic analysis grouped all SWIO haplotypes into one cluster, suggesting the hypothesis of genetic interconnectedness could not be rejected. Furthermore, the populations showed high haplotype diversity (0.96 ± 0.007) combined with low nucleotide diversity (1.09 ± 0.61%). They also showed significant negative Tajima's D values (−1.71, p &lt; 0.02) and a unimodal mismatch distribution, indicating the hypothesis of demographic expansion post a historical bottleneck could not be rejected. Nevertheless, the Bayesian skyline plot indicated an increase in population size post‐bottleneck, suggesting an ongoing recovery. These findings indicate that, despite varying management approaches in the SWIO, Barred mudskipper populations exhibit genetic exchange that transcends geographical and administrative borders. This implies that their populations should be considered a single stock within the region. These findings underscore the need for coordinated management approaches among the SWIO countries. This is crucial because inconsistencies in conservation measures across borders could disrupt the population's interconnectedness, potentially impacting their genetic integrity and long‐term survival.

Bayesian Phylogenetic Analysis on multi-core Compute Architectures: Implementation and evaluation of BEAGLE in RevBayes with MPI.

Phylogenies are central to many research areas in biology and commonly estimated using likelihood based methods. Unfortunately, any likelihood based method, including Bayesian inference, can be restrictively slow for large datasets -with many taxa and/or many sites in the sequence alignment- or complex substitutions models. The primary limiting factor when using large datasets and/or complex models in probabilistic phylogenetic analyses is the likelihood calculation, which dominates the total computation time. To address this bottleneck, we incorporated the high-performance phylogenetic library BEAGLE into RevBayes, which enables multi-threading on multi-core CPUs and GPUs, as well as hardware specific vectorized instructions for faster likelihood calculations. Our new implementation of RevBayes+BEAGLE retains the flexibility and dynamic nature that users expect from vanilla RevBayes. Additionally, we implemented native parallelization within RevBayes without an external library using the message passing interface (MPI); RevBayes+MPI. We evaluated our new implementation of RevBayes+BEAGLE using multi-threading on CPUs and two different powerful GPUs (NVidia Titan V and NVIDIA A100) against our native implementation of RevBayes+MPI. We found good improvements in speedup when multiple cores were used, with up to 20-fold speedup when using multiple CPU cores and over 90-fold speedup when using multiple GPU cores. The improvement depended on the data type used, DNA or amino acids, and the size of the alignment, but less on the size of the tree. We additionally investigated the cost of rescaling partial likelihoods to avoid numerical underflow and showed that unnecessarily frequent and inefficient rescaling can increase runtimes up to 4-fold. Finally, we presented and compared a new approach to store partial likelihoods on branches instead of nodes which can speed up computations up to 1.7-times but comes at twice the memory requirements.