Orthologous Sequences Research Articles

Testing Phylogenetic Placement Accuracy of DNA Barcode Sequences on a Fish Backbone Tree: Implications of Backbone Tree Completeness and Species Representation.

Advancements in DNA sequencing technology have facilitated the generation of a vast number of DNA sequences, posing opportunities and challenges for constructing large phylogenetic trees. DNA barcode sequences, particularly COI, represent extensive orthologous sequences suitable for phylogenetic analysis. Phylogenetic placement analysis offers a promising method to integrate COI data into tree-building efforts, yet the impacts of backbone tree completeness and species composition remain under-explored. Using a dataset comprising 27 genes and 4520 species of bony fishes, we assessed the accuracy of phylogenetic inference by "placing" COI sequences onto backbone trees. The backbone tree completeness was varied by subsampling 20%, 40%, 60%, 80%, and 99% of the total species separately, followed by placement of those missing species based on their COI sequences using software packages EPA-ng and APPLES. We also compared the effects of biased, random, and stratified sampling strategies; the latter ensured the representation of all major lineages (Family) of bony fish. Our findings indicate that the placement accuracy is consistently high across all levels of backbone tree completeness, where 70%-78% missing species are correctly placed (by EPA-ng) in the same locations as the reference tree derived from the complete data. High completeness produces slightly high placement accuracy, although in many cases the differences are nonsignificant. For example, at the 99% completeness level with stratified sampling, EPA-ng placed 78% missing species correctly, and when only considering placement with high confidence (LWR > 0.9), the percentage is 87%. Additionally, stratified sampling outperforms random sampling in most cases, and biased sampling has the worst performance. The likelihood-based EPA-ng consistently provide higher accurate placements than the distance-based APPLES. In conclusion, COI-based placement analysis represents a potential route of using the available vast barcoding data for building large phylogenetic trees.

The amniote-conserved DNA-binding domain of CGGBP1 restricts cytosine methylation of transcription factor binding sites in proximal promoters to regulate gene expression

CGGBP1 is a GC-rich DNA-binding protein which is important for genomic integrity, gene expression and epigenome maintenance through regulation of CTCF occupancy and cytosine methylation. It has remained unclear how CGGBP1 integrates multiple diverse functions with its simple architecture of only a DNA-binding domain tethered to a C-terminal tail with low structural rigidity. We have used truncated forms of CGGBP1 with or without the DNA-binding domain (DBD) to assay cytosine methylation and global gene expression. Proximal promoters of CGGBP1-repressed genes, although significantly GC-poor, contain GC-rich transcription factor binding motifs and exhibit base compositions indicative of low C-T transition rates due to prevention of cytosine methylation. Genome-wide analyses of cytosine methylation and binding of CGGBP1 DBD show that CGGBP1 restricts cytosine methylation in a manner that depends on its DBD and its DNA-binding. The CGGBP1-repressed genes show an increase in promoter cytosine methylation alongside a decrease in transcript abundance when the DBD-deficient CGGBP1 is expressed. Our findings suggest that CGGBP1 protects transcription factor binding sites (TFBS) from cytosine methylation-associated loss and thereby regulates gene expression. By analysing orthologous promoter sequences we show that restriction of cytosine methylation is a function of CGGBP1 progressively acquired during vertebrate evolution. A superimposition of our results and evolution of CGGBP1 suggests that mitigation of cytosine methylation is majorly achieved by its N-terminal DBD. Our results position CGGBP1 DNA-binding as a major evolutionarily acquired mechanism through which it keeps cytosine methylation under check and regulates TFBS retention and gene activity.

Molecular evolution and genetic diversity of defective chorion 1 in Anastrepha fraterculus and Anastrepha obliqua (Diptera, Tephritidae).

