Phenotypic divergence of Glossina morsitans (Diptera: Glossinidae) populations in Zambia: Application of landmark-based wing geometric morphometrics to discriminate population-level variation.

An important consequence of the discontinuous distribution of insect populations within their geographic range is phenotypic divergence. Detection of this divergence can be challenging when it occurs through subtle shifts in morphological traits with complex geometries, such as insect wing venation. Here, we used landmark-based wing geometric morphometrics to investigate the population-level phenotypic variation of the two subspecies of Glossina morsitans, G. m. centralis Machado and G. m. morsitans Westwood that occur in Zambia. Twelve homologous landmarks digitised on the right wings of 720 specimens collected from four and five sites (80 per site with 1:1 sex ratio) within the G. m. centralis and G. m. morsitans range respectively, were subjected to generalised Procrustes analysis to obtain wing centroid size (CS) and wing shape variables. Linear permutation models and redundancy analysis were then used to compare CS and wing shape between male and female G. morsitans, the two subspecies G. m. centralis and G. m. morsitans, the sexes of each subspecies and between sample locations within each subspecies range, respectively. Significant differences in CS and wing shape were observed between G. morsitans sexes, subspecies and sample locations within each subspecies range. A neighbour-joining cladogram derived from the analysis of Procrustes distances showed that tsetse within each subspecies range were highly divergent. We conclude that G. morsitans populations in Zambia exhibit significant population-level variation in fly size and wing shape which suggests high levels of population structuring. The main drivers of this structuring could be random genetic drift in G. m. centralis demes and local adaptation to environmental conditions in G. m. morsitans populations. We therefore recommend molecular studies to estimate the levels of gene flow between these populations and identify possible barriers to genetic flow.

Influence of straight and inclined baffles on enhancement of battery thermal management system performance.

Battery thermal management systems (BTMSs) are used in electric vehicles (EVs) to regulate the heat generated by batteries while in use. Existing research on the improvement of BTMSs revealed that more research can be done by exploring different strategies to extend the service life and capabilities of EV batteries beyond the current limitations. In this study, effects of orientation of baffles; straight ( 90 ° ), inclined ( 60 ° and 30 ° ), height of baffles, thickness of baffles, positioning of baffles and number of baffles on the performance of conventional Z - Type BTMS were studied using Computational Fluid Dynamics (CFD) approach. The CFD approach was validated with existing experimental result from literature. Findings from the study showed that for straight baffle, the maximum temperature ( T max ) reduces as the height of baffles decreases. Furthermore, increasing the baffle thickness from 1 mm to 2 mm, produced reduction in T max by 0.26 K. For inclined baffles, by comparing the BTMS without baffles and BTMS with 2 baffles, T max and maximum temperature difference ( Δ T max ) reduced by 0.82 K and 0.68 K, respectively at angle 60 ° , and reduced by 1.65 K and 1.43 K, respectively at angle 30 ° . The BTMS with 2 baffles, inclined at angle 60 ° , with height of 6 mm and thickness of 1 mm, yielded the optimum T max (Lowest) value of 333.15 K, a reduction by 2.65 K when compared to the conventional Z - Type BTMS. This was also accompanied with a slight increase in Δ P by 1.12 Pa. In conclusion, it can be said that findings from this study will be beneficial in enhancing the design of BTMSs through adequate selection and utilization of baffles orientation.

Chickadees sing different songs in sympatry versus allopatry.

Character displacement theory predicts that closely-related co-occurring species should diverge in relevant traits to reduce costly interspecific interactions such as competition or hybridization. While many studies document character shifts in sympatry, few provide corresponding evidence that these shifts are driven by the costs of co-occurrence. Black-capped (Poecile atricapillus) and mountain chickadees (P. gambeli) are closely-related, ecologically similar, and broadly distributed songbirds with both allopatric and sympatric populations. In sympatry, both species appear to suffer costs of their co-occurrence: 1) both species are in worse body condition compared to allopatry and 2) hybridization sometimes yields sterile offspring. Here, we explored character displacement in the songs of black-capped and mountain chickadees by characterizing variation in male songs from sympatric and allopatric populations. We found that mountain chickadees sing differently in sympatry versus allopatry. Specifically, they produced more notes per song, were more likely to include an extra introductory note, and produced a smaller glissando in their first notes compared to all other populations. Combined with previous research on social dominance and maladaptive hybridization between black-capped and mountain chickadees, we posit that differences in sympatric mountain chickadee song are population-wide shifts to reduce aggression from dominant black-capped chickadees and/or prevent maladaptive hybridization.

