Likelihood ratio tests for detecting positive selection and application to primate lysozyme evolution

Diversification of social complexity following a major evolutionary transition in bees

How social complexity evolved remains a long-standing enigma. In most animal groups, social complexity is typically classified into a few discrete classes. This approach is oversimplified and constrains our inference of social evolution to a narrow trajectory consisting of transitions between classes. Such categorical classifications also limit quantitative studies on the molecular and environmental drivers of social complexity. The recent accumulation of relevant quantitative data has set the stage to overcome these limitations. Here, we propose a data-driven, high-dimensional approach for studying the full diversity of social phenotypes. We curated and analyzed a comprehensive dataset encompassing 17 social traits across 80 species and studied the evolution of social complexity in bees. We found that honey bees, stingless bees, and bumble bees underwent a major evolutionary transition ∼80 mya, inconsistent with the stepwise progression of the social ladder conceptual framework. This major evolutionary transition was followed by a phase of substantial phenotypic diversification of social complexity. Other bee lineages display a continuum of social complexity, ranging from solitary to simple societies, but do not reach the levels of social complexity seen in honey bees, stingless bees, and bumble bees. Bee evolution, therefore, provides a remarkable demonstration of a macroevolutionary process in which a major transition removed biological constraints and opened novel evolutionary opportunities, driving the exploration of the landscape of social phenotypes. Our approach can be extended to incorporate additional data types and readily applied to illuminate the evolution of social complexity in other animal groups.

New Gain-of-Function Mutations Prioritize Mechanisms of HER2 Activation

ERBB2 (HER2) is a well-studied oncogene with several driver mutations apart from the well-known amplification defect in some breast cancers. We used the GigaAssay to test the functional effect of HER2 missense mutations on its receptor tyrosine kinase function. The GigaAssay is a modular high-throughput one-pot assay system for simultaneously measuring molecular function of thousands of genetic variants at very high accuracy. The activities of 5,886 mutations were classified, significantly more than mutants previously reported. These variants include 112 new in vitro, 10 known, and 9 new in vivo gain-of-function (GOF) mutations. Many of the GOFs spatially cluster in sequence and structure, supporting the activation mechanisms of heterodimerization with EGFR and release of kinase inhibition by the juxtamembrane domain. Retrospective analysis of patient outcomes from the Genomic Data Commons predicts increased survival with the newly identified HER2 GOF variants.

Development and patterning of a highly versatile visual system in spiders

Visual systems provide a key interface between organisms and their surroundings, and have evolved in many forms to perform diverse functions across the animal kingdom. Spiders exhibit a range of visual abilities and ecologies, the diversity of which is underpinned by a highly versatile, modular visual system architecture. This typically includes eight eyes of two developmentally distinct types, but the number, size, location and function of the eyes can vary dramatically between lineages. Previous studies of visual system development in spiders have confirmed that many components of the retinal determination gene (RDG) network are conserved with other arthropods, but so far, comparative studies among spiders are lacking. We characterized visual system development in seven species of spiders representing a range of morphologies, visual ecologies and phylogenetic positions, to determine how these diverse configurations are formed, and how they might evolve. Combining transcriptomics, in situ hybridization, and selection analyses, we characterize the repertoires and expression of key RDGs in relation to adult morphology. We identify key molecular players, timepoints and developmental events that may contribute to adult diversity, in particular the molecular and developmental underpinnings of eye size, number, position and identity across spiders.

Chloroplast genomic insights into adaptive evolution and rapid radiation in the genus Passiflora (Passifloraceae)

Chloroplasts are essential organelles in plants and eukaryotic algae, responsible for photosynthesis, fatty acid synthesis, amino acid production, and stress responses. The genus Passiflora, known for its species diversity and dynamic chloroplast (cp) genome evolution, serves as an excellent model for studying structural variations. This study investigates evolutionary relationships within Passiflora by sequencing 11 new chloroplast genomes, assessing selective pressures on cp genes, and comparing plastid and nuclear phylogenies. Passiflora cp genomes showed significant variations in size, gene content, and structure, ranging from 132,736 to 163,292 base pairs, especially in Decaloba. Structural rearrangements and species-specific repeat patterns were identified. Selective pressure tests revealed significant adaptive evolution in certain lineages, with several genes, including clpP and petL, under positive selection. Phylogenetic analyses confirmed the monophyly of subgenera Astrophea, Passiflora, and Decaloba, while Deidamioides appeared polyphyletic. Nuclear phylogenetic analysis based on 35S rDNA sequences supported the monophyly of Astrophea but showed inconsistencies within subgenus Passiflora compared to cp genome data. This study highlights the evolutionary complexity of Passiflora cp genomes, demonstrating significant structural variations and adaptive evolution. The findings underscore the effectiveness of plastid phylogenomics in resolving phylogenetic relationships and provide insights into adaptive mechanisms shaping cp genome diversity in angiosperms.

Uncovering selection pressures on the IRF gene family in bats' immune system

Unlike other mammals, bats serve as natural reservoirs for several highly pathogenic viruses without exhibiting symptoms of infection. Recent research has explored the complex mechanisms underlying the balance between bats' antiviral defenses and their pathological responses. However, the evolution of the molecular drivers behind bats' antiviral strategies remains largely unknown. Interferon regulatory factors (IRFs) are essential transcription factors that bind to DNA and regulate the expression of numerous genes involved in antiviral defense, inflammation, immune cell differentiation, apoptosis, and oncogenesis. Our research focused on members of the IRF family, using 17 bat species and four terrestrial mammals available in GenBank. We employed CodeML to detect signs of positive selection through three different models. Statistically significant results were obtained for the IRF-1, IRF-4, IRF-5, IRF-6, and IRF-9 genes, which are known to play pivotal roles in various regulation mechanisms. Specifically, IRF-4 and IRF-5 are key in modulating the inflammatory response, while IRF-1 is essential for antiviral defense in bats, and IRF-9 regulates genes activated by type I interferon. Although the role of IRF-6 in these mechanisms requires further investigation in bats, all these genes show signs of positive selection, suggesting an optimization of the processes they regulate. These findings highlight the adaptive role of IRF elements in enhancing, among other things, the bat immune system, potentially improving their resilience and efficacy. Our study not only provides new genetic insights into bats but also underscores the remarkable molecular evolution within this unique group of mammals.

Lineage classification and selective site identification of Orthoebolavirus zairense

As the high pathogenic species of Filoviridae virus family, Orthoebolavirus zairense (EBOV) shows frequent outbreaks in human in recently years since its first emerging in 1976 in Democratic Republic of the Congo (COD), bringing ongoing risks and burden on public health safety. Here, the phylogenetic relationship among major outbreaks was analyzed. The results showed that EBOV isolates could be divided into four lineages according to spatial and temporal epidemics. Then, the positive selection sites (PSSs) were detected on all proteins of the EBOV, exhibiting lineage characteristic. Particularly, sites in GP and VP24 were identified to be significantly under positive selection, and partial of which were maintained in the latest isolates in 2021. GP and L were found to have high variability between lineages. Substitutions including F443L and F443S in GP, as well as F1610L and I1951V in L could be characteristic of the two large outbreaks in COD (2018) and West Africa (2014), respectively. Further, substitutions of significant PSSs in VP24 and L proteins were visualized for analysis of structural changes, which may affect EBOV pathogenesis. In summary, our results gains insights in genetic characteristic and adaptive evolution of EBOV, which could facilitate gene functional research against EBOV.

Testing for heterogeneous rates of discrete character evolution on phylogenies

Many hypotheses in the field of phylogenetic comparative biology involve specific changes in the rate or process of trait evolution. This is particularly true of approaches designed to connect macroevolutionary pattern to microevolutionary process. We present a method to test whether the rate of evolution of a discrete character has changed in one or more clades, lineages, or time periods. This method differs from other related approaches (such as the "covarion" model) in that the "regimes" in which the rate or process is postulated to have changed are specified a priori by the user, rather than inferred from the data. Similarly, it differs from methods designed to model a correlation between two binary traits in that the regimes mapped onto the tree are fixed. We apply our method to investigate the rate of dewlap colour and/or caudal vertebra number evolution in Caribbean and mainland clades of the diverse lizard genus Anolis. We find little evidence to support any difference in the evolutionary process between mainland and island evolution for either character. We also examine the statistical properties of the method more generally and show that it has acceptable type I error, parameter estimation, and power. Finally, we discuss some general issues of frequentist hypothesis testing and model adequacy, as well as the relationship of our method to existing models of heterogeneity in the rate of discrete character evolution on phylogenies.

Positive Selection of TLR2 and MyD88 Genes Provides Insights Into the Molecular Basis of Immunological Adaptation in Amphibians

The transition from water to land of amphibians is evolutionarily significant in the history of vertebrates, and immunological adaptation is an important challenge for amphibians to respond to the dramatic changes of the environmental pathogens during their origin and diversification. Toll-like receptors (TLRs) are important pattern recognition receptors for the innate immune response and TLRs signaling pathway play essential roles in the immune responses to pathogens and inflammatory reaction. However, the evolutionary patterns and molecular mechanisms underlying their adaptation in amphibians are poorly documented to date. Here, we determined the coding regions, expression patterns of TLR2 and Myeloid differentiation factor 88 (MyD88) in the large treefrog (Zhangixalus dennysi), and explored the evolutionary patterns of these two genes in amphibians. Quantitative Real-time PCR analyses showed that the TLR2 and MyD88 mRNA were expressed in all the organs/tissues examined, both with the highest levels in the heart and the lowest levels in the body fat for TLR2 and lung for MyD88. The highly conservation and functional significance of these two genes in amphibians were supported based on the sequence characteristics and evolutionary analyses. Significantly positive selection was found to be acting on TLR2 and MyD88 in amphibians based on different site models. Strong signal of positive selection among different amphibian lineages for these two genes was also detected and a series of positively selected sites were identified based on the branch-site analysis. Our results suggest that amphibians have adapted to different pathogenic microorganisms during their transition from the aquatic to terrestrial environment and diversification into various habitats. The present study will provide new insights into the evolutionary process and molecular basis underlying the immunological adaptation in vertebrates.