The family Tephritidae comprises numerous fruit fly species, some of which are economically significant, such as several in the genus Anastrepha. Most pest species in this genus belong to the fraterculus group, characterized by closely related species that are difficult to differentiate due to recent divergence and gene flow. Identifying genetic markers for their study is paramount for understanding the group's evolution and eventual phytosanitary control. Because there is variation in eggshell morphology among species in the genus, the study of the rapidly evolving defective chorion 1 (dec-1) gene, which is crucial for chorion formation and reproduction, could provide relevant information for Anastrepha differentiation. We compared transcriptome sequences of dec-1 from two of the most important pest species in the genus, Anastrepha fraterculus and Anastrepha obliqua to dec-1 sequences from Anastrepha ludens, which was used for structure prediction. Furthermore, we amplified a conserved exon across populations of these species. These data revealed three alternative transcripts in A. fraterculus and A. obliqua, consistent with patterns found in other Tephritidae; we obtained orthologous sequences for these other tephritids from NCBI to investigate patterns of selection affecting this gene at different hierarchical levels using different methods. These analyses show a general pattern of purifying selection across the whole gene and throughout its history at different hierarchical levels, from populations to more distantly related species. That notwithstanding, we still found evidence of positive and episodic diversifying selection at different levels. Different parts of the gene have shown distinct evolutionary rates, which were associated with the diverse proproteins produced by posttranslational changes of DEC-1, with proproteins that are incorporated in the chorion earlier in egg formation being in general more conserved than others that are incorporated later. This correlation appears more evident in certain lineages, including the branch that separates Anastrepha, as well as other internal branches that differentiate species within the genus. Our data showed that this gene shows remarkable variation across its different exons, which has proven to be informative at different evolutionary levels. These changes hold promise not only for studying differentiation in Anastrepha but also for the eventual management of selected pest species.

Genomic characterization and comparative genomics of Chlorella sp. CH2018 from Musi River water, India

Microalgae, a diverse class of photosynthetic eukaryote, can provide food and energy sustainably. Selecting productive strains is crucial for commercial viability. Here, we isolated an axenic microalgal species from Musi River water, and based on its morphology, molecular makeup, and cell wall composition identified it as Chlorella sp. CH2018. Scanning electron microscopy was used to measure the size of cells, which were found to be 3-, 4-, or 5-μm in diameter with a discernible thick outer cell wall. The characterization of algal genomes is imperative for comprehending species, studying metabolic pathways, and modifying genetics. In this study, we sequenced the entire genome of Chlorella sp. CH2018, yielded genome size of 56.83 Mb with 11,143 functionally annotated protein-coding gene models. The Gene Ontology (GO) analysis revealed that Chlorella's predominant metabolic pathway is carbohydrate metabolism. Ortholog comparative analysis of species of phylum Chlorophyta with Chlorella sp. CH2018 showed that the isolated species possesses unique protein families with a maximum number of 6292 ortholog groups with Chlorella sorokiniana. The phylogenetic tree created by concatenating single-copy ortholog sequences demonstrates the uniqueness of Chlorella sp. CH2018, and its genome sequence serves as a genetic resource for future research.

A new wave of Mesoamerican bumblebees? Revising the weisi-complex to reject numts and pseudospecies (Apidae: Bombus).

COI-barcode-like sequences appear to show substantially more species diversity among Mesoamerican bumblebees than had been reported previously from morphological studies. Closer examination shows that some of this apparent diversity may be pseudospecies (groups falsely misinterpreted as separate species), often supported by paralogous 'numts' (nuclear copies of mitochondrial sequences). For the well-sampled weisi-complex, we seek to filter out pseudogenes in order to use the orthologous COI-barcode sequences for identifying estimates of evolutionary relationships and likely species' gene coalescents for candidate species. Even after this filtering, in contrast to recent purely morphological studies our results from an integrative assessment of species' gene coalescents together with skeletal morphology support that 'Bombus weisi' Friese in its recent broad sense consists of two species: B. weisi (which includes the taxon montezumae Cockerell); and B. nigrodorsalis Franklin. Our interpretation rejects likely numts-based pseudospecies and a candidate species that are unsupported by skeletal morphology. This shows that careful attention needs to be paid to both barcode analysis and to skeletal morphology, to avoid describing pseudospecies.

