Sirenian genomes illuminate the evolution of fully aquatic species within the mammalian superorder afrotheria.

Sirenians of the superorder Afrotheria were the first mammals to transition from land to water and are the only herbivorous marine mammals. Here, we generated a chromosome-level dugong (Dugong dugon) genome. A comparison of our assembly with other afrotherian genomes reveals possible molecular adaptations to aquatic life by sirenians, including a shift in daily activity patterns (circadian clock) and tolerance to a high-iodine plant diet mediated through changes in the iodide transporter NIS (SLC5A5) and its co-transporters. Functional in vitro assays confirm that sirenian amino acid substitutions alter the properties of the circadian clock protein PER2 and NIS. Sirenians show evidence of convergent regression of integumentary system (skin and its appendages) genes with cetaceans. Our analysis also uncovers gene losses that may be maladaptive in a modern environment, including a candidate gene (KCNK18) for sirenian cold stress syndrome likely lost during their evolutionary shift in daily activity patterns. Genomes from nine Australian locations and the functionally extinct Okinawan population confirm and date a genetic break ~10.7 thousand years ago on the Australian east coast and provide evidence of an associated ecotype, and highlight the need for whole-genome resequencing data from dugong populations worldwide for conservation and genetic management.

Functional vertical connectivity of microbial communities in the ocean.

Sinking particles are a critical conduit for the transport of surface microbes to the ocean's interior. Vertical connectivity of phylogenetic composition has been shown; however, the functional vertical connectivity of microbial communities has not yet been explored in detail. We investigated protein and taxa profiles of both free-living and particle-attached microbial communities from the surface to 3000 m depth using a combined metaproteomic and 16S rRNA amplicon sequencing approach. A clear compositional and functional vertical connectivity of microbial communities was observed throughout the water column with Oceanospirillales, Alteromonadales, and Rhodobacterales as key taxa. The surface-derived particle-associated microbes increased the expression of proteins involved in basic metabolism, organic matter processing, and environmental stress response in deep waters. This study highlights the functional vertical connectivity between surface and deep-sea microbial communities via sinking particles and reveals that a considerable proportion of the deep-sea microbes might originate from surface waters and have a major impact on the biogeochemical cycles in the deep sea.

A new invertebrate NPY-like polypeptide, ZoaNPY, from the Zoanthus sociatus, as a novel ligand of human NPY Y2 receptor rescues vascular insufficiency via PLC/PKC and Src- FAK-dependent signaling pathways.

Our recent multi-omics studies have revealed rich sources of novel bioactive proteins and polypeptides from marine organisms including cnidarians. In the present study, we initially conducted a transcriptomic analysis to review the composition profile of polypeptides from Zoanthus sociatus. Then, a newly discovered NPY-like polypeptide-ZoaNPY was selected for further in silico structural, binding and virtually pharmacological studies. To evaluate the pro-angiogenic effects of ZoaNPY, we employed an in vitro HUVECs model and an in vivo zebrafish model. Our results indicate that ZoaNPY, at 1-100 pmol, enhances cell survival, migration and tube formation in the endothelial cells. Besides, treatment with ZoaNPY could restore a chemically-induced vascular insufficiency in zebrafish embryos. Western blot results demonstrated the application of ZoaNPY could increase the phosphorylation of proteins related to angiogenesis signaling including PKC, PLC, FAK, Src, Akt, mTOR, MEK, and ERK1/2. Furthermore, through molecular docking and surface plasmon resonance (SPR) verification, ZoaNPY was shown to directly and physically interact with NPY Y2 receptor. In view of this, all evidence showed that the pro-angiogenic effects of ZoaNPY involve the activation of NPY Y2 receptor, thereby activating the Akt/mTOR, PLC/PKC, ERK/MEK and Src- FAK-dependent signaling pathways. Furthermore, in an excision wound model, the treatment with ZoaNPY was shown to accelerate the wound healing process in mice. Our findings provide new insights into the discovery and development of novel pro-angiogenic drugs derived from NPY-like polypeptides in the future.

A chromosome-level genome of electric catfish (Malapterurus electricus) provided new insights into order Siluriformes evolution.