Molecular evolution of toothed whale genes reveals adaptations to echolocating in different environments

BACKGROUND

Echolocation was a key development in toothed whale evolution, enabling their adaptation and diversification across various environments. Previous bioacoustic and morphological studies suggest that environmental pressures have influenced the evolution of echolocation in toothed whales. This hypothesis demands further investigation, especially regarding the molecular mechanisms involved in the adaptive radiation of toothed whales across multiple habitats. Here we show that the coding sequences of four hearing genes involved in echolocation (CDH23, prestin, TMC1, and CLDN14) have different signatures of molecular evolution among riverine, coastal, and oceanic dolphins, suggesting that the evolutionary constraints of these habitats shaped the underlying genetic diversity of the toothed whale sonar.

RESULTS

Our comparative analysis across 37 odontocete species revealed patterns of accelerated evolution within coastal and riverine lineages, supporting the hypothesis that shallow habitats pose specific selective pressures to sonar propagation, which are not found in the deep ocean. All toothed whales with genes evolving under positive selection are shallow coastal species, including three species that have recently diverged from freshwater lineages (Cephalorhynchus commersonii, Sotalia guianensis, and Orcaella heinsohni - CDH23), and three species that operate specialized Narrow Band High Frequency (NBHF) Sonars (Phocoena sinus - prestin, Neophocaena phocaenoides - TMC1 and Cephalorhynchus commersonii - CDH23). For river dolphins and deep-diving toothed whales, we found signatures of positive selection and molecular convergence affecting specific sites on CDH23, TMC1, and prestin. Positively selected sites (PSS) were different in number, identity, and substitution rates (dN/dS) across riverine, coastal, and oceanic toothed whales.

CONCLUSION

Here we shed light on potential molecular mechanisms underlying the diversification of toothed whale echolocation. Our results suggest that toothed whale hearing genes changed under different selective pressures in coastal, riverine, and oceanic environments.

Insight into the Phylogenetic Relationships of Phasmatodea and Selection Pressure Analysis of Phraortes liaoningensis Chen & He, 1991 (Phasmatodea: Lonchodidae) Using Mitogenomes

Stick and leaf insects are a group among the Insecta that are famous for their extraordinary mimicry ability. Since the establishment of the Phasmatodea, their internal classification has been constantly revised. Mitochondrial genes as molecular markers have been widely used for species classification, but the phylogenetic relationships within the Phasmatodea remain to be thoroughly discussed. In the present study, five mitogenomes of Phasmatodea ranging from 15,746 bp to 16,747 bp in length were sequenced. Bayesian inference (BI) and maximum likelihood (ML) analyses were carried out based on a 13 PCGs data matrix (nt123) and a combined matrix of 13 PCGs and two rRNA genes (nt123_rRNA). The present study supports the conclusion that Phylliidae was the basal group of Neophasmatodea and confirms the monophyly of Lonchodinae and Necrosciinae, but it shows that Lonchodidae was polyphyletic. A sister group of Bacillidae and Pseudophasmatidae was also recovered. The phylogenetic tree based on the nt_123 dataset showed higher node support values. The construction of a divergent time tree in this study supported the conclusion that extant Phasmatodea originated in the Jurassic (170 Mya) and most lineages diverged after the Cretaceous-Paleogene extinction event. To explore whether the mitochondrial genes of Phraortes liaoningensis collected from high latitudes where low temperatures occur for eight months of the year are under selection pressure, this study used the branch-site model and the branch model to analyze the selection pressure on the 13 mitochondria protein-coding genes (PCGs). We found that both ND2 and ND4L of Ph. liaoningensis exhibited positive selection sites using the branch-site model. This study shows that a low-temperature environment causes mitochondrial genes to be selected to meet the energy requirements for survival.

First record and molecular diversity of haemosporidian parasites in Sula leucogaster (Boddaert, 1783) from the southeastern Brazilian coast

Avian haemosporidian parasites are protozoans transmitted by vectors, belonging to the genera Plasmodium, Haemoproteus, and Leucocytozoon, capable of infecting various bird species. The prevalence of these parasites varies depending on factors such as the habitat and geographical location of the host and vector, as well as across different bird groups. Seabirds typically exhibit low infection rates due to the absence of vectors in marine environments. This study aimed to identify and characterize the molecular diversity of Haemosporida parasites in a population of Brown Boobies (Sula leucogaster) (n = 37) from the Santana Archipelago in Macaé, Rio de Janeiro state. We utilized PCR-based detection methods to confirm the presence of haemosporidian parasites and sequencing (ctyb gene) to delineate the lineages. The results uncovered a notable prevalence (48.6%) of haemosporidian parasites (Plasmodium and Haemoproteus), comprising 12 novel lineages and four previously documented lineages (establishing new host records). Our findings suggest that transmission and infection by these hemoparasites can occur both on and off the island. Overall, this study yields valuable insights into the genetic diversity of haemosporidian parasites in Brown Boobies, enriching our understanding of the ecology and epidemiology of these parasites in seabirds.

Parallel Spectral Tuning of a Cone Visual Pigment Provides Evidence for Ancient Deep-Sea Adaptations in Cetaceans

Dichromatic color vision is mediated by two cone visual pigments in many eutherian mammals. After reentry into the sea, early cetaceans lost their violet-sensitive visual pigment (short wavelength-sensitive 1) independently in the baleen and toothed whale ancestors and thus obtained only monochromatic cone vision. Subsequently, losses of the middle/long wavelength-sensitive (M/LWS) pigment have also been reported in multiple whale lineages, leading to rhodopsin (RH1)-mediated rod monochromatic vision. To further elucidate the phenotypic evolution of whale visual pigments, we assessed the spectral tuning of both M/LWS and RH1 from representative cetacean taxa. Interestingly, although the coding sequences for M/LWS are intact in both the pygmy right whale and the Baird's beaked whale, no spectral sensitivity was detected in vitro. Pseudogenization of other cone vision-related genes is observed in the pygmy right whale, suggesting a loss of cone-mediated vision. After ancestral sequence reconstructions, ancient M/LWS pigments from cetacean ancestors were resurrected and functionally measured. Spectral tuning of M/LWS from the baleen whale ancestor shows that it is green sensitive, with a 40-nm shift in sensitivity to a shorter wavelength. For the ancestor of sperm whales, although no spectral sensitivity could be recorded for its M/LWS pigment, a substantial sensitivity shift (20 to 30 nm) to a shorter wavelength may have also occurred before its functional inactivation. The parallel phenotypic evolution of M/LWS to shorter wavelength sensitivity might be visual adaptations in whales allowing more frequent deep-sea activities, although additional ecological differentiations may have led to their subsequent losses.

Relaxed selection in evolution of genes regulating limb development gives clue to variation in forelimb morphology of cetaceans and other mammals

Cetaceans have evolved unique limb structures, such as flippers, due to genetic changes during their transition to aquatic life. However, the full understanding of the genetic and evolutionary mechanisms behind these changes is still developing. By examining 25 limb-related protein-coding genes across various mammalian species, we compared genetic changes between aquatic mammals, like whales, and other mammals with unique limb structures such as bats, rodents and elephants. Our findings revealed significant modifications in limb-related genes, including variations in the Hox, GDF5 and Evx genes. Notably, a relaxed selection in several key genes was observed, suggesting a lifting of developmental constraints, which might have facilitated the emergence of morphological innovations in cetacean limb morphology. We also uncovered non-synonymous changes, insertions and deletions in these genes, particularly in the polyalanine tract of HOXD13, which are distinctive to cetaceans or convergent with other aquatic mammals. These genetic variations correlated with the diverse and specialized limb structures observed in cetaceans, indicating a complex interplay of relaxed selection and specific mutations in mammalian limb evolution.

Adaptive evolution of pancreatic ribonuclease gene (RNase1) in Cetartiodactyla

Pancreatic ribonuclease (RNase1), a digestive enzyme produced by the pancreas, is associated with the functional adaptation of dietary habits and is regarded as an attractive model system for studies of molecular evolution. In this study, we identified 218 functional genes and 48 pseudogenes from 114 species that span all four Cetartiodactyla lineages: two herbivorous lineages (Ruminantia and Tylopoda) and two non-herbivorous lineages (Cetancodonta and Suoidea). Multiple RNase1 genes were detected in all species of the two herbivorous lineages, and phylogenetic and genomic location analyses demonstrated that independent gene duplication events occurred in Ruminantia and Tylopoda. In Ruminantia, the gene duplication events occurred in the ancestral branches of the lineage in the Middle Eocene, a time of increasing climatic seasonality during which Ruminantia rapidly radiated. In contrast, only a single RNase1 gene was observed in the species of the two non-herbivorous lineages (Cetancodonta and Suoidea), suggesting that the previous Cetacea-specific loss hypothesis should be rejected. Moreover, the duplicated genes of RNase1 in the two herbivorous lineages (Ruminantia and Tylopoda) may have undergone functional divergence. In combination with the temporal coincidence between gene replication and the enhanced climatic seasonality during the Middle Eocene, this functional divergence suggests that RNase1 gene duplication was beneficial for Ruminantia to use the limited quantities of sparse fibrous vegetation and adapt to seasonal changes in climate. In summary, the findings indicate a complex and intriguing evolutionary pattern of RNase1 in Cetartiodactyla and demonstrate the molecular mechanisms by which organisms adapt to the environment.

