Drivers of genomic differentiation landscapes in populations of disparate ecological and geographical settings within mainland Apis cerana.

Elucidating the evolutionary processes that drive population divergence can enhance our understanding of the early stages of speciation and inform conservation management decisions. The honeybee Apis cerana displays extensive population divergence, providing an informative natural system for exploring these processes. The mainland lineage A. cerana includes several peripheral subspecies with disparate ecological and geographical settings radiated from a central ancestor. Under this evolutionary framework, we can explore the patterns of genome differentiation and the evolutionary models that explain them. We can also elucidate the contribution of non-genomic spatiotemporal mechanisms (extrinsic features) and genomic mechanisms (intrinsic features) that influence these genomic differentiation landscapes. Based on 293 whole genomes, a small part of the genome is highly differentiated between central-peripheral subspecies pairs, while low and partial parallelism partly reflects idiosyncratic responses to environmental differences. Combined elements of recurrent selection and speciation-with-gene-flow models generate the heterogeneous genome landscapes. These elements weight differently between central-island and other central-peripheral subspecies pairs, influenced by glacial cycles superimposed on different geomorphologies. Although local recombination rates exert a significant influence on patterns of genomic differentiation, it is unlikely that low-recombination rates regions were generated by structural variation. In conclusion, complex factors including geographical isolation, divergent ecological selection and non-uniform genome features have acted concertedly in the evolution of reproductive barriers that could reduce gene flow in part of the genome and facilitate the persistence of distinct populations within mainland lineage of A. cerana.

Microbial perspective on the skin-gut axis and atopic dermatitis.

Atopic dermatitis (AD) is a relapsing inflammatory skin condition that has become a global health issue with complex etiology and mounting prevalence. The association of AD with skin and gut microbiota has been revealed by virtue of the continuous development of sequencing technology and genomics analysis. Also, the gut-brain-skin axis and its mutual crosstalk mechanisms have been gradually verified. Accordingly, the microbiota-skin-gut axis also plays an important role in allergic skin inflammation. Herein, we reviewed the relationship between the microbiota-skin-gut axis and AD, explored the underlying signaling molecules and potential pathways, and focused on the potential mechanisms of probiotics, antimicrobial peptides (AMPs), coagulase-negative staphylococci transplantation, fecal microbiota transplantation, AMPs, and addition of essential fatty acids in alleviating AD, with the aim to provide a new perspective for targeting microbiota in the treatment of allergic skin inflammation.

Study on a Highly Thermostable Dy3+-Activated Borophosphate Phosphor.

Constructing a phosphor with multifunctional applications is an imperative challenge. Especially, highly thermostable luminescence of phosphor is indispensable for stable white-light-emitting diodes (LEDs). Nevertheless, good thermal quenching resistance behavior is unfavorable for a fluorescence intensity ratio (FIR)-based optical temperature sensor. Herein, a highly thermostable Ba3(ZnB5O10)PO4 (BZBP)-based phosphor is successfully achieved via replacing Ba2+ with Dy3+, demonstrating simultaneously promising lighting and thermometry utilizations. Under the excitation of 350 nm, the title phosphor only loses 12% of the initial intensity when the temperature is up to 473 K, ensuring sufficient luminescence thermostability for white-LED lighting. The white-LED device fabricated using the title phosphor emits high-quality white light with a high color rendering index (Ra = 93) and low correlated color temperature (CCT = 3996 K). Meanwhile, the yellow and blue emission intensities demonstrate a downtrend difference with rising temperature. Temperature sensing properties are assessed through FIR technology. The maximal relative sensitivity reaches as high as 0.0379 K-1 at 298 K. These results reveal that the title phosphor has a great potential for indoor lighting and thermometry applications.

Natural selection and genetic diversity maintenance in a parasitic wasp during continuous biological control application.

