Early Austronesians: Into and Out Of Taiwan

A Taiwan origin for the expansion of the Austronesian languages and their speakers is well supported by linguistic and archaeological evidence. However, human genetic evidence is more controversial. Until now, there had been no ancient skeletal evidence of a potential Austronesian-speaking ancestor prior to the Taiwan Neolithic ∼6,000 years ago, and genetic studies have largely ignored the role of genetic diversity within Taiwan as well as the origins of Formosans. We address these issues via analysis of a complete mitochondrial DNA genome sequence of an ∼8,000-year-old skeleton from Liang Island (located between China and Taiwan) and 550 mtDNA genome sequences from 8 aboriginal (highland) Formosan and 4 other Taiwanese groups. We show that the Liangdao Man mtDNA sequence is closest to Formosans, provides a link to southern China, and has the most ancestral haplogroup E sequence found among extant Austronesian speakers. Bayesian phylogenetic analysis allows us to reconstruct a history of early Austronesians arriving in Taiwan in the north ∼6,000 years ago, spreading rapidly to the south, and leaving Taiwan ∼4,000 years ago to spread throughout Island Southeast Asia, Madagascar, and Oceania.     

Genetic Origins and Sex-Biased Admixture of the Huis

The Hui people are unique among Chinese ethnic minorities in that they speak the same language as Han Chinese (HAN) but practice Islam. However, as the second-largest minority group in China numbering well over 10 million, the Huis are under-represented in both global and regional genomic studies. Here, we present the first whole-genome sequencing effort of 234 Hui individuals (NXH) aged over 60 who have been living in Ningxia, where the Huis are mostly concentrated. NXH are genetically more similar to East Asian than to any other global populations. In particular, the genetic differentiation between NXH and HAN (F_ST = 0.0015) is only slightly larger than that between northern and southern HAN (F_ST = 0.0010), largely attributed to the western ancestry in NXH (∼10%). Highly differentiated functional variants between NXH and HAN were identified in genes associated with skin pigmentation (e.g., SLC24A5), facial morphology (e.g., EDAR), and lipid metabolism (e.g., ABCG8). The Huis are also distinct from other Muslim groups such as the Uyghurs (F_ST = 0.0187), especially, NXH derived much less western ancestry (∼10%) compared with the Uyghurs (∼50%). Modeling admixture history indicated that NXH experienced an episode of two-wave admixture. An ancient admixture occurred ∼1,025 years ago, reflecting the intensive west–east contacts during the late Tang Dynasty, and the Five Dynasties and Ten Kingdoms period. A recent admixture occurred ∼500 years ago, corresponding to the Ming Dynasty. Notably, we identified considerable sex-biased admixture, that is, excess of western males and eastern females contributing to the NXH gene pool. The origins and the genomic diversity of the Hui people imply the complex history of contacts between western and eastern Eurasians.     

Oxygen Concentration Inside a Functioning Photosynthetic Cell

The excess oxygen concentration in the photosynthetic membranes of functioning oxygenic photosynthetic cells was estimated using classical diffusion theory combined with experimental data on oxygen production rates of cyanobacterial cells. The excess oxygen concentration within the plesiomorphic cyanobacterium Gloeobactor violaceus is only 0.025 μM, or four orders of magnitude lower than the oxygen concentration in air-saturated water. Such a low concentration suggests that the first oxygenic photosynthetic bacteria in solitary form could have evolved ∼2.8 billion years ago without special mechanisms to protect them against reactive oxygen species. These mechanisms instead could have been developed during the following ∼500 million years while the oxygen level in the Earth’s atmosphere was slowly rising. Excess oxygen concentrations within individual cells of the apomorphic cyanobacteria Synechocystis and Synechococcus are 0.064 and 0.25 μM, respectively. These numbers suggest that intramembrane and intracellular proteins in isolated oxygenic photosynthetic cells are not subjected to excessively high oxygen levels. The situation is different for closely packed colonies of photosynthetic cells. Calculations show that the excess concentration within colonies that are ∼40 μm or larger in diameter can be comparable to the oxygen concentration in air-saturated water, suggesting that species forming colonies require protection against reactive oxygen species even in the absence of oxygen in the surrounding atmosphere.     

Tyrant Dinosaur Evolution Tracks the Rise and Fall of Late Cretaceous Oceans