Comparison of organelle genomes between endangered mangrove plant Dolichandrone spathacea to terrestrial relative provides insights into its origin and adaptative evolution.

Dolichandrone spathacea is a mangrove associate with high medicinal and ecological values. However, due to the dual-pressure of climate change and human activities, D. spathacea has become endangered in China. Moreover, misidentification between D. spathacea and its terrestrial relative D. cauda-felina poses further challenges to field protection and proper medicinal usage of D. spathacea. Thus, to address these problems, we sequenced and assembled mitochondrial (mt) and chloroplast (cp) genomes for both D. spathacea and D. cauda-felina. Comparative analysis revealed apparently different size and scaffold number between the two mt genomes, but a high similarity between the cp genomes. Eight regions with high sequence divergence were identified between the two cp genomes, which might be used for developing candidate DNA markers for distinguishing the two species. The splitting between D. spathacea and D. cauda-felina was inferred to occur at ~6.8 - 7.7 million years ago (Mya), which may be driven by the environment fluctuations in late Miocene. In the cp genome, 12 genes related to the expression of photosynthesis-associated proteins were detected with signatures of positive selection, which may contribute to the origin and evolutionary adaptation of Dolichandrone mangrove species. These new findings do not only enrich organelle genomic resources of Dolichandrone species, but also provide important genetic clues for improving the conservation and proper usage of endangered mangrove associate D. spathacea.

Decrypting Corals: Does Regulatory Evolution Underlie Environmental Specialisation of Coral Cryptic Lineages?

A recent sequencing study has shown that two common Caribbean corals, Montastraea cavernosa and Siderastrea siderea, each consist of four genetically distinct lineages in the Florida Keys. These lineages are specialised to a certain depth and, to a lesser extent, to nearshore or offshore habitats. We hypothesised that the lineages' environmental specialisation is at least in part due to regulatory evolution, which would manifest as the emergence of groups of coregulated genes ('modules') demonstrating lineage-specific responses to different reef environments. Our hypothesis also predicted that genes belonging to such modules would show greater genetic divergence between lineages than other genes. Contrary to these expectations, responses of cryptic lineages to natural environmental variation were essentially the same at the genome-wide gene coexpression network level, with much fewer differences in gene expression between lineages compared to between habitats. Moreover, none of the identified coregulated gene expression modules exhibit elevated genetic divergence between lineages. Possible explanations of these unexpected results range from the leading role of algal symbionts and/or microbiome in adaptation to strong action of spatially varying selection equalising gene expression patterns despite different genetic background. We discuss how future studies could assist in discriminating between these possibilities.

Combining intraprocedural CT guided navigation with ventilatory strategy for atelectasis: A modified electromagnetic navigation bronchoscopy.

CT-body divergence limits the accuracy of electromagnetic navigation bronchoscopy(ENB) in peripheral lung lesions diagnosis. Finding one effective and safe method to eliminating CT-body divergence may improve ENB accuracy. Thus, we developed a modified ENB which combining intraprocedural CT guided Navigation with Ventilatory strategy for Atelectasis to eliminate CT-body divergence. We called it inCTNVA-ENB. We present the case of an 80-year-old female with peripheral pulmonary nodule (without bronchial direct connection). She underwent inCTNVA-ENB, and the navigation probe accurately reached 6mm next to the target lesion without complications. The operation time was 42 minutes, and rapid on-site evaluation showed adenocarcinoma cells. CT data revealed the CT-body divergence caused by atelectasis was reduced.

