Martínez-Castilla L, García-Sandoval R, Ballester P, Ferrándiz C, Sanchez MDLP, García-Ponce B, Garay-Arroyo A, R. Álvarez-Buylla E. HDACs MADS-domain protein interaction: a case study of HDA15 and XAL1 in Arabidopsis thaliana

Cellular behavior, cell differentiation and ontogenetic development in eukaryotes result from complex interactions between epigenetic and classic molecular genetic mechanisms, with many of these interactions still to be elucidated. Histone deacetylase enzymes (HDACs) promote the interaction of histones with DNA by compacting the nucleosome, thus causing transcriptional repression. MADS-domain transcription factors are highly conserved in eukaryotes and participate in controlling diverse developmental processes in animals and plants, as well as regulating stress responses in plants. In this work, we focused on finding out putative interactions of Arabidopsis thaliana HDACs and MADS-domain proteins using an evolutionary perspective combined with bioinformatics analyses and testing the more promising predicted interactions through classic molecular biology tools. Through bioinformatic analyses, we found similarities between HDACs proteins from different organisms, which allowed us to predict a putative protein-protein interaction between the Arabidopsis thaliana deacetylase HDA15 and the MADS-domain protein XAANTAL1 (XAL1). The results of two-hybrid and Bimolecular Fluorescence Complementation analysis demonstrated in vitro and in vivo HDA15-XAL1 interaction in the nucleus. Likely, this interaction might regulate developmental processes in plants as is the case for this type of interaction in animals.

A survey of k-mer methods and applications in bioinformatics

The rapid progression of genomics and proteomics has been driven by the advent of advanced sequencing technologies, large, diverse, and readily available omics datasets, and the evolution of computational data processing capabilities. The vast amount of data generated by these advancements necessitates efficient algorithms to extract meaningful information. K-mers serve as a valuable tool when working with large sequencing datasets, offering several advantages in computational speed and memory efficiency and carrying the potential for intrinsic biological functionality. This review provides an overview of the methods, applications, and significance of k-mers in genomic and proteomic data analyses, as well as the utility of absent sequences, including nullomers and nullpeptides, in disease detection, vaccine development, therapeutics, and forensic science. Therefore, the review highlights the pivotal role of k-mers in addressing current genomic and proteomic problems and underscores their potential for future breakthroughs in research.

Whole genome sequencing unravels cryptic circulation of divergent dengue virus lineages in the rainforest region of Nigeria

Dengue is often misclassified and underreported in Africa due to inaccurate differential diagnoses of nonspecific febrile illnesses such as malaria, sparsity of diagnostic testing and poor clinical and genomic surveillance. There are limited reports on the seroprevalence and genetic diversity of dengue virus (DENV) in humans and vectors in Nigeria. In this study, we investigated the epidemiology and genetic diversity of dengue in the rainforest region of Nigeria. We screened 515 febrile patients who tested negative for malaria and typhoid fever in three hospitals in Oyo and Ekiti States in southern Nigeria with a combination of anti-dengue IgG/IgM/NS1 rapid test kits and metagenomic sequencing. We found that approximately 28% of screened patients had previous DENV exposure, with the highest prevalence in persons over sixty. Approximately 8% of the patients showed evidence of recent or current infection, and 2.7% had acute infection. Following sequencing of sixty samples, we assembled twenty DENV-1 genomes (3 complete and 17 partial). We found that all assembled genomes belonged to DENV-1 genotype III. Our phylogenetic analyses showed evidence of prolonged cryptic circulation of divergent DENV lineages in Oyo state. We were unable to resolve the source of DENV in Nigeria owing to limited sequencing data from the region. However, our sequences clustered closely with sequences in Tanzania and sequences reported in Chinese with travel history to Tanzania in 2019. This may reflect the wider unsampled bidirectional transmission of DENV-1 in Africa, which strongly emphasizes the importance of genomic surveillance in monitoring ongoing DENV transmission in Africa.

Molecular characterization of the meq oncogene of Marek's disease virus in vaccinated Brazilian poultry farms reveals selective pressure on prevalent strains

Marek's disease virus (MDV) has become an increasingly virulent pathogen in the poultry industry despite vaccination efforts to control it. Brazil has experienced a significant rise of Marek's disease (MD) outbreaks in recent years. Our study aimed to analyze the complete meq gene sequences to understand the molecular epidemiological basis of MD outbreaks in Brazilian vaccinated layer farms. We detected a high incidence rate of visceral MD (67.74%) and multiple circulating MDV strains. The most prevalent and geographically widespread genotype presented several clinical and molecular characteristics of a highly virulent strain and evolving under positive selective pressure. Phylogenetic and phylogeographic analysis revealed a closer relationship with strains from the USA and Japan. This study sheds light on the circulation of MDV strains capable of infecting vaccinated birds. We emphasize the urgency of adopting preventive measures to manage MDV outbreaks threatening the poultry farming industry.

Viral diversity in wild and urban rodents of Yunnan Province, China

Rodents represent over 40% of known mammal species and are found in various terrestrial habitats. They are significant reservoirs for zoonotic viruses, including harmful pathogens such as arenaviruses and hantaviruses, yet knowledge of their hosts and distributions is limited. Therefore, characterizing the virome profile in these animals is invaluable for outbreak preparedness, especially in potential hotspots of mammal diversity. This study included 681 organs from 124 rodents and one Chinese tree shrew collected from Yunnan Province, China, during 2020-2021. Metagenomic analysis revealed unique features of mammalian viruses in rodent organs across habitats with varying human disturbances. R. tanezumi in locations with high anthropogenic disturbance exhibited the highest mammal viral diversity, with spleen and lung samples showing the highest diversities for these viruses at the organ level. Mammal viral diversity for both commensal and non-commensal rats was identified to positively correlate with landscape disturbance. Some virus families were associated with particular organs or host species, suggesting tropism for these pathogens. Notably, known and novel viral species that are likely to infect humans were identified. R. tanezumi was identified as a reservoir and carrier for various zoonotic viruses, including porcine bocavirus, hantavirus, cardiovirus, and lyssavirus. These findings highlight the influence of rodent community composition and anthropogenic activities on diverse virome profiles, with R. tanezumi as an important reservoir for zoonotic viruses.

Cloning, expression, and purification of an α-carbonic anhydrase from Toxoplasma gondii to unveil its kinetic parameters and anion inhibition profile

Toxoplasmosis, induced by the intracellular parasite Toxoplasma gondii, holds considerable implications for global health. While treatment options primarily focusing on folate pathway enzymes have notable limitations, current research endeavours concentrate on pinpointing specific metabolic pathways vital for parasite survival. Carbonic anhydrases (CAs, EC 4.2.1.1) have emerged as potential drug targets due to their role in fundamental reactions critical for various protozoan metabolic processes. Within T. gondii, the Carbonic Anhydrase-Related Protein (TgCA_RP) plays a pivotal role in rhoptry biogenesis. Notably, α-CA (TcCA) from another protozoan, Trypanosoma cruzi, exhibited considerable susceptibility to classical CA inhibitors (CAIs) such as anions, sulphonamides, thiols, and hydroxamates. Here, the recombinant DNA technology was employed to synthesise and clone the identified gene in the T. gondii genome, which encodes an α-CA protein (Tg_CA), with the purpose of heterologously overexpressing its corresponding protein. Tg_CA kinetic constants were determined, and its inhibition patterns explored with inorganic metal-complexing compounds, which are relevant for rational compound design. The significance of this study lies in the potential development of innovative therapeutic strategies that disrupt the vital metabolic pathways crucial for T. gondii survival and virulence. This research may lead to the development of targeted treatments, offering new approaches to manage toxoplasmosis.

Visualizing metagenomic and metatranscriptomic data: A comprehensive review

The fields of Metagenomics and Metatranscriptomics involve the examination of complete nucleotide sequences, gene identification, and analysis of potential biological functions within diverse organisms or environmental samples. Despite the vast opportunities for discovery in metagenomics, the sheer volume and complexity of sequence data often present challenges in processing analysis and visualization. This article highlights the critical role of advanced visualization tools in enabling effective exploration, querying, and analysis of these complex datasets. Emphasizing the importance of accessibility, the article categorizes various visualizers based on their intended applications and highlights their utility in empowering bioinformaticians and non-bioinformaticians to interpret and derive insights from meta-omics data effectively.

Parasitic findings on threatened pudu deer from Central Chile accounts first genetic characterization of lice parasitizing P. puda in Chile and the first molecular report of Taenia hydatigena metacestodes in this species