The Late Cretaceous (∼95–66 million years ago) western North American landmass of Laramidia displayed heightened non-marine vertebrate diversity and intracontinental regionalism relative to other latest Cretaceous Laurasian ecosystems. Processes generating these patterns during this interval remain poorly understood despite their presumed role in the diversification of many clades. Tyrannosauridae, a clade of large-bodied theropod dinosaurs restricted to the Late Cretaceous of Laramidia and Asia, represents an ideal group for investigating Laramidian patterns of evolution. We use new tyrannosaurid discoveries from Utah—including a new taxon which represents the geologically oldest member of the clade—to investigate the evolution and biogeography of Tyrannosauridae. These data suggest a Laramidian origin for Tyrannosauridae, and implicate sea-level related controls in the isolation, diversification, and dispersal of this and many other Late Cretaceous vertebrate clades.     

A New Malaria Agent in African Hominids

Plasmodium falciparum is the major human malaria agent responsible for 200 to 300 million infections and one to three million deaths annually, mainly among African infants. The origin and evolution of this pathogen within the human lineage is still unresolved. A single species, P. reichenowi, which infects chimpanzees, is known to be a close sister lineage of P. falciparum. Here we report the discovery of a new Plasmodium species infecting Hominids. This new species has been isolated in two chimpanzees (Pan troglodytes) kept as pets by villagers in Gabon (Africa). Analysis of its complete mitochondrial genome (5529 nucleotides including Cyt b, Cox I and Cox III genes) reveals an older divergence of this lineage from the clade that includes P. falciparum and P. reichenowi (∼21±9 Myrs ago using Bayesian methods and considering that the divergence between P. falciparum and P. reichenowi occurred 4 to 7 million years ago as generally considered in the literature). This time frame would be congruent with the radiation of hominoids, suggesting that this Plasmodium lineage might have been present in early hominoids and that they may both have experienced a simultaneous diversification. Investigation of the nuclear genome of this new species will further the understanding of the genetic adaptations of P. falciparum to humans. The risk of transfer and emergence of this new species in humans must be now seriously considered given that it was found in two chimpanzees living in contact with humans and its close relatedness to the most virulent agent of malaria.     

A RAD-Tag Genetic Map for the Platyfish (Xiphophorus maculatus) Reveals Mechanisms of Karyotype Evolution Among Teleost Fish

Mammalian genomes can vary substantially in haploid chromosome number even within a small taxon (e.g., 3–40 among deer alone); in contrast, teleost fish genomes are stable (24–25 in 58% of teleosts), but we do not yet understand the mechanisms that account for differences in karyotype stability. Among perciform teleosts, platyfish (Xiphophorus maculatus) and medaka (Oryzias latipes) both have 24 chromosome pairs, but threespine stickleback (Gasterosteus aculeatus) and green pufferfish (Tetraodon nigroviridis) have just 21 pairs. To understand the evolution of teleost genomes, we made a platyfish meiotic map containing 16,114 mapped markers scored on 267 backcross fish. We tiled genomic contigs along the map to create chromosome-length genome assemblies. Genome-wide comparisons of conserved synteny showed that platyfish and medaka karyotypes remained remarkably similar with few interchromosomal translocations but with numerous intrachromosomal rearrangements (transpositions and inversions) since their lineages diverged ∼120 million years ago. Comparative genomics with platyfish shows how reduced chromosome numbers in stickleback and green pufferfish arose by fusion of pairs of ancestral chromosomes after their lineages diverged from platyfish ∼195 million years ago. Zebrafish and human genomes provide outgroups to root observed changes. These studies identify likely genome assembly errors, characterize chromosome fusion events, distinguish lineage-independent chromosome fusions, show that the teleost genome duplication does not appear to have accelerated the rate of translocations, and reveal the stability of syntenies and gene orders in teleost chromosomes over hundreds of millions of years.     

Implications of Diet for the Extinction of Saber-Toothed Cats and American Lions

