Assessment of methods for amino acid matrix selection and their use on empirical data shows that ad hoc assumptions for choice of matrix are not justified

Phylogenomic analyses of all species of swordtail fishes (genus Xiphophorus) show that hybridization preceded speciation

Hybridization has been recognized to play important roles in evolution, however studies of the genetic consequence are still lagging behind in vertebrates due to the lack of appropriate experimental systems. Fish of the genus Xiphophorus are proposed to have evolved with multiple ancient and ongoing hybridization events. They have served as an informative research model in evolutionary biology and in biomedical research on human disease for more than a century. Here, we provide the complete genomic resource including annotations for all described 26 Xiphophorus species and three undescribed taxa and resolve all uncertain phylogenetic relationships. We investigate the molecular evolution of genes related to cancers such as melanoma and for the genetic control of puberty timing, focusing on genes that are predicted to be involved in pre-and postzygotic isolation and thus affect hybridization. We discovered dramatic size-variation of some gene families. These persisted despite reticulate evolution, rapid speciation and short divergence time. Finally, we clarify the hybridization history in the entire genus settling disputed hybridization history of two Southern swordtails. Our comparative genomic analyses revealed hybridization ancestries that are manifested in the mosaic fused genomes and show that hybridization often preceded speciation.

Ancient Mitochondrial Genomes Provide New Clues in the History of the Akhal-Teke Horse in China

This study analyzed ancient DNA from the remains of horses unearthed from the Shihuyao tombs. These were found to date from the Han and Tang Dynasties in Xinjiang (approximately 2200 to 1100 years ago). Two high-quality mitochondrial genomes were acquired and analyzed using next-generation sequencing. The genomes were split into two maternal haplogroups, B and D, according to a study that included ancient and contemporary samples from Eurasia. A close genetic affinity was observed between the horse of the Tang Dynasty and Akhal-Teke horses according to the primitive horse haplotype G1. Historical evidence suggests that the ancient Silk Road had a vital role in their dissemination. Additionally, the matrilineal history of the Akhal-Teke horse was accessed and suggested that the early domestication of the breed was for military purposes.

Astroviruses in terrestrial Malagasy mammals

Small terrestrial mammals are major hosts of infectious agents responsible for zoonotic diseases. Astroviruses (AstVs)-the cause of non-bacterial gastroenteritis mainly affecting young children-have been detected in a wide array of mammalian and avian host species. However, understanding the factors that influence AstV infection within and across hosts is limited. Here, we investigated the impact of land use changes on AstVs in terrestrial small mammals in rural northeastern Madagascar. We sampled 515 small mammals, representing seven endemic and four introduced species. Twenty-two positive samples were identified, all but one of which were found in the introduced species Mus musculus and Rattus rattus (family Muridae), with a positivity rate of 7.7% (6/78) and 5.6% (15/266), respectively. The non-introduced rodent case was from an endemic shrew-tenrec (family Tenrecidae). We found the highest positivity rate of AstVs infection in brushy regrowth (17.5%, 7/40) as compared to flooded rice fields (4.60%, 8/174), secondary forest (4.1%, 3/74), agroforest (3.6%, 1/28), village (2.61%, 3/115), and semi-intact forest (0%, 0/84). A phylogenetic analysis revealed an association between AstVs and their rodent host species. None of the viruses were phylogenetically related to AstVs previously described in Malagasy bats. This study supports AstV circulation in synanthropic animals in agricultural habitats of Madagascar and highlights the need to assess the spillover risk to human populations in rural areas.

High-throughput sequencing of Strongyloides stercoralis - a fatal disseminated infection in a dog

The rhabditid nematode Strongyloides stercoralis is known worldwide as the causative agent of strongyloidiasis in humans. In addition to public health concerns, S. stercoralis also infects dogs, which represent a possible reservoir for potentially zoonotic transmissions. We describe the first confirmed case of fatal disseminated infection in a dog in the Czech Republic. The microscopic and histological results were supported by a complex genotyping approach. Using high-throughput sequencing of the hypervariable region (HVR-IV) of 18S rDNA and Sanger sequencing of the partial cytochrome c oxidase subunit 1 gene (cox1), the potentially zoonotic haplotype/lineage A of S. stercoralis was confirmed, while the solely canine haplotype/lineage B was not found. The development of the disease is mainly associated with immunodeficiency, and in this case, it was triggered by inappropriate treatment, in particular the use of corticosteroids.

Lighting lantern above Psalteriomonadidae: Unveiling novel diversity within the genus Psalteriomonas (Discoba: Heterolobosea)

Psalteriomonadidae are a small family of anaerobic free-living protists belonging to Heterolobosea, Discoba. We cultured 74 new strains of mostly amoeboid Psalteriomonadidae obtained from mainly freshwater habitats and sequenced their 18S rRNA gene. Based on the phylogenetic analysis and genetic distances, we report multiple novel species, four of which we formally describe based on the light-microscopic morphology (Psalteriomonas minuta, P. australis, P. fimbriata, and P. parva). We also examined the ultrastructure of two Psalteriomonas species using transmission electron microscopy. We transfer Sawyeria marylandensis into the genus Psalteriomonas and synonymize Sawyeria with Psalteriomonas. In addition, we studied the flagellate stage of P. marylandensis comb. nov. for the first time, using light and scanning electron microscopy.

High-throughput analysis of the Trypanosoma cruzi minicirculome (mcDNA) unveils structural variation and functional diversity

Trypanosoma cruzi causes Chagas disease and has a unique extranuclear genome enclosed in a structure called the kinetoplast, which contains circular genomes known as maxi- and minicircles. While the structure and function of maxicircles are well-understood, many aspects of minicircles remain to be discovered. Here, we performed a high-throughput analysis of the minicirculome (mcDNA) in 50 clones isolated from Colombia's diverse T. cruzi I populations. Results indicate that mcDNA comprises four diverse subpopulations with different structures, lengths, and numbers of interspersed semi-conserved (previously termed ultra-conserved regions mHCV) and hypervariable (mHVPs) regions. Analysis of mcDNA ancestry and inter-clone differentiation indicates the interbreeding of minicircle sequence classes is placed along diverse strains and hosts. These results support evidence of the multiclonal dynamics and random bi-parental segregation. Finally, we disclosed the guide RNA repertoire encoded by mcDNA at a clonal scale, and several attributes of its abundance and function are discussed.

Selection among site-dependent structurally constrained substitution models of protein evolution by approximate Bayesian computation

MOTIVATION

The selection among substitution models of molecular evolution is fundamental for obtaining accurate phylogenetic inferences. At the protein level, evolutionary analyses are traditionally based on empirical substitution models but these models make unrealistic assumptions and are being surpassed by structurally constrained substitution (SCS) models. The SCS models often consider site-dependent evolution, a process that provides realism but complicates their implementation into likelihood functions that are commonly used for substitution model selection.

RESULTS

We present a method to perform selection among site-dependent SCS models, also among empirical and site-dependent SCS models, based on the approximate Bayesian computation (ABC) approach and its implementation into the computational framework ProteinModelerABC. The framework implements ABC with and without regression adjustments and includes diverse empirical and site-dependent SCS models of protein evolution. Using extensive simulated data, we found that it provides selection among SCS and empirical models with acceptable accuracy. As illustrative examples, we applied the framework to analyze a variety of protein families observing that SCS models fit them better than the corresponding best-fitting empirical substitution models.

AVAILABILITY AND IMPLEMENTATION

ProteinModelerABC is freely available from https://github.com/DavidFerreiro/ProteinModelerABC, can run in parallel and includes a graphical user interface. The framework is distributed with detailed documentation and ready-to-use examples.

PtoMYB031, the R2R3 MYB transcription factor involved in secondary cell wall biosynthesis in poplar

Introduction

The biosynthesis of the secondary cell wall (SCW) is orchestrated by an intricate hierarchical transcriptional regulatory network. This network is initiated by first-layer master switches, SCW-NAC transcription factors, which in turn activate the second-layer master switches MYBs. These switches play a crucial role in regulating xylem specification and differentiation during SCW formation. However, the roles of most MYBs in woody plants are yet to be fully understood.

Methods

In this study, we identified and isolated the R2R3-MYB transcription factor, PtoMYB031, from Populus tomentosa. We explored its expression, mainly in xylem tissues, and its role as a transcriptional repressor in the nucleus. We used overexpression and RNA interference techniques in poplar, along with Yeast two-hybrid (Y2H) and bimolecular fluorescence complementation (BiFC) assays, to analyze the regulatory effects of PtoMYB031.

Results