Aphidius gifuensis is a parasitoid wasp and primary endoparasitoid enemy of the peach potato aphid, Myzus persicae. Artificially reared, captive wasps of this species have been extensively and effectively used to control populations of aphids and limit crop loss. However, the consequences of large-scale releasing of captive A. gifuensis, such as genetic erosion and reduced fitness in wild populations of this species, remains unclear. Here, we sequence the genomes of 542 A. gifuensis individuals collected across China, including 265 wild and 277 human-intervened samples. Population genetic analyses on wild individuals recovered Yunnan populations as the ancestral group with the most complex genetic structure. We also find genetic signature of environmental adaptation during the dispersal of wild populations from Yunnan to other regions. While comparative genomic analyses of captive wasps revealed a decrease in genetic diversity during long-term rearing, population genomic analyses revealed signatures of natural selection by several biotic (host plants) or abiotic (climate) factors, which support maintenance of the gene pool of wild populations in spite of the introduction of captive wasps. Therefore, the impact of large-scale release is reduced. Our study suggests that A. gifuensis is a good system for exploring the genetic and evolutionary effects of mass rearing and release on species commonly used as biocontrol agents.

Chromosome-level genome of the poultry shaft louse Menopon gallinae provides insight into the host-switching and adaptive evolution of parasitic lice.

BACKGROUND: Lice (Psocodea: Phthiraptera) are one important group of parasites that infects birds and mammals. It is believed that the ancestor of parasitic lice originated on the ancient avian host, and ancient mammals acquired these parasites via host-switching from birds. Here we present the first chromosome-level genome of Menopon gallinae in Amblycera (earliest diverging lineage of parasitic lice). We explore the transition of louse host-switching from birds to mammals at the genomic level by identifying numerous idiosyncratic genomic variations. RESULTS: The assembled genome is 155 Mb in length, with a contig N50 of 27.42 Mb. Hi-C scaffolding assigned 97% of the bases to 5 chromosomes. The genome of M. gallinae retains a basal insect repertoire of 11,950 protein-coding genes. By comparing the genomes of lice to those of multiple representative insects in other orders, we discovered that gene families of digestion, detoxification, and immunity-related are generally conserved between bird lice and mammal lice, while mammal lice have undergone a significant reduction in genes related to chemosensory systems and temperature. This suggests that mammal lice have lost some of these genes through the adaption to environment and temperatures after host-switching. Furthermore, 7 genes related to hematophagy were positively selected in mammal lice, suggesting their involvement in the hematophagous behavior. CONCLUSIONS: Our high-quality genome of M. gallinae provides a valuable resource for comparative genomic research in Phthiraptera and facilitates further studies on adaptive evolution of host-switching within parasitic lice.

Population genomics of the muskox' resilience in the near absence of genetic variation.

Genomic studies of species threatened by extinction are providing crucial information about evolutionary mechanisms and genetic consequences of population declines and bottlenecks. However, to understand how species avoid the extinction vortex, insights can be drawn by studying species that thrive despite past declines. Here, we studied the population genomics of the muskox (Ovibos moschatus), an Ice Age relict that was at the brink of extinction for thousands of years at the end of the Pleistocene yet appears to be thriving today. We analysed 108 whole genomes, including present-day individuals representing the current native range of both muskox subspecies, the white-faced and the barren-ground muskox (O. moschatus wardi and O. moschatus moschatus) and a ~21,000-year-old ancient individual from Siberia. We found that the muskox' demographic history was profoundly shaped by past climate changes and post-glacial re-colonizations. In particular, the white-faced muskox has the lowest genome-wide heterozygosity recorded in an ungulate. Yet, there is no evidence of inbreeding depression in native muskox populations. We hypothesize that this can be explained by the effect of long-term gradual population declines that allowed for purging of strongly deleterious mutations. This study provides insights into how species with a history of population bottlenecks, small population sizes and low genetic diversity survive against all odds.

Genetic Modification of a Hox Locus Drives Mimetic Color Pattern Variation in a Highly Polymorphic Bumble Bee.