The electric catfish (Malapterurus electricus), belonging to the family Malapteruridae, order Siluriformes (Actinopterygii: Ostariophysi), is one of the six branches that has independently evolved electrical organs. We assembled a 796.75 Mb M. electricus genome and anchored 88.72% sequences into 28 chromosomes. Gene family analysis revealed 295 expanded gene families that were enriched on functions related to glutamate receptors. Convergent evolutionary analyses of electric organs among different lineage of electric fishes further revealed that the coding gene of rho guanine nucleotide exchange factor 4-like (arhgef4), which is associated with G-protein coupled receptor (GPCR) signaling pathway, underwent adaptive parallel evolution. Gene identification suggests visual degradation in catfishes, and an important role for taste in environmental adaptation. Our findings fill in the genomic data for a branch of electric fish and provide a relevant genetic basis for the adaptive evolution of Siluriformes. Supplementary Information: The online version contains supplementary material available at 10.1007/s42995-023-00197-8.

Spatially resolved single-cell atlas of ascidian endostyle provides insight into the origin of vertebrate pharyngeal organs.

The pharyngeal endoderm, an innovation of deuterostome ancestors, contributes to pharyngeal development by influencing the patterning and differentiation of pharyngeal structures in vertebrates; however, the evolutionary origin of the pharyngeal organs in vertebrates is largely unknown. The endostyle, a distinct pharyngeal organ exclusively present in basal chordates, represents a good model for understanding pharyngeal organ origins. Using Stereo-seq and single-cell RNA sequencing, we constructed aspatially resolved single-cell atlas for the endostyle of the ascidian Styela clava. We determined the cell composition of the hemolymphoid region, which illuminates a mixed ancestral structure for the blood and lymphoid system. In addition, we discovered a cluster of hair cell-like cells in zone 3, which has transcriptomic similarity with the hair cells of the vertebrate acoustico-lateralis system. These findings reshape our understanding of the pharynx of the basal chordate and provide insights into the evolutionary origin of multiplexed pharyngeal organs.

Single-cell transcriptome landscape of zebrafish liver reveals hepatocytes and immune cell interactions in understanding nonalcoholic fatty liver disease.

Nonalcoholic fatty liver disease (NAFLD) is becoming the most common chronic liver disease in the world. Immunity is the major contributing factor in NAFLD; however, the interaction of immune cells and hepatocytes in disease progression has not been fully elucidated. As a popular species for studying NAFLD, zebrafish, whose liver is a complex immune system mediated by immune cells and non-immune cells in maintaining immune tolerance and homeostasis. Understanding the cellular composition and immune environment of zebrafish liver is of great significance for its application in NAFLD. Here, we established a liver atlas that consists of 10 cell types using single-cell RNA sequencing (scRNA-seq). By examining the heterogeneity of hepatocytes and analyzing the expression of NAFLD-associated genes in the specific cluster, we provide a potential target cell model to study NAFLD. Additionally, our analysis identified two subtypes of distinct resident macrophages with inflammatory and non-inflammatory functions and characterized the successive stepwise development of T cell subclusters in the liver. Importantly, we uncovered the possible regulation of macrophages and T cells on target cells of fatty liver by analyzing the cellular interaction between hepatocytes and immune cells. Our data provide valuable information for an in-depth study of immune cells targeting hepatocytes to regulate the immune balance in NAFLD.

In situ community transcriptomics illuminates CO2-fixation potentials and supporting roles of phagotrophy and proton pump in plankton in a subtropical marginal sea.

Lineage-wise physiological activities of plankton communities in the ocean are important but challenging to characterize. Here, we conducted whole-assemblage metatranscriptomic profiling at continental shelf and slope sites in the South China Sea to investigate carbon fixation potential in different lineages. RuBisCO expression, the proxy of Calvin carbon fixation (CCF) potential, was mainly contributed by Bacillariophyta, Chlorophyta, Cyanobacteria, and Haptophyta, which was differentially affected by environmental factors among lineages. CCF potential exhibited positive or negative correlations with phagotrophy gene expression, suggesting phagotrophy possibly enhances or complements CCF. Our data also reveal significant non-Calvin carbon fixation (NCF) potential, as indicated by the active expression of genes in all five currently recognized NCF pathways, mainly contributed by Flavobacteriales, Alteromonadales, and Oceanospirillales. Furthermore, in Flavobacteriales, Alteromonadales, Pelagibacterales, and Rhodobacterales, NCF potential was positively correlated with proton-pump rhodopsin (PPR) expression, suggesting that NCF might be energetically supported by PPR. The novel insights into the lineage-differential potential of carbon fixation, widespread mixotrophy, and PPR as an energy source for NCF lay a methodological and informational foundation for further research to understand carbon fixation and the trophic landscape in the ocean.IMPORTANCEMarine plankton plays an important role in global carbon cycling and climate regulation. Phytoplankton and cyanobacteria fix CO2 to produce organic compounds using solar energy and mainly by the Calvin cycle, whereas autotrophic bacteria and archaea may fix CO2 by non-Calvin cycle carbon fixation pathways. How active individual lineages are in carbon fixation and mixotrophy, and what energy source bacteria may employ in non-Calvin carbon fixation, in a natural plankton assemblage are poorly understood and underexplored. Using metatranscriptomics, we studied carbon fixation in marine plankton with lineage resolution in tropical marginal shelf and slope areas. Based on the sequencing results, we characterized the carbon fixation potential of different lineages and assessed Calvin- and non-Calvin- carbon fixation activities and energy sources. Data revealed a high number of unigenes (4.4 million), lineage-dependent differential potentials of Calvin carbon fixation and responses to environmental conditions, major contributors of non-Calvin carbon fixation, and their potential energy source.