Overexpression of PtoMYB031 in poplar significantly reduced lignin, cellulose, and hemicellulose content, and inhibited vascular development in stems, resulting in decreased SCW thickness in xylem tissues. Gene expression analysis showed that structural genes involved in SCW biosynthesis were downregulated in PtoMYB031-OE lines. Conversely, RNA interference of PtoMYB031 increased these compounds. Additionally, PtoMYB031 was found to recruit the repressor PtoZAT11, forming a transcriptional inhibition complex.

Discussion

Our findings provide new insights into how PtoMYB031, through its interaction with PtoZAT11, forms a complex that can suppress the expression of key regulatory genes, PtoWND1A and PtoWND2B, in SCW biosynthesis. This study enhances our understanding of the transcriptional regulation involved in SCW formation in poplar, highlighting the significant role of PtoMYB031.

Characterisation of putative novel tick viruses and zoonotic risk prediction

Tick-associated viruses remain a substantial zoonotic risk worldwide, so knowledge of the diversity of tick viruses has potential health consequences. Despite their importance, large amounts of sequences in public data sets from tick meta-genomic and -transcriptomic projects remain unannotated, sequence data that could contain undocumented viruses. Through data mining and bioinformatic analysis of more than 37,800 public meta-genomic and -transcriptomic data sets, we found 83 unannotated contigs exhibiting high identity with known tick viruses. These putative viral contigs were classified into three RNA viral families (Alphatetraviridae, Orthomyxoviridae and Chuviridae) and one DNA viral family (Asfarviridae). After manual checking of quality and dissimilarity towards other sequences in the data set, these 83 contigs were reduced to five contigs in the Alphatetraviridae from four putative viruses, four in the Orthomyxoviridae from two putative viruses and one in the Chuviridae which clustered with known tick-associated viruses, forming a separate clade within the viral families. We further attempted to assess which previously known tick viruses likely represent zoonotic risks and thus deserve further investigation. We ranked the human infection potential of 133 known tick-associated viruses using a genome composition-based machine learning model. We found five high-risk tick-associated viruses (Langat virus, Lonestar tick chuvirus 1, Grotenhout virus, Taggert virus and Johnston Atoll virus) that have not been known to infect human and two viral families (Nairoviridae and Phenuiviridae) that contain a large proportion of potential zoonotic tick-associated viruses. This adds to the knowledge of tick virus diversity and highlights the importance of surveillance of newly emerging tick-associated diseases.

Hydrogenimonas cancrithermarum sp. nov., a hydrogen- and thiosulfate-oxidizing mesophilic chemolithoautotroph isolated from diffuse-flow fluids on the East Pacific Rise, and an emended description of the genus Hydrogenimonas

A novel mesophilic, hydrogen- and thiosulfate-oxidizing bacterium, strain ISO32T, was isolated from diffuse-flow hydrothermal fluids from the Crab Spa vent on the East Pacific Rise. Cells of ISO32T were rods, being motile by means of a single polar flagellum. The isolate grew at a temperature range between 30 and 55 °C (optimum, 43 °C), at a pH range between 5.3 and 7.6 (optimum, pH 5.8) and in the presence of 2.0-4.0 % NaCl (optimum, 2.5 %). The isolate was able to grow chemolithoautotrophically with molecular hydrogen, thiosulfate or elemental sulfur as the sole electron donor. Thiosulfate, elemental sulfur, nitrate and molecular oxygen were each used as a sole electron acceptor. Phylogenetic analysis of 16S rRNA gene sequences placed ISO32T in the genus Hydrogenimonas of the class Epsilonproteobacteria, with Hydrogenimonas thermophila EP1-55-1 %T as its closest relative (95.95 % similarity). On the basis of the phylogenetic, physiological and genomic characteristics, it is proposed that the organism represents a novel species within the genus Hydrogenimonas, Hydrogenimonas cancrithermarum sp. nov. The type strain is ISO32T (=JCM 39185T =KCTC 25252T). Furthermore, the genomic properties of members of the genus Hydrogenimonas are distinguished from those of members of other thermophilic genera in the orders Campylobacterales (Nitratiruptor and Nitrosophilus) and Nautiliales (Caminibacter, Nautilia and Lebetimonas), with larger genome sizes and lower 16S rRNA G+C content values. Comprehensive metabolic comparisons based on genomes revealed that genes responsible for the Pta-AckA pathway were observed exclusively in members of mesophilic genera in the order Campylobacterales and of the genus Hydrogenimonas. Our results indicate that the genus Hydrogenimonas contributes to elucidating the evolutionary history of Epsilonproteobacteria in terms of metabolism and transition from a thermophilic to a mesophilic lifestyle.

Study on the Interactions of Cyclins with CDKs Involved in Auxin Signal during Leaf Development by WGCNA in Populus alba

Cell division plays an indispensable role in leaf morphogenesis, which is regulated via the complexes formed by cyclin and cyclin-dependent kinase (CDK). In this study, gene family analysis, exogenous auxin stimulation, RNA-seq and WGCNA analysis were all used to investigate the molecular mechanisms by which cell-cycle-related factors participated in the auxin signaling pathway on leaf morphogenesis. Sixty-three cyclin members and seventeen CDK members in Populus alba were identified and systematically analyzed. During the evolution, WGD was the main reason that resulted in the expansion of cyclin and CDK genes. Firstly, after a short time treating with auxin to matured leaves of seedlings, genes related to cell division including GRF and ARGOS were both upregulated to restart the transition of cells from G1-to-S phase. Secondly, with three days of continuous auxin stimulation to leaves at different developmental stages, leaves area variation, transcriptomes and hormones were analyzed. By PCA, PCoA and WGCNA analyses, the turquoise module was both positively related to leaf development and auxin. Based on the co-expression analysis and Y2H experiment, PoalbCYCD1;4, PoalbCYCD3;3 and PoalbCYCD3;5 were supposed to interact with PoalbCDKA;1, which could be the trigger to promote the G1-to-S phase transition. The ARF transcription factor might play the key role of connecting the auxin signaling pathway and cell division in leaf morphogenesis by affecting CYC-CDK complexes.

The complete chloroplast genome of Salix lindleyana (salicaceae), a plateau plant species

Salix lindleyana Wallich ex Andersson 1851 is a species of genus Salix which mainly grows on mountains above 3000 m at sea level in Qinghai-Tibetan Plateau (including the Himalayas and Hengduan Mountains). To determine its phylogenetic position within Salix, we reconstructed S. lindleyana complete chloroplast (cp) genome sequence by de novo assembly using whole-genome sequencing data. The completed chloroplast genome was 155,304 bp, with a total GC content of 36.7%. It had a very typical tetrad structure, including a large single-copy (LSC) region of 84,539 bp, a small single-copy (SSC) region of 16,161 bp, and two inverted repeats (IR) regions of 27,302 bp. A total of 132 functional genes were distributed in the chloroplast genome, including 87 protein-coding genes, 37 tRNA genes, and 8 rRNA genes. Phylogenetic analysis showed that S. lindleyana was clustered with Salix dasyclados Wimmer 1849 and Salix variegata Franchet 1887. The complete chloroplast genome of S. lindleyana provides potential genetic resources for further phylogenetic studies.

Data-specific substitution models improve protein-based phylogenetics

Calculating amino-acid substitution models that are specific for individual protein data sets is often difficult due to the computational burden of estimating large numbers of rate parameters. In this study, we tested the computational efficiency and accuracy of five methods used to estimate substitution models, namely Codeml, FastMG, IQ-TREE, P4 (maximum likelihood), and P4 (Bayesian inference). Data-specific substitution models were estimated from simulated alignments (with different lengths) that were generated from a known simulation model and simulation tree. Each of the resulting data-specific substitution models was used to calculate the maximum likelihood score of the simulation tree and simulated data that was used to calculate the model, and compared with the maximum likelihood scores of the known simulation model and simulation tree on the same simulated data. Additionally, the commonly-used empirical models, cpREV and WAG, were assessed similarly. Data-specific models performed better than the empirical models, which under-fitted the simulated alignments, had the highest difference to the simulation model maximum-likelihood score, clustered further from the simulation model in principal component analysis ordination, and inferred less accurate trees. Data-specific models and the simulation model shared statistically indistinguishable maximum-likelihood scores, indicating that the five methods were reasonably accurate at estimating substitution models by this measure. Nevertheless, tree statistics showed differences between optimal maximum likelihood trees. Unlike other model estimating methods, trees inferred using data-specific models generated with IQ-TREE and P4 (maximum likelihood) were not significantly different from the trees derived from the simulation model in each analysis, indicating that these two methods alone were the most accurate at estimating data-specific models. To show the benefits of using data-specific protein models several published data sets were reanalysed using IQ-TREE-estimated models. These newly estimated models were a better fit to the data than the empirical models that were used by the original authors, often inferred longer trees, and resulted in different tree topologies in more than half of the re-analysed data sets. The results of this study show that software availability and high computation burden are not limitations to generating better-fitting data-specific amino-acid substitution models for phylogenetic analyses.