The saber-toothed cat, Smilodon fatalis, and American lion, Panthera atrox, were among the largest terrestrial carnivores that lived during the Pleistocene, going extinct along with other megafauna ∼12,000 years ago. Previous work suggests that times were difficult at La Brea (California) during the late Pleistocene, as nearly all carnivores have greater incidences of tooth breakage (used to infer greater carcass utilization) compared to today. As Dental Microwear Texture Analysis (DMTA) can differentiate between levels of bone consumption in extant carnivores, we use DMTA to clarify the dietary niches of extinct carnivorans from La Brea. Specifically, we test the hypothesis that times were tough at La Brea with carnivorous taxa utilizing more of the carcasses. Our results show no evidence of bone crushing by P. atrox, with DMTA attributes most similar to the extant cheetah, Acinonyx jubatus, which actively avoids bone. In contrast, S. fatalis has DMTA attributes most similar to the African lion Panthera leo, implying that S. fatalis did not avoid bone to the extent previously suggested by SEM microwear data. DMTA characters most indicative of bone consumption (i.e., complexity and textural fill volume) suggest that carcass utilization by the extinct carnivorans was not necessarily more complete during the Pleistocene at La Brea; thus, times may not have been “tougher” than the present. Additionally, minor to no significant differences in DMTA attributes from older (∼30–35 Ka) to younger (∼11.5 Ka) deposits offer little evidence that declining prey resources were a primary cause of extinction for these large cats.     

First Dinosaurs from Saudi Arabia

Dinosaur remains from the Arabian subcontinent are exceedingly rare, and those that have been documented manifest indeterminate affinities. Consequently the discovery of a small, but diagnostic, accumulation of elements from Campanian-Maastrichtian (∼75 Ma) deposits in northwestern Saudi Arabia is significant because it constitutes the first taxonomically identifiable dinosaur material described from the Arabian Peninsula. The fossils include a series of possible lithostrotian titanosaur caudal vertebrae, and some isolated theropod marginal teeth that share unique character states and metric parameters (analyzed using multivariate statistical methods) with derived abelisaurids – this is the first justifiable example of a non-avian carnivorous dinosaur clade from Arabia. The recognition of titanosaurians and abelisaurids from Saudi Arabia extends the palaeogeographical range of these groups along the entire northern Gondwanan margin during the latest Cretaceous. Moreover, given the extreme paucity of coeval occurrences elsewhere, the Saudi Arabian fossils provide a tantalizing glimpse into dinosaurian assemblage diversity within the region.     

Effects of Global Warming on Ancient Mammalian Communities and Their Environments

Background

Current global warming affects the composition and dynamics of mammalian communities and can increase extinction risk; however, long-term effects of warming on mammals are less understood. Dietary reconstructions inferred from stable isotopes of fossil herbivorous mammalian tooth enamel document environmental and climatic changes in ancient ecosystems, including C₃/C₄ transitions and relative seasonality.

Methodology/Principal Findings

Here, we use stable carbon and oxygen isotopes preserved in fossil teeth to document the magnitude of mammalian dietary shifts and ancient floral change during geologically documented glacial and interglacial periods during the Pliocene (∼1.9 million years ago) and Pleistocene (∼1.3 million years ago) in Florida. Stable isotope data demonstrate increased aridity, increased C₄ grass consumption, inter-faunal dietary partitioning, increased isotopic niche breadth of mixed feeders, niche partitioning of phylogenetically similar taxa, and differences in relative seasonality with warming.

Conclusion/Significance

Our data show that global warming resulted in dramatic vegetation and dietary changes even at lower latitudes (∼28°N). Our results also question the use of models that predict the long term decline and extinction of species based on the assumption that niches are conserved over time. These findings have immediate relevance to clarifying possible biotic responses to current global warming in modern ecosystems.     

Northern Hemisphere Glaciation during the Globally Warm Early Late Pliocene

The early Late Pliocene (3.6 to ∼3.0 million years ago) is the last extended interval in Earth''s history when atmospheric CO₂ concentrations were comparable to today''s and global climate was warmer. Yet a severe global glaciation during marine isotope stage (MIS) M2 interrupted this phase of global warmth ∼3.30 million years ago, and is seen as a premature attempt of the climate system to establish an ice-age world. Here we propose a conceptual model for the glaciation and deglaciation of MIS M2 based on geochemical and palynological records from five marine sediment cores along a Caribbean to eastern North Atlantic transect. Our records show that increased Pacific-to-Atlantic flow via the Central American Seaway weakened the North Atlantic Current and attendant northward heat transport prior to MIS M2. The consequent cooling of the northern high latitude oceans permitted expansion of the continental ice sheets during MIS M2, despite near-modern atmospheric CO₂ concentrations. Sea level drop during this glaciation halted the inflow of Pacific water to the Atlantic via the Central American Seaway, allowing the build-up of a Caribbean Warm Pool. Once this warm pool was large enough, the Gulf Stream–North Atlantic Current system was reinvigorated, leading to significant northward heat transport that terminated the glaciation. Before and after MIS M2, heat transport via the North Atlantic Current was crucial in maintaining warm climates comparable to those predicted for the end of this century.     