Deep genetic substructure within bonobos.

Establishing the genetic and geographic structure of populations is fundamental, both to understand their evolutionary past and preserve their future. Nevertheless, the patterns of genetic population structure are unknown for most endangered species. This is the case for bonobos (Pan paniscus), which, together with chimpanzees (Pan troglodytes), are humans' closest living relatives. Chimpanzees live across equatorial Africa and are classified into four subspecies,1 with some genetic population substructure even within subspecies. Conversely, bonobos live exclusively in the Democratic Republic of Congo and are considered a homogeneous group with low genetic diversity,2 despite some population structure inferred from mtDNA. Nevertheless, mtDNA aside, their genetic structure remains unknown, hampering our understanding of the species and conservation efforts. Mapping bonobo genetic diversity in space is, however, challenging because, being endangered, only non-invasive sampling is possible for wild individuals. Here, we jointly analyze the exomes and mtDNA from 20 wild-born bonobos, the whole genomes of 10 captive bonobos, and the mtDNA of 136 wild individuals. We identify three genetically distinct bonobo groups of inferred Central, Western, and Far-Western geographic origin within the bonobo range. We estimate the split time between the central and western populations to be ∼145,000 years ago and genetic differentiation to be in the order of that of the closest chimpanzee subspecies. Furthermore, our estimated long-term Ne for Far-West (∼3,000) is among the lowest estimated for any great ape lineage. Our results highlight the need to attend to the bonobo substructure, both in terms of research and conservation.

A time-calibrated phylogeny of the diversification of Holoadeninae frogs.

The phylogeny of the major lineages of Amphibia has received significant attention in recent years, although evolutionary relationships within families remain largely neglected. One such overlooked group is the subfamily Holoadeninae, comprising 73 species across nine genera and characterized by a disjunct geographical distribution. The lack of a fossil record for this subfamily hampers the formulation of a comprehensive evolutionary hypothesis for their diversification. Aiming to fill this gap, we inferred the phylogenetic relationships and divergence times for Holoadeninae using molecular data and calibration information derived from the fossil record of Neobatrachia. Our inferred phylogeny confirmed most genus-level associations, and molecular dating analysis placed the origin of Holoadeninae in the Eocene, with subsequent splits also occurring during this period. The climatic and geological events that occurred during the Oligocene-Miocene transition were crucial to the dynamic biogeographical history of the subfamily. However, the wide highest posterior density intervals in our divergence time estimates are primarily attributed to the absence of Holoadeninae fossil information and, secondarily, to the limited number of sampled nucleotide sites.

A new computational method to estimate adaptation time in Avicennia by using divergence time.

Evolutionary studies of plant groups which are distributed in vast geographical regions and face different ecological and environmental conditions are important as they through light on different mechanisms of local adaptation and species divergence through time. The genus Avicennia is one of these plant groups which inspire of few species show interesting geographical distribution with some degree of species-specific geographical isolations. In general, very limited molecular phylogenetic investigations have been carried out in the genus Avicennia, and therefore we conducted the present study with the following aims: 1. To estimate the species divergence time based on different nuclear and chloroplast DNA regions, separately. This will illustrate how different genetic regions evolved in this genus, 2. To identify the sequences with potential adaptive value against geographical variable by latent factor mixed models (LFMM) analysis, 3. To illustrate the phylogenetic signal of these DNA regions and their role in speciation within the genus and, 4. To introducing a new computational strategy for estimating adaptive time for the sequences. The results showed that different genetic regions may produce different species divergent time, both the nuclear ribosomal internal transcribed spacer (ITS) region and chloroplast DNA sequences, contained potentially adaptive single nucleotide polymorphisms (SNPs).We could present a suggestive time for these adaptive sequences for the first time. In conclusion both local adaptation and independent mutations seem to have played role in Avicennia speciation and evolution.

Sequence comparison of the mitochondrial genomes of five caridean shrimps of the infraorder Caridea: phylogenetic implications and divergence time estimation.