Anaerocyclidiidae fam. nov. (Oligohymenophorea, Scuticociliatia): A newly recognized major lineage of anaerobic ciliates hosting prokaryotic symbionts

The research on anaerobic ciliates, to date, has mainly been focused on representatives of obligately anaerobic classes such as Armophorea or Plagiopylea. In this study, we focus on the anaerobic representatives of the subclass Scuticociliatia, members of the class Oligohymenophorea, which is mainly composed of aerobic ciliates. Until now, only a single anaerobic species, Cyclidium porcatum (here transferred to the genus Anaerocyclidium gen. nov.), has been described both molecularly and morphologically. Our broad sampling of anoxic sediments together with cultivation and single cell sequencing approaches have shown that scuticociliates are common and diversified in anoxic environments. Our results show that anaerobic scuticociliates represent a distinctive evolutionary lineage not closely related to the family Cyclidiidae (order Pleuronematida), as previously suggested. However, the phylogenetic position of the newly recognized lineage within the subclass Scuticociliatia remains unresolved. Based on molecular and morphological data, we establish the family Anaerocyclidiidae fam. nov. to accommodate members of this clade. We further provide detailed morphological descriptions and 18S rRNA gene sequences for six new Anaerocyclidium species and significantly broaden the described diversity of anaerobic scuticociliates.

Diversity, phylogeny and intraspecific variability of Paradiplozoon species (Monogenea: Diplozoidae) parasitizing endemic cyprinoids in the Middle East

Diplozoidae are common monogenean ectoparasites of cyprinoid fish, with the genus Paradiplozoon being the most diversified. Despite recent studies on Diplozoidae from Europe, Africa and Asia, the diversity, distribution and phylogeny of this parasite group appears to be still underestimated in the Middle East. The objective of this study was to investigate the diversity, endemism and host specificity of diplozoids parasitizing cyprinoid fish from the Middle East, considering this region as an important historical interchange of fish fauna, and to elucidate the phylogenetic position of Middle Eastern Paradiplozoon species within Diplozoidae. Four Paradiplozoon species were collected from 48 out of 94 investigated cyprinoid species. Three known species, Paradiplozoon homoion, Paradiplozoon bliccae and Paradiplozoon bingolensis, were recorded on new cyprinoid host species, and a new species, Paradiplozoon koubkovae n. sp., was recorded on Luciobarbus capito and Capoeta capoeta from the Caspian Sea basin in Iran and Turkey. Paradiplozoon bliccae, exhibiting a wide host range in the Middle East, expressed both morphological and genetic intraspecific variabilities. The four Paradiplozoon species collected in the Middle East were placed in divergent clades, showing the rich evolutionary history of diplozoid parasites in the Middle East. Our study also revealed that two lineages of African diplozoids have a Middle Eastern origin. We stress the importance of applying an integrative approach combining morphological, ecological and molecular methods to reveal the real diversity of diplozoids.

Maternal Ancestry of First Parsi Settlers of India Using Ancient Mitogenome

The rich cultural and genetic diversity of South Asia emerged from multiple migrations and cultural assimilation of multiple waves of migrants. The Parsi community of North-western India were one of those who migrated from West Eurasia in the aftermath of 7^th century CE and assimilated into the local cultural framework. Earlier genetic studies further strengthened this notion with the finding that they harbour both Middle Eastern and South Asian genetic components. Although these studies covered both autosomal and uniparental markers, still maternal ancestry was not covered in depth and with good resolution of mitochondrial markers. Hence in our current study, we have first time generated a complete mitogenome of 19 ancient samples of the first Parsi settlers excavated from the archaeological site of Sanjan and performed detailed phylogenetic analysis to infer their maternal genetic affinity. In our analysis, we found that the Parsi mitogenome with mtDNA haplogroup M3a1+204 shares clade with both Middle Eastern and South Asian modern individuals in both the Maximum Likelihood tree and Bayesian phylogenetic tree. This haplogroup was also prevalent among the medieval Swat valley population of present-day Northern Pakistan and was also observed in two Roopkund A individuals. In the phylogenetic network this sample share haplotype with both South Asian and Middle Eastern samples. So conclusively, the first Parsi settlers' maternal ancestry encompasses both South Asian and Middle Eastern genetic composition.

Parallel evolution of two AIM24 protein subfamilies and their conserved functions in ER stress tolerance in land plants

Despite decades of efforts in genome sequencing and functional characterization, some important protein families remain poorly understood. In this study, we report the classification, evolution, and functions of the largely uncharacterized AIM24 protein family in plants, including the identification of a novel subfamily. We show that two AIM24 subfamilies (AIM24-A and AIM24-B) are commonly distributed in major plant groups. These two subfamilies not only have modest sequence similarities and different gene structures but also are of independent bacterial ancestry. We performed comparative functional investigations on the two AIM24 subfamilies using three model plants: the moss Physcomitrium patens, the liverwort Marchantia polymorpha, and the flowering plant Arabidopsis thaliana. Intriguingly, despite their significant differences in sequence and gene structure, both AIM24 subfamilies are involved in ER stress tolerance and the unfolded protein response (UPR). In addition, transformation of the AIM24-A gene from P. patens into the AIM24-B null mutant of A. thaliana could at least partially rescue ER stress tolerance and the UPR. We also discuss the role of AIM24 genes in plant development and other cellular activities. This study provides a unique example of parallel evolution in molecular functions and can serve as a foundation for further investigation of the AIM24 family in plants.

Identification of potentially zoonotic parasites in captive orangutans and semi-captive mandrills: Phylogeny and morphological comparison

Cysts and trophozoites of vestibuliferid ciliates and larvae of Strongyloides were found in fecal samples from captive orangutans Pongo pygmaeus and P. abelii from Czech and Slovak zoological gardens. As comparative material, ciliates from semi-captive mandrills Mandrillus sphinx from Gabon were included in the study. Phylogenetic analysis of the detected vestibuliferid ciliates using ITS1-5.8s-rRNA-ITS2 and partial 18S ribosomal deoxyribonucleic acid (rDNA) revealed that the ciliates from orangutans are conspecific with Balantioides coli lineage A, while the ciliates from mandrills clustered with Buxtonella-like ciliates from other primates. Morphological examination of the cysts and trophozoites using light microscopy did not reveal differences robust enough to identify the genera of the ciliates. Phylogenetic analysis of detected L1 larvae of Strongyloides using partial cox1 revealed Strongyloides stercoralis clustering within the cox1 lineage A infecting dogs, humans, and other primates. The sequences of 18S rDNA support these results. As both B. coli and S. stercoralis are zoonotic parasites and the conditions in captive and semi-captive settings may facilitate transmission to humans, prophylactic measures should reflect the findings.

Lost characters and lost taxonomists: Coelocarteriidae fam. nov., Poecilosclerida with and without chelae, to accommodate supposed phloeodictyids (Haplosclerida)

The taxonomic study of Great Amazon Reef System sponges yielded three Oceanapia-like (Phloeodictyidae, Haplosclerida) species of similar gross morphology, two preliminarily identified as O. bartschi and another as Coelocarteria (Poecilosclerida), due to the presence of common palmate isochelae. The remarkable overall similarity of all three species in gross morphology necessitated an integrative assessment of the phylogenetic affinities. A selection of haplosclerid and poecilosclerid sequences (18S and 28S rRNA) was gathered from Genbank and compared to sequences mapped to reference from metagenome of two Oceanapia-like species from the Amazon River mouth, one of which matched O. bartschi. Both Brazilian species clustered with Coelocarteria singaporensis (Singapore). These species nest in the Poecilosclerida, far from Oceanapia (sp. and isodictyiformis) and other haplosclerids (Amphimedon, Petrosia, Siphonodictyon and Xestospongia) but also far from the poecilosclerid Isodictya that is currently classified in the same family as Coelocarteria, the Isodictyidae. Specimens with chelae are named Coelocarteria alcoladoi sp. nov. herein, while those without chelae represent the other two species. One matched Inflatella bartschi (O. bartchi’s holotype, here re-examined), thereby supporting the transfer of this species to Coelocarteria. The other is proposed as C. amadoi sp. nov. and is the second known lipochelous species in this genus. The 28S phylogenies recovered Coelocarteria bartschi comb. nov. (formerly Oceanapia bartschi) in the Poecilosclerida clade that clustered with Coelocarteria spp., including the type species of this genus, C. singaporensis, with 100% support. Coelocarteria alcoladoi sp. nov., also without chelae, grouped in the same clade, thereby corroborating the classification in this genus. The clade composed of Coelocarteria spp. grouped with Cornulella sp., suggesting an affinity between these genera. Coelocarteria is currently situated within Isodictyidae and Cornulella in Acarnidae. Isodictya (Isodictyidae) grouped with Amphilectus (Esperiopsidae) with high support and is only distantly related to Coelocarteria. Acarnus (Acarnidae) grouped with Clathria, also with high support, far from Coelocarteria and Cornulella. These results suggest the polyphyletic nature of Isodictyidae and Acarnidae, for which reason we preferred to propose a new, currently monotypic family for Coelocarteria spp., Coelocarteriidae fam. nov. ZooBank: urn:lsid:zoobank.org:pub:71FDB6FD-4A5F-4180-8DA7-79EA4CB615D1