Comparative analysis of the complete chloroplast genomes of thirteen Bougainvillea cultivars from South China with implications for their genome structures and phylogenetic relationships

Bougainvillea spp., belonging to the Nyctaginaceae family, have high economic and horticultural value in South China. Despite the high similarity in terms of leaf appearance and hybridization among Bougainvillea species, especially Bougainvillea × buttiana, their phylogenetic relationships are very complicated and controversial. In this study, we sequenced, assembled and analyzed thirteen complete chloroplast genomes of Bougainvillea cultivars from South China, including ten B. × buttiana cultivars and three other Bougainvillea cultivars, and identified their phylogenetic relationships within the Bougainvillea genus and other species of the Nyctaginaceae family for the first time. These 13 chloroplast genomes had typical quadripartite structures, comprising a large single-copy (LSC) region (85,169-85,695 bp), a small single-copy (SSC) region (18,050-21,789 bp), and a pair of inverted-repeat (IR) regions (25,377-25,426 bp). These genomes each contained 112 different genes, including 79 protein-coding genes, 29 tRNAs and 4 rRNAs. The gene content, codon usage, simple sequence repeats (SSRs), and long repeats were essentially conserved among these 13 genomes. Single-nucleotide polymorphisms (SNPs) and insertions/deletions (indels) were detected among these 13 genomes. Four divergent regions, namely, trnH-GUG_psbA, trnS-GCU_trnG-UCC-exon1, trnS-GGA_rps4, and ccsA_ndhD, were identified from the comparative analysis of 16 Bougainvillea cultivar genomes. Among the 46 chloroplast genomes of the Nyctaginaceae family, nine genes, namely, rps12, rbcL, ndhF, rpoB, rpoC2, ndhI, psbT, ycf2, and ycf3, were found to be under positive selection at the amino acid site level. Phylogenetic relationships within the Bougainvillea genus and other species of the Nyctaginaceae family based on complete chloroplast genomes and protein-coding genes revealed that the Bougainvillea genus was a sister to the Belemia genus with strong support and that 35 Bougainvillea individuals were divided into 4 strongly supported clades, namely, Clades Ⅰ, Ⅱ, Ⅲ and Ⅳ. Clade Ⅰ included 6 individuals, which contained 2 cultivars, namely, B. × buttiana 'Gautama's Red' and B. spectabilis 'Flame'. Clades Ⅱ only contained Bougainvillea spinosa. Clade Ⅲ comprised 7 individuals of wild species. Clade Ⅳ included 21 individuals and contained 11 cultivars, namely, B. × buttiana 'Mahara', B. × buttiana 'California Gold', B. × buttiana 'Double Salmon', B. × buttiana 'Double Yellow', B. × buttiana 'Los Banos Beauty', B. × buttiana 'Big Chitra', B. × buttiana 'San Diego Red', B. × buttiana 'Barbara Karst', B. glabra 'White Stripe', B. spectabilis 'Splendens' and B. × buttiana 'Miss Manila' sp. 1. In conclusion, this study not only provided valuable genome resources but also helped to identify Bougainvillea cultivars and understand the chloroplast genome evolution of the Nyctaginaceae family.

Genome-wide identification of the phenylalanine ammonia-lyase gene from Epimedium Pubescens Maxim. (Berberidaceae): novel insight into the evolution of the PAL gene family

BACKGROUND

Phenylalanine ammonia-lyase (PAL) serves as a key gateway enzyme, bridging primary metabolism and the phenylpropanoid pathway, and thus playing an indispensable role in flavonoid, anthocyanin and lignin biosynthesis. PAL gene families have been extensively studied across species using public genomes. However, a comprehensive exploration of PAL genes in Epimedium species, especially those involved in prenylated flavonol glycoside, anthocyanin, or lignin biosynthesis, is still lacking. Moreover, an in-depth investigation into PAL gene family evolution is warranted.

RESULTS

Seven PAL genes (EpPAL1-EpPAL7) were identified. EpPAL2 and EpPAL3 exhibit low sequence identity to other EpPALs (ranging from 61.09 to 64.38%) and contain two unique introns, indicating distinct evolutionary origins. They evolve at a rate ~ 10 to ~ 54 times slower compared to EpPAL1 and EpPAL4-7, suggesting strong purifying selection. EpPAL1 evolved independently and is another ancestral gene. EpPAL1 formed EpPAL4 through segmental duplication, which lead to EpPAL5 and EpPAL6 through tandem duplications, and EpPAL7 through transposed duplication, shaping modern EpPALs. Correlation analysis suggests EpPAL1, EpPAL2 and EpPAL3 play important roles in prenylated flavonol glycosides biosynthesis, with EpPAL2 and EpPAL3 strongly correlated with both Epimedin C and total prenylated flavonol glycosides. EpPAL1, EpPAL2 and EpPAL3 may play a role in anthocyanin biosynthesis in leaves. EpPAL2, EpPAL3, EpPAL6, and EpPAL7 might be engaged in anthocyanin production in petals, and EpPAL2 and EpPAL3 might also contribute to anthocyanin synthesis in sepals. Further experiments are needed to confirm these hypotheses. Novel insights into the evolution of PAL gene family suggest that it might have evolved from a monophyletic group in bryophytes to large-scale sequence differentiation in gymnosperms, basal angiosperms, and Magnoliidae. Ancestral gene duplications and vertical inheritance from gymnosperms to angiosperms likely occurred during PAL evolution. Most early-diverging eudicotyledons and monocotyledons have distinct histories, while modern angiosperm PAL gene families share similar patterns and lack distant gene types.

CONCLUSIONS

EpPAL2 and EpPAL3 may play crucial roles in biosynthesis of prenylated flavonol glycosides and anthocyanins in leaves and flowers. This study provides novel insights into PAL gene family evolution. The findings on PAL genes in E. pubescens will aid in synthetic biology research on prenylated flavonol glycosides production.

Expression of defensin genes across house fly (Musca domestica) life history gives insight into immune system subfunctionalization

Animals encounter diverse microbial communities throughout their lifetime, which exert varying selection pressures. Antimicrobial peptides (AMPs), which lyse or inhibit microbial growth, are a first line of defense against some of these microbes. Here we examine how developmental variation in microbial exposure has affected the evolution of expression and amino acid sequences of Defensins (an ancient class of AMPs) in the house fly (Musca domestica). The house fly is a well-suited model for this work because it trophically associates with varying microbial communities throughout its life history and its genome contains expanded families of AMPs, including Defensins. We identified two subsets of house fly Defensins: one expressed in larvae or pupae, and the other expressed in adults. The amino acid sequences of these two Defensin subsets form distinct monophyletic clades, and they are located in separate gene clusters in the genome. The adult-expressed Defensins evolve faster than larval/pupal Defensins, consistent with different selection pressures across developmental stages. Our results therefore suggest that varied microbial communities encountered across life history can shape the evolutionary trajectories of immune genes.

Whole genome sequencing and genome characterization of Aichivirus isolated from Korean adults

The whole-genome sequence (WGS) analysis of Aichivirus (AiV) identified in Korea was performed in this study. Using Sanger and Nanopore sequencing, the 8228-nucleotide-long genomic sequence of AiV (OQ121963) was determined and confirmed to belong to genotype A. The full-length genome of OQ121963 consisted of a 7296 nt open reading frame (ORF) that encodes a single polyprotein, and 5' UTR (676 nt) and 3' UTR (256 nt) at 5' and 3' ends, respectively. The ORF consisted of leader protein (L), structural protein P1 (VP0, VP1, and VP3), and nonstructural protein P2 (2A, 2B, and 2C) and P3 (3A, 3B, 3C, and 3D). The secondary structure analysis of the 5' UTR identified only stem-loop C (SL-C) and not SL-A and SL-B. The variable region of the AiV genome was analyzed by MegAlign Pro and reconfirmed by SimPlot analysis using 16 AiV whole genomes known to date. Among the entire regions, structural protein region P1 showed the lowest amino acid identity (96.07%) with reference sequence AB040749 (originated in Japan; genotype A), while the highest amino acid identity (98.26%) was confirmed in the 3D region among nonstructural protein region P2 and P3. Moreover, phylogenetic analysis of the WGS of OQ121963 showed the highest homology (96.96%) with JX564249 (originated in Taiwan; genotype A) and lowest homology (90.14%) with DQ028632 (originated in Brazil; genotype B). Therefore, the complete genome characterization of OQ121963 and phylogenetic analysis of the AiV conducted in this study provide useful information allowing to improve diagnostic tools and epidemiological studies of AiVs.

Draft Genome of Akame (Lates Japonicus) Reveals Possible Genetic Mechanisms for Long-Term Persistence and Adaptive Evolution with Low Genetic Diversity

