Exploring the evolutionary dynamics of myrmecophytism: Perspectives from the Southeast Asian Macaranga ant-plant symbiosis

Myrmecophytic plants utilise defensive services offered by obligate ant partners nesting in their domatia in a novel means of survival in tropical habitats. Although much is known about the ecology of myrmecophytism, there aren't enough empirical examples to demonstrate whether it substantially influences evolutionary patterns in host plant lineages. In this study, we make use of the species-rich Macaranga (Euphorbiaceae) ant-plant symbiosis distributed in the Southeast Asian Sundaland to delve into the evolutionary dynamics of myrmecophytism in host plants. We generated the most comprehensive dated phylogeny of myrmecophytic Macaranga till date using genotyping-by-sequencing (GBS). With this in hand, we traced the evolutionary history of myrmecophytism in Macaranga using parametric biogeography and ancestral state reconstruction. Diversification rate analysis methods were employed to determine if myrmecophytism enhanced diversification rates in the genus. Our results demonstrate that myrmecophytism is labile and easily lost. Ancestral state reconstruction supported a single origin of myrmecophytism in Macaranga ∼18 mya on Borneo followed by multiple losses. Diversification rate analysis methods did not yield sufficient evidence to support the hypothesis that myrmecophytism enhanced diversification rates in Macaranga; we found that topographical features on Borneo may have played a more direct role in the divergence of clades instead. Our study provides evidence that while the acquisition of domatia clearly functions as a key innovation that has enabled host plants to exploit the environment in novel ways, it may not necessarily enhance diversification rates. In fact, we hypothesise that overly specialised cases of myrmecophytism may even be an evolutionary dead end.

Net-zero energy transition in ASEAN countries: The evolutionary model brings novel perspectives to the cooperative mechanism of climate governance

In an era marked by escalating climate change, the fragile ecological balance faces increasing strain. Whilst significant knowledge exists regarding the accumulation of carbon emission within the Association of Southeast Asian Nations, little is known about when and how countries could reach net-zero emission goal as agreed in Paris Agreement. For this purpose, our study examines the primary driving factors of carbon emission from 1990 to 2020 across the ten countries using the Logarithmic Mean Divisia Index model. We leverage the random forest model to explore the net-zero scenarios and the Autoregressive Integrated Moving Average approach to identify the evolutionary trajectories of carbon emission trends. Our findings underscore the imperative need for expediting decarbonization efforts, emphasizing the urgency for widespread adoption of clean technologies and substantial investment in green initiatives. Countries at similar stages of progress might establish a cooperation mechanism of clean energy base construction, energy storage allocation and policy formulation. These insights can help us better estimate future demand of clean energy, explore strategies for decarbonization, and inform historical commonalities of carbon emission growth.

A novel framework for investigating the mechanisms of climate change and anthropogenic activities on the evolution of hydrological drought

Climate change and anthropogenic activity are the primary drivers of water cycle changes. Hydrological droughts are caused by a shortage of surface and/or groundwater resources caused by climate change and/or anthropogenic activity. Existing hydrological models have primarily focused on simulating natural water cycle processes, while limited research has investigated the influence of anthropogenic activities on water cycle processes. This study proposes a novel framework that integrates a distributed hydrological model and an attribution analysis method to assess the impacts of climate change and anthropogenic activities on hydrological drought The distributed dualistic water cycle model was applied to the Fuhe River Basin (FRB), and it generated a Nash-Sutcliffe efficiency coefficient > 0.85 with a relative error of <5 %. Excluding the year with extreme drought conditions, our analysis revealed that climate change negatively impacted the average drought duration (-105.5 %) and intensity (-23.6 %) because of increasing precipitation. However, anthropogenic activities continued to contribute positively to the drought, accounting for 5.5 % and 123.6 % of the average drought duration and intensity, respectively, because of increased water consumption. When accounting for extreme drought years, our results suggested that climate change has contributed negatively to the average duration of drought (-113.2 %) but positively to its intensity (7.8 %). Further, we found that anthropogenic activities contributed positively to both the average drought duration and intensity (13.2 % and 92.2 %, respectively). While climate change can potentially mitigate hydrological drought in the FRB by boosting precipitation levels, its overall effect may exacerbate drought through the amplification of extreme climate events resulting from global climate change. Therefore, greater attention should be paid to the effects of extreme drought.