Genomics of post-bottleneck recovery in the northern elephant seal.

Populations and species are threatened by human pressure, but their fate is variable. Some depleted populations, such as that of the northern elephant seal (Mirounga angustirostris), recover rapidly even when the surviving population was small. The northern elephant seal was hunted extensively and taken by collectors between the early 1800s and 1892, suffering an extreme population bottleneck as a consequence. Recovery was rapid and now there are over 200,000 individuals. We sequenced 260 modern and 8 historical northern elephant seal nuclear genomes to assess the impact of the population bottleneck on individual northern elephant seals and to better understand their recovery. Here we show that inbreeding, an increase in the frequency of alleles compromised by lost function, and allele frequency distortion, reduced the fitness of breeding males and females, as well as the performance of adult females on foraging migrations. We provide a detailed investigation of the impact of a severe bottleneck on fitness at the genomic level and report on the role of specific gene systems.

Combined ATAC-seq, RNA-seq, and GWAS analysis reveals glycogen metabolism regulatory network in Jinjiang oyster ( Crassostrea ariakensis).

Glycogen serves as the principal energy reserve for metabolic processes in aquatic shellfish and substantially contributes to the flavor and quality of oysters. The Jinjiang oyster ( Crassostrea ariakensis) is an economically and ecologically important species in China. In the present study, RNA sequencing (RNA-seq) and assay for transposase-accessible chromatin using sequencing (ATAC-seq) were performed to investigate gene expression and chromatin accessibility variations in oysters with different glycogen contents. Analysis identified 9 483 differentially expressed genes (DEGs) and 7 215 genes with significantly differential chromatin accessibility (DCAGs) were obtained, with an overlap of 2 600 genes between them. Notably, a significant proportion of these genes were enriched in pathways related to glycogen metabolism, including "Glycogen metabolic process" and "Starch and sucrose metabolism". In addition, genome-wide association study (GWAS) identified 526 single nucleotide polymorphism (SNP) loci associated with glycogen content. These loci corresponded to 241 genes, 63 of which were categorized as both DEGs and DCAGs. This study enriches basic research data and provides insights into the molecular mechanisms underlying the regulation of glycogen metabolism in C. ariakensis.

Unravelling the genomic features, phylogeny and genetic basis of tooth ontogenesis in Characiformes through analysis of four genomes.

Characiformes is a diverse and evolutionarily significant order of freshwater fish encompassing over 2,300 species. Despite its diversity, our understanding of Characiformes' evolutionary relationships and adaptive mechanisms is limited due to insufficient genome sequences. In this study, we sequenced and assembled the genomes of four Characiformes species, three of which were chromosome-level assemblies. Our analyses revealed dynamic changes in gene family evolution, repeat sequences and variations in chromosomal collinearity within these genomes. With the assembled genomes, we were not only able to elucidate the evolutionary relationship of the four main orders in Otophysi but also indicated Characiformes as the paraphyletic group. Comparative genomic analysis with other available fish genomes shed light on the evolution of genes related to tooth development in Characiformes. Notably, variations in the copy number of secretory calcium-binding phosphoproteins (SCPP) genes were observed among different orders of Otophysi, indicating their potential contribution to the diversity of tooth types. Our study offers invaluable genome sequences and novel insights into Characiformes' evolution, paving the way for further genomic and evolutionary research in fish.