It is known that some endangered species have persisted for thousands of years despite their very small effective population sizes and low levels of genetic polymorphisms. To understand the genetic mechanisms of long-term persistence in threatened species, we determined the whole genome sequences of akame (Lates japonicus), which has survived for a long time with extremely low genetic variations. Genome-wide heterozygosity in akame was estimated to be 3.3 to 3.4 × 10-4/bp, one of the smallest values in teleost fishes. Analysis of demographic history revealed that the effective population size in akame was around 1,000 from 30,000 years ago to the recent past. The relatively high ratio of nonsynonymous to synonymous heterozygosity in akame indicated an increased genetic load. However, a detailed analysis of genetic diversity in the akame genome revealed that multiple genomic regions, including genes involved in immunity, synaptic development, and olfactory sensory systems, have retained relatively high nucleotide polymorphisms. This implies that the akame genome has preserved the functional genetic variations by balancing selection, to avoid a reduction in viability and loss of adaptive potential. Analysis of synonymous and nonsynonymous nucleotide substitution rates has detected signs of positive selection in many akame genes, suggesting adaptive evolution to temperate waters after the speciation of akame and its close relative, barramundi (Lates calcarifer). Our results indicate that the functional genetic diversity likely contributed to the long-term persistence of this species by avoiding the harmful effects of the population size reduction.

Unique myoglobin adaptation to endothermy and flight since the origin of birds

Myoglobin (Mb) mediates oxygen diffusion and storage in muscle tissue and thus is important for the energy utilization and activity of animals. Birds generally have a high body temperature, and most species also possess the capability of powered flight. Both of these require high levels of aerobic metabolism. Within endothermic mammals, bats also independently evolved flight. Although the functional evolution of myoglobins in deep-diving amniote vertebrates has been well-studied, the functional evolution of myoglobin since the origins of both birds and bats is unclear. Here, with Mb-coding sequences from >200 extant amniote species, we reconstructed ancestral sequences to estimate the functional properties of myoglobin through amniote evolution. A dramatic change in net surface charge on myoglobin occurred during the origin of Aves, which might have been driven by positively selected amino acid substitutions that occurred on the lineage leading to all birds. However, in bats, no change in net surface charge occurred and instead, the Mb genes show evidence of strong purifying selection. The increased net surface charge on bird myoglobins implies an adaptation to flight-related endothermic and higher body temperatures, possibly by reducing harmful protein aggregations. Different from the findings of net surface charge, myoglobins of extant birds show lower stability compared with other amniotes, which probably accelerates the rate of oxygen utilization in muscles. In bats and other mammals, higher stability of Mb may be an alternative pathway for adaptation to endothermy, indicating divergent evolution of myoglobin in birds and bats.

Identification and Characterization of an Alphacoronavirus in Rhinolophus sinicus and a Betacoronavirus in Apodemus ilex in Yunnan, China

Coronaviruses (CoVs), the largest positive-sense RNA viruses, have caused infections in both humans and animals. The cross-species transmission of CoVs poses a serious threat to public health. Rodents and bats, the two largest orders of mammals, serve as significant natural reservoirs for CoVs. It is important to monitor the CoVs carried by bats and rodents. In this study, we collected 410 fecal samples from bats and 74 intestinal samples from rats in Yunnan Province, China. Using RT-PCR, we identified one positive sample for alphacoronavirus (TC-14) from Rhinolophus sinicus (Chinese rufous horseshoe bat) and two positive samples for betacoronavirus (GS-53, GS-56) from Apodemus ilex (Rodentia: Muridae). We successfully characterized the complete genomes of TC-14 and GS-56. Phylogenetic analysis revealed that TC-14 clustered with bat CoV HKU2 and SADS-CoV, while GS-56 was closely related to rat CoV HKU24. The identification of positive selection sites and estimation of divergence dates further helped characterize the genetic evolution of TC-14 and GS-56. In summary, this research reveals the genetic evolution characteristics of TC-14 and GS-56, providing valuable references for the study of CoVs carried by bats and rodents in Yunnan Province.

Purifying selection drove the adaptation of mitochondrial genes along with correlation of gene rearrangements and evolutionary rates in two subfamilies of Whitefly (Insecta: Hemiptera)

Insect mitochondrial genomes (mitogenomes) are usually represented by a conserved gene order. Whiteflies exhibit gene rearrangement in their mitogenomes; however, understanding how nucleotide substitution rates shape gene rearrangement in whiteflies is unclear due to the limited number of mitogenomes. Additionally, the mechanisms by which selection pressure drives adaptations in mitochondrial genes in the two subfamilies of whiteflies are not yet known. Here, we analyzed 18 whitefly mitogenomes, including one newly generated mitogenome, to compare nucleotide substitution rates, selection pressure, and gene arrangements. The newly generated mitogenome is reported along with reannotation of Pealius mori and comparisons to other whitefly mitogenomes. Comparative studies on nucleotide composition of 18 whiteflies revealed the positive GC skewness, confirming the reversal of strand asymmetry. We found 11 rearranged gene orders within two subfamilies of whiteflies with 8-18 breakpoints of gene rearrangements. Members of the subfamily Aleyrodinae exhibit more complex pathways in the evolution of gene order as compared to the subfamily Aleurodicinae. Our findings also revealed that the increase or reduction of nucleotide substitution rates does not have an impact on any of the gene rearrangement scenarios depicting neutral correlation. Selection pressure analysis revealed that the mitogenomes from members of both the subfamilies Aleurodicinae and Aleyrodinae are characterized by intense purifying selection pressure.

Gene Protein Sequence Evolution Can Predict the Rapid Divergence of Ovariole Numbers in the Drosophila melanogaster Subgroup

Ovaries play key roles in fitness and evolution: they are essential female reproductive structures that develop and house the eggs in sexually reproducing animals. In Drosophila, the mature ovary contains multiple tubular egg-producing structures known as ovarioles. Ovarioles arise from somatic cellular structures in the larval ovary called terminal filaments (TFs), formed by TF cells and subsequently enclosed by sheath (SH) cells. As in many other insects, ovariole number per female varies extensively in Drosophila. At present, however, there is a striking gap of information on genetic mechanisms and evolutionary forces that shape the well-documented rapid interspecies divergence of ovariole numbers. To address this gap, here we studied genes associated with Drosophila melanogaster ovariole number or functions based on recent experimental and transcriptional datasets from larval ovaries, including TFs and SH cells, and assessed their rates and patterns of molecular evolution in five closely related species of the melanogaster subgroup that exhibit species-specific differences in ovariole numbers. From comprehensive analyses of protein sequence evolution (dN/dS), branch-site positive selection, expression specificity (tau), and phylogenetic regressions (phylogenetic generalized least squares), we report evidence of 42 genes that showed signs of playing roles in the genetic basis of interspecies evolutionary change of Drosophila ovariole number. These included the signaling genes upd2 and Ilp5 and extracellular matrix genes vkg and Col4a1, whose dN/dS predicted ovariole numbers among species. Together, we propose a model whereby a set of ovariole-involved gene proteins have an enhanced evolvability, including adaptive evolution, facilitating rapid shifts in ovariole number among Drosophila species.

Molecular diversity and polyparasitism of avian trypanosomes in the Brazilian Atlantic Rainforest

The current study proposes to investigate the diversity and phylogeny of trypanosomes parasitizing wild birds from the Brazilian Atlantic Forest. Cytological examination was carried out by light microscopy of blood smears and positive birds were selected for amplification of the 18S rDNA sequence through PCR. The resulting amplicons were subjected to purification, cloning, and sequencing analysis. Phylogenetic reconstruction was conducted, including all avian trypanosomes representative's lineages. A total of ten bird samples from species of Turdus flavipes (N=1/12), T. albicollis (N=1/8), Tachyphonus coronatus (N=6/121), Thamnophilus caerulescens (N=1/22) and Synallaxis spixi (N=1/8) were positive for Trypanosoma spp. In the six specimens of T. coronatus, five distinct lineages of Trypanosoma spp. 18S-rRNA were observed in ninety sequences obtained, and using the strategy of cloning independent PCR, it was possible to observe that two of them were related to T. avium (JB01/JB02), and three were closed related to T. bennetti (JB03/ JB04/JB05). Addionaly, all fifteen sequences obtained from T. caerulescens/ S. spixi/T. flavipes/T. albicollis were identical. The present research is the first study to access molecular diversity and polyparasitism by avian trypanosomes in Brazil. The current research exhibits the wide genetic variability in avian trypanosomes and its non-specific relationship with its avian hosts.

Binding Evolution of the Dengue Virus Envelope Against DC-SIGN: A Combined Approach of Phylogenetics and Molecular Dynamics Analyses Over 30 Years of Dengue Virus in Brazil

The Red Queen Hypothesis (RQH), derived from Lewis Carroll's "Through the Looking-Glass", postulates that organisms must continually adapt in response to each other to maintain relative fitness. Within the context of host-pathogen interactions, the RQH implies an evolutionary arms race, wherein viruses evolve to exploit hosts and hosts evolve to resist viral invasion. This study delves into the dynamics of the RQH in the context of virus-cell interactions, specifically focusing on virus receptors and cell receptors. We observed multiple virus-host systems and noted patterns of co-evolution. As viruses evolved receptor-binding proteins to effectively engage with cell receptors, cells countered by altering their receptor genes. This ongoing mutual adaptation cycle has influenced the molecular intricacies of receptor-ligand interactions. Our data supports the RQH as a driving force behind the diversification and specialization of both viral and host cell receptors. Understanding this co-evolutionary dance offers insights into the unpredictability of emerging viral diseases and potential therapeutic interventions. Future research is crucial to dissect the nuanced molecular changes and the broader ecological consequences of this ever-evolving battle. Here, we combine phylogenetic inferences, structural modeling, and molecular dynamics analyses to describe the epidemiological characteristics of major Brazilian DENV strains that circulated from 1990 to 2022 from a combined perspective, thus providing us with a more detailed picture on the dynamics of such interactions over time.