Genome-Wide Identification of 2-Oxoglutarate and Fe (II)-Dependent Dioxygenase (2ODD-C) Family Genes and Expression Profiles under Different Abiotic Stresses in Camellia sinensis (L.)

The 2-oxoglutarate and Fe (II)-dependent dioxygenase (2ODD-C) family of 2-oxoglutarate-dependent dioxygenases potentially participates in the biosynthesis of various metabolites under various abiotic stresses. However, there is scarce information on the expression profiles and roles of 2ODD-C genes in Camellia sinensis. We identified 153 Cs2ODD-C genes from C. sinensis, and they were distributed unevenly on 15 chromosomes. According to the phylogenetic tree topology, these genes were divided into 21 groups distinguished by conserved motifs and an intron/exon structure. Gene-duplication analyses revealed that 75 Cs2ODD-C genes were expanded and retained after WGD/segmental and tandem duplications. The expression profiles of Cs2ODD-C genes were explored under methyl jasmonate (MeJA), polyethylene glycol (PEG), and salt (NaCl) stress treatments. The expression analysis showed that 14, 13, and 49 Cs2ODD-C genes displayed the same expression pattern under MeJA and PEG treatments, MeJA and NaCl treatments, and PEG and NaCl treatments, respectively. A further analysis showed that two genes, Cs2ODD-C36 and Cs2ODD-C21, were significantly upregulated and downregulated after MeJA, PEG, and NaCl treatments, indicating that these two genes played positive and negative roles in enhancing the multi-stress tolerance. These results provide candidate genes for the use of genetic engineering technology to modify plants by enhancing multi-stress tolerance to promote phytoremediation efficiency.

Genome-wide identification, expression profile and evolutionary relationships of TPS genes in the neotropical fruit tree species Psidium cattleyanum

Terpenoids are essential for plant growth, development, defense, and adaptation mechanisms. Psidium cattleyanum (Myrtaceae) is a fleshy fruit tree species endemics from Atlantic Forest, known for its pleasant fragrance and sweet taste, attributed to terpenoids in its leaves and fruits. In this study, we conducted genome-wide identification, evolutionary and expression analyses of the terpene synthase gene (TPS) family in P. cattleyanum red guava (var. cattleyanum), and yellow guava (var. lucidum Hort.) morphotypes. We identified 32 full-length TPS in red guava (RedTPS) and 30 in yellow guava (YlwTPS). We showed different expression patterns of TPS paralogous in the two morphotypes, suggesting the existence of distinct gene regulation mechanisms and their influence on the final essential oil content in both morphotypes. Moreover, the oil profile of red guava was dominated by 1,8-cineole and linalool and yellow guava was enriched in α-pinene, coincident in proportion to TPS-b1 genes, which encode enzymes that produce cyclic monoterpenes, suggesting a lineage-specific subfamily expansion of this family. Finally, we identified amino acid residues near the catalytic center and functional areas under positive selection. Our findings provide valuable insights into the terpene biosynthesis in a Neotropical Myrtaceae species and their potential involvement in adaptation mechanisms.

Morphology and phylogenetic position of three anaerobic ciliates from the classes Odontostomatea and Muranotrichea (Ciliophora)

The diversity of the classes Odontostomatea and Muranotrichea, which contain solely obligate anaerobes, is poorly understood. We studied two populations of Mylestoma sp., one of Saprodinium dentatum (Odontostomatea), two of Muranothrix felix sp. nov., and one of Muranothrix sp. (Muranotrichea) employing live observation, protargol impregnation, scanning electron microscopy, and 18S rRNA gene sequencing. Conspecificity of Mylestoma sp., described here, with a previously described species of this genus cannot be excluded since no species have been studied with modern methods. Phylogenetically, the genus Mylestoma is closely related to the odontostomatid Discomorphella pedroeneasi, although the phylogenetic position of class Odontostomatea itself remains unresolved. The newly described muranotrichean species, Muranothrix felix sp. nov., is morphologically similar to M. gubernata but can be distinguished by its fewer macronuclear nodules and fewer adoral membranelles; moreover, it is clearly distinguished from M. gubernata by its 18S rRNA gene sequence. Another population, designated here as Muranothrix sp., most likely represents a separate species.

Canonical or noncanonical? Structural plasticity of serine protease-binding loops in Kunitz-STI protease inhibitors

The Kunitz-Soybean Trypsin Inhibitor (Kunitz-STI) family is a large family of proteins with most of its members being protease inhibitors. The versatility of the inhibitory profile and the structural plasticity of these proteins, make this family a promising scaffold for designing new multifunctional proteins. Historically, Kunitz-STI inhibitors have been classified as canonical serine protease inhibitors, but new inhibitors with novel inhibition mechanisms have been described in recent years. Different inhibition mechanisms could be the result of different evolutionary pathways. In the present work, we performed a structural analysis of all the crystallographic structures available for Kunitz-STI inhibitors to characterize serine protease-binding loop structural features and locations. Our study suggests a relationship between the conformation of serine protease-binding loops and the inhibition mechanism, their location in the β-trefoil fold, and the plant source of the inhibitors. The classical canonical inhibitors of this family are restricted to plants from the Fabales order and bind their targets via the β4-β5 loop, whereas serine protease-binding loops in inhibitors from other plants lie mainly in the β5-β6 and β9-β10 loops. In addition, we found that the β5-β6 loop is used to inhibit two different families of serine proteases through a steric blockade inhibition mechanism. This work will help to change the general perception that all Kunitz-STI inhibitors are canonical inhibitors and proteins with protease-binding loops adopting noncanonical conformations are exceptions. Additionally, our results will help in the identification of protease-binding loops in uncharacterized or newly discovered inhibitors, and in the design of multifunctional proteins.

Pratylenchus Araucensis (Rhabditida: Pratylenchidae) a Widely Distributed Nematode in Musa Spp. from Colombia

Pratylenchus is one of the most limiting nematodes of Musaceae production in the world. Knowledge of the nematode species is one of the requirements for its management in the field. This study aimed to identify up to the species level Pratylenchus populations associated with plantain and banana crops in the states of Caldas, Quindío, and Valle del Cauca in Colombia. In these regions, Pratylenchus has been reported to affect these crops in the past, but with records of the nematode only up to the genus level. For this purpose, five populations of Pratylenchus extracted from samples composed of roots and rhizospheric soil, four from plantain and one from banana, were identified through morphological, morphometric, and molecular analysis (sequencing of the D2-D3 of rRNA and cytochrome oxidase subunit I of mtDNA). All populations were identified as P. araucensis, a species reported previously in eastern Colombia, and one that the present study found in the center and southwest of the country, indicating that this species of nematode is widely distributed in the Musaceae-producing areas of Colombia. The present study reports the first COI mtDNA sequences for this species of nematode.

Hidden Diversity within Tetralophozia filiformis (Marchantiophyta, Anastrophyllaceae) in East Asia

Tetralophozia filiformis s.l. is known from a number of localities mostly in amphi-oceanic areas in Northern Hemisphere, including Atlantic Europe, amphi-Pacific Asia, South Siberia, and western North America. The newly obtained collections of this 'species' show strong variation in morphology of the taxon across amphi-Pacific Asia although connected by some 'intergrading' modifications. This implies the genetic diversity within this unit earlier recognized as a single taxon. Authors used molecular-genetic, morphological, and chorological methods to understand if the geographically correlated morphological variation also correlates with genetic differences and if it is possible to distinguish some additional taxa within the series of specimens originating from the various areas in amphi-Pacific Asia. It was found that Tetralophozia filiformis is a complex of at least three morphologically similar species, including one long forgotten name (Chandonanthus pusillus) that should be reinstated as separate species and one taxon (Tetralophozia sibirica) that should be described as new. Tetralophozia filiformis and Chandonanthus pusillus are lectotypified, and the new combination is provided for the latter. The three accepted taxa distinctly differ one from another in distribution patterns, preferable climate characteristics, and genetic distances, besides minor differentiations in morphology. The main morphological distinguishing features are the leaf cell size, height of undivided part in leaf lamina, and leaf dentation characteristics. Taking into account the robust correlation between the climate-based and molecular-genetic-based clusters, one more (fourth) taxon could be probably segregated from Tetralophozia filiformis.