Southern pudu (Pudu puda) is a threatened endemic deer of the temperate forests of Chile. In recent years pudu populations rates have decreased mainly due to anthropogenic causes including forest loss and landscape fragmentation. In this context, the parasitic fauna of Chilean pudu has been scarcely investigated. The aim of this study was to determine the parasitic status of rescued pudu n = 13 from its natural habitat in Central Chile (Maule region) during March 2022 and June 2023 by applying morphological, histopathological, and molecular analyses. As result, we report the presence of transmission of parasites from dogs to pudus as showed by the presence of metacestodes of the parasite Taenia hydatigena on omentum, liver, and pleura of pudus during postmortem examinations, being the first molecular report on the presence of this parasite on Chilean pudu. Meanwhile, ectoparasite examinations determined the presence of chewing and sucking lice on pudu exemplars here analysed. Molecular and phylogenetic analysis of lice revealed new insights on Bovicola and Anoplura lice parasitizing P. puda in Chile, equally being the first genetic characterization of lice parasitizing pudu exemplars in Chile. In addition, parasite loads of lice and metacestodes were analysed. However, no statistically significance was observed when comparing environmental and individual traits influence on parasite load variation. Overall, the study area is the northern limit of habitat distribution of this specie in Chile and we here provide novel information on pudu deer parasites, thus making a useful and valuable contribution to the parasitological knowledge on this threatened species.

Tracing histoplasmosis genomic epidemiology and species occurrence across the USA

ABSTRACTHistoplasmosis is an endemic mycosis in North America frequently reported along the Ohio and Mississippi River Valleys, although autochthonous cases occur in non-endemic areas. In the United States, the disease is provoked by two genetically distinct clades of Histoplasma capsulatum sensu lato, Histoplasma mississippiense (Nam1) and H. ohiense (Nam2). To bridge the molecular epidemiological gap, we genotyped 93 Histoplasma isolates (62 novel genomes) including clinical, environmental, and veterinarian samples from a broader geographical range by whole-genome sequencing, followed by evolutionary and species niche modelling analyses. We show that histoplasmosis is caused by two major lineages, H. ohiense and H. mississippiense; with sporadic cases caused by H. suramericanum in California and Texas. While H. ohiense is prevalent in eastern states, H. mississipiense was found to be prevalent in the central and western portions of the United States, but also geographically overlapping in some areas suggesting that these species might co-occur. Species Niche Modelling revealed that H. ohiense thrives in places with warmer and drier conditions, while H. mississippiense is endemic to areas with cooler temperatures and more precipitation. In addition, we predicted multiple areas of secondary contact zones where the two species co-occur, potentially facilitating gene exchange and hybridization. This study provides the most comprehensive understanding of the genomic epidemiology of histoplasmosis in the USA and lays a blueprint for the study of invasive fungal diseases.

A lineage-specific protein network at the trypanosome nuclear envelope

The nuclear envelope (NE) separates translation and transcription and is the location of multiple functions, including chromatin organization and nucleocytoplasmic transport. The molecular basis for many of these functions have diverged between eukaryotic lineages. Trypanosoma brucei, a member of the early branching eukaryotic lineage Discoba, highlights many of these, including a distinct lamina and kinetochore composition. Here, we describe a cohort of proteins interacting with both the lamina and NPC, which we term lamina-associated proteins (LAPs). LAPs represent a diverse group of proteins, including two candidate NPC-anchoring pore membrane proteins (POMs) with architecture conserved with S. cerevisiae and H. sapiens, and additional peripheral components of the NPC. While many of the LAPs are Kinetoplastid specific, we also identified broadly conserved proteins, indicating an amalgam of divergence and conservation within the trypanosome NE proteome, highlighting the diversity of nuclear biology across the eukaryotes, increasing our understanding of eukaryotic and NPC evolution.

The complete catalog of antimicrobial resistance secondary active transporters in Clostridioides difficile: evolution and drug resistance perspective

Secondary active transporters shuttle substrates across eukaryotic and prokaryotic membranes, utilizing different electrochemical gradients. They are recognized as one of the antimicrobial efflux pumps among pathogens. While primary active transporters within the genome of C. difficile 630 have been completely cataloged, the systematical study of secondary active transporters remains incomplete. Here, we not only identify secondary active transporters but also disclose their evolution and role in drug resistance in C. difficile 630. Our analysis reveals that C. difficile 630 carries 147 secondary active transporters belonging to 27 (super)families. Notably, 50 (34%) of them potentially contribute to antimicrobial resistance (AMR). AMR-secondary active transporters are structurally classified into five (super)families: the p-aminobenzoyl-glutamate transporter (AbgT), drug/metabolite transporter (DMT) superfamily, major facilitator (MFS) superfamily, multidrug and toxic compound extrusion (MATE) family, and resistance-nodulation-division (RND) family. Surprisingly, complete RND genes found in C. difficile 630 are likely an evolutionary leftover from the common ancestor with the diderm. Through protein structure comparisons, we have potentially identified six novel AMR-secondary active transporters from DMT, MATE, and MFS (super)families. Pangenome analysis revealed that half of the AMR-secondary transporters are accessory genes, which indicates an important role in adaptive AMR function rather than innate physiological homeostasis. Gene expression profile firmly supports their ability to respond to a wide spectrum of antibiotics. Our findings highlight the evolution of AMR-secondary active transporters and their integral role in antibiotic responses. This marks AMR-secondary active transporters as interesting therapeutic targets to synergize with other antibiotic activity.

Detection and phylogenetic analysis of contemporary H14N2 Avian influenza A virus in domestic ducks in Southeast Asia (Cambodia)

Avian influenza virus (AIV) in Asia is a complex system with numerous subtypes and a highly porous wild birds-poultry interface. Certain AIV subtypes, such as H14, are underrepresented in current surveillance efforts, leaving gaps in our understanding of their ecology and evolution. The detection of rare subtype H14 in domestic ducks in Southeast Asia comprises a geographic region and domestic bird population previously unassociated with this subtype. These H14 viruses have a complex evolutionary history involving gene reassortment events. They share sequence similarity to AIVs endemic in Cambodian ducks, and Eurasian low pathogenicity and high pathogenicity H5Nx AIVs. The detection of these H14 viruses in Southeast Asian domestic poultry further advances our knowledge of the ecology and evolution of this subtype and reinforces the need for continued, longitudinal, active surveillance in domestic and wild birds. Additionally, in vivo and in vitro risk assessment should encompass rare AIV subtypes, as they have the potential to establish in poultry systems.

Metabarcoding using nanopore sequencing enables identification of diverse and zoonotic vector-borne pathogens from neglected regions: A case study investigating dogs from Bhutan

The diversity and prevalence of canine vector-borne pathogens (VBPs) in Bhutan have to date remained unexplored, whilst recent epidemiological surveys in other South Asian nations have found diseases caused by VBPs to be rife in local dog populations. Importantly, many of such VBPs can infect people as well, with a building body of evidence identifying potentially zoonotic rickettsial organisms infecting humans in Bhutan. Given the lack of data on canine pathogens in Bhutan we employed a suite of deep-sequencing metabarcoding methods using Oxford Nanopore Technologies' MinION™ device to holistically characterise the bacterial, apicomplexan and filarial worm blood-borne pathogens of dogs in the country's south. Of the 95 stray, owned and community dogs sampled 78% (95% CI = 69%-85%) were infected with at least one VBP. Pathogen species detected were highly diverse including the bacteria Mycoplasma haemocanis in 16% (95% CI: 10-24%), Ehrlichia canis in 4% (95% CI: 2-10%), Anaplasma platys in 2% (95% CI: 0.5-7%) of dogs as well as the zoonotic species Bartonella clarridgeiae in 1% (95% CI: 0.1-6%), a potentially novel Bartonella spp. and an Ehrlichia chaffeensis-like bacterium, both in 1% (95% CI: 0.1-6%) of dogs. The apicomplexan haemoparasites Hepatozoon canis in 62% (95% CI: 52-71%), Babesia gibsoni in 45% (95% CI: 36-55%) and Babesia vogeli in 3% (95% CI: 1-9%) of dogs were also detected. Finally, 5% (95% CI: 2-12%) of dogs were found to be infected with the filarioid Acanthocheilonema reconditum and 1% (95% CI: 0.1-6%) with zoonotic Dirofilaria sp. hongkongensis. One canine was found positive to the filarioid Setaria tundra, a species normally found infecting cervids. The elucidated diversity of VBP communities highlights the strength of assumption-free diagnostics, such as metabarcoding, in detecting rare, novel, and unexpected pathogens. This approach to identifying pathogen diversity is of critical importance when investigating regions and populations that have thus far been neglected, with the findings aiding the development of future One Health informed strategies for disease control.

Genetic characteristics and potential pathogenic agents in Campylobacter upsaliensis based on genomic analysis