BACKGROUND: The Caridea, affiliated with Malacostraca, Decapoda, and Pleocyemata, constitute one of the most significant shrimp groups. They are widely distributed across diverse aquatic habitats worldwide, enriching their evolutionary history. In recent years, considerable attention has been focused on the classification and systematic evolution of Caridea, yet controversies still exist regarding the phylogenetic relationships among families. METHODS: Here, the complete mitochondrial genome (mitogenome) sequences of five caridean species, namely Heterocarpus sibogae, Procletes levicarina, Macrobrachium sp., Latreutes anoplonyx, and Atya gabonensis, were determined using second-generation high-throughput sequencing technology. The basic structural characteristics, nucleotide composition, amino acid content, and codon usage bias of their mitogenomes were analyzed. Selection pressure values of protein-coding genes (PCGs) in species within the families Pandalidae, Palaemonidae, and Atyidae were also computed. Phylogenetic trees based on the nucleotide and amino acid sequences of 13 PCGs from 103 caridean species were constructed, and divergence times for various families within Caridea were estimated. RESULTS: The mitogenome of these five caridean species vary in length from 15,782 to 16,420 base pairs, encoding a total of 37 or 38 genes, including 13 PCGs, 2 rRNA genes, and 22 or 23 tRNA genes. Specifically, L. anoplonyx encodes an additional tRNA gene, bringing its total gene count to 38. The base composition of the mitogenomes of these five species exhibited a higher proportion of adenine-thymine (AT) bases. Six start codons and four stop codons were identified across the five species. Analysis of amino acid content and codon usage revealed variations among the five species. Analysis of selective pressure in Pandalidae, Palaemonidae, and Atyidae showed that the Ka/Ks values of PCGs in all three families were less than 1, indicating that purifying selection is influencing on their evolution. Phylogenetic analysis revealed that each family within Caridea is monophyletic. The results of gene rearrangement and phylogenetic analysis demonstrated correlations between these two aspects. Divergence time estimation, supported by fossil records, indicated that the divergence of Caridea species occurred in the Triassic period of the Mesozoic era, with subsequent differentiation into two major lineages during the Jurassic period. CONCLUSIONS: This study explored the fundamental characteristics and phylogenetic relationships of mitogenomes within the infraorder Caridea, providing valuable insights into their classification, interspecific evolutionary patterns, and the evolutionary status of various Caridea families. The findings provide essential references for identifying shrimp species and detecting significant gene rearrangements within the Caridea infraorder.

Local inconsistency detection using the Kullback-Leibler divergence measure.

BACKGROUND: The standard approach to local inconsistency assessment typically relies on testing the conflict between the direct and indirect evidence in selected treatment comparisons. However, statistical tests for inconsistency have low power and are subject to misinterpreting a p-value above the significance threshold as evidence of consistency. METHODS: We propose a simple framework to interpret local inconsistency based on the average Kullback-Leibler divergence (KLD) from approximating the direct with the corresponding indirect estimate and vice versa. Our framework uses directly the mean and standard error (or posterior mean and standard deviation) of the direct and indirect estimates obtained from a local inconsistency method to calculate the average KLD measure for selected comparisons. The average KLD values are compared with a semi-objective threshold to judge the inconsistency as acceptably low or material. We exemplify our novel interpretation approach using three networks with multiple treatments and multi-arm studies. RESULTS: Almost all selected comparisons in the networks were not associated with statistically significant inconsistency at a significance level of 5%. The proposed interpretation framework indicated 14%, 66%, and 75% of the selected comparisons with an acceptably low inconsistency in the corresponding networks. Overall, information loss was more notable when approximating the posterior density of the indirect estimates with that of the direct estimates, attributed to indirect estimates being more imprecise. CONCLUSIONS: Using the concept of information loss between two distributions alongside a semi-objectively defined threshold helped distinguish target comparisons with acceptably low inconsistency from those with material inconsistency when statistical tests for inconsistency were inconclusive.

A comparison of basal and activity-dependent exon splicing in cortical-patterned neurons of human and mouse origin.