Breathing Air and Living Underwater: Molecular Evolution of Genes Related to Antioxidant Response in Cetaceans and Pinnipeds

Cetaceans and pinnipeds are lineages of mammals that have independently returned to the aquatic environment, acquiring varying degrees of dependence on it while sharing adaptations for underwater living. Here, we focused on one critical adaptation from both groups, their ability to withstand the ischemia and reperfusion experienced during apnea diving, which can lead to the production of reactive oxygen species (ROS) and subsequent oxidative damage. Previous studies have shown that cetaceans and pinnipeds possess efficient antioxidant enzymes that protect against ROS. In this study, we investigated the molecular evolution of key antioxidant enzyme genes (CAT, GPX3, GSR, PRDX1, PRDX3, and SOD1) and the ROS-producing gene XDH, in cetaceans and pinnipeds lineages. We used the ratio of non-synonymous (dN) to synonymous (dS) substitutions as a measure to identify signatures of adaptive molecular evolution in these genes within and between the two lineages. Additionally, we performed protein modeling and variant impact analyzes to assess the functional consequences of observed mutations. Our findings revealed distinct selective regimes between aquatic and terrestrial mammals in five of the examined genes, including divergences within cetacean and pinniped lineages, between ancestral and recent lineages and between crowns groups. We identified specific sites under positive selection unique to Cetacea and Pinnipedia, with one site showing evidence of convergent evolution in species known for their long and deep-diving capacities. Notably, many sites under adaptive selection exhibited radical changes in amino acid properties, with some being damaging mutations in human variations, but with no apparent detrimental impacts on aquatic mammals. In conclusion, our study provides insights into the adaptive changes that have occurred in the antioxidant systems of aquatic mammals throughout their evolutionary history. We observed both distinctive features within each group of Cetacea and Pinnipedia and instances of convergence. These findings highlight the dynamic nature of the antioxidant system in response to challenges of the aquatic environment and provide a foundation for further investigations into the molecular mechanisms underlying these adaptations.

Comparative chloroplast genomics, phylogenetic relationships and molecular markers development of Aglaonema commutatum and seven green cultivars of Aglaonema

Aglaonema commutatum is a famous species in the Aglaonema genus, which has important ornamental and economic value. However, its chloroplast genome information and phylogenetic relationships among popular green cultivars of Aglaonema in southern China have not been reported. Herein, chloroplast genomes of one variety of A. commutatum and seven green cultivars of Aglaonema, namely, A. commutatum 'San Remo', 'Kai Sa', 'Pattaya Beauty', 'Sapphire', 'Silver Queen', 'Snow White', 'White Gem', and 'White Horse Prince', were sequenced and assembled for comparative analysis and phylogeny. These eight genomes possessed a typical quadripartite structure that consisted of a LSC region (90,799-91,486 bp), an SSC region (20,508-21,137 bp) and a pair of IR regions (26,661-26,750 bp). Each genome contained 112 different genes, comprising 79 protein-coding genes, 29 tRNA genes and 4 rRNA genes. The gene orders, GC contents, codon usage frequency, and IR/SC boundaries were highly conserved among these eight genomes. Long repeats, SSRs, SNPs and indels were analyzed among these eight genomes. Comparative analysis of 15 Aglaonema chloroplast genomes identified 7 highly variable regions, including trnH-GUG-exon1-psbA, trnS-GCU-trnG-UCC-exon1, trnY-GUA-trnE-UUC, psbC-trnS-UGA, trnF-GAA-ndhJ, ccsA-ndhD, and rps15-ycf1-D2. Reconstruction of the phylogenetic trees based on chloroplast genomes, strongly supported that Aglaonema was a sister to Anchomanes, and that the Aglaonema genus was classified into two sister clades including clade I and clade II, which corresponded to two sections, Aglaonema and Chamaecaulon, respectively. One variety and five cultivars, including A. commutatum 'San Remo', 'Kai Sa', 'Pattaya Beauty', 'Silver Queen', 'Snow White', and 'White Horse Prince', were classified into clade I; and the rest of the two cultivars, including 'Sapphire' and 'White Gem', were classified into clade II. Positive selection was observed in 34 protein-coding genes at the level of the amino acid sites among 77 chloroplast genomes of the Araceae family. Based on the highly variable regions and SSRs, 4 DNA markers were developed to differentiate the clade I and clade II in Aglaonema. In conclusion, this study provided chloroplast genomic resources for Aglaonema, which were useful for its classification and phylogeny.

Ascorbic acid biosynthesis in Pacific abalone Haliotis discus hannai Ino and L-gulonolactone oxidase gene loss as an independent event

Up to now, it has been believed that invertebrates are unable to synthesize ascorbic acid (AA) in vivo. However, in the present study, the full-length CDs (Coding sequence) of L-gulonolactone oxidase (GLO) from Pacific abalone (Haliotis discus hannai Ino) were obtained through molecular cloning. The Pacific abalone GLO contained a FAD-binding domain in the N-termination, and ALO domain and conserved HWAK motif in the C-termination. The GLO gene possesses 12 exons and 11 introns. The Pacific abalone GLO was expressed in various tissues, including the kidney, digestive gland, gill, intestine, muscle and mantle. The GLO activity assay revealed that GLO activity was only detected in the kidney of Pacific abalone. After a 100-day feeding trial, dietary AA levels did not significantly affect the survival, weight gain, daily increment in shell length, and feed conversion ratio of Pacific abalone. The expression of GLO in the kidney was downregulated by dietary AA. These results implied that the ability to synthesize AA in abalone had not been lost. From the evolutionary perspective, the loss of GLO occurred independently as an independent event by matching with the genomes of various species. The positive selection analysis revealed that the GLO gene underwent purifying selective pressure during its evolution. In conclusion, the present study provided direct evidence to prove that the GLO activity and the ability to synthesize AA exist in abalone. The AA synthesis ability in vertebrates might have originated from invertebrates dating back 930.31 million years.

Adaptive evolution of antioxidase-related genes in hypoxia-tolerant mammals

To cope with the damage from oxidative stress caused by hypoxia, mammals have evolved a series of physiological and biochemical traits, including antioxidant ability. Although numerous research studies about the mechanisms of hypoxia evolution have been reported, the molecular mechanisms of antioxidase-related genes in mammals living in different environments are yet to be completely understood. In this study, we constructed a dataset comprising 7 antioxidase-related genes (CAT, SOD1, SOD2, SOD3, GPX1, GPX2, and GPX3) from 43 mammalian species to implement evolutionary analysis. The results showed that six genes (CAT, SOD1, SOD2, SOD3, GPX1, and GPX3) have undergone divergent evolution based on the free-ratio (M1) model. Furthermore, multi-ratio model analyses uncovered the divergent evolution between hypoxic and non-hypoxic lineages, as well as various hypoxic lineages. In addition, the branch-site model identified 9 positively selected branches in 6 genes (CAT, SOD1, SOD2, SOD3, GPX2, and GPX3) that contained 35 positively selected sites, among which 31 positively selected sites were identified in hypoxia-tolerant branches, accounting for 89% of the total number of positively selected sites. Interestingly, 65 parallel/convergent sites were identified in the 7 genes. In summary, antioxidase-related genes are subjected to different selective pressures among hypoxia-tolerant species living in different habitats. This study provides a valuable insight into the molecular evolution of antioxidase-related genes in hypoxia evolution in mammals.

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.

Genomic insights into Yak (Bos grunniens) adaptations for nutrient assimilation in high-altitudes

High-altitude environments present formidable challenges for survival and reproduction, with organisms facing limited oxygen availability and scarce nutrient resources. The yak (Bos grunniens), indigenous to the Tibetan Plateau, has notably adapted to these extreme conditions. This study delves into the genomic basis of the yak's adaptation, focusing on the positive selection acting on genes involved in nutrient assimilation pathways. Employing techniques in comparative genomics and molecular evolutionary analyses, we selected genes in the yak that show signs of positive selection associated with nutrient metabolism, absorption, and transport. Our findings reveal specific genetic adaptations related to nutrient metabolism in harsh climatic conditions. Notably, genes involved in energy metabolism, oxygen transport, and thermoregulation exhibited signs of positive selection, suggesting their crucial role in the yak's successful colonization of high-altitude regions. The study also sheds light on the yak's immune system adaptations, emphasizing genes involved in response to various stresses prevalent at elevated altitudes. Insights into the yak's genomic makeup provide valuable information for understanding the broader implications of high-altitude adaptations in mammalian evolution. They may contribute to efforts in enhancing livestock resilience to environmental challenges.

A comprehensive evolutionary scenario for the origin and neofunctionalization of the Drosophila speciation gene Odysseus (OdsH)

Odysseus (OdsH) was the first speciation gene described in Drosophila related to hybrid sterility in offspring of mating between Drosophila mauritiana and Drosophila simulans. Its origin is attributed to the duplication of the gene unc-4 in the subgenus Sophophora. By using a much larger sample of Drosophilidae species, we showed that contrary to what has been previously proposed, OdsH origin occurred 62 MYA. Evolutionary rates, expression, and transcription factor-binding sites of OdsH evidence that it may have rapidly experienced neofunctionalization in male sexual functions. Furthermore, the analysis of the OdsH peptide allowed the identification of mutations of D. mauritiana that could result in incompatibility in hybrids. In order to find if OdsH could be related to hybrid sterility, beyond Sophophora, we explored the expression of OdsH in Drosophila arizonae and Drosophila mojavensis, a pair of sister species with incomplete reproductive isolation. Our data indicated that OdsH expression is not atypical in their male-sterile hybrids. In conclusion, we have proposed that the origin of OdsH occurred earlier than previously proposed, followed by neofunctionalization. Our results also suggested that its role as a speciation gene might be restricted to D. mauritiana and D. simulans.