Earliest Archaeological Evidence of Persistent Hominin Carnivory

The emergence of lithic technology by ∼2.6 million years ago (Ma) is often interpreted as a correlate of increasingly recurrent hominin acquisition and consumption of animal remains. Associated faunal evidence, however, is poorly preserved prior to ∼1.8 Ma, limiting our understanding of early archaeological (Oldowan) hominin carnivory. Here, we detail three large well-preserved zooarchaeological assemblages from Kanjera South, Kenya. The assemblages date to ∼2.0 Ma, pre-dating all previously published archaeofaunas of appreciable size. At Kanjera, there is clear evidence that Oldowan hominins acquired and processed numerous, relatively complete, small ungulate carcasses. Moreover, they had at least occasional access to the fleshed remains of larger, wildebeest-sized animals. The overall record of hominin activities is consistent through the stratified sequence – spanning hundreds to thousands of years – and provides the earliest archaeological evidence of sustained hominin involvement with fleshed animal remains (i.e., persistent carnivory), a foraging adaptation central to many models of hominin evolution.     

The Oxytricha trifallax Macronuclear Genome: A Complex Eukaryotic Genome with 16,000 Tiny Chromosomes

The macronuclear genome of the ciliate Oxytricha trifallax displays an extreme and unique eukaryotic genome architecture with extensive genomic variation. During sexual genome development, the expressed, somatic macronuclear genome is whittled down to the genic portion of a small fraction (∼5%) of its precursor “silent” germline micronuclear genome by a process of “unscrambling” and fragmentation. The tiny macronuclear “nanochromosomes” typically encode single, protein-coding genes (a small portion, 10%, encode 2–8 genes), have minimal noncoding regions, and are differentially amplified to an average of ∼2,000 copies. We report the high-quality genome assembly of ∼16,000 complete nanochromosomes (∼50 Mb haploid genome size) that vary from 469 bp to 66 kb long (mean ∼3.2 kb) and encode ∼18,500 genes. Alternative DNA fragmentation processes ∼10% of the nanochromosomes into multiple isoforms that usually encode complete genes. Nucleotide diversity in the macronucleus is very high (SNP heterozygosity is ∼4.0%), suggesting that Oxytricha trifallax may have one of the largest known effective population sizes of eukaryotes. Comparison to other ciliates with nonscrambled genomes and long macronuclear chromosomes (on the order of 100 kb) suggests several candidate proteins that could be involved in genome rearrangement, including domesticated MULE and IS1595-like DDE transposases. The assembly of the highly fragmented Oxytricha macronuclear genome is the first completed genome with such an unusual architecture. This genome sequence provides tantalizing glimpses into novel molecular biology and evolution. For example, Oxytricha maintains tens of millions of telomeres per cell and has also evolved an intriguing expansion of telomere end-binding proteins. In conjunction with the micronuclear genome in progress, the O. trifallax macronuclear genome will provide an invaluable resource for investigating programmed genome rearrangements, complementing studies of rearrangements arising during evolution and disease.     

Autosomal Resequence Data Reveal Late Stone Age Signals of Population Expansion in Sub-Saharan African Foraging and Farming Populations

Background

A major unanswered question in the evolution of Homo sapiens is when anatomically modern human populations began to expand: was demographic growth associated with the invention of particular technologies or behavioral innovations by hunter-gatherers in the Late Pleistocene, or with the acquisition of farming in the Neolithic?

Methodology/Principal Findings

We investigate the timing of human population expansion by performing a multilocus analysis of≥20 unlinked autosomal noncoding regions, each consisting of ∼6 kilobases, resequenced in ∼184 individuals from 7 human populations. We test the hypothesis that the autosomal polymorphism data fit a simple two-phase growth model, and when the hypothesis is not rejected, we fit parameters of this model to our data using approximate Bayesian computation.

Conclusions/Significance

The data from the three surveyed non-African populations (French Basque, Chinese Han, and Melanesians) are inconsistent with the simple growth model, presumably because they reflect more complex demographic histories. In contrast, data from all four sub-Saharan African populations fit the two-phase growth model, and a range of onset times and growth rates is inferred for each population. Interestingly, both hunter-gatherers (San and Biaka) and food-producers (Mandenka and Yorubans) best fit models with population growth beginning in the Late Pleistocene. Moreover, our hunter-gatherer populations show a tendency towards slightly older and stronger growth (∼41 thousand years ago, ∼13-fold) than our food-producing populations (∼31 thousand years ago, ∼7-fold). These dates are concurrent with the appearance of the Late Stone Age in Africa, supporting the hypothesis that population growth played a significant role in the evolution of Late Pleistocene human cultures.     