Hypo-Osmoregulatory Roles of Vasotocinergic and Isotocinergic Systems in the Intestines of Two European Sea Bass Lineages

European sea bass (Dicentrarchus labrax) are a major aquaculture species that live in habitats with fluctuating salinities that are sometimes higher than in seawater (SW). Atlantic and West-Mediterranean genetic lineages were compared regarding intestinal neuropeptide receptor expression in SW (36%) and following a two-week transfer to hypersalinity (HW, 55%). Phylogenetic analysis revealed seven neuropeptide receptors belonging to the arginine vasotocine (AVTR) family and two isotocin receptors (ITR). Among AVTR paralogs, the highest mRNA levels were recorded for v1a2, with a two- to fourfold upregulation in the European sea bass intestinal sections after transfer of fish to HW. Principal component analysis in posterior intestines showed that v1a2 expression grouped together with the expression and activity of main ion transporters and channels involved in solute-coupled water uptake, indicating a possible role of this receptor in triggering water absorption. v1a1 expression, however, was decreased or did not change after transfer to hypersaline water. Among ITR paralogs, itr1 was the most expressed paralog in the intestine and opposite expression patterns were observed following salinity transfer, comparing intestinal sections. Overall, different expression profiles were observed between genetic lineages for several analyzed genes which could contribute to different osmotic stress-related responses in D. labrax lineages.

Study on the interaction preference between CYCD subclass and CDK family members at the poplar genome level

Cyclin-dependent kinases (CDKs) control the progression of the cell cycle. D-type cyclin (CYCD) is generally believed to form a complex with CDK and control the G1/S transition. In plants, CYCD and CDK gene families can be divided into 6 (D1-D7) and 7 (CDKA-CDKG) subclasses, respectively. Different subclasses in the CYCD and CDK families have different numbers, structures and functions. In some heterologous woody plants, the functions of these subclass family members remain unclear. In this study, 43 CYCD and 27 CDK gene family members were identified in the allodiploid Populus tomentosa Carr. Phylogenetic analysis suggested that these CYCDs and CDKs were divided into 6 and 7 subclasses, respectively, which were the same as other species. The analysis of protein properties, gene structure, motifs, domains, cis-acting elements and tissue-specific expression of all members of these CYCDs and CDKs showed that the differences between members of different subclasses varied widely, but members of the same subclass especially in the CDK gene family were very similar. These findings also demonstrated a strong correlation between CYCD and CDK gene family members in response to hormones and specific expression. The collinear analysis of P. tomentosa, Populus trichocarpa and Arabidopsis thaliana showed that the expansion patterns of CYCD and CDK gene families were predominantly whole genome duplications (WGD). The protein interaction prediction results of different subclasses of CYCD and CDKs showed that the interaction between different subclasses of CYCD and CDKs was significantly different. Our previous study found that transgenic PtoCYCD2;1 and PtoCYCD3;3 poplars exhibited opposite phenotypes. Y2H and BIFC results showed that the interaction between PtoCYCD2;1 and PtoCYCD3;3 was significantly different with CDKs. This finding might suggest that the functional differences of different CYCD subclasses in plant growth and development were closely related to the different interactions between CYCD and CDK. Our results provide a good idea and direction for the functional study of CYCD and CDK proteins in woody plants.

A Comparative Study of the Dynamics and Diversity of Bdellovibrio and Like Organisms in Lakes Annecy and Geneva

Bdellovibrio and like organisms (BALOs) are obligate bacterial predators of other Gram-negative bacteria. Here, we used quantitative PCR (qPCR) and recently developed specific primers which target the 16S rRNA gene to explore the abundance and distribution of three families of BALO belonging to the Oligoflexia class (i.e., Bdellovibrionaceae, Peredibacteraceae and Bacteriovoracaceae) over one year in the epilimnion and hypolimnion of Lakes Annecy and Geneva. Peredibacteraceae was the dominant group at all sampling points except at the bottom of Lake Geneva, where Bdellovibrionaceae was found in higher number. In addition, the abundance of BALOs increased significantly during the warmer months. Using high-throughput sequencing (Illumina Miseq), hundreds of OTUs were identified for Bdellovibrionaceae and Peredibacteraceae. Phylogenetic analysis suggests that Bdellovibrionaceae are more diverse than Peredibacteraceae and that some OTUs belong to new species of Bdellovibrionaceae. We also found that dominant OTUs were present simultaneously in the two lakes, while some others were specific to each lake, suggesting an adaptive pattern. Finally, both abundance and diversity of BALOs were poorly associated with abiotic factors except temperature, suggesting the importance of studying biotic relationships, assumed to play a greater role than physico-chemical variables in BALOs' dynamics and distribution.

Comparative chloroplast genome and phylogenetic analyses of Chinese Polyspora

Polyspora Sweet (Theaceae) are winter ornamental landscape plants native to southern and southeastern Asia, some of which have medicinal value. The chloroplast (cp) genome data of Polyspora are scarce, and the gene evolution and interspecific relationship are still unclear. In this study, we sequenced and annotated Polyspora chrysandra cp genome and combined it with previously published genomes for other Chinese Polyspora species. The results showed that cp genomes of six Chinese Polyspora varied in length between 156,452 bp (P. chrysandra) and 157,066 bp (P. speciosa), but all contained 132 genes, with GC content of 37.3%, and highly similar genes distribution and codon usage. A total of eleven intergenic spacer regions were found having the highest levels of divergence, and eight divergence hotspots were identified as molecular markers for Phylogeography and genetic diversity studies in Polyspora. Gene selection pressure suggested that five genes were subjected to positive selection. Phylogenetic relationships among Polyspora species based on the complete cp genomes were supported strongly, indicating that the cp genomes have the potential to be used as super barcodes for further analysis of the phylogeny of the entire genus. The cp genomes of Chinese Polyspora species will provide valuable information for species identification, molecular breeding and evolutionary analysis of genus Polyspora.

Analysis of lineage-specific protein family variability in prokaryotes combined with evolutionary reconstructions

Background

Evolutionary rate is a key characteristic of gene families that is linked to the functional importance of the respective genes as well as specific biological functions of the proteins they encode. Accurate estimation of evolutionary rates is a challenging task that requires precise phylogenetic analysis. Here we present an easy to estimate protein family level measure of sequence variability based on alignment column homogeneity in multiple alignments of protein sequences from Clade-Specific Clusters of Orthologous Genes (csCOGs).

Results

We report genome-wide estimates of variability for 8 diverse groups of bacteria and archaea and investigate the connection between variability and various genomic and biological features. The variability estimates are based on homogeneity distributions across amino acid sequence alignments and can be obtained for multiple groups of genomes at minimal computational expense. About half of the variance in variability values can be explained by the analyzed features, with the greatest contribution coming from the extent of gene paralogy in the given csCOG. The correlation between variability and paralogy appears to originate, primarily, not from gene duplication, but from acquisition of distant paralogs and xenologs, introducing sequence variants that are more divergent than those that could have evolved in situ during the lifetime of the given group of organisms. Both high-variability and low-variability csCOGs were identified in all functional categories, but as expected, proteins encoded by integrated mobile elements as well as proteins involved in defense functions and cell motility are, on average, more variable than proteins with housekeeping functions. Additionally, using linear discriminant analysis, we found that variability and fraction of genomes carrying a given gene are the two variables that provide the best prediction of gene essentiality as compared to the results of transposon mutagenesis in Sulfolobus islandicus.

Conclusions

Variability, a measure of sequence diversity within an alignment relative to the overall diversity within a group of organisms, offers a convenient proxy for evolutionary rate estimates and is informative with respect to prediction of functional properties of proteins. In particular, variability is a strong predictor of gene essentiality for the respective organisms and indicative of sub- or neofunctionalization of paralogs.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13062-022-00337-7.

Evidence for the involvement of branchial Vacuolar-type H+-ATPase in the acidification of the external medium by the West African lungfish, Protopterus annectens, exposed to ammonia-loading conditions