Campylobacter upsaliensis was the most common Campylobacter species in pets' gastrointestinal tracts and has been isolated from patients with bacteremia, hemolytic-uremic syndrome, spontaneous abortion, and Guillain-Barré syndrome. However, the genetic characteristics and the full extent of its significance as a human pathogen remain to be fully understood. This study involved an investigation for genomic analysis of 154 strains from different sources and additional antimicrobial resistance profiles of 26 strains for this species. The genomes contained 1,558-1,971 CDS and the genome sizes were estimated to vary from 1.53 Mb to 1.86 Mb, with an average GC content of  34.71%. The entire analyzed genomes could be divided into three clades (A, B, and C) based on ANI and phylogenomic analysis. Significantly, nearly all strains in Clade B were isolated from patient samples, and the virulence-related sequences FlgD, GmhA, and CdtC might serve as determining factors for the classification of Clade B. Half of the tested isolates had MIC values over 64 μg mL-1 for nalidixic acid, gentamicin, and streptomycin. Isolates from pets in China carried more resistant elements in the genomes. This study both provided a comprehensive profile of C. upsaliensis for its genomic features and suggested some pathogenic agents for human infection with this species.

Spatio-temporal spread of Lassa virus and a new rodent host in the Mano River Union area, West Africa

The spread of Lassa virus (LASV) in Guinea, Liberia and Sierra Leone, which together are named the Mano River Union (MRU) area, was examined phylogeographically. To provide a reliable evolutionary scenario, new rodent-derived, whole LASV sequences were included. These were generated by metatranscriptomic next-generation sequencing from rodents sampled between 2003 and 2020 in 21 localities of Guinea and Sierra Leone. An analysis was performed using BEAST to perform continuous phylogeographic inference and EvoLaps v36 to visualize spatio-temporal spread. LASV was identified as expected in its primary host reservoir, the Natal multimammate mouse (Mastomys natalensis), and also in two Guinean multimammate mice (Mastomys erythroleucus) in northern Sierra Leone and two rusty-bellied brush-furred mice (Lophuromys sikapusi) in southern Sierra Leone. This finding is consistent with the latter two species being secondary host reservoirs. The strains in these three species were very closely related in LASV lineage IV. Phylogenetic analysis indicated that the most recent common ancestor of lineage IV existed 316-374 years ago and revealed distinct, well-supported clades from Sierra Leone (Bo, Kabala and Kenema), Guinea (Faranah, Kissidougou-Guekedou and Macenta) and Liberia (Phebe-Ganta). The phylogeographic scenario suggests southern Guinea as the point of origin of LASV in the MRU area, with subsequent spread to towards Mali, Liberia and Sierra Leone at a mean speed of 1.6 to 1.1 km/year.

Genome-wide analysis reveals the contributors to fast molecular evolution of the Chinese hook snout carp (Opsariichthys bidens)

Variations in molecular evolutionary rate have been widely investigated among lineages and genes. However, it remains an open question whether fast rate of molecular evolution is driven by natural selection or random drift, and how the fast rate is linked to metabolic rate. Additionally, previous studies on fast molecular evolution have been largely restricted to concatenated matrix of genes or a few specifically selected genes, but less is known for individual genes at the genome-wide level. Here we addressed these questions using more than 5000 single-copy orthologous (SCO) genes through comparative genomic and phylogenetic analyses among fishes, with a special focus on a newly-sequenced clupeocephalan fish the Chinese hook snout carp Opsariichthys bidens. We showed O. bidens displays significantly higher mean substitution rate and more fast-evolving SCO genes (2172 genes) than most fishes studied here. The rapidly evolving genes are enriched in highly conserved and very basic functions such as translation and ribosome that are critical for biological fitness. We further revealed that ∼25 % of these fast-evolving genes exhibit a constant increase of substitution rate from the common ancestor down to the present, suggesting a neglected but important contribution from ancestral states. Model fitting showed that ∼85 % of fast-evolving genes exclusive to O. bidens and related species follow the adaptive evolutionary model rather than random-drift model, and 7.6 % of fast-evolving genes identified in O. bidens have experienced positive selection, both indicating the reflection of adaptive selection. Finally, metabolic rate was observed to be linked with substitution rate in a gene-specific manner. Overall, our findings reveal fast molecular evolution of SCO genes at genome-wide level in O. bidens, and uncover the evolutionary and ecological contributors to it.

A comprehensive investigation of the anion inhibition profile of a β-carbonic anhydrase from Acinetobacter baumannii for crafting innovative antimicrobial treatments

This study refers to the intricate world of Acinetobacter baumannii, a resilient pathogenic bacterium notorious for its propensity at antibiotic resistance in nosocomial infections. Expanding upon previous findings that emphasised the bifunctional enzyme PaaY, revealing unexpected γ-carbonic anhydrase (CA) activity, our research focuses on a different class of CA identified within the A. baumannii genome, the β-CA, designated as 𝛽-AbauCA (also indicated as CanB), which plays a crucial role in the resistance mechanism mediated by AmpC beta-lactamase. Here, we cloned, expressed, and purified the recombinant 𝛽-AbauCA, unveiling its distinctive kinetic properties and inhibition profile with inorganic anions (classical CA inhibitors). The exploration of 𝛽-AbauCA not only enhances our understanding of the CA repertoire of A. baumannii but also establishes a foundation for targeted therapeutic interventions against this resilient pathogen, promising advancements in combating its adaptability and antibiotic resistance.

Characteristics and phylogenetic distribution of megaplasmids and prediction of a putative chromid in Pseudomonas aeruginosa

Research on megaplasmids that contribute to the spread of antimicrobial resistance (AMR) in Pseudomonas aeruginosa strains has grown in recent years due to the now widely used technologies allowing long-read sequencing. Here, we systematically analyzed distinct and consistent genetic characteristics of megaplasmids found in P. aeruginosa. Our data provide information on their phylogenetic distribution and hypotheses tracing the potential evolutionary paths of megaplasmids. Most of the megaplasmids we found belong to the IncP-2-type, with conserved and syntenic genetic backbones carrying modules of genes associated with chemotaxis apparatus, tellurite resistance and plasmid replication, segregation, and transmission. Extensively variable regions harbor abundant AMR genes, especially those encoding β-lactamases such as VIM-2, IMP-45, and KPC variants, which are high-risk elements in nosocomial infection. IncP-2 megaplasmids act as effective vehicles transmitting AMR genes to diverse regions. One evolutionary model of the origin of megaplasmids claims that chromids can develop from megaplasmids. These chromids have been characterized as an intermediate between a megaplasmid and a chromosome, also containing core genes that can be found on the chromosome but not on the megaplasmid. Using in silico prediction, we identified the "PABCH45 unnamed replicon" as a putative chromid in P. aeruginosa, which shows a much higher similarity and closer phylogenetic relationship to chromosomes than to megaplasmids while also encoding plasmid-like partition genes. We propose that such a chromid could facilitate genome expansion, allowing for more rapid adaptations to novel ecological niches or selective conditions, in comparison to megaplasmids.

Genomic characterization of highly pathogenic avian influenza A H5N1 virus newly emerged in dairy cattle

In March 2024, the emergence of highly pathogenic avian influenza (HPAI) A (H5N1) infections in dairy cattle was detected in the United Sates for the first time. We genetically characterize HPAI viruses from dairy cattle showing an abrupt drop in milk production, as well as from two cats, six wild birds, and one skunk. They share nearly identical genome sequences, forming a new genotype B3.13 within the 2.3.4.4b clade. B3.13 viruses underwent two reassortment events since 2023 and exhibit critical mutations in HA, M1, and NS genes but lack critical mutations in PB2 and PB1 genes, which enhance virulence or adaptation to mammals. The PB2 E627 K mutation in a human case associated with cattle underscores the potential for rapid evolution post infection, highlighting the need for continued surveillance to monitor public health threats.

Low depth sequencing reveals the critically endangered Batagur kachuga (Red-crowned roofed turtle) mitochondrial genome and its evolutionary implications

The Batagur kachuga (B. kachuga), commonly known as the Red-crowned roofed turtle, is a critically endangered species native to India and its neighboring countries like Bangladesh, and Nepal. The present study is the first report of the complete mitochondrial genome of B. kachuga (16,517 bp) construed via the next-generation sequencing (NGS) approach from eggshell DNA. There are 22 transfer RNAs (tRNAs), 2 ribosomal RNAs (rRNAs), 13 protein-coding genes (PCGs), and one putative control region (CR/D-loop) in the mitogenome. The CR region from the current study reveals conserved TAS, CD, and CSB domains and two AT-rich tandem repeat regions. Most genes are encoded in the heavy strand except the NADH dehydrogenase subunit 6 (ND6) gene and seven tRNA genes. Most PCGs start with the initiation codon ATG, except the COI (Cytochrome Oxidase Subunit-I) gene, which starts with the GTG codon. The present investigation also predicts the distinctive cloverleaf structures of tRNAs except for tRNA-Ser1 and tRNA-Ser-2, which lack a DHU arm. The comparative analysis of Ka/Ks with other 33 species from Order Testudines, in relation to B. kachuga, revealed negative selection in most PCGs, indicating a process of preservation and purification that aids in eliminating undesirable or detrimental substitutes. Phylogenetic analysis of this species has been analysed using the complete mitogenome of 33 turtle species. The maximum likelihood phylogenetic tree strongly supports each family in different clades and also reveals a close relationship between the Pangashura and Batagur genera. Our study suggests the generation of genome-wide molecular data, in terms of mitogenomes, SNPs, and SSRs, is needed to improve the understanding of this species and their phylogenetics and evolutionary relationships, which will help to improve the conservation efforts of this species.