Intraspecific differentiation of Lindera obtusiloba as revealed by comparative plastomic and evolutionary analyses

Lindera obtusiloba Blume is the northernmost tree species in the family Lauraceae, and it is a key species in understanding the evolutionary history of this family. The species of L. obtusiloba in East Asia has diverged into the Northern and Southern populations, which are geographically separated by an arid belt. Though the morphological differences between populations have been observed and well documented, intraspecific variations at the plastomic level have not been systematically investigated to date. Here, ten chloroplast genomes of L. obtusiloba individuals were sequenced and analyzed along with three publicly available plastomes. Comparative plastomic analysis suggests that both the Northern and the Southern populations share similar overall structure, gene order, and GC content in their plastomes although the size of the plasome and the level of intraspecific variability do vary between the two populations. The Northern have relatively larger plastomes while the Southern population possesses higher intraspecific variability, which could be attributed to the complexity of the geological environments in the South. Phylogenomic analyses also support the split of the Northern and Southern clades among L. obtusiloba individuals. However, there is no obvious species boundary between var. obtusiloba and var. heterophylla in the Southern population, indicating that gene flow could still occur between these two varieties, and this could be used as a good example of reticulate evolution. It is also found that a few photosynthesis-related genes are under positive selection, which is mainly related to the geological and environmental differences between the Northern and the Southern regions. Our results provide a reference for phylogenetic analysis within species and suggest that phylogenomic analyses with a sufficient number of nuclear and chloroplast genomic target loci from widely distributed individuals could provide a deeper understanding of the population evolution of the widespread species.

Rapid evolution of genes with anti-cancer functions during the origins of large bodies and cancer resistance in elephants

Elephants have emerged as a model system to study the evolution of body size and cancer resistance because, despite their immense size, they have a very low prevalence of cancer. Previous studies have found that duplication of tumor suppressors at least partly contributes to the evolution of anti-cancer cellular phenotypes in elephants. Still, many other mechanisms must have contributed to their augmented cancer resistance. Here, we use a suite of codon-based maximum-likelihood methods and a dataset of 13,310 protein-coding gene alignments from 261 Eutherian mammals to identify positively selected and rapidly evolving elephant genes. We found 496 genes (3.73% of alignments tested) with statistically significant evidence for positive selection and 660 genes (4.96% of alignments tested) that likely evolved rapidly in elephants. Positively selected and rapidly evolving genes are statistically enriched in gene ontology terms and biological pathways related to regulated cell death mechanisms, DNA damage repair, cell cycle regulation, epidermal growth factor receptor (EGFR) signaling, and immune functions, particularly neutrophil granules and degranulation. All of these biological factors are plausibly related to the evolution of cancer resistance. Thus, these positively selected and rapidly evolving genes are promising candidates for genes contributing to elephant-specific traits, including the evolution of molecular and cellular characteristics that enhance cancer resistance.

Thirteen complete chloroplast genomes of the costaceae family: insights into genome structure, selective pressure and phylogenetic relationships

BACKGROUND

Costaceae, commonly known as the spiral ginger family, consists of approximately 120 species distributed in the tropical regions of South America, Africa, and Southeast Asia, of which some species have important ornamental, medicinal and ecological values. Previous studies on the phylogenetic and taxonomic of Costaceae by using nuclear internal transcribed spacer (ITS) and chloroplast genome fragments data had low resolutions. Additionally, the structures, variations and molecular evolution of complete chloroplast genomes in Costaceae still remain unclear. Herein, a total of 13 complete chloroplast genomes of Costaceae including 8 newly sequenced and 5 from the NCBI GenBank database, representing all three distribution regions of this family, were comprehensively analyzed for comparative genomics and phylogenetic relationships.

RESULT

The 13 complete chloroplast genomes of Costaceae possessed typical quadripartite structures with lengths from 166,360 to 168,966 bp, comprising a large single copy (LSC, 90,802 - 92,189 bp), a small single copy (SSC, 18,363 - 20,124 bp) and a pair of inverted repeats (IRs, 27,982 - 29,203 bp). These genomes coded 111 - 113 different genes, including 79 protein-coding genes, 4 rRNA genes and 28 - 30 tRNAs genes. The gene orders, gene contents, amino acid frequencies and codon usage within Costaceae were highly conservative, but several variations in intron loss, long repeats, simple sequence repeats (SSRs) and gene expansion on the IR/SC boundaries were also found among these 13 genomes. Comparative genomics within Costaceae identified five highly divergent regions including ndhF, ycf1-D2, ccsA-ndhD, rps15-ycf1-D2 and rpl16-exon2-rpl16-exon1. Five combined DNA regions (ycf1-D2 + ndhF, ccsA-ndhD + rps15-ycf1-D2, rps15-ycf1-D2 + rpl16-exon2-rpl16-exon1, ccsA-ndhD + rpl16-exon2-rpl16-exon1, and ccsA-ndhD + rps15-ycf1-D2 + rpl16-exon2-rpl16-exon1) could be used as potential markers for future phylogenetic analyses and species identification in Costaceae. Positive selection was found in eight protein-coding genes, including cemA, clpP, ndhA, ndhF, petB, psbD, rps12 and ycf1. Maximum likelihood and Bayesian phylogenetic trees using chloroplast genome sequences consistently revealed identical tree topologies with high supports between species of Costaceae. Three clades were divided within Costaceae, including the Asian clade, Costus clade and South American clade. Tapeinochilos was a sister of Hellenia, and Parahellenia was a sister to the cluster of Tapeinochilos + Hellenia with strong support in the Asian clade. The results of molecular dating showed that the crown age of Costaceae was about 30.5 Mya (95% HPD: 14.9 - 49.3 Mya), and then started to diverge into the Costus clade and Asian clade around 23.8 Mya (95% HPD: 10.1 - 41.5 Mya). The Asian clade diverged into Hellenia and Parahellenia at approximately 10.7 Mya (95% HPD: 3.5 - 25.1 Mya).

CONCLUSION

The complete chloroplast genomes can resolve the phylogenetic relationships of Costaceae and provide new insights into genome structures, variations and evolution. The identified DNA divergent regions would be useful for species identification and phylogenetic inference in Costaceae.

Multiple mechanisms explain loss of anthocyanins from betalain-pigmented Caryophyllales, including repeated wholesale loss of a key anthocyanidin synthesis enzyme

In this study, we investigate the genetic mechanisms responsible for the loss of anthocyanins in betalain-pigmented Caryophyllales, considering our hypothesis of multiple transitions to betalain pigmentation. Utilizing transcriptomic and genomic datasets across 357 species and 31 families, we scrutinize 18 flavonoid pathway genes and six regulatory genes spanning four transitions to betalain pigmentation. We examined evidence for hypotheses of wholesale gene loss, modified gene function, altered gene expression, and degeneration of the MBW (MYB-bHLH-WD40) trasnscription factor complex, within betalain-pigmented lineages. Our analyses reveal that most flavonoid synthesis genes remain conserved in betalain-pigmented lineages, with the notable exception of TT19 orthologs, essential for the final step in anthocyanidin synthesis, which appear to have been repeatedly and entirely lost. Additional late-stage flavonoid pathway genes upstream of TT19 also manifest strikingly reduced expression in betalain-pigmented species. Additionally, we find repeated loss and alteration in the MBW transcription complex essential for canonical anthocyanin synthesis. Consequently, the loss and exclusion of anthocyanins in betalain-pigmented species appear to be orchestrated through several mechanisms: loss of a key enzyme, downregulation of synthesis genes, and degeneration of regulatory complexes. These changes have occurred iteratively in Caryophyllales, often coinciding with evolutionary transitions to betalain pigmentation.

Evolution towards simplicity in bacterial small heat shock protein system

Evolution can tinker with multi-protein machines and replace them with simpler single-protein systems performing equivalent functions in an equally efficient manner. It is unclear how, on a molecular level, such simplification can arise. With ancestral reconstruction and biochemical analysis, we have traced the evolution of bacterial small heat shock proteins (sHsp), which help to refold proteins from aggregates using either two proteins with different functions (IbpA and IbpB) or a secondarily single sHsp that performs both functions in an equally efficient way. Secondarily single sHsp evolved from IbpA, an ancestor specialized in strong substrate binding. Evolution of an intermolecular binding site drove the alteration of substrate binding properties, as well as the formation of higher-order oligomers. Upon two mutations in the α-crystallin domain, secondarily single sHsp interacts with aggregated substrates less tightly. Paradoxically, less efficient binding positively influences the ability of sHsp to stimulate substrate refolding, since the dissociation of sHps from aggregates is required to initiate Hsp70-Hsp100-dependent substrate refolding. After the loss of a partner, IbpA took over its role in facilitating the sHsp dissociation from an aggregate by weakening the interaction with the substrate, which became beneficial for the refolding process. We show that the same two amino acids introduced in modern-day systems define whether the IbpA acts as a single sHsp or obligatorily cooperates with an IbpB partner. Our discoveries illuminate how one sequence has evolved to encode functions previously performed by two distinct proteins.

Developmental regulation of the male urogenital papilla in the male marine teleost black rockfish, Sebastes schlegelii (Hilgendorf, 1880)†