Rodent studies have shown that alternative splicing in neurons plays important roles in development and maturity, and is regulatable by signals such as electrical activity. However, rodent-human similarities are less well explored. We compared basal and activity-dependent exon splicing in cortical-patterned human ESC-derived neurons with that in cortical mouse ESC-derived neurons, primary mouse cortical neurons at two developmental stages, and mouse hippocampal neurons, focussing on conserved orthologous exons. Both basal exon inclusion levels and activity-dependent changes in splicing showed human-mouse correlation. Conserved activity regulated exons are enriched in RBFOX, SAM68, NOVA and PTBP targets, and centered on cytoskeletal organization, mRNA processing, and synaptic signaling genes. However, human-mouse correlations were weaker than inter-mouse comparisons of neurons from different brain regions, developmental stages and origin (ESC vs. primary), suggestive of some inter-species divergence. The set of genes where activity-dependent splicing was observed only in human neurons were dominated by those involved in lipid biosynthesis, signaling and trafficking. Study of human exon splicing in mouse Tc1 neurons carrying human chromosome-21 showed that neuronal basal exon inclusion was influenced by cis-acting sequences, although may not be sufficient to confer activity-responsiveness in an allospecific environment. Overall, these comparisons suggest that neuronal alternative splicing should be confirmed in a human-relevant system even when exon structure is evolutionarily conserved.

Fine-scale genetic structure and phenotypic divergence of a passerine bird population inhabiting a continuous Mediterranean woodland.

Genetic differentiation between populations inhabiting ecologically different habitats might appear because of limited dispersal and gene flow, which may lead to patterns of phenotypic divergence and local adaptation. In this study, we use dispersal, genotypic (24 microsatellite loci) and phenotypic (body size and clutch size) data to analyse patterns of genetic structuring and phenotypic divergence in a blue tit (Cyanistes caeruleus) population inhabiting a continuous and heterogeneous woodland along a valley. The two slopes of the valley differ in their forest formations and environmental conditions. Findings showed that most blue tits reproduced within their natal slope. Accordingly, microsatellite analyses revealed that populations of blue tits established in the two slopes show subtle genetic differentiation. The two genetic populations diverged in clutch size, exceeding the level of differentiation expected based on genetic drift, hence suggesting divergent selection (or other processes promoting divergence) on this life-history trait. Our findings reveal that restricted dispersal and spatial heterogeneity may lead to genetic differentiation among bird populations at a surprisingly small scale. In this respect, it is worth highlighting that such differentiation occurs for an organism with high dispersal capacity and within a continuous woodland. Moreover, we show that small-scale ecological differences, together with limited gene flow, can result in selection favouring different phenotypes even within the same continuum population.

Long genetic and social isolation in Neanderthals before their extinction.

Neanderthal genomes have been recovered from sites across Eurasia, painting an increasingly complex picture of their populations' structure that mostly indicates that late European Neanderthals belonged to a single metapopulation with no significant evidence of population structure. Here, we report the discovery of a late Neanderthal individual, nicknamed "Thorin," from Grotte Mandrin in Mediterranean France, and his genome. These dentognathic fossils, including a rare example of distomolars, are associated with a rich archeological record of Neanderthal final technological traditions in this region ∼50-42 thousand years ago. Thorin's genome reveals a relatively early divergence of ∼105 ka with other late Neanderthals. Thorin belonged to a population with a small group size that showed no genetic introgression with other known late European Neanderthals, revealing some 50 ka of genetic isolation of his lineage despite them living in neighboring regions. These results have important implications for resolving competing hypotheses about causes of the disappearance of the Neanderthals.

The role of mate competition in speciation and divergence: a systematic review.