A single-cell atlas of West African lungfish respiratory system reveals evolutionary adaptations to terrestrialization.

The six species of lungfish possess both lungs and gills and are the closest extant relatives of tetrapods. Here, we report a single-cell transcriptome atlas of the West African lungfish (Protopterus annectens). This species manifests the most extreme form of terrestrialization, a life history strategy to survive dry periods that can last for years, characterized by dormancy and reversible adaptive changes of the gills and lungs. Our atlas highlights the cell type diversity of the West African lungfish, including gene expression consistent with phenotype changes of terrestrialization. Comparison with terrestrial tetrapods and ray-finned fishes reveals broad homology between the swim bladder and lung cell types as well as shared and idiosyncratic changes of the external gills of the West African lungfish and the internal gills of Atlantic salmon. The single-cell atlas presented here provides a valuable resource for further exploration of the respiratory system evolution in vertebrates and the diversity of lungfish terrestrialization.

Chromosome-level genome of butterflyfish unveils genomic features of unique colour patterns and morphological traits.

Chaetodontidae, known as butterflyfishes, are typical fish in coral ecosystems, exhibiting remarkable interspecific differences including body colour patterns and feeding ecology. In this study, we report genomes of three butterflyfish species (Chelmon rostratus, Chaetodon trifasciatus and Chaetodon auriga) and a closely related species from the Pomacanthidae family, Centropyge bicolour, with an average genome size of 65,611 Mb. Chelmon rostratus, comprising 24 chromosomes assembled to the chromosome level, could be served as a reference genome for butterflyfish. By conducting a collinearity analysis between butterflyfishes and several fishes, we elucidated the specific and conserved genomic features of butterflyfish, with particular emphasis on novel genes arising from tandem duplications and their potential functions. In addition to the two melanocyte-specific tyr genes commonly found in fish, we found the gene tyrp3, a new tyrosinase-related proteins gene in the reef fish, including butterflyfish and clownfish, implicating their involvement in the pigmentation diversity of fish. Additionally, we observed a tandem duplication expansion of three copies of nell1 gene in C. rostratus genome, which likely contribute to its unique jaw development and distinctive morphology of its sharp mouth. These results provided valuable genomic resources for further investigations into the genetic diversity and evolutionary adaptations of reef fish.

The complete and fully-phased diploid genome of a male Han Chinese.

Since the release of the complete human genome, the priority of human genomic study has now been shifting towards closing gaps in ethnic diversity. Here, we present a fully phased and well-annotated diploid human genome from a Han Chinese male individual (CN1), in which the assemblies of both haploids achieve the telomere-to-telomere (T2T) level. Comparison of this diploid genome with the CHM13 haploid T2T genome revealed significant variations in the centromere. Outside the centromere, we discovered 11,413 structural variations, including numerous novel ones. We also detected thousands of CN1 alleles that have accumulated high substitution rates and a few that have been under positive selection in the East Asian population. Further, we found that CN1 outperforms CHM13 as a reference genome in mapping and variant calling for the East Asian population owing to the distinct structural variants of the two references. Comparison of SNP calling for a large cohort of 8869 Chinese genomes using CN1 and CHM13 as reference respectively showed that the reference bias profoundly impacts rare SNP calling, with nearly 2 million rare SNPs miss-called with different reference genomes. Finally, applying the CN1 as a reference, we discovered 5.80 Mb and 4.21 Mb putative introgression sequences from Neanderthal and Denisovan, respectively, including many East Asian specific ones undetected using CHM13 as the reference. Our analyses reveal the advances of using CN1 as a reference for population genomic studies and paleo-genomic studies. This complete genome will serve as an alternative reference for future genomic studies on the East Asian population.

VT3D: a visualization toolbox for 3D transcriptomic data.

Data visualization empowers researchers to communicate their results that support scientific reasoning in an intuitive way. Three-dimension (3D) spatially resolved transcriptomic atlases constructed from multi-view and high-dimensional data have rapidly emerged as a powerful tool to unravel spatial gene expression patterns and cell type distribution in biological samples, revolutionizing the understanding of gene regulatory interactions and cell niches. However, limited accessible tools for data visualization impede the potential impact and application of this technology. Here we introduce VT3D, a visualization toolbox that allows users to explore 3D transcriptomic data, enabling gene expression projection to any 2D plane of interest, 2D virtual slice creation and visualization, and interactive 3D data browsing with surface model plots. In addition, it can either work on personal devices in standalone mode or be hosted as a web-based server. We apply VT3D to multiple datasets produced by the most popular techniques, including both sequencing-based approaches (Stereo-seq, spatial transcriptomics, and Slide-seq) and imaging-based approaches (MERFISH and STARMap), and successfully build a 3D atlas database that allows interactive data browsing. We demonstrate that VT3D bridges the gap between researchers and spatially resolved transcriptomics, thus accelerating related studies such as embryogenesis and organogenesis processes. The source code of VT3D is available at https://github.com/BGI-Qingdao/VT3D, and the modeled atlas database is available at http://www.bgiocean.com/vt3d_example.