Müllerian mimicry provides natural replicates ideal for exploring mechanisms underlying adaptive phenotypic divergence and convergence, yet the genetic mechanisms underlying mimetic variation remain largely unknown. The current study investigates the genetic basis of mimetic color pattern variation in a highly polymorphic bumble bee, Bombus breviceps (Hymenoptera, Apidae). In South Asia, this species and multiple comimetic species converge onto local Müllerian mimicry patterns by shifting the abdominal setal color from orange to black. Genetic crossing between the orange and black phenotypes suggested the color dimorphism being controlled by a single Mendelian locus, with the orange allele being dominant over black. Genome-wide association suggests that a locus at the intergenic region between 2 abdominal fate-determining Hox genes, abd-A and Abd-B, is associated with the color change. This locus is therefore in the same intergenic region but not the same exact locus as found to drive red black midabdominal variation in a distantly related bumble bee species, Bombus melanopygus. Gene expression analysis and RNA interferences suggest that differential expression of an intergenic long noncoding RNA between abd-A and Abd-B at the onset setal color differentiation may drive the orange black color variation by causing a homeotic shift late in development. Analysis of this same color locus in comimetic species reveals no sequence association with the same color shift, suggesting that mimetic convergence is achieved through distinct genetic routes. Our study establishes Hox regions as genomic hotspots for color pattern evolution in bumble bees and demonstrates how pleiotropic developmental loci can drive adaptive radiations in nature.

Convergent molecular evolution of thermogenesis and circadian rhythm in Arctic ruminants.

The muskox and reindeer are the only ruminants that have evolved to survive in harsh Arctic environments. However, the genetic basis of this Arctic adaptation remains largely unclear. Here, we compared a de novo assembled muskox genome with reindeer and other ruminant genomes to identify convergent amino acid substitutions, rapidly evolving genes and positively selected genes among the two Arctic ruminants. We found these candidate genes were mainly involved in brown adipose tissue (BAT) thermogenesis and circadian rhythm. Furthermore, by integrating transcriptomic data from goat adipose tissues (white and brown), we demonstrated that muskox and reindeer may have evolved modulating mitochondrion, lipid metabolism and angiogenesis pathways to enhance BAT thermogenesis. In addition, results from co-immunoprecipitation experiments prove that convergent amino acid substitution of the angiogenesis-related gene hypoxia-inducible factor 2alpha (HIF2A), resulting in weakening of its interaction with prolyl hydroxylase domain-containing protein 2 (PHD2), may increase angiogenesis of BAT. Altogether, our work provides new insights into the molecular mechanisms involved in Arctic adaptation.

Chromosome-level genome assembly of bean flower thrips Megalurothrips usitatus (Thysanoptera: Thripidae).

Bean flower thrips Megalurothrips usitatus is a staple pest of cowpea and other legumes and causes dramatic economic losses. Its small size allows for easy concealment, and large reproductive capacity easily leads to infestations. Despite the importance of a genome in developing novel management strategies, genetic studies on M. usitatus remain limited. Thus, we generated a chromosome-level M. usitatus genome using a combination of PacBio long read and Hi-C technologies. The assembled genome was 238.14 Mb with a scaffold N50 of 13.85 Mb. The final genome was anchored into 16 pseudo-chromosomes containing 14,000 genes, of which 91.74% were functionally annotated. Comparative genomic analyses revealed that expanded gene families were enriched in fatty acid metabolism and detoxification metabolism (ABC transporters), and contracted gene families were strongly associated with chitin-based cuticle development and sensory perception of taste. In conclusion, this high-quality genome provides an invaluable resource for us to understand the thrips' ecology and genetics, contributing to pest management.

Defining honeybee subspecies in an evolutionary context warrants strategized conservation.