Unraveling the Drivers of Tumorigenesis in the Context of Evolution: Theoretical Models and Bioinformatics Tools

Cancer originates from somatic cells that have accumulated mutations. These mutations alter the phenotype of the cells, allowing them to escape homeostatic regulation that maintains normal cell numbers. The emergence of malignancies is an evolutionary process in which the random accumulation of somatic mutations and sequential selection of dominant clones cause cancer cells to proliferate. The development of technologies such as high-throughput sequencing has provided a powerful means to measure subclonal evolutionary dynamics across space and time. Here, we review the patterns that may be observed in cancer evolution and the methods available for quantifying the evolutionary dynamics of cancer. An improved understanding of the evolutionary trajectories of cancer will enable us to explore the molecular mechanism of tumorigenesis and to design tailored treatment strategies.

Gene fusion as an important mechanism to generate new genes in the genus Oryza

Background

Events of gene fusion have been reported in several organisms. However, the general role of gene fusion as part of new gene origination remains unknown.

Results

We conduct genome-wide interrogations of four Oryza genomes by designing and implementing novel pipelines to detect fusion genes. Based on the phylogeny of ten plant species, we detect 310 fusion genes across four Oryza species. The estimated rate of origination of fusion genes in the Oryza genus is as high as 63 fusion genes per species per million years, which is fixed at 16 fusion genes per species per million years and much higher than that in flies. By RNA sequencing analysis, we find more than 44% of the fusion genes are expressed and 90% of gene pairs show strong signals of purifying selection. Further analysis of CRISPR/Cas9 knockout lines indicates that newly formed fusion genes regulate phenotype traits including seed germination, shoot length and root length, suggesting the functional significance of these genes.

Conclusions

We detect new fusion genes that may drive phenotype evolution in Oryza. This study provides novel insights into the genome evolution of Oryza.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13059-022-02696-w.

Genome-wide identification and characterization of UBP gene family in wheat (Triticum aestivum L.)

Ubiquitination is essential for plant growth and development. Deubiquitination cooperates with ubiquitination to regulate the ubiquitination levels of target proteins. The ubiquitin-specific protease (UBP) family is the largest group of deubiquitinases (DUBs), which perform extensive and significant roles in eukaryotic organisms. However, the UBP genes in wheat (TaUBPs) are not identified, and the functions of TaUBPs are unknown. The present study identified 97 UBP genes in the whole genome of T. aestivum. These genes were divided into 15 groups and non-randomly distributed on chromosomes of T. aestivum. Analyses of evolutionary patterns revealed that TaUBPs mainly underwent purification selection. The studies of cis-acting regulatory elements indicated that they might be involved in response to hormones. Quantitative real-time PCR (qRT-PCR) results showed that TaUBPs were differentially expressed in different tissues. Besides, several TaUBPs were significantly up-regulated when plants were treated with salicylic acid (SA), implying that these DUBs may play a role in abiotic stress responses in plants and few TaUBPs displayed differential expression after viral infection. Furthermore, TaUBP1A.1 (TraesCS1A02G432600.1) silenced by virus-induced gene silencing (VIGS) facilitates Chinese wheat mosaic virus (CWMV) infection in wheat, indicating that TaUBP1A.1 may be involved in a defense mechanism against viruses. This study comprehensively analyzed the UBP gene family in wheat and provided a basis for further research of TaUBPs functions in wheat plant response to viral infection.

Transcriptome-wide identification and characterization of toll-like receptors response to Vibrio anguillarum infection in Manila clam (Ruditapes philippinarum)