African lungfishes are obligatory air-breathers with exceptionally high environmental ammonia tolerance. They can lower the pH of the external medium during exposure to ammonia-loading conditions. This study aimed to demonstrate the possible involvement of branchial vacuolar-type H⁺-ATPase (Vha) in the ammonia-induced acidification of the external medium by the West African lungfish, Protopterus annectens, and to examine whether its capacity to acidify the medium could be augmented after exposure to 100 mmol l^-1 NH₄Cl for six days. Two full coding cDNA sequences of Vha subunit B (atp6v1b), atp6v1b1 and atp6v1b2, were obtained from the internal gills of P. annectens. The sequence of atp6v1b1 comprised 1548 bp, encoding 515 amino acids (57.4 kDa), while that of atp6v1b2 comprised 1536 bp, encoding 511 amino acids (56.6 kDa). Using a custom-made antibody reactive to both isoforms, immunofluorescence microscopy revealed the collective localization of Atp6v1b (atp6v1b1 and atp6v1b2) at the apical or the basolateral membrane of two different types of branchial Na⁺/K⁺-ATPase-immunoreactive ionocyte. The ionocytes labelled apically with Atp6v1b presumably expressed Atp6v1b1 containing a PDZ-binding domain, indicating that the apical Vha was positioned to transport H⁺ to the external medium. The expression of Atp6v1b was regulated post-transcriptionally, as the protein abundance of Atp6v1b and Vha activity increased significantly in the gills of fish exposed to 100 mmol l^-1 NH₄Cl for six days. Correspondingly, the fish exposed to ammonia had a greater capacity to acidify the external medium, presumably to decrease the ratio of [NH₃] to [NH₄⁺] in order to reduce the influx of exogenous NH₃.

Evolution of the Noncoding Features of Sea Snake Mitochondrial Genomes within Elapidae

Mitochondrial genomes of four elapid snakes (three marine species [Emydocephalus ijimae, Hydrophis ornatus, and Hydrophis melanocephalus], and one terrestrial species [Sinomicrurus japonicus]) were completely sequenced by a combination of Sanger sequencing, next-generation sequencing and Nanopore sequencing. Nanopore sequencing was especially effective in accurately reading through long tandem repeats in these genomes. This led us to show that major noncoding regions in the mitochondrial genomes of those three sea snakes contain considerably long tandem duplications, unlike the mitochondrial genomes previously reported for same and other sea snake species. We also found a transposition of the light-strand replication origin within a tRNA gene cluster for the three sea snakes. This change can be explained by the Tandem Duplication—Random Loss model, which was further supported by remnant intervening sequences between tRNA genes. Mitochondrial genomes of true snakes (Alethinophidia) have been shown to contain duplicate major noncoding regions, each of which includes the control region necessary for regulating the heavy-strand replication and transcription from both strands. However, the control region completely disappeared from one of the two major noncoding regions for two Hydrophis sea snakes, posing evolutionary questions on the roles of duplicate control regions in snake mitochondrial genomes. The timing and molecular mechanisms for these changes are discussed based on the elapid phylogeny.

Stolonicaulon: A Section-Puzzle within Marsupella (Gymnomitriaceae, Marchantiophyta)

Marsupella sect. Stolonicaulon is not speciose and is a commonly neglected section within the genus, which currently includes three species with somewhat similar morphologies (wiry shoots with distanced leaves) and distributions in the mountains of tropical and subtropical regions (SE (Southeast) Asia, the Venezuelan Andes, and the high mountains of SE Brazil). After studying materials that were found to be dissimilar to the “traditional” Marsupella that were collected in the last decade by the authors of this article, it was found that these plants belong to three new-for-science species, and all of these species should be included in Marsupella sect. Stolonicaulon. The newly described species have expanded the boundaries of morphological variability, not only for the section itself, but also for the genus based on two findings: (1) the leaves of Marsupella sect. Stolonicaulon can be either appressed and entire or spaced and deeply divided (thus, the plants could occasionally be similar to Cephaloziella or Anastrophyllum); (2) some species of the section possess regular underleaf production. The first discovery of regular underleaves in Marsupella, as noted in two of the three newly described taxa, is the main morphological novelty described in this paper. The development of regular underleaves is a presumable relict character that brings Marsupella closer to Nardia, which was recently transferred to the Gymnomitriaceae and occupies an isolated position within its own subfamily, Nardioideae.

Genetic differentiation and phylogeography of rotifer Polyarthra dolichoptera and P. vulgaris populations between Southeastern China and eastern North America: High intercontinental differences

Genetic differentiations and phylogeographical patterns of small organisms may be shaped by spatial isolation, environmental gradients, and gene flow. However, knowledge about genetic differentiation of rotifers at the intercontinental scale is still limited. Polyarthra dolichoptera and P. vulgaris are cosmopolitan rotifers that are tolerant to environmental changes, offering an excellent model to address the research gap. Here, we investigated the populations in Southeastern China and eastern North America and evaluated the phylogeographical patterns from their geographical range sizes, geographic–genetic distance relationships and their responses to spatial‐environmental factors. Using the mitochondrial cytochrome c oxidase subunit I gene as the DNA marker, we analyzed a total of 170 individuals. Our results showed that some putative cryptic species, also known as entities were widely distributed, but most of them were limited to single areas. The divergence of P. dolichoptera and P. vulgaris indicated that gene flow between continents was limited while that within each continent was stronger. Oceanographic barriers do affect the phylogeographic pattern of rotifers in continental waters and serve to maintain genetic diversity in nature. The genetic distance of P. dolichoptera and P. vulgaris populations showed significant positive correlation with geographic distance. This might be due to the combined effects of habitat heterogeneity, long‐distance colonization, and oceanographic barriers. Furthermore, at the intercontinental scale, spatial distance had a stronger influence than environmental variables on the genetic differentiations of both populations. Wind‐ and animal‐mediated transport and even historical events of continental plate tectonics are potential factors for phylogeography of cosmopolitan rotifers.

Radiocarbon and genomic evidence for the survival of Equus Sussemionus until the late Holocene

The exceptionally rich fossil record available for the equid family has provided textbook examples of macroevolutionary changes. Horses, asses, and zebras represent three extant subgenera of Equus lineage, while the Sussemionus subgenus is another remarkable Equus lineage ranging from North America to Ethiopia in the Pleistocene. We sequenced 26 archaeological specimens from Northern China in the Holocene that could be assigned morphologically and genetically to Equus ovodovi, a species representative of Sussemionus. We present the first high-quality complete genome of the Sussemionus lineage, which was sequenced to 13.4× depth of coverage. Radiocarbon dating demonstrates that this lineage survived until ~3500 years ago, despite continued demographic collapse during the Last Glacial Maximum and the great human expansion in East Asia. We also confirmed the Equus phylogenetic tree and found that Sussemionus diverged from the ancestor of non-caballine equids ~2.3–2.7 million years ago and possibly remained affected by secondary gene flow post-divergence. We found that the small genetic diversity, rather than enhanced inbreeding, limited the species’ chances of survival. Our work adds to the growing literature illustrating how ancient DNA can inform on extinction dynamics and the long-term resilience of species surviving in cryptic population pockets.

Describing novel mitochondrial genomes of Antarctic amphipods

To date, only one mitogenome from an Antarctic amphipod has been published. Here, novel complete mitochondrial genomes (mitogenomes) of two morphospecies are assembled, namely, Charcotia amundseni and Eusirus giganteus. For the latter species, we have assembled two mitogenomes from different genetic clades of this species. The lengths of Eusirus and Charcotia mitogenomes range from 15,534 to 15,619 base pairs and their mitogenomes are composed of 13 protein coding genes, 22 transfer RNAs, 2 ribosomal RNAs, and 1 putative control region CR. Some tRNAs display aberrant structures suggesting that minimalization is also ongoing in amphipod mitogenomes. The novel mitogenomes of the two Antarctic species have features distinguishing them from other amphipod mitogenomes such as a lower AT-richness in the whole mitogenomes and a negative GC- skew in both strands of protein coding genes. The genetically most variable mitochondrial regions of amphipods are nad6 and atp8, while cox1 shows low nucleotide diversity among closely and more distantly related species. In comparison to the pancrustacean mitochondrial ground pattern, E. giganteus shows a translocation of the nad1 gene, while cytb and nad6 genes are translocated in C. amundseni. Phylogenetic analysis based on mitogenomes illustrates that Eusirus and Charcotia cluster together with other species belonging to the same amphipod superfamilies. In the absence of reference nuclear genomes, mitogenomes can be useful to develop markers for studying population genetics or evolutionary relationships at higher taxonomic levels.

Divergence of active site motifs among different classes of Populus glutaredoxins results in substrate switches

Enzymes are essential components of all biological systems. The key characteristics of proteins functioning as enzymes are their substrate specificities and catalytic efficiencies. In plants, most genes encoding enzymes are members of large gene families. Within such families, the contributions of active site motifs to the functional divergence of duplicate genes have not been well elucidated. In this study, we identified 41 glutaredoxin (GRX) genes in the Populus trichocarpa genome. GRXs are ubiquitous enzymes in plants that play important roles in developmental and stress tolerance processes. In poplar, GRX genes were divided into four classes based on clear differences in gene structure and expression pattern, subcellular localization, enzymatic activity, and substrate specificity of the encoded proteins. Using site-directed mutagenesis, this study revealed that the divergence of the active site motif among different classes of GRX proteins resulted in substrate switches and thus provided new insights into the molecular evolution of these important plant enzymes.