Species-specific variation in mitochondrial genome tandem repeat polymorphisms in hares (Lepus spp., Lagomorpha, Leporidae) provides insight into their evolution

The non-coding regions of the mitochondrial DNAs (mtDNAs) of hares, rabbits, and pikas (Lagomorpha) contain short (∼20 bp) and long (130-160 bp) tandem repeats, absent in related mammalian orders. In the presented study, we provide in-depth analysis for mountain hare (Lepus timidus) and brown hare (L. europaeus) mtDNA non-coding regions, together with a species- and population-level analysis of tandem repeat variation. Mountain hare short tandem repeats (SRs) as well as other analyzed hare species consist of two conserved 10 bp motifs, with only brown hares exhibiting a single, more variable motif. Long tandem repeats (LRs) also differ in sequence and copy number between species. Mountain hares have four to seven LRs, median value five, while brown hares exhibit five to nine LRs, median value six. Interestingly, introgressed mountain hare mtDNA in brown hares obtained an intermediate LR length distribution, with median copy number being the same as with conspecific brown hare mtDNA. In contrast, transfer of brown hare mtDNA into cultured mtDNA-less mountain hare cells maintained the original LR number, whereas the reciprocal transfer caused copy number instability, suggesting that cellular environment rather than the nuclear genomic background plays a role in the LR maintenance. Due to their dynamic nature and separation from other known conserved sequence elements on the non-coding region of hare mitochondrial genomes, the tandem repeat elements likely to represent signatures of ancient genetic rearrangements. clarifying the nature and dynamics of these rearrangements may shed light on the possible role of NCR repeated elements in mitochondria and in species evolution.

Mitogenome selection shaped the terrestrial adaptation of Grapsidae (Decapoda: Brachyura)

The colonization of aquatic to terrestrial habitats by brachyuran crabs requires genetic innovations as well as morphological adaptations to adapt to terrestrial environments. The genetic basis of such adaptive evolution, however, is largely unknown. This study focuses on terrestrialization in Geograpsus (Grapsidae) the only highly terrestrial genus in this family, which represents a notable example of terrestrial adaptive radiation. Here, we sequenced the mitogenomes of two Geograpsus species and used the mitogenomes of 215 representative crabs to construct phylogenetic and time frameworks that we used to infer terrestrial origins and evolution. Using mitochondrial genomic data, we demonstrated that marine crab ancestors began to settle on land during the early Eocene. Ocean acidification, the Paleocene-Eocene Thermal Maximum (PETM), and mangrove expansion at that time may have driven the diversification and ecological expansion of these terrestrial crabs. Evolutionary analyses reveal strong positive selection signals on monophyletic lineages of Grapsidae, especially the terrestrial species of Geograpsus. Positively selected sites in functionally important regions of ND5 and ND4 may imply enhanced energy metabolism in Grapsidae compared to other crabs, and may have played an important role in their terrestrial adaptation. Overall, our work provides valuable resources and opportunities to reveal the adaptation of crabs to complex terrestrial environments.

Geographic conditions impact the relationship between plant phenology and phylogeny

Plant phenology is influenced by a combined effect of phylogeny and climate, although it is yet unclear how these two variables work together to change phenology. We synthesized 107 previously published studies to examine whether phenological changes were impacted by both phylogeny and climate changes in various geographical settings globally. Phenological observation data from 52,463 plant species at 71 sites worldwide revealed that 90 % of phenological records showed phylogenetic conservation. i.e., closely related species exhibited similar phenology. To explore the significant and non-significant phylogenetic conservation between plant phenophases, our dataset comprises 5,47,000 observation records from the four main phenophases (leaf bud, leaf, flower, and fruit). Three-dimensional geographical distribution (altitude, latitude, and longitude) data analysis revealed that plant phenology may exhibit phylogenetic signals at finer special scales (optimal environmental conditions) that vanish in high altitude and latitude regions. Additionally, climatic sensitivity analysis suggested that phylogenetic signals were associated with plant phenophases and were stronger in the regions of ideal temperature (7-18 °C) and photoperiod (10-14 h) and weaker in harsh climatic conditions. These results show that phylogenetic conservation in plant phenological traits is frequently influenced by the interaction of harsh climatic conditions and geographical ranges. This meta-analysis enhances our knowledge of predicting species responses over geographic gradients under varied climatic conditions.

A novel synthesized Vanillin-Based Deep Eutectic Agent (V-DEA) mitigates postharvest fungal decay and improve shelf life and quality of cherry tomatoes

Fusarium oxysporum and Botrytis cinerea are the main pathogens that cause fruit decay and reduce the postharvest shelf life of cherry tomatoes. Boosting the potency of natural products requires implementing structural modification to combat postharvest pathogens. Herein, we developed a novel Vanillin-Deep Eutectic Agent (V-DEA) from natural compounds and evaluated its effectiveness against tomato fruit rot pathogens. The results demonstrated that V-DEA suppressed mycelium growth and spore germination of F. oxysporum and B. cinerea by enhancing cell membrane permeability, increasing lipid peroxidation, and inhibiting enzyme activities. Importantly, using 8-mM V-DEA successfully prevented postharvest decay in cherry tomatoes, while 4-mM significantly extended their shelf life by reducing weight loss and shriveling, and enhancing key fruit qualities such as total soluble solids, ascorbic acid, tartaric acid, and lycopene. In conclusion, V-DEA exhibits dual properties as a potent pathogen inhibitor and antioxidant activity, thus prolonging the shelf life of cherry tomatoes.

Sea anemone (Anthozoa, Actiniaria) diversity in Mo'orea (French Polynesia)

Sea anemones (Order Actiniaria) are a diverse group of marine invertebrates ubiquitous across marine ecosystems. Despite their wide distribution and success, a knowledge gap persists in our understanding of their diversity within tropical systems, owed to sampling bias of larger and more charismatic species overshadowing cryptic lineages. This study aims to delineate the sea anemone diversity in Mo'orea (French Polynesia) with the use of a dataset from the Mo'orea Biocode's "BioBlitz" initiative, which prioritized the sampling of more cryptic and understudied taxa. Implementing a target enrichment approach, we integrate 71 newly sequenced samples into an expansive phylogenetic framework and contextualize Mo'orea's diversity within global distribution patterns of sea anemones. Our analysis corroborates the presence of several previously documented sea anemones in French Polynesia and identifies for the first time the occurrence of members of genera Andvakia and Aiptasiomorpha. This research unveils the diverse sea anemone ecosystem in Mo'orea, spotlighting the area's ecological significance and emphasizing the need for continued exploration. Our methodology, encompassing a broad BLAST search coupled with phylogenetic analysis, proved to be a practical and effective approach for overcoming the limitations posed by the lack of comprehensive sequence data for sea anemones. We discuss the merits and limitations of current molecular methodologies and stress the importance of further research into lesser-studied marine organisms like sea anemones. Our work sets a precedent for future phylogenetic studies stemming from BioBlitz endeavors.

Utility of various genome lengths in diversity and evolution analyses of Hepatitis E virus

The aim of this study was to investigate to what extent fragments of the HEV genome could be used for accurate diagnostics and inference of viral population-scale processes. For this, we selected all the published whole genome sequences from the NCBI GenBank and trimmed them to various fragment lengths (ORF1,2,3, ORF1, ORF2, ORF3, 493 nt in ORF2 and 148 nt in ORF2). Each of the fragment lengths was used to infer the richness and diversity of the viral sequence types, typing accuracy, and potential use in phylodynamics. The results obtained from the different fragments were compared. We observed that, generally, the longer the nucleic acid fragment used in typing, the better the accuracy in predicting the viral subtype. However, the dominant HEV subtypes circulating in Europe were relatively well classified even by the 493 nt fragment, with false negative rates as low as 8 in 1000 typed sequences. Most fragments also give comparable results in analyses of population size, albeit with shorter fragments showing a broader 95 % highest posterior density interval and less obvious increase of the viral effective population size. The reconstructed phylogenies of a heterochronous subset indicated a good concordance between all the fragments, with the major clades following similar branching patterns. Furthermore, we have used the HEV sequence data from the Netherlands available in the HEVnet database as a case study for reconstruction of population size changes in the past decades. This data showed that molecular and epidemiological results are concordant and point to an increase in the viral effective population size underlying the observed increase in incidence of acute HEV infection cases. In the absence of whole genome sequencing data, the 493 bp fragment can be used for analyzing HEV strains currently circulating in Europe, as it is informative for describing short term population-scale processes.