Despite the urgent need for conservation consideration, strategic action plans for the preservation of the Asian honeybee, Apis cerana Fabricius, 1793, remain lacking. Both the convergent and divergent adaptations of this widespread insect have led to confusing phenotypical traits and inconsistent infraspecific taxonomy. Unclear subspecies boundaries pose a significant challenge to honeybee conservation efforts, as it is difficult to effectively prioritize conservation targets without a clear understanding of subspecies identities. Here, we investigated genome variations in 362 worker bees representing almost all populations of mainland A. cerana to understand how evolution has shaped its population structure. Whole-genome single nucleotide polymorphisms (SNPs) based on nuclear sequences revealed eight putative subspecies, with all seven peripheral subspecies exhibiting mutually exclusive monophyly and distinct genetic divergence from the widespread central subspecies. Our results demonstrated that most classic morphological traits, including body size, were related to the climatic variables of the local habitats and did not reflect the true evolutionary history of the organism. Thus, such morphological traits were not suitable for subspecific delineation. Conversely, wing vein characters showed relative independence to the environment and supported the subspecies boundaries inferred from nuclear genomes. Mitochondrial phylogeny further indicated that the present subspecies structure was a result of multiple waves of population divergence from a common ancestor. Based on our findings, we propose that criteria for subspecies delineation should be based on evolutionary independence, trait distinction, and geographic isolation. We formally defined and described eight subspecies of mainland A. cerana. Elucidation of the evolutionary history and subspecies boundaries enables a customized conservation strategy for both widespread and endemic honeybee conservation units, guiding colony introduction and breeding.

Evolutionary analysis of a complete chicken genome.

Microchromosomes are prevalent in nonmammalian vertebrates [P. D. Waters et al., Proc. Natl. Acad. Sci. U.S.A. 118 (2021)], but a few of them are missing in bird genome assemblies. Here, we present a new chicken reference genome containing all autosomes, a Z and a W chromosome, with all gaps closed except for the W. We identified ten small microchromosomes (termed dot chromosomes) with distinct sequence and epigenetic features, among which six were newly assembled. Those dot chromosomes exhibit extremely high GC content and a high level of DNA methylation and are enriched for housekeeping genes. The pericentromeric heterochromatin of dot chromosomes is disproportionately large and continues to expand with the proliferation of satellite DNA and testis-expressed genes. Our analyses revealed that the 41-bp CNM repeat frequently forms higher-order repeats (HORs) at the centromeres of acrocentric chromosomes. The centromere core regions where the kinetochore attaches often encompass telomeric sequence (TTAGGG)n, and in a one of the dot chromosomes, the centromere core recruits an endogenous retrovirus (ERV). We further demonstrate that the W chromosome shares some common features with dot chromosomes, having large arrays of hypermethylated tandem repeats. Finally, using the complete chicken chromosome models, we reconstructed a fine picture of chordate karyotype evolution, revealing frequent chromosomal fusions before and after vertebrate whole-genome duplications. Our sequence and epigenetic characterization of chicken chromosomes shed insights into the understanding of vertebrate genome evolution and chromosome biology.

Genomes of the cosmopolitan fruit pest Bactrocera dorsalis (Diptera: Tephritidae) reveal its global invasion history and thermal adaptation.

INTRODUCTION: The oriental fruit fly Bactrocera dorsalis is one of the most destructive agricultural pests worldwide, with highly debated species delimitation, origin, and global spread routes. OBJECTIVES: Our study intended to (i) resolve the taxonomic uncertainties between B. dorsalis and B. carambolae, (ii) reveal the population structure and global invasion routes of B. dorsalis across Asia, Africa, and Oceania, and (iii) identify genomic regions that are responsible for the thermal adaptation of B. dorsalis. METHODS: Based on a high-quality chromosome-level reference genome assembly, we explored the population relationship using a genome-scale single nucleotide polymorphism dataset generated from the resequencing data of 487 B. dorsalis genomes and 25 B. carambolae genomes. Genome-wide association studies and silencing using RNA interference were used to identify and verify the candidate genes associated with extreme thermal stress. RESULTS: We showed that B. dorsalis originates from the Southern India region with three independent invasion and spread routes worldwide: (i) from Northern India to Northern Southeast Asia, then to Southern Southeast Asia; (ii) from Northern India to Northern Southeast Asian, then to China and Hawaii; and (iii) from Southern India toward the African mainland, then to Madagascar, which is mainly facilitated by human activities including trade and immigration. Twenty-seven genes were identified by a genome-wide association study to be associated with 11 temperature bioclimatic variables. The Cyp6a9 gene may enhance the thermal adaptation of B. dorsalis and thus boost its invasion, which tended to be upregulated at a hardening temperature of 38 °C. Functional verification using RNA interference silencing against Cyp6a9, led to the specific decrease in Cyp6a9 expression, reducing the survival rate of dsRNA-feeding larvae exposed to extreme thermal stress of 45 °C after heat hardening treatments in B. dorsalis. CONCLUSION: This study provides insights into the evolutionary history and genetic basis of temperature adaptation in B. dorsalis.