Early development of the facial region in a non-avian theropod dinosaur

An isolated maxilla of the theropod dinosaur Allosaurus from the Late Jurassic (the Kimmeridgian, 153 million years ago) of Portugal is the first cranial remain of a non-coelurosaurian theropod hatchling reported so far, and sheds new light on the early cranial development of non-avian theropods. Allosaurus hatchlings seem to have been one-seventh or less of the adult length and are thus comparable in relative size to hatchlings of large extant crocodile species, but are unlike the relatively larger hatchlings in coelurosaurs. The snout experienced considerable positive allometry and an increase in tooth count during early development. The element is especially noteworthy for the abundant and well-developed features associated with the paranasal pneumatic system. Pneumatic structures present include all those found in adult allosaurids and most are even more developed than in adult skulls. Together with evidence on the ontogeny of the tympanic pneumatic system in allosaurids, these findings demonstrate that cranial pneumaticity developed early in theropod ontogeny. The strong development of pneumatic features in early ontogenetic stages of non-avian theropods supports the hypothesis that pneumatization of cranial bones was opportunistic and indicates that heterochrony played an important role in the evolution of craniofacial pneumaticity in this group.     

New dinosaur tracksites from the Middle Jurassic of Madagascar: ichnotaxonomical,behavioural and palaeoenvironmental implications

Abstract: New dinosaur tracksites are described from the Bajocian–Bathonian Bemaraha Formation of western Madagascar. Two track‐bearing surfaces can be followed over a distance of at least 4 km, suggesting the existence of a hitherto unrecognized megatracksite. The track assemblage is theropod dominated, but sauropod tracks also occur at one site. Qualitative and quantitative analysis of the abundant theropod track material suggests that most, if not all, theropod footprints are attributable to a single trackmaker and are referred to Kayentapus isp. Although this ichnogenus, originally described from the Lower Jurassic of North America, has never been recorded from Gondwana nor from the Middle Jurassic, track morphology strongly suggests this attribution. Palaeogeographical, sedimentological and ichnological data suggest that the dinosaur tracks formed in an intertidal to supratidal setting where the coastline influenced the preferred walking direction of the animals.     

Bifurcation and Enhancement of Autonomous-Nonautonomous Retrotransposon Partnership through LTR Swapping in Soybean

Long terminal repeat (LTR) retrotransposons, the most abundant genomic components in flowering plants, are classifiable into autonomous and nonautonomous elements based on their structural completeness and transposition capacity. It has been proposed that selection is the major force for maintaining sequence (e.g., LTR) conservation between nonautonomous elements and their autonomous counterparts. Here, we report the structural, evolutionary, and expression characterization of a giant retrovirus-like soybean (Glycine max) LTR retrotransposon family, SNARE. This family contains two autonomous subfamilies, SARE^A and SARE^B, that appear to have evolved independently since the soybean genome tetraploidization event ∼13 million years ago, and a nonautonomous subfamily, SNRE, that originated from SARE^A. Unexpectedly, a subset of the SNRE elements, which amplified from a single founding SNRE element within the last ∼3 million years, have been dramatically homogenized with either SARE^A or SARE^B primarily in the LTR regions and bifurcated into distinct subgroups corresponding to the two autonomous subfamilies. We uncovered evidence of region-specific swapping of nonautonomous elements with autonomous elements that primarily generated various nonautonomous recombinants with LTR sequences from autonomous elements of different evolutionary lineages, thus revealing a molecular mechanism for the enhancement of preexisting partnership and the establishment of new partnership between autonomous and nonautonomous elements.     

Demographic factors shaped diversity in the two gene pools of wild common bean Phaseolus vulgaris L.