High-throughput amino acid-level characterization of the interactions of plasminogen activator inhibitor-1 with variably divergent proteases.

While members of large paralogous protein families share structural features, their functional niches often diverge significantly. Serine protease inhibitors (SERPINs), whose members typically function as covalent inhibitors of serine proteases, are one such family. Plasminogen activator inhibitor-1 (PAI-1) is a prototypic SERPIN, which canonically inhibits tissue-and urokinase-type plasminogen activators (tPA and uPA) to regulate fibrinolysis. PAI-1 has been shown to also inhibit other serine proteases, including coagulation factor XIIa (FXIIa) and transmembrane serine protease 2 (TMPRSS2). The structural determinants of PAI-1 inhibitory function toward these non-canonical protease targets, and the biological significance of these functions, are unknown. We applied deep mutational scanning (DMS) to assess the effects of ∼80% of all possible single amino acid substitutions in PAI-1 on its ability to inhibit three putative serine protease targets (uPA, FXIIa, and TMPRSS2). Selection with each target protease generated a unique PAI-1 mutational landscape, with the determinants of protease specificity distributed throughout PAI-1's primary sequence. Next, we conducted a comparative analysis of extant orthologous sequences, demonstrating that key residues modulating PAI-1 inhibition of uPA and FXIIa, but not TMPRSS2, are maintained by purifying selection. PAI-1's activity toward FXIIa may reflect how protease evolutionary relationships predict SERPIN functional divergence, which we support via a cophylogenetic analysis of all secreted SERPINs and their cognate serine proteases. This work provides insight into the functional diversification of SERPINs and lays the framework for extending these studies to other proteases and their regulators.

The Human Accelerated Region HAR202 Controls NPAS3 Expression in the Developing Forebrain Displaying Differential Enhancer Activity Between Modern and Archaic Human Sequences.

It has been proposed that the phenotypic differences in cognitive abilities between humans and our closest living relatives, chimpanzees, are largely due to changes in the regulation of neurodevelopmental genes. We have previously found that the neurodevelopmental transcription factor gene NPAS3 accumulates the largest number of human accelerated regions (HARs), suggesting it may play some role in the phenotypic evolution of the human nervous system. In this work, we performed a comparative functional analysis of NPAS3-HAR202 using enhancer reporter assays in transgenic zebrafish and mice. We found that the Homo sapiens HAR202 ortholog failed to drive reporter expression to the zebrafish nervous system, in high contrast to the strong expression displayed by the rest of the vertebrate ortholog sequences tested. Remarkably, the HAR202 ortholog from archaic humans (Neanderthals/Denisovans) also displayed a pan-vertebrate expression pattern, despite the fact that archaic and modern humans have only one nucleotide substitution. Moreover, similar results were found when comparing enhancer activity in transgenic mice, where we observed a loss of activity of the modern human version in the mouse developing brain. To investigate the functional importance of HAR202, we generated mice lacking HAR202 and found a remarkable decrease of Npas3 expression in the forebrain during development. Our results place HAR202 as one of the very few examples of a neurodevelopmental transcriptional enhancer displaying functional evolution in the brain as a result of a fast molecular evolutionary process that specifically occurred in the human lineage.

Antennal Transcriptome Screening and Identification of Chemosensory Proteins in the Double-Spine European Spruce Bark Beetle, Ips duplicatus (Coleoptera: Scolytinae).