Identification and characterization of two DMD pedigrees with large inversion mutations based on a long-read sequencing pipeline.

Pathogenic large inversions are rarely reported on DMD gene due to the lack of effective detection methods. Here we report two DMD pedigrees and proposed a reliable pipeline to define large inversions in DMD patients. In the first pedigree, conventional approaches including multiplex ligation-dependent probe amplification, and whole-exome sequencing by next generation sequencing were failed to detect any pathologic variant. Then an advanced analysis pipeline which consists of RNA-seq, cDNA array capture sequencing, optical mapping, long-read sequencing was built. RNA-seq and cDNA capture sequencing showed a complete absence of transcripts of exons 3-55. Optical mapping identified a 55 Mb pericentric inversion between Xp21 and Xq21. Subsequently, long-read sequencing and Sanger sequencing determined the inversion breakpoints at 32,915,769 and 87,989,324 of the X chromosomes. In the second pedigree, long-read sequencing was directly conducted and Sanger sequencing was performed to verify the mutation. Long-read sequencing and Sanger sequencing found breakpoints at 32,581,576 and 127,797,236 on DMD gene directly. In conclusion, large inversion might be a rare but important mutation type in DMD gene. An effective pipeline was built in detecting large inversion mutations based on long-read sequencing platforms.

Evolution and diversification of Mountain voles (Rodentia: Cricetidae).

The systematics of the Cricetid genus Neodon have long been fraught with uncertainty due to sampling issues and a lack of comprehensive datasets. To gain better insights into the phylogeny and evolution of Neodon, we systematically sampled Neodon across the Hengduan and Himalayan Mountains, which cover most of its range in China. Analyses of skulls, teeth, and bacular structures revealed 15 distinct patterns corresponding to 15 species of Neodon. In addition to morphological analyses, we generated a high-quality reference genome for the mountain vole and generated whole-genome sequencing data for 47 samples. Phylogenomic analyses supported the recognition of six new species, revealing a long-term underestimation of Neodon diversity. We further identified positively selected genes potentially related to high-elevation adaptation. Together, our results illuminate how climate change caused the plateau to become the centre of Neodon origin and diversification and how mountain voles have adapted to the hypoxic high-altitude plateau environment.

The population genomic legacy of the second plague pandemic.

Human populations have been shaped by catastrophes that may have left long-lasting signatures in their genomes. One notable example is the second plague pandemic that entered Europe in ca. 1,347 CE and repeatedly returned for over 300 years, with typical village and town mortality estimated at 10%-40%.1 It is assumed that this high mortality affected the gene pools of these populations. First, local population crashes reduced genetic diversity. Second, a change in frequency is expected for sequence variants that may have affected survival or susceptibility to the etiologic agent (Yersinia pestis).2 Third, mass mortality might alter the local gene pools through its impact on subsequent migration patterns. We explored these factors using the Norwegian city of Trondheim as a model, by sequencing 54 genomes spanning three time periods: (1) prior to the plague striking Trondheim in 1,349 CE, (2) the 17th-19th century, and (3) the present. We find that the pandemic period shaped the gene pool by reducing long distance immigration, in particular from the British Isles, and inducing a bottleneck that reduced genetic diversity. Although we also observe an excess of large FST values at multiple loci in the genome, these are shaped by reference biases introduced by mapping our relatively low genome coverage degraded DNA to the reference genome. This implies that attempts to detect selection using ancient DNA (aDNA) datasets that vary by read length and depth of sequencing coverage may be particularly challenging until methods have been developed to account for the impact of differential reference bias on test statistics.