Exploring a blue-light-sensing transcription factor to double the peak productivity of oil in Nannochloropsis oceanica

Oleaginous microalgae can produce triacylglycerol (TAG) under stress, yet the underlying mechanism remains largely unknown. Here, we show that, in Nannochloropsis oceanica, a bZIP-family regulator NobZIP77 represses the transcription of a type-2 diacylgycerol acyltransferase encoding gene NoDGAT2B under nitrogen-repletion (N+), while nitrogen-depletion (N−) relieves such inhibition and activates NoDGAT2B expression and synthesis of TAG preferably from C16:1. Intriguingly, NobZIP77 is a sensor of blue light (BL), which reduces binding of NobZIP77 to the NoDGAT2B-promoter, unleashes NoDGAT2B and elevates TAG under N−. Under N+ and white light, NobZIP77 knockout fully preserves cell growth rate and nearly triples TAG productivity. Moreover, exposing the NobZIP77-knockout line to BL under N− can double the peak productivity of TAG. These results underscore the potential of coupling light quality to oil synthesis in feedstock or bioprocess development.

Microalgae are promising feedstock for oil production. The authors report that a transcription factor NobZIP77 can regulate oil synthesis by sensing the blue light, and explore these findings to greatly enhance oil productivity via genetic and process engineering in Nannochloropsis oceanica.

Ectopic expression of Populus MYB10 promotes secondary cell wall thickening and inhibits anthocyanin accumulation

Secondary cell wall (SCW) formation is regulated by a multilevel transcriptional regulatory network, in which MYB transcription factors (TFs) play key roles. In woody plants, hundreds of MYB TFs have been identified, most of which have unknown functions in wood SCW biosynthesis. Here, we characterized the function of a Populus MYB gene, PtoMYB10. PtoMYB10 was found to encode an R2R3-MYB TF and exhibit dominant expression in xylem tissues. PtoMYB10 was determined to be located in the nucleus with the ability to activate transcription. Overexpression of PtoMYB10 in Populus resulted in a drastic increase in SCW thickening in xylem fiber cells as well as ectopic deposition of lignin in cortex cells. The expression of genes associated with lignin biosynthesis was induced in PtoMYB10 overexpressing plants, whereas repressed gene expression was found with the anthocyanin biosynthesis pathway. Lignin and anthocyanin are both produced from metabolites of the phenylpropanoid pathway. Accordingly, the anthocyanin content of Populus overexpressing PtoMYB10 decreased by more than 68%. These results indicate that PtoMYB10 can positively regulate xylary fiber SCW thickening, accompanied by the reprogramming of phenylpropanoid metabolism, which redirects metabolic flux from anthocyanin biosynthesis to monolignol biosynthesis.

Mechanism for the genomic and functional evolution of the MIR2118 family in the grass lineage

The MIR2118 family has undergone tremendous expansion in the grass lineage, in which the miRNA targets numerous noncoding PHAS loci to produce 21-nt phased small interfering RNAs (phasiRNAs) involved in male fertility. However, the evolutionary trajectory of the grass MIR2118 genes and the functions of phasiRNAs have not yet been fully elucidated. We conducted comparative genomic, molecular evolution, expression and parallel analysis of RNA ends (PARE) analyses of MIR2118 and the miR2118-mediated regulatory pathway in grasses, focusing on Oryza sativa. In total, 617 MIR2118 and eight MIR1859 novel members were identified. Phylogenetic analyses showed that grass MIR2118 genes form a distinct clade from the MIR482/2118 genes of nongrass species. We reconstructed hypothetical evolutionary histories of the grass MIR2118 clusters and its MIR1859 variants, and examined the polycistronic composition and the differential expression of the osa-MIR2118 clusters. PARE data showed that osa-miR2118 might also direct the cleavage of some protein-coding gene transcripts. Importantly, we found that PARE analysis is inherently prone to false-positive target predictions when a large number of small RNAs, such as phasiRNAs, are analysed. Our results revealed the evolution and diversification of the MIR2118 family, and provide new insights into the functions of phasiRNAs in the grasses.

Genome-Wide Comprehensive Analysis of PtLACs: Prediction and Verification of the Functional Divergence of Tandem-Duplicated Genes

Laccases (EC 1.10.3.2) have been widely considered to participate in the metabolic processes of lignin synthesis, osmotic stress response, and flavonoid oxidation in higher plants. The research into Populus trichocarpa laccase focused on the synthesis of lignin in the past few years. In this study, for the first time, a comprehensive analysis of 53 laccase copies in the P. trichocarpa genome was conducted. Positive selection analysis using the branch-site model indicated that LAC genes in terrestrial plants have undergone selective pressure for adaptive evolution. On the basis of the phylogenetic relationship, we reconstructed the evolutionary process of terrestrial plant laccase and found that this gene family began to expand during the evolution of angiosperms. Tandem duplication is the main form of expansion of the PtLAC gene family. The analysis of the sequence characteristics, gene structure, expression pattern, and gene synonymous mutation rate of PtLACs provided a theoretical basis for the functional divergence of tandem duplicated genes. The synonymous mutation rate was used to quantify the divergence time of 11 tandem duplicated gene clusters. Cluster 2, with the earliest divergence time and lower share of sequence similarity, and cluster 5, with the latest divergence time and higher share of similarity, were selected in this study to explore the functional divergence of tandem-duplicated gene clusters. Tobacco subcellular localization and Arabidopsis transgenes verified the functional differentiation of PtLAC genes in cluster 2 and the functional non-differentiation of PtLAC genes in cluster 5. The results of this study provide a reference for the functional differentiation of tandem-duplicated PtLAC.

The Importance of Glycerophospholipid Production to the Mutualist Symbiosis of Trypanosomatids

The symbiosis in trypanosomatids is a mutualistic relationship characterized by extensive metabolic exchanges between the bacterium and the protozoan. The symbiotic bacterium can complete host essential metabolic pathways, such as those for heme, amino acid, and vitamin production. Experimental assays indicate that the symbiont acquires phospholipids from the host trypanosomatid, especially phosphatidylcholine, which is often present in bacteria that have a close association with eukaryotic cells. In this work, an in-silico study was performed to find genes involved in the glycerophospholipid (GPL) production of Symbiont Harboring Trypanosomatids (SHTs) and their respective bacteria, also extending the search for trypanosomatids that naturally do not have symbionts. Results showed that most genes for GPL synthesis are only present in the SHT. The bacterium has an exclusive sequence related to phosphatidylglycerol production and contains genes for phosphatidic acid production, which may enhance SHT phosphatidic acid production. Phylogenetic data did not indicate gene transfers from the bacterium to the SHT nucleus, proposing that enzymes participating in GPL route have eukaryotic characteristics. Taken together, our data indicate that, differently from other metabolic pathways described so far, the symbiont contributes little to the production of GPLs and acquires most of these molecules from the SHT.

HIV Protease and Integrase Empirical Substitution Models of Evolution: Protein-Specific Models Outperform Generalist Models

Diverse phylogenetic methods require a substitution model of evolution that should mimic, as accurately as possible, the real substitution process. At the protein level, empirical substitution models have traditionally been based on a large number of different proteins from particular taxonomic levels. However, these models assume that all of the proteins of a taxonomic level evolve under the same substitution patterns. We believe that this assumption is highly unrealistic and should be relaxed by considering protein-specific substitution models that account for protein-specific selection processes. In order to test this hypothesis, we inferred and evaluated four new empirical substitution models for the protease and integrase of HIV and other viruses. We found that these models more accurately fit, compared with any of the currently available empirical substitution models, the evolutionary process of these proteins. We conclude that evolutionary inferences from protein sequences are more accurate if they are based on protein-specific substitution models rather than taxonomic-specific (generalist) substitution models. We also present four new empirical substitution models of protein evolution that could be useful for phylogenetic inferences of viral protease and integrase.

Characterization of Pterotylenchus cecidogenus in Desmodium ovalifolium cover crop from oil palm plantations in central Colombia

Until recently, the stem gall nematode Pterotylenchus cecidogenus was only registered in eastern Colombia. However, the disease has recently been observed in central Colombian oil palm plantations that use Desmodium ovalifolium as a cover crop. Soil, root, stem, and leaf samples were collected from D. ovalifolium. Plants showed foliar yellowing, leaf drying, and galls within stem nodes. Nematodes were identified, and the distribution, population density, and relative importance of different genera were determined. We performed morphometric and molecular identification of nematodes associated with gall symptoms. The D2-D3 segment of the large subunit-28S of ribosomal ribonucleic acid (RNA) and internal transcribed spacer (ITS) was sequenced, and phylogenetic analysis was performed. P. cecidogenus mainly occurred in the galls and to a lesser extent in the roots and soil. Nematodes were not found in leaf or inflorescence tissue. Morphological and morphometric data confirm the presence of P. cecidogenus in the stems of D. ovalifolium with gall symptoms. This study is the first to report deoxyribonucleic acid (DNA) sequences of P. cecidogenus. Based on D2-D3 and ITS partial sequences, P. cecidogenus is a sister species of the leaf-galling nematode Ditylenchus phyllobius (Sinm. Orrina phyllobia).