The northern bark beetle, Ips duplicatus, is an emerging economic pest, reportedly infesting various species of spruce (Picea spp.), pine (Pinus spp.), and larch (Larix spp.) in Central Europe. Recent climate changes and inconsistent forest management practices have led to the rapid spread of this species, leaving the current monitoring strategies inefficient. As understanding the molecular components of pheromone detection is key to developing novel control strategies, we generated antennal transcriptomes from males and females of this species and annotated the chemosensory proteins. We identified putative candidates for 69 odorant receptors (ORs), 50 ionotropic receptors (IRs), 25 gustatory receptors (GRs), 27 odorant-binding proteins (OBPs), including a tetramer-OBP, 9 chemosensory proteins (CSPs), and 6 sensory neuron membrane proteins (SNMPs). However, no sex-specific chemosensory genes were detected. The phylogenetic analysis revealed conserved orthology in bark beetle chemosensory proteins, especially with a major forest pest and co-habitant, Ips typographus. Recent large-scale functional studies in I. typographus chemoreceptors add greater significance to the orthologous sequences reported here. Nevertheless, identifying chemosensory genes in I. duplicatus is valuable to understanding the chemosensory system and its evolution in bark beetles (Coleoptera) and, generally, insects.

Comparative genome analysis and the genome-shaping role of long terminal repeat retrotransposons in the evolutionary divergence of fungal pathogens Blastomyces dermatitidis and Blastomyces gilchristii.

Blastomyces dermatitidis and Blastomyces gilchristii are cryptic species of fungi that cause blastomycosis, an often severe disease involving pulmonary infection capable of systemic dissemination. While these species appear morphologically identical, differences exist in the genetic makeup, geographical range, and possibly the clinical presentation of infection. Here, we show genetic divergence between the cryptic species through both a Blastomyces species tree constructed from orthologous protein sequences and whole genome single-nucleotide variant phylogenomic analysis. Following linked-read sequencing and de novo genome assembly, we characterized and compared the genomes of 3 B. dermatitidis and 3 B. gilchristii isolates. The B. gilchristii genomes (73.25-75.4 Mb) were ∼8 Mb larger than the B. dermatitidis genomes (64.88-66.61 Mb). Average nucleotide identity was lower between genomes of different species than genomes of the same species, yet functional classification of genes suggested similar proteomes. The most striking difference involved long terminal repeat retrotransposons. Although the same retrotransposon elements were detected in the genomes, the quantity of elements differed between the 2 species. Gypsy retrotransposon content was significantly higher in B. gilchristii (38.04-39.26 Mb) than in B. dermatitidis (30.85-32.40 Mb), accounting for the majority of genome size difference between species. Age estimation and phylogenetic analysis of the reverse transcriptase domains suggested that these retrotransposons are relatively ancient, with genome insertion predating the speciation of B. dermatitidis and B. gilchristii. We postulate that different trajectories of genome contraction led to genetic incompatibility, reproductive isolation, and speciation, highlighting the role of transposable elements in fungal evolution.

Evolution of neuropeptide Y/RFamide-like receptors in nematodes

The Neuropeptide Y/RFamide-like receptors belong to the Rhodopsin-like G protein-coupled receptors G protein-coupled receptors (GPCRs) and are involved in functions such as locomotion, feeding and reproduction. With 41 described receptors they form the best-studied group of neuropeptide GPCRs in Caenorhabditis elegans. In order to understand the expansion of the Neuropeptide Y/RFamide-like receptor family in nematodes, we started from the sequences of selected receptor paralogs in C. elegans as query and surveyed the corresponding orthologous sequences in another 159 representative nematode target genomes. To this end we employed a automated pipeline based on ExonMatchSolver, a tool that solves the paralog-to-contig assignment problem. Utilizing subclass-specific HMMs we were able to detect a total of 1557 Neuropeptide Y/RFamide-like receptor sequences (1100 NPRs, 375 FRPRs and 82 C09F12.3) in the 159 target nematode genomes investigated here. These sequences demonstrate a good conservation of the Neuropeptide Y/RFamide-like receptors across the Nematoda and highlight the diversification of the family in nematode evolution. No other genus shares all Neuropeptide Y/RFamide-like receptors with the genus Caenorhabditis. At the same time, we observe large numbers of clade specific duplications and losses of family members across the phylum Nematoda.