The Manila clam (Ruditapes philippinarum), one of the major marine aquaculture species in China, is susceptible to infection with the pathogen Vibrio, which results in massive mortality and economic losses. Toll-like receptors (TLRs) are significant pattern recognition receptors (PRRs) of innate immunity that are involved in immune regulation against pathogenic invasion. Molecular characterization of Manila clam TLRs and investigations of their immune functions are essential to prevent and control Vibrio infection. In the present research, eight cDNA sequences of R. philippinarum TLRs (RpTLRs) were identified from previous transcriptome libraries and then classified into four groups, namely, P-TLR (one sequence), V-TLR (one sequence), Ls-TLR (two sequences) and sP-TLR (four sequences), based on the corresponding LRR domain arrangement of their protein structures within the typical TLR motifs. A selective pressure test firstly suggested that the molluscan P-TLR, V-TLR, Ls-TLR and sP-TLR families underwent positive selection, and different numbers of positive selection sites (PSSs) were identified in different domains of the four types of RpTLRs, as determined by PAML and analysis of website data. These findings indicated that the evolution of RpTLRs may be associated with their immune recognition and function. Furthermore, tissue-specific expression analysis showed that all RpTLRs were ubiquitously expressed in all test tissues and were dominant in hemocytes. Quantitative real-time PCR revealed that the cDNA expression of all eight RpTLRs was upregulated after injection with Vibrio anguillarum (P < 0.01) in R. philippinarum hemocytes, revealing that these RpTLRs play important roles in responding to pathogenic stimulation. In summary, these findings provide a foundation for future investigations of the molecular classification and evolutionary patterns of Toll-like receptors in invertebrates, and the innate immune responses of TLR signaling pathways in Mollusca.

Evidence of genomic information and structural restrictions of HIV-1 PR and RT gene regions from individuals experiencing antiretroviral virologic failure

OBJECTIVES

This study analyzed Protease-PR and Reverse Transcriptase-RT HIV-1 genomic information entropy metrics among patients under antiretroviral virologic failure, according to the numbers of virologic failures or resistance mutations.

METHODS

For this purpose, we used genomic sequences from PR and RT of HIV-1 from a cohort of chronic patients followed up at São Paulo Hospital.

RESULTS

Informational entropy proportionally increases with the number of antiretroviral virologic failures in PR and RT (p < .001). Affected regions of PR were related to catalytic and structural functions, such as Fulcrum (K20) Flap (M46) and Cantilever (A71). In RT, this occurred at Fingers (E44) and Palm (K219). Informational entropy increases according to the number of resistance mutations in PR and RT (p < .001). Higher PR entropy was proportional to the resistance mutation numbers in Fulcrum (L10), Active site (L24) Flap (M46), Cantilever (L63) and near Interface (L90). In RT, they related to regions responsible for protein stability such as Fingers (T39) and Palm (L100).

CONCLUSIONS

The antiretroviral selective pressure affects HIV genomic informational entropy at the PR and RT regions, leading to the emergence of more unstable virions. Mapping the three-dimensional structure in these HIV-1 proteins is relevant to designing new antiretroviral targeting resistant strains.

Evolutionary fates of universal stress protein paralogs in Platyhelminthes

Background

Universal stress proteins (USPs) are present in all domains of life. Their expression is upregulated in response to a large variety of stress conditions. The functional diversity found in this protein family, paired with the sequence degeneration of the characteristic ATP-binding motif, suggests a complex evolutionary pattern for the paralogous USP-encoding genes. In this work, we investigated the origin, genomic organization, expression patterns and evolutionary history of the USP gene family in species of the phylum Platyhelminthes.

Results

Our data showed a cluster organization, a lineage-specific distribution, and the presence of several pseudogenes among the USP gene copies identified. The absence of a well conserved -CCAATCA- motif in the promoter region was positively correlated with low or null levels of gene expression, and with amino acid changes within the ligand binding motifs. Despite evidence of the pseudogenization of various USP genes, we detected an important functional divergence at several residues, mostly located near sites that are critical for ligand interaction.

Conclusions

Our results provide a broad framework for the evolution of the USP gene family, based on the emergence of new paralogs that face very contrasting fates, including pseudogenization, subfunctionalization or neofunctionalization. This framework aims to explain the sequence and functional diversity of this gene family, providing a foundation for future studies in other taxa in which USPs occur.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1129-x) contains supplementary material, which is available to authorized users.

Comparative and evolutionary analysis of an adapter molecule MyD88 in invertebrate metazoans