The enormous repetitive Antarctic krill genome reveals environmental adaptations and population insights.

Antarctic krill (Euphausia superba) is Earth's most abundant wild animal, and its enormous biomass is vital to the Southern Ocean ecosystem. Here, we report a 48.01-Gb chromosome-level Antarctic krill genome, whose large genome size appears to have resulted from inter-genic transposable element expansions. Our assembly reveals the molecular architecture of the Antarctic krill circadian clock and uncovers expanded gene families associated with molting and energy metabolism, providing insights into adaptations to the cold and highly seasonal Antarctic environment. Population-level genome re-sequencing from four geographical sites around the Antarctic continent reveals no clear population structure but highlights natural selection associated with environmental variables. An apparent drastic reduction in krill population size 10 mya and a subsequent rebound 100 thousand years ago coincides with climate change events. Our findings uncover the genomic basis of Antarctic krill adaptations to the Southern Ocean and provide valuable resources for future Antarctic research.

Hologenome analysis reveals independent evolution to chemosymbiosis by deep-sea bivalves.

BACKGROUND: Bivalves have independently evolved a variety of symbiotic relationships with chemosynthetic bacteria. These relationships range from endo- to extracellular interactions, making them ideal for studies on symbiosis-related evolution. It is still unclear whether there are universal patterns to symbiosis across bivalves. Here, we investigate the hologenome of an extracellular symbiotic thyasirid clam that represents the early stages of symbiosis evolution. RESULTS: We present a hologenome of Conchocele bisecta (Bivalvia: Thyasiridae) collected from deep-sea hydrothermal vents with extracellular symbionts, along with related ultrastructural evidence and expression data. Based on ultrastructural and sequencing evidence, only one dominant Thioglobaceae bacteria was densely aggregated in the large bacterial chambers of C. bisecta, and the bacterial genome shows nutritional complementarity and immune interactions with the host. Overall, gene family expansions may contribute to the symbiosis-related phenotypic variations in different bivalves. For instance, convergent expansions of gaseous substrate transport families in the endosymbiotic bivalves are absent in C. bisecta. Compared to endosymbiotic relatives, the thyasirid genome exhibits large-scale expansion in phagocytosis, which may facilitate symbiont digestion and account for extracellular symbiotic phenotypes. We also reveal that distinct immune system evolution, including expansion in lipopolysaccharide scavenging and contraction of IAP (inhibitor of apoptosis protein), may contribute to the different manners of bacterial virulence resistance in C. bisecta. CONCLUSIONS: Thus, bivalves employ different pathways to adapt to the long-term co-existence with their bacterial symbionts, further highlighting the contribution of stochastic evolution to the independent gain of a symbiotic lifestyle in the lineage.

Inbreeding depression explains killer whale population dynamics.

Understanding the factors that cause endangered populations to either grow or decline is crucial for preserving biodiversity. Conservation efforts often address extrinsic threats, such as environmental degradation and overexploitation, that can limit the recovery of endangered populations. Genetic factors such as inbreeding depression can also affect population dynamics but these effects are rarely measured in the wild and thus often neglected in conservation efforts. Here we show that inbreeding depression strongly influences the population dynamics of an endangered killer whale population, despite genomic signatures of purging of deleterious alleles via natural selection. We find that the 'Southern Residents', which are currently endangered despite nearly 50 years of conservation efforts, exhibit strong inbreeding depression for survival. Our population models suggest that this inbreeding depression limits population growth and predict further decline if the population remains genetically isolated and typical environmental conditions continue. The Southern Residents also had more inferred homozygous deleterious alleles than three other, growing, populations, further suggesting that inbreeding depression affects population fitness. These results demonstrate that inbreeding depression can substantially limit the recovery of endangered populations. Conservation actions focused only on extrinsic threats may therefore fail to account for key intrinsic genetic factors that also limit population growth.