The male external genitalia of the black rockfish (Sebastes schlegelii Hilgendorf, 1880) is a fleshy protrusion known as the urogenital papilla (UGP), which functions to deliver sperm into the female reproductive tract for internal fertilization. It is not known which genes regulate the development of the UGP. The aim of this study was to identify key genes that regulate the development of the UGP in black rockfish and to determine the distribution of androgen receptor gene (ar) in the UGP. A total of 26 adult males and 560 juvenile fish were used in the experiment, in which we divided all normally developing juveniles into normal development and androgen groups. We added methyltestosterone solution (100 μg/l) to the androgen group-treated fish tank, soaked for 2 h per day for 38 days, and sampled 5~10 samples each time every 5 days during the culture process. Gene expression changes related to UGP were analyzed with tissue specificity between control and androgen groups during sex differentiation, adult male maturation, and the copulation stage (September to December) using real-time quantitative polymerase chain reaction. The expression of ar was also localized by two-color in situ hybridization in the UGP region of juvenile fish. Androgen treatment enhanced ar expression levels and the ar signal was stronger in the UGP region of both adult breeding fish and androgen-treated juvenile fish. This study provides insights into the regulation of the external genitalia of black rockfish and presents vital information for the artificial breeding of viviparous fish.

Comparative Chloroplast Genomics of 21 Species in Zingiberales with Implications for Their Phylogenetic Relationships and Molecular Dating

Zingiberales includes eight families and more than 2600 species, with many species having important economic and ecological value. However, the backbone phylogenetic relationships of Zingiberales still remain controversial, as demonstrated in previous studies, and molecular dating based on chloroplast genomes has not been comprehensively studied for the whole order. Herein, 22 complete chloroplast genomes from 21 species in Zingiberales were sequenced, assembled, and analyzed. These 22 genomes displayed typical quadripartite structures, which ranged from 161,303 bp to 163,979 bp in length and contained 111-112 different genes. The genome structures, gene contents, simple sequence repeats, long repeats, and codon usage were highly conserved, with slight differences among these genomes. Further comparative analysis of the 111 complete chloroplast genomes of Zingiberales, including 22 newly sequenced ones and the remaining ones from the national center for biotechnology information (NCBI) database, identified three highly divergent regions comprising ccsA, psaC, and psaC-ndhE. Maximum likelihood and Bayesian inference phylogenetic analyses based on chloroplast genome sequences found identical topological structures and identified a strongly supported backbone of phylogenetic relationships. Cannaceae was sister to Marantaceae, forming a clade that was collectively sister to the clade of (Costaceae, Zingiberaceae) with strong support (bootstrap (BS) = 100%, and posterior probability (PP) = 0.99-1.0); Heliconiaceae was sister to the clade of (Lowiaceae, Strelitziaceae), then collectively sister to Musaceae with strong support (BS = 94-100%, and PP = 0.93-1.0); the clade of ((Cannaceae, Marantaceae), (Costaceae, Zingiberaceae)) was sister to the clade of (Musaceae, (Heliconiaceae, (Lowiaceae, Strelitziaceae))) with robust support (BS = 100%, and PP = 1.0). The results of divergence time estimation of Zingiberales indicated that the crown node of Zingiberales occurred approximately 85.0 Mya (95% highest posterior density (HPD) = 81.6-89.3 million years ago (Mya)), with major family-level lineages becoming from 46.8 to 80.5 Mya. These findings proved that chloroplast genomes could contribute to the study of phylogenetic relationships and molecular dating in Zingiberales, as well as provide potential molecular markers for further taxonomic and phylogenetic studies of Zingiberales.

The mitochondrial genome of the mountain wooly tapir, Tapirus pinchaque and a formal test of the effect of altitude on the adaptive evolution of mitochondrial protein coding genes in odd-toed ungulates

BACKGROUND

The harsh conditions of high-altitude environments are known to drive the evolution of physiological and morphological traits in endothermic animals. These conditions are expected to result in the adaptive evolution of protein coding genes encoded in mitochondrial genomes that are vital for the oxidative phosphorylation pathway. In this study, we formally tested for signatures of adaptive evolution on mitochondrial protein coding genes in Tapirus pinchaque and other odd-toed ungulates inhabiting high-elevation environments.

RESULTS

The AT-rich mitochondrial genome of T. pinchaque is 16,750 bp long. A phylomitogenomic analysis supports the monophyly of the genus Tapirus and families in the Perissodactyla. The ratio of non-synonymous to synonymous substitutions demonstrated that all mitochondrial genes undergo purifying selection in T. pinchaque and other odd ungulates living at high elevations. Over this negative background selection, Branch Models suggested that cox3 and nad6 might be undergoing stronger purifying selection than other mitochondrial protein coding genes. Furthermore, Site Models suggested that one and four sites in nad2 and nad5, respectively, could be experiencing positive selection. However, these results were supported by Likelihood Ratio Tests but not Bayesian Empirical Bayes posterior probabilities. Additional analyses (in DataMonkey) indicated a relaxation of selection strength in nad6, evidence of episodic diversifying selection in cob, and revealed episodic positive/diversifying selection signatures for two sites in nad1, and one site each in nad2 and nad4.

CONCLUSION

The mitochondrial genome of T. pinchaque is an important genomic resource for conservation of this species and this study contributes to the understanding of adaptive evolution of mitochondrial protein coding genes in odd-toed ungulates inhabiting high-altitude environments.

An Insight of Betula platyphylla SWEET Gene Family through Genome-Wide Identification, Expression Profiling and Function Analysis of BpSWEET1c under Cold Stress

SWEET proteins play important roles in plant growth and development, sugar loading in phloem and resistance to abiotic stress through sugar transport. In this study, 13 BpSWEET genes were identified from birch genome. Collinearity analysis showed that there were one tandem repeating gene pair (BpSWEET1b/BpSWEET1c) and two duplicative gene pairs (BpSWEET17a/BpSWEET17b) in the BpSWEET gene family. The BpSWEET gene promoter regions contained several cis-acting elements related to stress resistance, for example: hormone-responsive and low-temperature-responsive cis-elements. Analysis of transcriptome data showed that BpSWEET genes were highly expressed in several sink organs, and the most BpSWEET genes were rapidly up-regulated under cold stress. BpSWEET1c, which was highly expressed in cold stress, was selected for further analysis. It was found that BpSWEET1c was located on the cell membrane. After 6 h of 4 °C stress, sucrose content in the leaves and roots of transient overexpressed BpSWEET1c was significantly higher than that of the control. MDA content in roots was significantly lower than that of the control. These results indicate that BpSWEET1c may play a positive role in the response to cold stress by promoting the metabolism and transport of sucrose. In conclusion, 13 BpSWEET genes were identified from the whole genome level. Most of the SWEET genes of birch were expressed in the sink organs and could respond to cold stress. Transient overexpression of BpSWEET1c changed the soluble sugar content and improved the cold tolerance of birch.

Relaxation of Natural Selection in the Evolution of the Giant Lungfish Genomes

Nonadaptive hypotheses on the evolution of eukaryotic genome size predict an expansion when the process of purifying selection becomes weak. Accordingly, species with huge genomes, such as lungfish, are expected to show a genome-wide relaxation signature of selection compared with other organisms. However, few studies have empirically tested this prediction using genomic data in a comparative framework. Here, we show that 1) the newly assembled transcriptome of the Australian lungfish, Neoceratodus forsteri, is characterized by an excess of pervasive transcription, or transcriptional leakage, possibly due to suboptimal transcriptional control, and 2) a significant relaxation signature in coding genes in lungfish species compared with other vertebrates. Based on these observations, we propose that the largest known animal genomes evolved in a nearly neutral scenario where genome expansion is less efficiently constrained.

Adaptive evolution in virulence effectors of the rice blast fungus Pyricularia oryzae

Plant pathogens secrete proteins called effectors that target host cellular processes to promote disease. Recently, structural genomics has identified several families of fungal effectors that share a similar three-dimensional structure despite remarkably variable amino-acid sequences and surface properties. To explore the selective forces that underlie the sequence variability of structurally-analogous effectors, we focused on MAX effectors, a structural family of effectors that are major determinants of virulence in the rice blast fungus Pyricularia oryzae. Using structure-informed gene annotation, we identified 58 to 78 MAX effector genes per genome in a set of 120 isolates representing seven host-associated lineages. The expression of MAX effector genes was primarily restricted to the early biotrophic phase of infection and strongly influenced by the host plant. Pangenome analyses of MAX effectors demonstrated extensive presence/absence polymorphism and identified gene loss events possibly involved in host range adaptation. However, gene knock-in experiments did not reveal a strong effect on virulence phenotypes suggesting that other evolutionary mechanisms are the main drivers of MAX effector losses. MAX effectors displayed high levels of standing variation and high rates of non-synonymous substitutions, pointing to widespread positive selection shaping the molecular diversity of MAX effectors. The combination of these analyses with structural data revealed that positive selection acts mostly on residues located in particular structural elements and at specific positions. By providing a comprehensive catalog of amino acid polymorphism, and by identifying the structural determinants of the sequence diversity, our work will inform future studies aimed at elucidating the function and mode of action of MAX effectors.

Evolution of chemosensory and detoxification gene families across herbivorous Drosophilidae