The myeloid differentiation factor 88 (MyD88) is an essential adapter in Toll-like receptor (TLR) signalling pathways, with TLR the first pattern-recognition receptor (PRR) that was discovered in Drosophila. In the present study, a MyD88 gene was identified and characterized from a commercially important shellfish, Scapharca subcrenata, including a DEATH domain and TIR domain conserved within other molluscs. Furthermore, comparative genomic evidence revealed that MyD88 was of different lengths and contained quantitative exon and intron regions, which might be involved in specific mechanisms. To further explore the phylogenetic relationships of invertebrate metazoan MyD88, we applied MrBayes and PhyML software to construct phylogenetic trees using Bayesian and maximum likelihood approaches, respectively, which suggested that the MyD88 of Arthropoda was closely related to lower invertebrates, in contrast to morphological taxonomy. Finally, we investigated the evolutionary patterns and location of positive selection sites (PSSs) in the MyD88 gene from Arthropoda, Mollusca and Insecta using PAML software with the maximum likelihood method. The data showed that positive selection sites were detected in these groups, and partial sites were located in the TIR domain but were not found in the DEATH domain. To summarize, in this study, we report on the diversification of MyD88 in invertebrate metazoans, the specific evolutionary position of Arthropoda MyD88, and the positive selection pressures on MyD88 of Arthropoda, Mollusca and Insecta. These results are a valuable contribution to understand and clarify the evolutionary pattern of TLR/MyD88 signalling pathways in invertebrate and vertebrate taxa.

Defining the RNA-binding glycine-rich (RBG) gene superfamily: new insights into nomenclature, phylogeny, and evolutionary trends obtained by genome-wide comparative analysis of Arabidopsis, Chinese cabbage, rice and maize genomes

RNA-binding glycine-rich (RBG) proteins play diverse roles in plant growth, development, protection and genome organization. An overly broad definition for class IV glycine-rich proteins (GRPs), namely RNA-binding activity and a glycine-rich C-terminus, has resulted in many distantly related and/or non-related proteins being grouped into this class of RBGs. This definition has hampered the study of RBG evolution. In this study, we used a comparative genomic approach consisting of ortholog, homolog, synteny and phylogenetic analyses to legitimately exclude all distantly/non-related proteins from class IV GRPs and to identify 15, 22, 12 and 18 RBG proteins in Arabidopsis, Chinese cabbage, rice and maize genomes, respectively. All identified RBGs could be divided into three subclasses, namely RBGA, RBGB and RBGD, which may be derived from a common ancestor. We assigned RBGs excluded from class IV GRPs to a separate RBG superfamily. RBGs have evolved and diversified in different species via different mechanisms; segmental duplication and recombination have had major effects, with tandem duplication, intron addition/deletion and domain recombination/deletion playing minor roles. Loss and retention of duplicated RBGs after polyploidization has been species and subclass specific. For example, following recent whole-genome duplication and triplication in maize and Chinese cabbage, respectively, most duplicated copies of RBGA have been lost in maize while RBGD duplicates have been retained; in Chinese cabbage, in contrast, RBGA duplicates have been retained while RBGD duplicates have been lost. Our findings reveal fundamental information and shed new light on the structural characteristics and evolutionary dynamics of RBGs.

Evolutionary analysis of RB/Rpi-blb1 locus in the Solanaceae family

Late blight caused by the oomycete Phytophthora infestans is one of the most severe threats to potato production worldwide. Numerous studies suggest that the most effective protective strategy against the disease would be to provide potato cultivars with durable resistance (R) genes. However, little is known about the origin and evolutional history of these durable R-genes in potato. Addressing this might foster better understanding of the dynamics of these genes in nature and provide clues for identifying potential candidate R-genes. Here, a systematic survey was executed at RB/Rpi-blb1 locus, an exclusive broad-spectrum R-gene locus in potato. As indicated by synteny analysis, RB/Rpi-blb1 homologs were identified in all tested genomes, including potato, tomato, pepper, and Nicotiana, suggesting that the RB/Rpi-blb1 locus has an ancient origin. Two evolutionary patterns, similar to those reported on RGC2 in Lactuca, and Pi2/9 in rice, were detected at this locus. Type I RB/Rpi-blb1 homologs have frequent copy number variations and sequence exchanges, obscured orthologous relationships, considerable nucleotide divergence, and high non-synonymous to synonymous substitutions (Ka/Ks) between or within species, suggesting rapid diversification and balancing selection in response to rapid changes in the oomycete pathogen genomes. These characteristics may serve as signatures for cloning of late blight resistance genes.