Wild common bean (Phaseolus vulgaris L.) is distributed throughout the Americas from Mexico to northern Argentina. Within this range, the species is divided into two gene pools (Andean and Middle American) along a latitudinal gradient. The diversity of 24 wild common bean genotypes from throughout the geographic range of the species was described by using sequence data from 13 loci. An isolation–migration model was evaluated using a coalescent analysis to estimate multiple demographic parameters. Using a Bayesian approach, Andean and Middle American subpopulations with high percentage of parentages were observed. Over all loci, the Middle American gene pool was more diverse than the Andean gene pool (π_sil=0.0089 vs 0.0068). The two subpopulations were strongly genetically differentiated over all loci (F_st=0.29). It is estimated that the two current wild gene pools diverged from a common ancestor ∼111 000 years ago. Subsequently, each gene pool underwent a bottleneck immediately after divergence and lasted ∼40 000 years. The Middle American bottleneck population size was ∼46% of the ancestral population size, whereas the Andean was 26%. Continuous asymmetric gene flow was detected between the two gene pools with a larger number of migrants entering Middle American gene pool from the Andean gene pool. These results suggest that because of the complex population structure associated with the ancestral divergence, subsequent bottlenecks in each gene pool, gene pool-specific domestication and intense selection within each gene pool by breeders; association mapping would best be practised within each common bean gene pool.     

Small Theropod Teeth from the Late Cretaceous of the San Juan Basin,Northwestern New Mexico and Their Implications for Understanding Latest Cretaceous Dinosaur Evolution

Studying the evolution and biogeographic distribution of dinosaurs during the latest Cretaceous is critical for better understanding the end-Cretaceous extinction event that killed off all non-avian dinosaurs. Western North America contains among the best records of Late Cretaceous terrestrial vertebrates in the world, but is biased against small-bodied dinosaurs. Isolated teeth are the primary evidence for understanding the diversity and evolution of small-bodied theropod dinosaurs during the Late Cretaceous, but few such specimens have been well documented from outside of the northern Rockies, making it difficult to assess Late Cretaceous dinosaur diversity and biogeographic patterns. We describe small theropod teeth from the San Juan Basin of northwestern New Mexico. These specimens were collected from strata spanning Santonian – Maastrichtian. We grouped isolated theropod teeth into several morphotypes, which we assigned to higher-level theropod clades based on possession of phylogenetic synapomorphies. We then used principal components analysis and discriminant function analyses to gauge whether the San Juan Basin teeth overlap with, or are quantitatively distinct from, similar tooth morphotypes from other geographic areas. The San Juan Basin contains a diverse record of small theropods. Late Campanian assemblages differ from approximately co-eval assemblages of the northern Rockies in being less diverse with only rare representatives of troodontids and a Dromaeosaurus-like taxon. We also provide evidence that erect and recurved morphs of a Richardoestesia-like taxon represent a single heterodont species. A late Maastrichtian assemblage is dominated by a distinct troodontid. The differences between northern and southern faunas based on isolated theropod teeth provide evidence for provinciality in the late Campanian and the late Maastrichtian of North America. However, there is no indication that major components of small-bodied theropod diversity were lost during the Maastrichtian in New Mexico. The same pattern seen in northern faunas, which may provide evidence for an abrupt dinosaur extinction.     

Chromosome-level genome assembly and characterization of Sophora Japonica

Sophora japonica is a medium-size deciduous tree belonging to Leguminosae family and famous for its high ecological, economic and medicinal value. Here, we reveal a draft genome of S. japonica, which was ∼511.49 Mb long (contig N50 size of 17.34 Mb) based on Illumina, Nanopore and Hi-C data. We reliably assembled 110 contigs into 14 chromosomes, representing 91.62% of the total genome, with an improved N50 size of 31.32 Mb based on Hi-C data. Further investigation identified 271.76 Mb (53.13%) of repetitive sequences and 31,000 protein-coding genes, of which 30,721 (99.1%) were functionally annotated. Phylogenetic analysis indicates that S. japonica separated from Arabidopsis thaliana and Glycine max ∼107.53 and 61.24 million years ago, respectively. We detected evidence of species-specific and common-legume whole-genome duplication events in S. japonica. We further found that multiple TF families (e.g. BBX and PAL) have expanded in S. japonica, which might have led to its enhanced tolerance to abiotic stress. In addition, S. japonica harbours more genes involved in the lignin and cellulose biosynthesis pathways than the other two species. Finally, population genomic analyses revealed no obvious differentiation among geographical groups and the effective population size continuously declined since 2 Ma. Our genomic data provide a powerful comparative framework to study the adaptation, evolution and active ingredients biosynthesis in S. japonica. More importantly, our high-quality S. japonica genome is important for elucidating the biosynthesis of its main bioactive components, and improving its production and/or processing.