COATi: Statistical Pairwise Alignment of Protein-Coding Sequences.

Sequence alignment is an essential method in bioinformatics and the basis of many analyses, including phylogenetic inference, ancestral sequence reconstruction, and gene annotation. Sequencing artifacts and errors made during genome assembly, such as abiological frameshifts and incorrect early stop codons, can impact downstream analyses leading to erroneous conclusions in comparative and functional genomic studies. More significantly, while indels can occur both within and between codons in natural sequences, most amino-acid- and codon-based aligners assume that indels only occur between codons. This mismatch between biology and alignment algorithms produces suboptimal alignments and errors in downstream analyses. To address these issues, we present COATi, a statistical, codon-aware pairwise aligner that supports complex insertion-deletion models and can handle artifacts present in genomic data. COATi allows users to reduce the amount of discarded data while generating more accurate sequence alignments. COATi can infer indels both within and between codons, leading to improved sequence alignments. We applied COATi to a dataset containing orthologous protein-coding sequences from humans and gorillas and conclude that 41% of indels occurred between codons, agreeing with previous work in other species. We also applied COATi to semiempirical benchmark alignments and find that it outperforms several popular alignment programs on several measures of alignment quality and accuracy.

Improved selection of canonical proteins for reference proteomes.

The 'canonical' protein sets distributed by UniProt are widely used for similarity searching, and functional and structural annotation. For many investigators, canonical sequences are the only version of a protein examined. However, higher eukaryotes often encode multiple isoforms of a protein from a single gene. For unreviewed (UniProtKB/TrEMBL) protein sequences, the longest sequence in a Gene-Centric group is chosen as canonical. This choice can create inconsistencies, selecting >95% identical orthologs with dramatically different lengths, which is biologically unlikely. We describe the ortho2tree pipeline, which examines Reference Proteome canonical and isoform sequences from sets of orthologous proteins, builds multiple alignments, constructs gap-distance trees, and identifies low-cost clades of isoforms with similar lengths. After examining140 000 proteins from eight mammals in UniProtKB release 2022_05, ortho2tree proposed 7804 canonical changes for release 2023_01, while confirming 53 434 canonicals. Gap distributions for isoforms selected by ortho2tree are similar to those in bacterial and yeast alignments, organisms unaffected by isoform selection, suggesting ortho2tree canonicals more accurately reflect genuine biological variation. 82% of ortho2tree proposed-changes agreed with MANE; for confirmed canonicals, 92% agreed with MANE. Ortho2tree can improve canonical assignment among orthologous sequences that are >60% identical, a group that includes vertebrates and higher plants.

Refining taxonomic identification of microalgae through molecular and genetic evolution: a case study of Prorocentrum lima and Prorocentrum arenarium.

Species delimitation based on lineage definition has become increasingly popular. However, these methods have been limited, especially for species that lack genomic data and are morphologically similar. The trickiest part for the species identification is that the interspecific and intraspecific boundaries are vague. Taking Prorocentrum (Dinophyta) as an example, analysis of cell morphology, growth, and toxin synthesis in both species of P. lima and P. arenarium does not provide a reliable basis for species delineation. However, through phylogenetic and genetic distance analyses of their ITS and LSU sequences, establishment of evolutionary tree based on orthologous gene sequences, and combining the results of automatic barcode gap discovery and Poisson tree processes models, it was sustained that P. arenarium does not belong to the P. lima complex and should be considered as an independent species. Interspecies genetic evolution analysis revealed that P. lima and P. arenarium may contribute to evolutionary direction that favors combating reverse environmental factors. In P. lima, viral invasion may be one of the reasons for its large genome size. In the study, P. lima complex has been selected as an example to enhance the taxonomic identification of microalgae through molecular and genetic evolution, offering valuable insights into refining taxonomic identification and promoting microbial biodiversity research in other species.IMPORTANCEMicroalgae, especially the species known as Prorocentrum, have received significant attention due to their ability to trigger harmful algal blooms and produce toxins. However, the boundaries between species and within species are ambiguous. Clear and comprehensive species delineation indicates that Prorocentrum arenarium should be considered as an independent species, separate from the Prorocentrum lima complex. Improving the classification and identification of microalgae through molecular and genetic evolution will provide reference points for other cryptic species. Prorocentrum occupy multiple ecological niches in marine environments, and studying their evolutionary direction contributes to understanding their ecological adaptations and community succession.