Instructional Influencers: Teaching Professors as Potential Departmental Change Agents in Diversity, Equity, and Inclusion

At many research-intensive universities in North America, there is a disproportionate loss of minoritized undergraduate students from Science, Technology, Engineering, and Mathematics (STEM) majors. Efforts to confront this diversity, equity, and inclusion (DEI) challenge, such as faculty adoption of evidenced-based instructional approaches that promote student success, have been slow. Instructional and pedagogical change efforts at the academic department level have been demonstrated to be effective at enacting reform. One potential strategy is to embed change agent individuals within STEM departments that can drive change efforts. This study seeks to assess whether tenure-track, teaching-focused faculty housed in STEM departments are perceived as influential on the instructional and pedagogical domains of their colleagues. To answer this, individuals across five STEM departments at large, research-intensive campuses identified faculty who were influential upon six domains of their instruction and pedagogy. Social network analysis of individuals in these departments revealed heterogeneity across the instructional domains. Some, like the teaching strategies network, are highly connected and involve the majority of the department; while others, like the DEI influence network, comprise a significantly smaller population of faculty. Importantly, we demonstrate that tenure-track, teaching-focused faculty are influential across all domains of instruction, but are disproportionately so in the sparsely populated DEI influence networks.

Genomic SNPs resolve the phylogeny of an ancient amphibian island radiation from the Seychelles

Unusually for oceanic islands, the granitic Seychelles host multiple lineages of endemic amphibians. This includes an ancient (likely ca. 60 million years) radiation of eight caecilian species, most of which occur on multiple islands.These caecilians have a complicated taxonomic history and their phylogenetic inter-species relationships have been difficult to resolve. Double-digest RAD sequencing (ddRADseq) has been applied extensively to phylogeography and increasingly to phylogenetics but its utility for resolving ancient divergences is less well established. To address this, we applied ddRADseq to generate a genome-wide SNP panel for phylogenomic analyses of the Seychelles caecilians, whose phylogeny has so far not been satisfactorily resolved with traditional DNA markers. Based on 129,154 SNPs, we resolved deep and shallow splits, with strong support. Our findings demonstrate the capability of genome-wide SNPs for evolutionary inference at multiple taxonomic levels and support the recently proposed synonymy of Grandisonia Taylor, 1968 with Hypogeophis Peters, 1879. We revealed three clades of Hypogeophis (large-, medium- and short-bodied) and identify a single origin of the diminutive, stocky-bodied and pointy-snouted phenotype.

Museomics help resolving the phylogeny of snowfinches (Aves, Passeridae, Montifringilla and allies)

Historical specimens from museum collections provide a valuable source of material also from remote areas or regions of conflict that are not easily accessible to scientists today. With this study, we are providing a taxon-complete phylogeny of snowfinches using historical DNA from whole skins of an endemic species from Afghanistan, the Afghan snowfinch, Pyrgilauda theresae. To resolve the strong conflict between previous phylogenetic hypotheses, we generated novel mitogenome sequences for selected taxa and genome-wide SNP data using ddRAD sequencing for all extant snowfinch species endemic to the Qinghai-Tibet Plateau (QTP) and for an extended intraspecific sampling of the sole Central and Western Palearctic snowfinch species (Montifringilla nivalis). Our phylogenetic reconstructions unanimously refuted the previously suggested paraphyly of genus Pyrgilauda. Misplacement of one species-level taxon (Onychostruthus tazcanowskii) in previous snowfinch phylogenies was undoubtedly inferred from chimeric mitogenomes that included heterospecific sequence information. Furthermore, comparison of novel and previously generated sequence data showed that the presumed sister-group relationship between M. nivalis and the QTP endemic M. henrici was suggested based on flawed taxonomy. Our phylogenetic reconstructions based on genome-wide SNP data and on mitogenomes were largely congruent and supported reciprocal monophyly of genera Montifringilla and Pyrgilauda with monotypic Onychostruthus being sister to the latter. The Afghan endemic P. theresae likely originated from a rather ancient Pliocene out-of-Tibet dispersal probably from a common ancestor with P. ruficollis. Our extended trans-Palearctic sampling for the white-winged snowfinch, M. nivalis, confirmed strong lineage divergence between an Asian and a European clade dated to 1.5 - 2.7 million years ago (mya). Genome-wide SNP data suggested subtle divergence among European samples from the Alps and from the Cantabrian mountains.

Tracing the genealogy origin of geographic populations based on genomic variation and deep learning

Assigning a query individual animal or plant to its derived population is a prime task in diverse applications related to organismal genealogy. Such endeavors have conventionally relied on short DNA sequences under a phylogenetic framework. These methods naturally show constraints when the inferred population sources are ambiguously phylogenetically structured, a scenario demanding substantially more informative genetic signals. Recent advances in cost-effective production of whole-genome sequences and artificial intelligence have created an unprecedented opportunity to trace the population origin for essentially any given individual, as long as the genome reference data are comprehensive and standardized. Here, we developed a convolutional neural network method to identify population origins using genomic SNPs. Three empirical datasets (an Asian honeybee, a red fire ant, and a chicken datasets) and two simulated populations are used for the proof of concepts. The performance tests indicate that our method can accurately identify the genealogy origin of query individuals, with success rates ranging from  93 % to 100 %. We further showed that the accuracy of the model can be significantly increased by refining the informative sites through FST filtering. Our method is robust to configurations related to batch sizes and epochs, whereas model learning benefits from the setting of a proper preset learning rate. Moreover, we explained the importance score of key sites for algorithm interpretability and credibility, which has been largely ignored. We anticipate that by coupling genomics and deep learning, our method will see broad potential in conservation and management applications that involve natural resources, invasive pests and weeds, and illegal trades of wildlife products.

Several secondary metabolite gene clusters in the genomes of ten Penicillium spp. raise the risk of multiple mycotoxin occurrence in chestnuts

Penicillium spp. produce a great variety of secondary metabolites, including several mycotoxins, on food substrates. Chestnuts represent a favorable substrate for Penicillium spp. development. In this study, the genomes of ten Penicillium species, virulent on chestnuts, were sequenced and annotated: P. bialowiezense. P. pancosmium, P. manginii, P. discolor, P. crustosum, P. palitans, P. viridicatum, P. glandicola, P. taurinense and P. terrarumae. Assembly size ranges from 27.5 to 36.8 Mb and the number of encoded genes ranges from 9,867 to 12,520. The total number of predicted biosynthetic gene clusters (BGCs) in the ten species is 551. The most represented families of BGCs are non ribosomal peptide synthase (191) and polyketide synthase (175), followed by terpene synthases (87). Genome-wide collections of gene phylogenies (phylomes) were reconstructed for each of the newly sequenced Penicillium species allowing for the prediction of orthologous relationships among our species, as well as other 20 annotated Penicillium species available in the public domain. We investigated in silico the presence of BGCs for 10 secondary metabolites, including 5 mycotoxins, whose production was validated in vivo through chemical analyses. Among the clusters present in this set of species we found andrastin A and its related cluster atlantinone A, mycophenolic acid, patulin, penitrem A and the cluster responsible for the synthesis of roquefortine C/glandicoline A/glandicoline B/meleagrin. We confirmed the presence of these clusters in several of the Penicillium species conforming our dataset and verified their capacity to synthesize them in a chestnut-based medium with chemical analysis. Interestingly, we identified mycotoxin clusters in some species for the first time, such as the andrastin A cluster in P. flavigenum and P. taurinense, and the roquefortine C cluster in P. nalgiovense and P. taurinense. Chestnuts proved to be an optimal substrate for species of Penicillium with different mycotoxigenic potential, opening the door to risks related to the occurrence of multiple mycotoxins in the same food matrix.

Artificial intelligence enables unified analysis of historical and landscape influences on genetic diversity

While genetic variation in any species is potentially shaped by a range of processes, phylogeography and landscape genetics are largely concerned with inferring how environmental conditions and landscape features impact neutral intraspecific diversity. However, even as both disciplines have come to utilize SNP data over the last decades, analytical approaches have remained for the most part focused on either broad-scale inferences of historical processes (phylogeography) or on more localized inferences about environmental and/or landscape features (landscape genetics). Here we demonstrate that an artificial intelligence model-based analytical framework can consider both deeper historical factors and landscape-level processes in an integrated analysis. We implement this framework using data collected from two Brazilian anurans, the Brazilian sibilator frog (Leptodactylus troglodytes) and granular toad (Rhinella granulosa). Our results indicate that historical demographic processes shape most the genetic variation in the sibulator frog, while landscape processes primarily influence variation in the granular toad. The machine learning framework used here allows both historical and landscape processes to be considered equally, rather than requiring researchers to make an a priori decision about which factors are important.