Competition for mates can play a critical role in determining reproductive success, shaping phenotypic variation within populations, and influencing divergence. Yet, studies of the role of sexual selection in divergence and speciation have focused disproportionately on mate choice. Here, we synthesize the literature on how mate competition may contribute to speciation and integrate concepts from work on sexual selection within populations - mating systems, ecology, and mate choice. Using this synthesis, we generate testable predictions for how mate competition may contribute to divergence. Then, we identify the extent of existing support for these predictions in the literature with a systematic review of the consequences of mate competition for population divergence across a range of evolutionary, ecological, and geographic contexts. We broadly evaluate current evidence, identify gaps in available data and hypotheses that need testing, and outline promising directions for future work. A major finding is that mate competition may commonly facilitate further divergence after initial divergence has occurred, e.g., upon secondary contact and between allopatric populations. Importantly, current hypotheses for how mate competition contributes to divergence do not fully explain observed patterns. While results from many studies fit predictions of negative frequency dependent selection, agonistic character displacement, and ecological selection, results from ~30% studies did not fit existing conceptual models. This review identifies future research aims for scenarios in which mate competition is likely important but has been understudied, including how ecological context and interactions between mate choice and mate competition can facilitate or hinder divergence and speciation.

Human leukocyte antigen evolutionary divergence as a novel risk factor for donor selection in acute lymphoblastic leukemia patients undergoing haploidentical hematopoietic stem cell transplantation.

Introduction: The human leukocyte antigen (HLA) evolutionary divergence (HED) reflects immunopeptidome diversity and has been shown to predict the response of tumors to immunotherapy. Its impact on allogeneic hematopoietic stem cell transplantation (HSCT) is controversial in different studies. Methods: In this study, we retrospectively analyzed the clinical impact of class I and II HED in 225 acute lymphoblastic leukemia patients undergoing HSCT from related haploidentical donors. The HED for recipient, donor, and donor-recipient pair was calculated based on Grantham distance, which accounts for variations in the composition, polarity, and volume of each amino acid within the peptide-binding groove of two HLA alleles. The median value of HED scores was used as a cut-off to stratify patients with high or low HED. Results: The class I HED for recipient (R_HEDclass I) showed the strongest association with cumulative incidence of relapse (12.2 vs. 25.0%, P = 0.00814) but not with acute graft-versus-host disease. The patients with high class II HED for donor-recipient (D/R_HEDclass II) showed a significantly higher cumulative incidence of severe aGVHD than those with low D/R_HEDclass II (24.0% vs. 6.1%, P = 0.0027). Multivariate analysis indicated that a high D/R_HEDclass II was an independent risk factor for the development of severe aGVHD (P = 0.007), and a high R_HEDclass I had a more than two-fold reduced risk of relapse (P = 0.028). However, there was no discernible difference in overall survival (OS) or disease-free survival (DFS) for patients with high or low HED, which was inconsistent with the previous investigation. Discussion: While the observation are limited by the presented single center retrospective cohort, the results show that HED has poor prognostic value in OS or DFS, as well as the associations with relapse and aGVHD. In haploidentical setting, class II HED for donor-recipient pair (D/R_HEDclass II) is an independent and novel risk factor for finding the best haploidentical donor, which could potentially influence clinical practice if verified in larger cohorts.

Multivariate analysis of phenotypic diversity elite bread wheat (Triticum aestivum L.) genotypes from ICARDA in Ethiopia.

Wheat is an important crop for food security, providing a source of protein and energy for the growing population in Ethiopia. However, both biotic and abiotic factors limit national wheat productivity. The availability of genetically diverse wheat genotypes is crucial for developing new wheat varieties that are both high-yielding and resilient to stress. Therefore, this field trial aimed to assess phenotypic variation and relationship among ICARDA-derived bread wheat genotypes using multivariate analysis techniques. The trial was conducted at three locations: Enewari, Wogere, and Kulumsa using an alpha lattice design with two replications during the main cropping seasons of 2022 and 2023. Phenotypic data on eight agronomic traits and the severity of yellow rust were collected and R programming was used for data analysis. Individual and combined location data analysis of variance showed significant differences (p ≤ 0.05) among genotypes for most of the studied traits. The highest heritability and genetic advance as a percentage of the mean were observed in days to heading (90.8, 21.29), plant height (72.4, 28.6), seeds per spike (61.7, 28), thousand kernel weight (61.9, 12), and area under the disease progress curve (67, 39.8), suggesting a predominance of additive gene action. Grain yield showed a strong positive correlation with days to maturity, plant height, spike length, spikelet per spike, and thousand kernel weight for each location. Dendrogram and phylogenetic tree methods were used to group genotypes into four genetically distinct clusters. Cluster II and III had the greatest inter-cluster distance, indicating higher diversity among their genotypes. This study identified new candidate genotypes with superior agronomic performance, high grain yield traits, and robust resistance to yellow rust, making them valuable for both current and future wheat breeding programs. Additionally, the comprehensive dataset produced in this study could facilitate the identification of genetic variations influencing desirable traits through genome-wide association analysis.