Cloning, characterization, and evolutionary patterns of KCNQ4 genes in anurans.

Acoustic communication plays important roles in the survival and reproduction of anurans. The perception and discrimination of conspecific sound signals of anurans were always affected by masking background noise. Previous studies suggested that some frogs evolved the high-frequency hearing to minimize the low-frequency noise. However, the molecular mechanisms underlying the high-frequency hearing in anurans have not been well explored. Here, we cloned and obtained the coding regions of a high-frequency hearing-related gene (KCNQ4) from 11 representative anuran species and compared them with orthologous sequences from other four anurans. The sequence characteristics and evolutionary analyses suggested the highly conservation of the KCNQ4 gene in anurans, which supported their functional importance. Branch-specific analysis showed that KCNQ4 genes were under different evolutionary forces in anurans and most anuran lineages showed a generally strong purifying selection. Intriguingly, one significantly positively selected site was identified in the anuran KCNQ4 gene based on FEL model. Positive selection was also found along the common ancestor of Ranidae and Rhacophoridae as well as the ancestral O. tianmuii based on the branch-site analysis, and the positively selected sites identified were involved in or near the N-terminal ion transport and the potassium ion channel functional domain of the KCNQ4 genes. The present study revealed valuable information regarding the KCNQ4 genes in anurans and provided some new insights for the underpinnings of the high-frequency hearing in frogs.

Comparative Genomics Reveal Phylogenetic Relationship and Chromosomal Evolutionary Events of Eight Cervidae Species.

Cervidae represents a family that is not only rich in species diversity but also exhibits a wide range of karyotypes. The controversies regarding the phylogeny and classification of Cervidae still persist. The flourishing development of the genomic era has made it possible to address these issues at the genomic level. Here, the genomes of nine species were used to explore the phylogeny and chromosomal evolutionary events of Cervidae. By conducting whole-genome comparisons, we identified single-copy orthologous genes across the nine species and constructed a phylogenetic tree based on the single-copy orthologous genes sequences, providing new insights into the phylogeny of Cervidae, particularly the phylogenetic relationship among sika deer, red deer, wapiti and Tarim red deer. Gene family analysis revealed contractions in the olfactory receptor gene family and expansions in the histone gene family across eight Cervidae species. Furthermore, synteny analysis was used to explore the chromosomal evolutionary events of Cervidae species, revealing six chromosomal fissions during the evolutionary process from Bovidae to Cervidae. Notably, specific chromosomal fusion events were found in four species of Cervus, and a unique chromosomal fusion event was identified in Muntiacus reevesi. Our study further completed the phylogenetic relationship within the Cervidae and demonstrated the feasibility of inferring species phylogeny at the whole-genome level. Additionally, our findings on gene family evolution and the chromosomal evolutionary events in eight Cervidae species lay a foundation for comprehensive research of the evolution of Cervidae.

Unraveling the genetics of heat tolerance in chickpea landraces (Cicer arietinum L.) using genome-wide association studies.