Incomplete lineage sorting and hybridization underlie tree discordance in Petunia and related genera (Petunieae, Solanaceae)

Despite the overarching history of species divergence, phylogenetic studies often reveal distinct topologies across regions of the genome. The sources of these gene tree discordances are variable, but incomplete lineage sorting (ILS) and hybridization are among those with the most biological importance. Petunia serves as a classic system for studying hybridization in the wild. While field studies suggest that hybridization is frequent, the extent of reticulation within Petunia and its closely related genera has never been examined from a phylogenetic perspective. In this study, we used transcriptomic data from 11 Petunia, 16 Calibrachoa, and 10 Fabiana species to illuminate the relationships between these species and investigate whether hybridization played a significant role in the diversification of the clade. We inferred that gene tree discordance within genera is linked to hybridization events along with high levels of ILS due to their rapid diversification. Moreover, network analyses estimated deeper hybridization events between Petunia and Calibrachoa, genera that have different chromosome numbers. Although these genera cannot hybridize at the present time, ancestral hybridization could have played a role in their parallel radiations, as they share the same habitat and life history.

Horizontally transferred glycoside hydrolase 26 may aid hemipteran insects in plant tissue digestion

Glycoside hydrolases are enzymes that break down complex carbohydrates into simple sugars by catalyzing the hydrolysis of glycosidic bonds. There have been multiple instances of adaptive horizontal gene transfer of genes belonging to various glycoside hydrolase families from microbes to insects, as glycoside hydrolases can metabolize constituents of the carbohydrate-rich plant cell wall. In this study, we characterize the horizontal transfer of a gene from the glycoside hydrolase family 26 (GH26) from bacteria to insects of the order Hemiptera. Our phylogenies trace the horizontal gene transfer to the common ancestor of the superfamilies Pentatomoidea and Lygaeoidea, which include stink bugs and seed bugs. After horizontal transfer, the gene was assimilated into the insect genome as indicated by the gain of an intron, and a eukaryotic signal peptide. Subsequently, the gene has undergone independent losses and expansions in copy number in multiple lineages, suggesting an adaptive role of GH26s in some insects. Finally, we measured tissue-level gene expression of multiple stink bugs and the large milkweed bug using publicly available RNA-seq datasets. We found that the GH26 genes are highly expressed in tissues associated with plant digestion, especially in the principal salivary glands of the stink bugs. Our results are consistent with the hypothesis that this horizontally transferred GH26 was co-opted by the insect to aid in plant tissue digestion and that this HGT event was likely adaptive.

Phylogenomic analyses of ochrophytes (stramenopiles) with an emphasis on neglected lineages

Ochrophyta is a photosynthetic lineage that crowns the phylogenetic tree of stramenopiles, one of the major eukaryotic supergroups. Due to their ecological impact as a major primary producer, ochrophytes are relatively well-studied compared to the rest of the stramenopiles, yet their evolutionary relationships remain poorly understood. This is in part due to a number of missing lineages in large-scale multigene analyses, and an apparently rapid radiation leading to many short internodes between ochrophyte subgroups in the tree. These short internodes are also found across deep-branching lineages of stramenopiles with limited phylogenetic signal, leaving many relationships controversial overall. We have addressed this issue with other deep-branching stramenopiles recently, and now examine whether contentious relationships within the ochrophytes may be resolved with the help of filling in missing lineages in an updated phylogenomic dataset of ochrophytes, along with exploring various gene filtering criteria to identify the most phylogenetically informative genes. We generated ten new transcriptomes from various culture collections and a single-cell isolation from an environmental sample, added these to an existing phylogenomic dataset, and examined the effects of selecting genes with high phylogenetic signal or low phylogenetic noise. For some previously contentious relationships, we find a variety of analyses and gene filtering criteria consistently unite previously unstable groupings with strong statistical support. For example, we recovered a robust grouping of Eustigmatophyceae with Raphidophyceae-Phaeophyceae-Xanthophyceae while Olisthodiscophyceae formed a sister-lineage to Pinguiophyceae. Selecting genes with high phylogenetic signal or data quality recovered more stable topologies. Overall, we find that adding under-represented groups across different lineages is still crucial in resolving phylogenetic relationships, and discrete gene properties affect lineages of stramenopiles differently. This is something which may be explored to further our understanding of the molecular evolution of stramenopiles.

Congruent patterns of cryptic cladogenesis revealed using RADseq and Sanger sequencing in a velvet worm species complex (Onychophora: Peripatopsidae: Peripatopsis sedgwicki)

In the present study, first generation DNA sequencing (mitochondrial cytochrome c oxidase subunit one, COI) and reduced-representative genomic RADseq data were used to understand the patterns and processes of diversification of the velvet worm, Peripatopsis sedgwicki species complex across its distribution range in South Africa. For the RADseq data, three datasets (two primary and one supplementary) were generated corresponding to 1,259-11,468 SNPs, in order to assess the diversity and phylogeography of the species complex. Tree topologies for the two primary datasets were inferred using maximum likelihood and Bayesian inferences methods. Phylogenetic analyses using the COI datasets retrieved four distinct, well-supported clades within the species complex. Five species delimitation methods applied to the COI data (ASAP, bPTP, bGMYC, STACEY and iBPP) all showed support for the distinction of the Fort Fordyce Nature Reserve specimens. In the main P. sedgwicki species complex, the species delimitation methods revealed a variable number of operational taxonomic units and overestimated the number of putative taxa. Divergence time estimates coupled with the geographic exclusivity of species and phylogeographic results suggest recent cladogenesis during the Plio/Pleistocene. The RADseq data were subjected to a principal components analysis and a discriminant analysis of principal components, under a maximum-likelihood framework. The latter results corroborate the four main clades observed using the COI data, however, applying additional filtering revealed additional diversity. The high overall congruence observed between the RADseq data and COI data suggest that first generation sequence data remain a cheap and effective method for evolutionary studies, although RADseq does provide a far greater resolution of contemporary temporo-spatial patterns.

Phylogenomic analyses revealed widely occurring hybridization events across Elsholtzieae (Lamiaceae)

Obtaining a robust phylogeny proves challenging due to the intricate evolutionary history of species, where processes such as hybridization and incomplete lineage sorting can introduce conflicting signals, thereby complicating phylogenetic inference. In this study, we conducted comprehensive sampling of Elsholtzieae, with a particular focus on its largest genus, Elsholtzia. We utilized 503 nuclear loci and complete plastome sequences obtained from 99 whole-genome sequencing datasets to elucidate the interspecific relationships within the Elsholtzieae. Additionally, we explored various sources of conflicts between gene trees and species trees. Fully supported backbone phylogenies were recovered, and the monophyly of Elsholtzia and Keiskea was not supported. Significant gene tree heterogeneity was observed at numerous nodes, particularly regarding the placement of Vuhuangia and the E. densa clade. Further investigations into potential causes of this discordance revealed that incomplete lineage sorting (ILS), coupled with hybridization events, has given rise to substantial gene tree discordance. Several species, represented by multiple samples, exhibited a closer association with geographical distribution rather than following a strictly monophyletic pattern in plastid trees, suggesting chloroplast capture within Elsholtzieae and providing evidence of hybridization. In conclusion, this study provides phylogenomic insights to untangle taxonomic problems in the tribe Elsholtzieae, especially the genus Elsholtzia.

Optimizing reproductive performance in pangasius catfish broodstock: A review of dietary and molecular strategies

Pangasius catfish, a significant player in the global whitefish market, encounters challenges in aquaculture production sustainability. Quality broodstock maintenance and seed production are impeded by growth, maturation, and fecundity issues. This review investigates the efficacy of strategic nutrient composition and molecular strategies in enhancing broodstock conditions and reproductive performance across various fish species. A notable knowledge gap for Pangasius catfish hampers aquaculture progress. The review assesses nutrient manipulation's impact on reproductive physiology, emphasizing pangasius broodstock. A systematic review analysis following PRISMA guidelines was conducted to identify research trends and hotspots quantitatively, revealing a focus on P. bocourti and fertilization techniques. Addressing this gap, the review offers insights into dietary nutrients manipulation and genetic tool utilization for improved seed production, contributing to pangasius catfish aquaculture sustainability.

A novel genus of Pectobacterium bacteriophages display broad host range by targeting several species of Danish soft rot isolates