Epi-Clock: A sensitive platform to help understand pathogenic disease outbreaks and facilitate the response to future outbreaks of concern.

To predict potential epidemic outbreaks, we tested our strategy, Epi-Clock, which applies the novel ZHU algorithm to different SARS-CoV-2 datasets before outbreaks to search for significant mutational accumulation patterns correlated with outbreak events. Surprisingly, some inter-species genetic distances in Coronaviridae may represent intermediate states of different species or subspecies in the evolutionary history of Coronaviridae. The insertions and deletions in whole-genome sequences between different hosts were separately associated with important roles in host transmission and shifts in Coronaviridae. Furthermore, we believe that non-nucleosomal DNA may play a dominant role in the divergence of different lineages of SARS-CoV-2 in different regions of the world owing to the lack of nucleosome protection. We suggest that strong selective variation among different lineages of SARS-CoV-2 is required to produce strong codon usage bias, which appears in B.1.640.2 and B.1.617.2 (Delta). Notably, we found that an increasing number of other types of substitutions, such as those resulting from the hitchhiking effect, accumulated, especially in the pre-breakout phase, although some of the previous substitutions were replaced by other dominant genotypes. From most validations, we could accurately predict the potential pre-phase of outbreaks with a median interval of 5 days.

Mitochondrial phylogeography of grassland caterpillars (Lepidoptera: Lymantriinae: Gynaephora) endemic to the Qinghai-Tibetan plateau.

Grassland caterpillars (Lepidoptera: Lymantriinae: Gynaephora) are the most damaging pests to alpine meadows in the Qinghai-Tibetan Plateau (QTP). Here, we conducted extensive sampling from 39 geographic populations covering almost the entire distribution of the eight QTP Gynaephora (Hübner) species to investigate phylogeographic patterns and speciation based on two mitochondrial genes (COI and ND5). A total of 40 haplotypes were detected in the 39 populations, with >70% of all haplotypes not shared between populations. The monophyletic QTP Gynaephora migrated from non-QTP regions during the Pliocene, corresponding to the uplift of the QTP, suggesting a mode of transport into the QTP. Among the eight QTP Gynaephora species described by morphological characteristics, two species (G. alpherakii and G. menyuanensis) were recovered as monophyletic groups (Clades B and C), while the remaining six formed two monophyletic clades: Clade A (G. qinghaiensis, G. jiuzhiensis, and G. qumalaiensis) and Clade D (G. aureata, G. ruoergensis, and G. minora). These results suggested that the number of the QTP Gynaephora species may be overestimated and further studies based on both morphological and nuclear gene data are needed. Genetic differentiation and speciation of the QTP Gynaephora were likely driven by the QTP uplifts and associated climate fluctuations during the Pleistocene, indicated by divergence time estimation, suggesting that isolation and subsequent divergence was the dominant mode of speciation. The Sanjiangyuan region (i.e., Clade A, characterized by high genetic diversity) may have been a glacial refugium of the QTP Gynaephora, as supported by analyses of gene flow and biogeography. High levels of genetic diversity were found in QTP Gynaephora, without population expansion, which may explain the high-altitude adaptation and outbreaks of grassland caterpillars in alpine meadows of the QTP. This study provides the largest phylogeographic analysis of QTP Gynaephora and improves our understanding of the diversity and speciation of QTP insects.

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.