Hepatic AMPK signaling dynamic activation in response to REDOX balance are sentinel biomarkers of exercise and antioxidant intervention to improve blood glucose control.

Antioxidant intervention is considered to inhibit reactive oxygen species (ROS) and alleviate hyperglycemia. Paradoxically, moderate exercise can produce ROS to improve diabetes. The exact redox mechanism of these two different approaches remains largely unclear. Here, by comparing exercise and antioxidant intervention on type 2 diabetic rats, we found moderate exercise upregulated compensatory antioxidant capability and reached a higher level of redox balance in the liver. In contrast, antioxidant intervention achieved a low-level redox balance by inhibiting oxidative stress. Both of these two interventions could promote glucose catabolism and inhibit gluconeogenesis through activation of hepatic AMP-activated protein kinase (AMPK) signaling; therefore, ameliorating diabetes. During exercise, different levels of ROS generated by exercise have differential regulations on the activity and expression of hepatic AMPK. Moderate exercise-derived ROS promoted hepatic AMPK glutathionylation activation. However, excessive exercise increased oxidative damage and inhibited the activity and expression of AMPK. Overall, our results illustrate that both exercise and antioxidant intervention improve blood glucose control in diabetes by promoting redox balance, despite different levels of redox state(s). These results indicate that the AMPK signaling activation, combined with oxidative damage markers, could act as sentinel biomarkers, reflecting the threshold of redox balance that is linked to effective glucose control in diabetes. These findings provide theoretical evidence for the precise management of diabetes by antioxidants and exercise.

Molecules known as reactive oxygen species or ROS play vital roles in healthy cells. However, ROS can act as a double-edged sword: if their levels become too high, they can be harmful and interfere with many physiological processes. Indeed, diabetes, high blood pressure and many other chronic diseases are associated with imbalances in the levels of ROS in the body. To counter high ROS levels, cells have antioxidant mechanisms that reduce the excess ROS in the cell and keep the 'redox' (from reduction and oxidation) balance of the cell. Exercise and antioxidant nutritional supplements have attracted much attention as drug-free interventions for diabetes. Both strategies alter the levels of ROS in the body, with exercise increasing the levels of ROS, and antioxidant supplements reducing them. Individuals with diabetes and other metabolic health issues have different ROS levels depending on the severity of the disease, age, genetics and other factors, leading to different redox states in their cells. Thus, approaches that can accurately evaluate the redox balance status of individuals are necessary for clinicians to identify what types of exercise and antioxidant supplements are beneficial and which treatments are most appropriate for each patient. Wu, Zhao, Yan, Gao et al. examined the effects of exercise and antioxidant supplements on rats with diabetes, with the aim of identifying molecules ­ also known as biomarkers ­ that reflect the bodies' redox balance. They found that moderate exercise increased the levels of ROS in the liver, which, in turn, compensated by increasing the production of antioxidants to protect against the higher levels of ROS. This resulted in a healthy 'high-level' redox balance, in which both ROS and antioxidants levels were high in the rats. On the other hand, giving the rats antioxidant supplements decreased their levels of ROS, leading to a healthy low-level redox balance with low levels of ROS. These findings indicate that regular moderate exercise may be appropriate for people with pre-diabetes symptoms to restore a healthy redox balance. This is because the compensatory antioxidant mechanisms that kick in during exercise may be enough to counteract the excessive levels of ROS in these people. For patients with mild diabetes, exercise, antioxidant supplements, or a combination of both may be appropriate treatment, depending on their levels of ROS. Finally, patients with severe diabetes, who already have high levels of ROS, may benefit from antioxidant supplements to help reduce their excessive levels of ROS. In the future, the biomarkers identified by Wu, Zhao, Yan, Gao et al. may be used to monitor and assess the change in the redox balance status of various populations and guide personalized interventions to maintain health. Additionally, these findings provide a new strategy for precision prevention and treatment of diabetes and other metabolic diseases.