The bacterial diseases black leg and soft rot in potatoes cause heavy losses of potatoes worldwide. Bacteria within the genus Pectobacteriaceae are the causative agents of black leg and soft rot. The use of antibiotics in agriculture is heavily regulated and no other effective treatment currently exists, but bacteriophages (phages) have shown promise as potential biocontrol agents. In this study we isolated soft rot bacteria from potato tubers and plant tissue displaying soft rot or black leg symptoms collected in Danish fields. We then used the isolated bacterial strains as hosts for phage isolation. Using organic waste, we isolated phages targeting different species within Pectobacterium. Here we focus on seven of these phages representing a new genus primarily targeting P. brasiliense; phage Ymer, Amona, Sabo, Abuela, Koroua, Taid and Pappous. TEM image of phage Ymer showed siphovirus morphotype, and the proposed Ymer genus belongs to the class Caudoviricetes, with double-stranded DNA genomes varying from 39 kb to 43 kb. In silico host range prediction using a CRISPR-Cas spacer database suggested both P. brasiliense, P. polaris and P. versatile as natural hosts for phages within the proposed Ymer genus. A following host range experiment, using 47 bacterial isolates from Danish tubers and plants symptomatic with soft rot or black leg disease verified the in silico host range prediction, as the genus as a group were able to infect all three Pectobacterium species. Phages did, however, primarily target P. brasiliense isolates and displayed differences in host range even within the species level. Two of the phages were able to infect two or more Pectobacterium species. Despite no nucleotide similarity with any phages in the NCBI database, the proposed Ymer genus did share some similarity at the protein level, as well as gene synteny, with currently known phages. None of the phages encoded integrases or other genes typically associated with lysogeny. Similarly, no virulence factors nor antimicrobial resistance genes were found, and combined with their ability to infect several soft rot-causing Pectobacterium species from Danish fields, demonstrates their potential as biocontrol agents against soft rot and black leg diseases in potatoes.

Phylogenomics and biogeography of the small carpenter bees (Apidae: Xylocopinae: Ceratina)

Small carpenter bees in the genus Ceratina are behaviourally diverse, species-rich, and cosmopolitan, with over 370 species and a range including all continents except Antarctica. Here, we present the first comprehensive phylogeny of the genus based on ultraconserved element (UCE) phylogenomic data, covering a total of 185 ingroup specimens representing 22 of the 25 current subgenera. Our results support most recognized subgenera as natural groups, but we also highlight several groups in need of taxonomic revision - particularly the nominate subgenus Ceratina sensu stricto - and several clades that likely need to be described as new subgenera. In addition to phylogeny, we explore the evolutionary history of Ceratina through divergence time estimation and biogeographic reconstruction. Our findings suggest that Ceratinini split from its sister tribe Allodapini about 72 million years ago. The common ancestor of Ceratina emerged in the Afrotropical realm approximately 42 million years ago, near the Middle Eocene Climatic Optimum. Multiple subsequent dispersal events led to the present cosmopolitan distribution of Ceratina, with the majority of transitions occurring between the Afrotropics, Indomalaya, and the Palearctic. Additional movements also led to the arrival of Ceratina in Madagascar, Australasia, and a single colonization of the Americas. Dispersal events were asymmetrical overall, with temperate regions primarily acting as destinations for migrations from tropical source regions.

Molecular phylogeny of Hesperiidae (Lepidoptera) with an emphasis on Asian and African genera

Despite considerable research efforts in recent years, the deeper phylogenetic relationships among skipper butterflies (Hesperiidae) remain unresolved. This is primarily because of limited sampling, especially within Asian and African lineages. In this study, we consolidated previous data and extensively sampled Asian and African taxa to elucidate the phylogenetic relationships within Hesperiidae. The molecular dataset comprised sequences from two mitochondrial and two nuclear gene regions from 563 species that represented 353 genera. Our analyses revealed seven subfamilies within Hesperiidae: Coeliadinae, Euschemoninae, Eudaminae, Pyrginae, Heteropterinae, Trapezitinae, and Hesperiinae. The systematics of most tribes and genera aligned with those of prior studies. However, notable differences were observed in several tribes and genera. Overall, the position of taxa assigned to incertae sedis in Hesperiinae is largely clarified in this study. Our results strongly support the monophyly of the tribe Tagiadini (Pyrginae), and the systematics of some genera are clarified with comprehensive discussion. We recognize 15 tribes within the subfamily Hesperiinae. Of these, nine tribes are discussed in detail: Aeromachini, Astictopterini, Erionotini, Unkanini (new status), Ancistroidini, Ismini (confirmed status), Plastingini (new status), Gretnini (confirmed status), and Eetionini (confirmed status). We propose four subtribes within Astictopterini: Hypoleucina subtrib.n., Aclerosina, Cupithina, and Astictopterina. Furthermore, we describe a new genus (Hyarotoidesgen.n.) and reinstate two genera (Zeareinst.stat. and Separeinst.stat.) as valid. Additionally, we propose several new combinations: Zea mythecacomb.n.,Sepa bononiacomb.n. & reinst.stat., and Sepa umbrosacomb.n. This study, with extensive sampling of Asian and African taxa, greatly enhances the understanding of the knowledge of the skipper tree of life.

C2H2 zinc finger protein PagIDD15A regulates secondary wall thickening and lignin biosynthesis in poplar

Wood production is largely determined by the activity of cambial cell proliferation, and the secondary cell wall (SCW) thickening of xylem cells determines the wood property. In this study, we identified an INDETERMINATE DOMAIN (IDD) type C2H2 zinc finger transcription factor PagIDD15A as a regulator of wood formation in Populus alba × Populus glandulosa. Downregulation of PagIDD15A expression by RNA interference (RNAi) inhibited xylem development and xylem cell secondary wall thickening. RNA-seq analysis showed that PagPAL1, PagCCR2 and PagCCoAOMT1 were downregulated in the differentiating xylem of the PagIDD15A-RNAi transgenic plants, showing that PagIDD15A may regulate SCW biosynthesis through inhibiting lignin biosynthesis. The downregulation of PagVND6-B2, PagMYB10 and PagMYC4 and upregulation of PagWRKY12 in the differentiating xylem of RNAi transgenic plants suggest that PagIDD15A may also regulate these transcription factor (TF) genes to affect SCW thickening. RT-qPCR analysis in the phloem-cambium of RNAi transgenic demonstrates that PagIDD15A may regulate the expression of the genes associated with cell proliferation, including, PagSHR (SHORTROOT), PagSCR (SCARECROW), PagCYCD3;1 (CYCLIN D3;1) and PagSMR4 (SIAMESE-RELATED4), to affect the cambial activity. This study provides the knowledge of the IDD-type C2H2 zinc finger protein in regulating wood formation.

Testing ultraconserved elements (UCEs) for phylogenetic inference across bivalves (Mollusca: Bivalvia)

Bivalves constitute an important resource for fisheries and as cultural objects. Bivalve phylogenetics has had a long tradition using both morphological and molecular characters, and genomic resources are available for a good number of commercially important species. However, relationships among bivalve families have been unstable and major conflicting results exist between mitogenomics and results based on Sanger-based amplicon sequencing or phylotranscriptomics. Here we design and test an ultraconserved elements probe set for the class Bivalvia with the aim to use hundreds of loci without the need to sequence full genomes or transcriptomes, which are expensive and complex to analyze, and to open bivalve phylogenetics to museum specimens. Our probe set successfully captured 1,513 UCEs for a total of 263,800 bp with an average length of 174.59 ± 3.44 per UCE (ranging from 28 to 842 bp). Phylogenetic testing of this UCE probe set across Bivalvia and within the family Donacidae using different data matrices and methods for phylogenetic inference shows promising results at multiple taxonomic levels. In addition, our probe set was able to capture large numbers of UCEs for museum specimens collected before 1900 and from DNAs properly stored, of which many museums and laboratories are well stocked. Overall, this constitutes a novel and useful resource for bivalve phylogenetics.

Molecular evolution of Hokkaido virus, a genotype of Orthohantavirus puumalaense, among Myodes rodents

Viruses in the genus Orthohantavirus within the family Hantaviridae cause human hantavirus infections and represent a threat to public health. Hokkaido virus (HOKV), a genotype of Orthohantavirus puumalaense (Puumala virus; PUUV), was first identified in Tobetsu, Hokkaido, Japan. Although it is genetically related to the prototype of PUUV, the evolutionary pathway of HOKV is unclear. We conducted a field survey in a forest in Tobetsu in 2022 and captured 44 rodents. Complete coding genome sequences of HOKVs were obtained from five viral-RNA-positive rodents (four Myodes rufocanus bedfordiae and one Apodemus speciosus). Phylogenetic analysis revealed a close relationship between the phylogenies and geographical origins of M. rufocanus-related orthohantaviruses. Comparison of the phylogenetic trees of the S segments of orthohantaviruses and the cytochrome b genes of Myodes species suggested that Myodes-related orthohantaviruses evolved in Myodes rodent species as a result of genetic isolation and host switching.