Herbivorous insects are exceptionally diverse, accounting for a quarter of all known eukaryotic species, but the genomic basis of adaptations that enabled this dietary transition remains poorly understood. Many studies have suggested that expansions and contractions of chemosensory and detoxification gene families-genes directly mediating interactions with plant chemical defenses-underlie successful plant colonization. However, this hypothesis has been challenging to test because the origins of herbivory in many insect lineages are ancient (>150 million years ago (mya)), obscuring genomic evolutionary patterns. Here, we characterized chemosensory and detoxification gene family evolution across Scaptomyza, a genus nested within Drosophila that includes a recently derived (<15 mya) herbivore lineage of mustard (Brassicales) specialists and carnation (Caryophyllaceae) specialists, and several nonherbivorous species. Comparative genomic analyses revealed that herbivorous Scaptomyza has among the smallest chemosensory and detoxification gene repertoires across 12 drosophilid species surveyed. Rates of gene turnover averaged across the herbivore clade were significantly higher than background rates in over half of the surveyed gene families. However, gene turnover was more limited along the ancestral herbivore branch, with only gustatory receptors and odorant-binding proteins experiencing strong losses. The genes most significantly impacted by gene loss, duplication, or changes in selective constraint were those involved in detecting compounds associated with feeding on living plants (bitter or electrophilic phytotoxins) or their ancestral diet (fermenting plant volatiles). These results provide insight into the molecular and evolutionary mechanisms of plant-feeding adaptations and highlight gene candidates that have also been linked to other dietary transitions in Drosophila.

Natural selection shapes the evolution of SARS-CoV-2 Omicron in Bangladesh

Severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has evolved to give rise to a highly transmissive and immune-escaping variant of concern, known as Omicron. Many aspects of the evolution of SARS-CoV-2 and the driving forces behind the ongoing Omicron outbreaks remain unclear. Substitution at the receptor-binding domain (RBD) in the spike protein is one of the primary strategies of SARS-CoV-2 Omicron to hinder recognition by the host angiotensin-converting enzyme 2 (ACE2) receptor and avoid antibody-dependent defense activation. Here, we scanned for adaptive evolution within the SARS-CoV-2 Omicron genomes reported from Bangladesh in the public database GISAID (www.gisaid.org; dated 2 April 2023). The ratio of the non-synonymous (Ka) to synonymous (Ks) nucleotide substitution rate, denoted as ω, is an indicator of the selection pressure acting on protein-coding genes. A higher proportion of non-synonymous to synonymous substitutions (Ka/Ks or ω > 1) indicates positive selection, while Ka/Ks or ω near zero indicates purifying selection. An equal amount of non-synonymous and synonymous substitutions (Ka/Ks or ω = 1) refers to neutrally evolving sites. We found evidence of adaptive evolution within the spike (S) gene of SARS-CoV-2 Omicron isolated from Bangladesh. In total, 22 codon sites of the S gene displayed a signature of positive selection. The data also highlighted that the receptor-binding motif within the RBD of the spike glycoprotein is a hotspot of adaptive evolution, where many of the codons had ω > 1. Some of these adaptive sites at the RBD of the spike protein are known to be associated with increased viral fitness. The M gene and ORF6 have also experienced positive selection. These results suggest that although purifying selection is the dominant evolutionary force, positive Darwinian selection also plays a vital role in shaping the evolution of SARS-CoV-2 Omicron in Bangladesh.

E7 Oncogene HPV58 Variants Detected in Northeast Brazil: Genetic and Functional Analysis

Cervical cancer is associated with persistent infections by high-risk Human Papillomavirus (HPV) types that may have nucleotide polymorphisms and, consequently, different oncogenic potentials. Therefore, this study aimed to evaluate the genetic variability and structural effects of the E7 oncogene of HPV58 in cervical scraping samples from Brazilian women. The study was developed with patients from hospitals in the metropolitan area of Recife, PE, Brazil. The most frequent HPV types were, in descending order of abundance, HPV16, 31, and 58. Phylogenetic analysis demonstrated that the isolates were classified into sublineages A2, C1, and D2. Two positively selected mutations were found in E7: 63G and 64T. The mutations G41R, G63D, and T64A in the E7 protein reduced the stability of the protein structure. Utilizing an NF-kB reporter assay, we observed a decrease in the NK-kB pathway activity with the HPV58-E7 variant 54S compared to the WT E7. The other detected E7 HPV58 variants presented similar NF-kB pathway activity compared to the WT E7. In this study, it was possible to identify mutations that may interfere with the molecular interaction between the viral oncoproteins and host proteins.

Repeated evolution of similar phenotypes: Integrating comparative methods with developmental pathways

Repeated phenotypes, often referred to as 'homoplasies' in cladistic analyses, may evolve through changes in developmental processes. Genetic bases of recurrent evolution gained attention and have been studied in the past years using approaches that combine modern analytical phylogenetic tools with the stunning assemblage of new information on developmental mechanisms. In this review, we evaluated the topic under an integrated perspective, revisiting the classical definitions of convergence and parallelism and detailing comparative methods used to evaluate evolution of repeated phenotypes, which include phylogenetic inference, estimates of evolutionary rates and reconstruction of ancestral states. We provide examples to illustrate how a given methodological approach can be used to identify evolutionary patterns and evaluate developmental mechanisms associated with the intermittent expression of a given trait along the phylogeny. Finally, we address why repeated trait loss challenges strict definitions of convergence and parallelism, discussing how changes in developmental pathways might explain the high frequency of repeated trait loss in specific lineages.

Evolutionary Shortcuts via Multinucleotide Substitutions and Their Impact on Natural Selection Analyses

Inference and interpretation of evolutionary processes, in particular of the types and targets of natural selection affecting coding sequences, are critically influenced by the assumptions built into statistical models and tests. If certain aspects of the substitution process (even when they are not of direct interest) are presumed absent or are modeled with too crude of a simplification, estimates of key model parameters can become biased, often systematically, and lead to poor statistical performance. Previous work established that failing to accommodate multinucleotide (or multihit, MH) substitutions strongly biases dN/dS-based inference towards false-positive inferences of diversifying episodic selection, as does failing to model variation in the rate of synonymous substitution (SRV) among sites. Here, we develop an integrated analytical framework and software tools to simultaneously incorporate these sources of evolutionary complexity into selection analyses. We found that both MH and SRV are ubiquitous in empirical alignments, and incorporating them has a strong effect on whether or not positive selection is detected (1.4-fold reduction) and on the distributions of inferred evolutionary rates. With simulation studies, we show that this effect is not attributable to reduced statistical power caused by using a more complex model. After a detailed examination of 21 benchmark alignments and a new high-resolution analysis showing which parts of the alignment provide support for positive selection, we show that MH substitutions occurring along shorter branches in the tree explain a significant fraction of discrepant results in selection detection. Our results add to the growing body of literature which examines decades-old modeling assumptions (including MH) and finds them to be problematic for comparative genomic data analysis. Because multinucleotide substitutions have a significant impact on natural selection detection even at the level of an entire gene, we recommend that selection analyses of this type consider their inclusion as a matter of routine. To facilitate this procedure, we developed, implemented, and benchmarked a simple and well-performing model testing selection detection framework able to screen an alignment for positive selection with two biologically important confounding processes: site-to-site synonymous rate variation, and multinucleotide instantaneous substitutions.

Whole-genome and dispersed duplication, including transposed duplication, jointly advance the evolution of TLP genes in seven representative Poaceae lineages

BACKGROUND

In the evolutionary study of gene families, exploring the duplication mechanisms of gene families helps researchers understand their evolutionary history. The tubby-like protein (TLP) family is essential for growth and development in plants and animals. Much research has been done on its function; however, limited information is available with regard to the evolution of the TLP gene family. Herein, we systematically investigated the evolution of TLP genes in seven representative Poaceae lineages.

RESULTS

Our research showed that the evolution of TLP genes was influenced not only by whole-genome duplication (WGD) and dispersed duplication (DSD) but also by transposed duplication (TRD), which has been neglected in previous research. For TLP family size, we found an evolutionary pattern of progressive shrinking in the grass family. Furthermore, the evolution of the TLP gene family was at least affected by evolutionary driving forces such as duplication, purifying selection, and base mutations.

CONCLUSIONS

This study presents the first comprehensive evolutionary analysis of the TLP gene family in grasses. We demonstrated that the TLP gene family is also influenced by a transposed duplication mechanism. Several new insights into the evolution of the TLP gene family are presented. This work provides a good reference for studying gene evolution and the origin of duplication.

Convergent evolutionary shifts in rhodopsin retinal release explain shared opsin repertoires in monotremes and crocodilians

The visual ecology of early mammals remains poorly resolved. Studies of ancestral photopigments suggest an ancient transition from nocturnal to more crepuscular conditions. By contrast, the phenotypic shifts following the split of monotremes and therians-which lost their SWS1 and SWS2 opsins, respectively-are less clear. To address this, we obtained new phenotypic data on the photopigments of extant and ancestral monotremes. We then generated functional data for another vertebrate group that shares the same photopigment repertoire as monotremes: the crocodilians. By characterizing resurrected ancient pigments, we show that the ancestral monotreme underwent a dramatic acceleration in its rhodopsin retinal release rate. Moreover, this change was likely mediated by three residue replacements, two of which also arose on the ancestral branch of crocodilians, which exhibit similarly accelerated retinal release. Despite this parallelism in retinal release, we detected minimal to moderate changes in the spectral tuning of cone visual pigments in these groups. Our results imply that ancestral forms of monotremes and crocodilians independently underwent niche expansion to encompass quickly changing light conditions. This scenario-which accords with reported crepuscular activity in extant monotremes-may help account for their loss of the ultraviolet-sensitive SWS1 pigment but retention of the blue-sensitive SWS2.