A genome and single-nucleus cerebral cortex transcriptome atlas of the short-finned pilot whale Globicephala macrorhynchus.

Cetaceans (dolphins, whales, and porpoises) have large and anatomically sophisticated brains. To expand our understanding of the cellular makeup of cetacean brains and the similarities and divergence between the brains of cetaceans and terrestrial mammals, we report a short-finned pilot whale (Globicephala macrorhynchus) single-nucleus transcriptome atlas. To achieve this goal, we assembled a chromosome-scale reference genome spanning 2.25 Gb on 22 chromosomes and profiled the gene expression of five major anatomical cortical regions of the short-finned pilot whale by single-nucleus RNA-sequencing (snRNA-seq). We identified six major cell lineages in the cerebral cortex (excitatory neurons, inhibitory neurons, oligodendrocytes, oligodendrocyte precursor cells, astrocytes, and endothelial cells), eight molecularly distinct subclusters of excitatory neurons, and four subclusters of inhibitory neurons. Finally, a comparison of snRNA-seq data from the short-finned pilot whale, human, and rhesus macaque revealed a broadly conserved cellular makeup of brain cell types. Our study provides genomic resources and molecular insights into cetacean brain evolution.

Single-cell-resolution transcriptome map revealed novel genes involved in testicular germ cell progression and somatic cells specification in Chinese tongue sole with sex reversal.

Female-to-male sex reversals (pseudomales) are common in lower vertebrates and have been found in natural populations, which is a concern under rapid changes in environmental conditions. Pseudomales can exhibit altered spermatogenesis. However, the regulatory mechanisms underlying pseudomale spermatogenesis remain unclear. Here, we characterized spermatogenesis in Chinese tongue sole (Cynoglossus semilaevis), a species with genetic and environmental sex determination, based on a high-resolution single-cell RNA-seq atlas of cells derived from the testes of genotypic males and pseudomales. We identified five germ cell types and six somatic cell types and obtained a single-cell atlas of dynamic changes in gene expression during spermatogenesis in Chinese tongue sole, including alterations in pseudomales. We detected decreased levels of Ca2+ signaling pathway-related genes in spermatogonia, insufficient meiotic initiation in spermatocytes, and a malfunction of somatic niche cells in pseudomales. However, a cluster of CaSR genes and MAPK signaling factors were upregulated in undifferentiated spermatogonia of pseudomales. Additionally, we revealed that Z chromosome-specific genes, such as piwil2, dhx37, and ehmt1, were important for spermatogenesis. These results improve our understanding of reproduction after female-to-male sex-reversal and provide new insights into the adaptability of reproductive strategies in lower vertebrates.

Genomic insights into the genetic signatures of selection and seed trait loci in cultivated peanut.

INTRODUCTION: Cultivated peanut (Arachis hypogaea L.) is an important oil crop for human nutrition and is cultivated in >100 countries. However, the present knowledge of its genomic diversity, evolution, and loci related to the seed traits is limited. OBJECTIVES: Our study intended to (1) uncover the population structure and the demographic history of peanuts, (2) identify signatures of selection that occurred during peanut improvement breeding, and (3) detect and verify the functions of candidate genes associated with seed traits. METHODS: We explored the population relationship and the evolution of peanuts using a largescale single nucleotide polymorphism dataset generated from the genome-wide resequencing of 203 cultivated peanuts. Genetic diversity and genomic scan analyses were applied to identify selective loci for genomic-selection breeding. Genome-wide association studies, transgenic experiments, and RNA-seq were employed to identify the candidate genes associated with seed traits. RESULTS: Our study revealed that the 203 resequenced accessions were divided into four genetic groups, consistent with their botanical classification. Moreover, the var. peruviana and var. fastigiata subpopulations have diverged to a greater extent than the others, and var. peruviana may be the earliest variant in the evolution from tetraploid ancestors. A recent dramatic expansion in the effective population size of the cultivated peanuts ca. 300-500 years ago was also noted. Selective sweeps underlying quantitative trait loci and genes of seed size, plant architecture, and disease resistance coincide with the major goals of improved peanut breeding compared with the landrace and cultivar populations. Genome-wide association testing with functional analysis led to the identification of two genes involved in seed weight and seed length regulation. CONCLUSION: Our study provides valuable information for understanding the genomic diversity and the evolution of peanuts and serves as a genomic basis for improving peanut cultivars.