Expression and functional analysis of Fushi Tarazu transcription factor 1 (FTZ-F1) in the regulation of steroid hormones during the gonad development of Fujian Oyster, Crassostrea angulata

Crassostrea angulata, a major shellfish cultivated in Southern China, has experienced a notable surge in commercial value in recent years. Understanding the molecular mechanisms governing their reproductive processes holds significant implications for advancing aquaculture practices. In this study, we cloned the orphan nuclear receptor gene, Fushi Tarazu transcription factor 1 (FTZ-F1), of C. angulata and investigated its functional role in the gonadal development. The full-length cDNA of FTZ-F1 spans 2357 bp and encodes a protein sequence of 530 amino acids. Notably, the amino acid sequence of FTZ-F1 in C. angulata shares remarkable similarity with its homologues in other species, particularly in the DNA-binding region (>90%) and ligand-binding region (>44%). In C. angulata, the highest expression level of FTZ-F1 was observed in the ovary, exhibiting more than a 200-fold increase during the maturation stage compared to the initiation stage (P < 0.001). Specifically, FTZ-F1 was mainly expressed in the follicular cells surrounding the oocytes of C. angulata. Upon inhibiting FTZ-F1 gene expression in C. angulata through RNA interference (RNAi), a substantial reduction in the expression of genes involved in the synthesis of sex steroids in the gonads, including 3β-HSD, Cyp17, and follistatin, was observed. In addition, estradiol (E2) and testosterone (T) levels also showed a decrease upon FTZ-F1 silencing, resulting in a delayed gonadal development. These results indicate that FTZ-F1 acts as a steroidogenic factor, participating in the synthesis and regulation of steroid hormones and thus playing an important role in the reproductive and endocrine systems within oysters.

Ice phenology interactions with water and air temperatures in high mountain lakes

Ice phenology is of great importance for the thermal structure of lakes and ponds and the biology of lake species. Under the current climate change conditions, ice-cover duration has been reduced by an advance in ice-off, and a delay in ice-on, and future projections foresee this trend as continuing. Here, we describe the current ice phenology of Pyrenean high mountain lakes and ponds, including ice-cover duration and ice-on and ice-off dates. We used mixed models to identify the variables that explained the observed patterns, extrapolated them across all water bodies in the mountain range, and related the seasonality of air and water temperatures with ice phenology using structural equation models. Ice phenology was obtained from the temperature series of 85 lakes and ponds for fourteen years, including 2001 to 2004 and 2009 to 2019. We discovered that high autumn precipitation was related to earlier ice-on dates, and that earlier ice-off dates were associated with higher following-summer water temperatures. We found a greater predictability of ice-off dates and ice-cover duration than ice-on dates. Altitude was the most important variable explaining the variation in ice phenology, followed by latitude, which was related to climatic differences among the northern and southern slopes of the mountain range. The lake area was significant for ice-on dates and ice-cover duration. The interannual variability in air temperature and radiation was remarkable for the ice-off date and ice-cover duration but not for the ice-on date. In contrast, wind speed was related to an earlier ice-off date and shorter ice-cover duration. All the measured lakes and ponds froze in winter during the studied period, a feature maintained in the extrapolation to the whole set of water bodies.

Integrative taxonomy, phylogenetics and historical biogeography of subgenus Aeschyntelus Stål, 1872 (Hemiptera: Heteroptera: Rhopalidae)

The subgenus Aeschyntelus includes six species that show variations in body color and shape, thus making it difficult to identify them based on morphological identification alone. To date, no genetic study has evaluated species within this genus. Herein, we collected 171 individuals from 90 localities of Rhopalus and employed an integrative taxonomic approach that incorporated morphological data, mitochondrial genomic data (COI, whole mitochondrial data) and nuclear genomic data (18S + 28S rRNAs, nuclear genome-wide SNPs) to delineate species boundaries. Our analyses confirmed the status of nine described species of Rhopalus and proposed the recognition of one new species known as Rhopalus qinlinganus sp. nov., which is classified within the subgenus Aeschyntelus. Discrepancies arising from nuclear and mitochondrial data suggest the presence of mito-nuclear discordance. Specifically, mitochondrial data indicated admixture within Clade A, comprising R. kerzhneri and R. latus, whereas genome-wide SNPs unambiguously identified two separate species, aligning with morphological classification. Conversely, mitochondrial data clearly distinguished Clade B- consisting of R. sapporensis into two lineages, whereas genome-wide SNPs unequivocally identified a single species. Our study also provides insights into the evolutionary history of Aeschyntelus, thus indicating that it likely originated in East Asia during the middle Miocene. The development of Aeschyntelus biodiversity in the southwestern mountains of China occurred via an uplift-driven diversification process. Our findings highlight the necessity of integrating both morphological and multiple molecular datasets for precise species identification, particularly when delineating closely related species. Additionally, it reveals the important role of mountain orogenesis on speciation within the southwestern mountains of China.

Complex patterns of gene flow and convergence in the evolutionary history of the spiral-horned antelopes (Tragelaphini)

The Tragelaphini, also known as spiral-horned antelope, is a phenotypically diverse mammalian tribe comprising a single genus, Tragelaphus. The evolutionary history of this tribe has attracted the attention of taxonomists and molecular geneticists for decades because its diversity is characterised by conflicts between morphological and molecular data as well as between mitochondrial, nuclear and chromosomal DNA. These inconsistencies point to a complex history of ecological diversification, coupled by either phenotypic convergence or introgression. Therefore, to unravel the phylogenetic relationships among spiral-horned antelopes, and to further investigate the role of divergence and gene flow in trait evolution, we sequenced genomes for all nine accepted species of the genus Tragelaphus, including a genome each for the highly divergent bushbuck lineages (T. s. scriptus and T. s. sylvaticus). We successfully reconstructed the Tragelaphus species tree, providing genome-level support for the early Pliocene divergence and monophyly of the nyala (T. angasii) and lesser kudu (T. imberbis), the monophyly of the two eland species (T. oryx and T. derbianus) and, importantly, the monophyly of kéwel (T. s. scriptus) and imbabala (T. s. sylvaticus) bushbuck. We found strong evidence for gene flow in at least four of eight nodes on the species tree. Among the six phenotypic traits assessed here, only habitat type mapped onto the species tree without homoplasy, showing that trait evolution was the result of complex patterns of divergence, introgression and convergent evolution.

Form, function, mind: What doesn't compute (and what might)

The applicability of computational and dynamical systems models to organisms is scrutinized, using examples from developmental biology and cognition. Developmental morphogenesis is dependent on the inherent material properties of developing animal (metazoan) tissues, a non-computational modality, but cell differentiation, which utilizes chromatin-based revisable memory banks and program-like function-calling, via the developmental gene co-expression system unique to the metazoans, has a quasi-computational basis. Multi-attractor dynamical models are argued to be misapplied to global properties of development, and it is suggested that along with computationalism, classic forms of dynamicism are similarly unsuitable to accounting for cognitive phenomena. Proposals are made for treating brains and other nervous tissues as novel forms of excitable matter with inherent properties which enable the intensification of cell-based basal cognition capabilities present throughout the tree of life. Finally, some connections are drawn between the viewpoint described here and active inference models of cognition, such as the Free Energy Principle.

Insights into the panorama of multiple DNA viruses in municipal wastewater and recycled sludge in Tianjin, China

Environmental viruses in wastewater and sludge are widely recognized for their roles in waterborne diseases. However, previous studies mainly focused on RNA viruses, and little is known about the diversity of DNA viral communities and their driving factors in municipal wastewater treatment environments. Herein, we conducted a pilot study to explore DNA virus profiles in municipal wastewater and recycled sludge by metagenomics method, and track their temporal changes in northern China. Results showed that 467 viral species were co-shared among all the samples. We identified six families of human viruses with a prevalence of 0.1%, which were rare but relatively stable in wastewater and sludge for six months. Adenoviridae, Parvoviridae, and Herpersviridae were the most dominant human viral families in municipal wastewater and recycled sludge. A time series of samples revealed that the dynamic changes of human DNA viruses were stable based on qPCR results, particularly for high-risk fecal-oral transmission viruses of adenovirus, bocavirus, polyomavirus, human gamma herpesvirus, human papillomavirus, and hepatitis B virus. Concentrations of Adenovirus (5.39-7.48 log10 copies/L) and bocavirus (4.36-7.48 log10 copies/L) were observed to be the highest in these samples compared to other viruses. Our findings demonstrated the DNA viruses' high prevalence and persistence in municipal wastewater treatment environments, highlighting the value of enhancing public health responses based on wastewater-based epidemiology.