DNA Barcoding unveils cryptic lineages of Hoplias malabaricus from Northeastern Brazil

The trahira or wolf fish - Hoplias malabaricus- is a valid species, although recent cytogenetic and molecular studies have indicated the existence of a species complex. In this context, the present study analyzed the mitochondrial COI marker to determine the levels of genetic diversity of specimens from the Brazilian state of Maranhão, and verify the occurrence of distinct lineages within the study area. Samples were collected from the basins of the Turiaçu, Pindaré, Mearim, Itapecuru, and Parnaíba rivers. A 630-bp fragment was obtained from 211 specimens, with 484 conserved and 108 variable sites, and 60 haplotypes (Hd = 0,947; π = 0,033). The phylogenetic analyses indicated the existence of three distinct lineages of H. malabaricus from Maranhão. Genetic distances of 1.5-8.2% were found between all the populations analyzed, while the variation between haplogroups ranged from 2.1% to 7.7%. The AMOVA indicated that most of the molecular variation was found among groups, with high FST values. The high levels of genetic variability found in the present study are supported by the available cytogenetic data. These findings reinforce the need for the development of effective programs of conservation and management independently for each river basin, in order to preserve the genetic variability found in this taxon.

Identification of Glutathione Peroxidase Gene Family in Ricinus communis and Functional Characterization of RcGPX4 in Cold Tolerance

Glutathione peroxidases (GPXs) protect cells against damage caused by reactive oxygen species (ROS) and play key roles in regulating many biological processes. Here, five GPXs were identified in the Ricinus communis genome. Phylogenetic analysis displayed that the GPXs were categorized into five groups. Conserved domain and gene structure analyses showed that the GPXs from different plant species harbored four highly similar motifs and conserved exon-intron arrangement patterns, indicating that their structure and function may have been conserved during evolution. Several abiotic stresses and hormone-responsive cis-acting elements existed in the promoters of the RcGPXs. The expression profiles indicated that the RcGPXs varied substantially, and some RcGPXs were coordinately regulated under abiotic stresses. Overexpression of RcGPX4 in Arabidopsis enhanced cold tolerance at seed germination but reduced freezing tolerance at seedlings. The expression of abscisic acid (ABA) signaling genes (AtABI4 and AtABI5), ABA catabolism genes (AtCYP707A1 and AtCYP707A2), gibberellin acid (GA) catabolism gene (AtGA2ox7), and cytokinin (CTK)-inducible gene (AtARR6) was regulated in the seeds of transgenic lines under cold stress. Overexpression of RcGPX4 can disturb the hydrogen peroxide (H2O2) homeostasis through the modulation of some antioxidant enzymes and compounds involved in the GSH-ascorbate cycle in transgenic plants. Additionally, RcGPX4 depended on the MAPK3-ICE1-C-repeat-binding factor (CBF)-COR signal transduction pathway and ABA-dependent pathway to negatively regulate the freezing tolerance of transgenic plants. This study provides valuable information for understanding the potential function of RcGPXs in regulating the abiotic stress responses of castor beans.

Population genetics of Sambar (Rusa unicolor) from the Western Himalayas: preliminary findings

BACKGROUND

The sambar (Rusa unicolor) is one of the largest deer distributed across diverse habitats of India and threatened due to habitat loss, changes in the land-use patterns, illegal poaching, and anthropogenic disturbances. Here, we report the first population genetics account of sambar deer from the Western Himalayas.

METHODS AND RESULTS

We observed relatively compromised genetic diversity (π = 0.0008 ± 0.0006 at mtDNA and Ho = 0.499 ± 0.056 at nDNA). We identified 60 unique individuals using a select panel of seven loci (P_ID sib cum 1.60E-03). Bayesian skyline plot showed a stable demographic history since the past 8 kyr with a decline in recent years. The population lacked genetic structuring, likely due to the contiguous distribution and large dispersal patterns of sambar.

CONCLUSION

The preliminary findings are valuable in exploring the utility of genetic diversity in monitoring the sambar population, estimating density following capture-recapture analysis, and aid to the conservation planning of sambar in large landscapes.

Duplication of spiralian-specific TALE genes and evolution of the blastomere specification mechanism in the bivalve lineage

Background

Despite the conserved pattern of the cell-fate map among spiralians, bivalves display several modified characteristics during their early development, including early specification of the D blastomere by the cytoplasmic content, as well as the distinctive fate of the 2d blastomere. However, it is unclear what changes in gene regulatory mechanisms led to such changes in cell specification patterns. Spiralian-TALE (SPILE) genes are a group of spiralian-specific transcription factors that play a role in specifying blastomere cell fates during early development in limpets. We hypothesised that the expansion of SPILE gene repertoires influenced the evolution of the specification pattern of blastomere cell fates.

Results

We performed a transcriptome analysis of early development in the purplish bifurcate mussel and identified 13 SPILE genes. Phylogenetic analysis of the SPILE gene in molluscs suggested that duplications of SPILE genes occurred in the bivalve lineage. We examined the expression patterns of the SPILE gene in mussels and found that some SPILE genes were expressed in quartet-specific patterns, as observed in limpets. Furthermore, we found that several SPILE genes that had undergone gene duplication were specifically expressed in the D quadrant, C and D quadrants or the 2d blastomere. These expression patterns were distinct from the expression patterns of SPILE in their limpet counterparts.

Conclusions

These results suggest that, in addition to their ancestral role in quartet specification, certain SPILE genes in mussels contribute to the specification of the C and D quadrants. We suggest that the expansion of SPILE genes in the bivalve lineage contributed to the evolution of a unique cell fate specification pattern in bivalves.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13227-021-00181-2.

Advances, challenges and opportunities of phylogenetic and social network analysis using COVID-19 data

Coronavirus disease 2019 (COVID-19) has attracted research interests from all fields. Phylogenetic and social network analyses based on connectivity between either COVID-19 patients or geographic regions and similarity between syndrome coronavirus 2 (SARS-CoV-2) sequences provide unique angles to answer public health and pharmaco-biological questions such as relationships between various SARS-CoV-2 mutants, the transmission pathways in a community and the effectiveness of prevention policies. This paper serves as a systematic review of current phylogenetic and social network analyses with applications in COVID-19 research. Challenges in current phylogenetic network analysis on SARS-CoV-2 such as unreliable inferences, sampling bias and batch effects are discussed as well as potential solutions. Social network analysis combined with epidemiology models helps to identify key transmission characteristics and measure the effectiveness of prevention and control strategies. Finally, future new directions of network analysis motivated by COVID-19 data are summarized.

Multi-tissue transcriptome analysis of two Begonia species reveals dynamic patterns of evolution in the chalcone synthase gene family

Begonia is an important horticultural plant group, as well as one of the most speciose Angiosperm genera, with over 2000 described species. Genus wide studies of genome size have shown that Begonia has a highly variable genome size, and analysis of paralog pairs has previously suggested that Begonia underwent a whole genome duplication. We address the contribution of gene duplication to the generation of diversity in Begonia using a multi-tissue RNA-seq approach. We chose to focus on chalcone synthase (CHS), a gene family having been shown to be involved in biotic and abiotic stress responses in other plant species, in particular its importance in maximising the use of variable light levels in tropical plants. We used RNA-seq to sample six tissues across two closely related but ecologically and morphologically divergent species, Begonia conchifolia and B. plebeja, yielding 17,012 and 19,969 annotated unigenes respectively. We identified the chalcone synthase gene family members in our Begonia study species, as well as in Hillebrandia sandwicensis, the monotypic sister genus to Begonia, Cucumis sativus, Arabidopsis thaliana, and Zea mays. Phylogenetic analysis suggested the CHS gene family has high duplicate turnover, all members of CHS identified in Begonia arising recently, after the divergence of Begonia and Cucumis. Expression profiles were similar within orthologous pairs, but we saw high inter-ortholog expression variation. Sequence analysis showed relaxed selective constraints on some ortholog pairs, with substitutions at conserved sites. Evidence of pseudogenisation and species specific duplication indicate that lineage specific differences are already beginning to accumulate since the divergence of our study species. We conclude that there is evidence for a role of gene duplication in generating diversity through sequence and expression divergence in Begonia.