The IMD pathway in Hemipteran: A comparative analysis and discussion.

The evolutionary patterns of the genes in the IMD pathway in hemipterans were characterized and compared. The hemipteran insects were clustered into two groups. One group that encompasses whitefly, plant lice, and scale insect partially lacks the IMD pathway and all antimicrobial peptide (AMP) genes, with the vast majority of IMD pathway and all AMP genes being absent in aphids. The reasons for the absence of the IMD pathway and AMP genes in hemipterans were analyzed based on aphids, in terms of fitness costs. In case of limited resources, aphids have to make a trade-off between the necessary costs such as clean food sources, the essential amino acids supplied by primary bacterial symbionts for survival, nutrients and/or protections against stress provided by secondary symbionts, and the high reproductive capacity, and the costs that do not increase the fitness. Obviously, aphids have to abandon the strong immune system, especially the AMPs and IMD pathway which is mainly against Gram-negative bacteria. The common ground shared with aphids may be the reason for the absence of the IMD pathway and AMP genes in other hemipteran insects.

Genomics and morphometrics reveal the adaptive evolution of pikas.

Pikas (Lagomorpha: Ochotonidae) are small mouse-like lagomorphs. To investigate their adaptation to different ecological environments during their dispersal from the Qinghai-Xizang (Tibet) Plateau (QTP), we collected 226 pikas and measured 20 morphological characteristics and recorded habitat information. We also sequenced the genome of 81 specimens, representing 27 putative pika species. The genome-wide tree based on 4 090 coding genes identified five subgenera, i.e., Alienauroa, Conothoa, Lagotona, Ochotona, and Pika, consistent with morphometric data. Morphologically, Alienauroa and Ochotona had similar traits, including smaller size and earlier divergence time compared to other pikas. Consistently, the habitats of Alienauroa and Ochotona differed from those of the remaining subgenera. Phylogenetic signal analysis detected 83 genes significantly related to morphological characteristics, including several visual and hearing-related genes. Analysis of shared amino acid substitutions and positively selected genes (PSGs) in Alienauroa and Ochotona identified two genes, i.e., mitochondrial function-related TSFM (p.Q155E) and low-light visual sensitivity-related PROM1 (p.H419Y). Functional experiments demonstrated that TSFM-155E significantly enhanced mitochondrial function compared to TSFM-155Q in other pikas, and PROM1-419Y decreased the modeling of dynamic intracellular chloride efflux upon calcium uptake. Alienauroa and Ochotona individuals mostly inhabit different environments (e.g., subtropical forests) than other pikas, suggesting that a shift from the larger ancestral type and changes in sensory acuity and energy enhancement may have been required in their new environments. This study increases our understanding of the evolutionary history of pikas.

A super pan-genomic landscape of rice.

Pan-genomes from large natural populations can capture genetic diversity and reveal genomic complexity. Using de novo long-read assembly, we generated a graph-based super pan-genome of rice consisting of a 251-accession panel comprising both cultivated and wild species of Asian and African rice. Our pan-genome reveals extensive structural variations (SVs) and gene presence/absence variations. Additionally, our pan-genome enables the accurate identification of nucleotide-binding leucine-rich repeat genes and characterization of their inter- and intraspecific diversity. Moreover, we uncovered grain weight-associated SVs which specify traits by affecting the expression of their nearby genes. We characterized genetic variants associated with submergence tolerance, seed shattering and plant architecture and found independent selection for a common set of genes that drove adaptation and domestication in Asian and African rice. This super pan-genome facilitates pinpointing of lineage-specific haplotypes for trait-associated genes and provides insights into the evolutionary events that have shaped the genomic architecture of various rice species.