Chickpea, being an important grain legume crop, is often confronted with the adverse effects of high temperatures at the reproductive stage of crop growth, drastically affecting yield and overall productivity. The current study deals with an extensive evaluation of chickpea genotypes, focusing on the traits associated with yield and their response to heat stress. Notably, we observed significant variations for these traits under both normal and high-temperature conditions, forming a robust basis for genetic research and breeding initiatives. Furthermore, the study revealed that yield-related traits exhibited high heritability, suggesting their potential suitability for marker-assisted selection. We carried out single-nucleotide polymorphism (SNP) genotyping using the genotyping-by-sequencing (GBS) method for a genome-wide association study (GWAS). Overall, 27 marker-trait associations (MTAs) linked to yield-related traits, among which we identified five common MTAs displaying pleiotropic effects after applying a stringent Bonferroni-corrected p-value threshold of <0.05 [-log10(p) > 4.95] using the BLINK (Bayesian-information and linkage-disequilibrium iteratively nested keyway) model. Through an in-depth in silico analysis of these markers against the CDC Frontier v1 reference genome, we discovered that the majority of the SNPs were located at or in proximity to gene-coding regions. We further explored candidate genes situated near these MTAs, shedding light on the molecular mechanisms governing heat stress tolerance and yield enhancement in chickpeas such as indole-3-acetic acid-amido synthetase GH3.1 with GH3 auxin-responsive promoter and pentatricopeptide repeat-containing protein, etc. The harvest index (HI) trait was associated with marker Ca3:37444451 encoding aspartic proteinase ortholog sequence of Oryza sativa subsp. japonica and Medicago truncatula, which is known for contributing to heat stress tolerance. These identified MTAs and associated candidate genes may serve as valuable assets for breeding programs dedicated to tailoring chickpea varieties resilient to heat stress and climate change.

Secondary metabolite gene clusters from the phytopathogenic fungus Gaeumannomyces tritici

The take-all disease is one of the most important maladies in cereals and grasses, being caused by the fungus Gaeumannomyces tritici. Secondary metabolites are known to perform critical functions during the infection process of various phytopathogens. However, the current understanding of the biosynthesis of secondary metabolites in G. tritici is limited. Similarly, comprehensive analyses of the expression, conservation, and evolution of these biosynthesis-related genes are crucial for enhancing our knowledge of the molecular mechanisms that drive the development of the take-all disease. Here we have performed a deep survey and description of secondary metabolite biosynthetic gene clusters in G. tritici, analyzed a previously published RNA-seq of a mimicked infection condition, and assessed the conservation among 10 different Magnaporthales order members. Notably, the majority of the 35 putative gene clusters identified were conserved among these species, with GtPKS1, GtPKS3, and GtTERP4 uniquely identified in G. tritici. In the mimicked infection condition, seven gene clusters, including the GtPKS1 cluster, exhibited upregulated expression. Through comparative genomic analysis, GtPKS1 was associated with the production of dichlorodiaporthin, a metabolite with cytotoxic and antifungal activity. In addition, GtPKS10 and GtPKSNRPS3 showed similarities to already characterized biosynthetic pathways involved in the synthesis of ACR-toxin (phytotoxic) and trichosetin (phytotoxic and antibiotic), respectively. These three gene clusters were further scrutinized through phylogenetic inference, which revealed the distribution of orthologous sequences across various plant-associated fungi. Finally, the detailed identification of several genes enrolled in secondary metabolite biosynthesis provides the foundation for future in-depth research, supporting the potential impact of several small molecules on G. tritici lifecycle and host interactions.

Evolutionary Analysis of Cnidaria Small Cysteine-Rich Proteins (SCRiPs), an Enigmatic Neurotoxin Family from Stony Corals and Sea Anemones (Anthozoa: Hexacorallia).

Cnidarians (corals, sea anemones, and jellyfish) produce toxins that play central roles in key ecological processes, including predation, defense, and competition, being the oldest extant venomous animal lineage. Cnidaria small cysteine-rich proteins (SCRiPs) were the first family of neurotoxins detected in stony corals, one of the ocean's most crucial foundation species. Yet, their molecular evolution remains poorly understood. Moreover, the lack of a clear classification system has hindered the establishment of an accurate and phylogenetically informed nomenclature. In this study, we extensively surveyed 117 genomes and 103 transcriptomes of cnidarians to identify orthologous SCRiP gene sequences. We annotated a total of 168 novel putative SCRiPs from over 36 species of stony corals and 12 species of sea anemones. Phylogenetic reconstruction identified four distinct SCRiP subfamilies, according to strict discrimination criteria based on well-supported monophyly with a high percentage of nucleotide and amino acids' identity. Although there is a high prevalence of purifying selection for most SCRiP subfamilies, with few positively selected sites detected, a subset of Acroporidae sequences is influenced by diversifying positive selection, suggesting potential neofunctionalizations related to the fine-tuning of toxin potency. We propose a new nomenclature classification system relying on the phylogenetic distribution and evolution of SCRiPs across Anthozoa, which will further assist future proteomic and functional research efforts.

Cross-species immunoprotective antigens (subolesin, ferritin 2 and P0) provide protection against Rhipicephalus sanguineus sensu lato

BackgroundTick control is mostly hampered by the rise of acaricide-resistant tick populations. Significant efforts have focused on developing alternative control methods, including cross-species protective and/or cocktail-based anti-tick vaccines, to achieve protection against various tick species.MethodsIn this study, full-length open reading frames encoding subolesin (SUB) from Rhipicephalus microplus and ferritin 2 (FER2) from Hyalomma anatolicum as well as the partial 60S acidic ribosomal protein (P0) from R. microplus were cloned, expressed in Escherichia coli and used as vaccine antigens against Rhipicephalus sanguineus sensu lato (R. sanguineus s.l.) infestation in rabbits.ResultsIn silico analyses revealed that the SUB, P0 and FER2 proteins were antigenic and displayed limited similarity to the host's homologous proteins. The proteins shared identities of 97.5%, 100% and 89.5% with their SUB, P0 and FER2 R. sanguineus s.l. orthologous sequences, respectively. Antibodies against each recombinant protein cross-recognized the native proteins in the different tissues and developmental stages of R. sanguineus s.l. Overall efficacy of the SUB, FER2 and cocktail (SUB+FER2+P0) vaccines against R. sanguineus s.l. infestation was 86.3%, 95.9% and 90.9%, respectively.ConclusionsBoth mono-antigen and the cocktail anti-tick vaccines affected the biological parameters of R. sanguineus s.l. infestation in the rabbit model, which could be extrapolated to its infested host under natural conditions. These findings support the possibility of using mono-antigenic and cocktail-based vaccines for large-scale anti-tick vaccine development against multiple tick species.Graphical

Mitogenomic variation in the Black-throated Tit (Aegithalos concinnus): Conserved structure, concerted evolution of duplicate control regions and multiple distinct evolutionary lineages

The mitochondrial genome is a prominent research topic due to its indispensable role in organisms and its application in many research disciplines. However, few studies have investigated intraspecies mitogenomic variation. In this study, 69 mitogenomes of the Black-throated Tit (Aegithalos concinnus) were assembled and annotated from a large number of short reads generated using high-throughput sequencing technology. Comparative analyses revealed that mitogenomic characteristics such as length, gene and nucleotide composition, codon usage, and duplicated control regions were relatively conserved despite substantial intraspecies morphological changes. Yet, all the individuals from the subspecies A. c. iredalei had one more nucleotide in the 12S rRNA than the other studied subspecies. Phylogenetic analyses showed five distinct lineages based on the complete mitogenomes and the 13 combined protein-coding genes, whereas only four lineages were observed when using the duplicate control regions. Most interestingly, each lineage had both copies of the control regions of the comprising individuals, indicating that the paralogous control regions were more similar than the orthologous sequences from the distinct lineages. This suggested the control regions had undergone concerted evolution. The Black-throated Tit has complex evolutionary history and needs further investigating the taxonomic status of these lineages, as well as the underlying evolutionary processes. Our findings call for more research on intraspecies mitogenomic variation.