Original Genome Research Articles

Insights from the genomes of 4 diploid Camelina spp.

Plant evolution has been a complex process involving hybridization and polyploidization making understanding the origin and evolution of a plant's genome challenging even once a published genome is available. The oilseed crop, Camelina sativa (Brassicaceae), has a fully sequenced allohexaploid genome with 3 unknown ancestors. To better understand which extant species best represent the ancestral genomes that contributed to C. sativa's formation, we sequenced and assembled chromosome level draft genomes for 4 diploid members of Camelina: C. neglecta C. hispida var. hispida, C. hispida var. grandiflora, and C. laxa using long and short read data scaffolded with proximity data. We then conducted phylogenetic analyses on regions of synteny and on genes described for Arabidopsis thaliana, from across each nuclear genome and the chloroplasts to examine evolutionary relationships within Camelina and Camelineae. We conclude that C. neglecta is closely related to C. sativa's sub-genome 1 and that C. hispida var. hispida and C. hispida var. grandiflora are most closely related to C. sativa's sub-genome 3. Further, the abundance and density of transposable elements, specifically Helitrons, suggest that the progenitor genome that contributed C. sativa's sub-genome 3 maybe more similar to the genome of C. hispida var. hispida than that of C. hispida var. grandiflora. These diploid genomes show few structural differences when compared to C. sativa's genome indicating little change to chromosome structure following allopolyploidization. This work also indicates that C. neglecta and C. hispida are important resources for understanding the genetics of C. sativa and potential resources for crop improvement.

Genome origin and evolution of common oat.

Genome origin and evolution of common oat.

Invader at the edge - Genomic origins and physiological differences of round gobies across a steep urban salinity gradient.

Species invasions are a global problem of increasing concern, especially in highly connected aquatic environments. Despite this, salinity conditions can pose physiological barriers to their spread, and understanding them is important for management. In Scandinavia's largest cargo port, the invasive round goby (Neogobius melanostomus) is established across a steep salinity gradient. We used 12,937 SNPs to identify the genetic origin and diversity of three sites along the salinity gradient and round goby from western, central and northern Baltic Sea, as well as north European rivers. Fish from two sites from the extreme ends of the gradient were also acclimated to freshwater and seawater, and tested for respiratory and osmoregulatory physiology. Fish from the high-salinity environment in the outer port showed higher genetic diversity, and closer relatedness to the other regions, compared to fish from lower salinity upstream the river. Fish from the high-salinity site also had higher maximum metabolic rate, fewer blood cells and lower blood Ca2+. Despite these genotypic and phenotypic differences, salinity acclimation affected fish from both sites in the same way: seawater increased the blood osmolality and Na+ levels, and freshwater increased the levels of the stress hormone cortisol. Our results show genotypic and phenotypic differences over short spatial scales across this steep salinity gradient. These patterns of the physiologically robust round goby are likely driven by multiple introductions into the high-salinity site, and a process of sorting, likely based on behaviour or selection, along the gradient. This euryhaline fish risks spreading from this area, and seascape genomics and phenotypic characterization can inform management strategies even within an area as small as a coastal harbour inlet.

Origins, genomic structure and copy number variation of snake venom myotoxins

Crotamine, myotoxin a and homologs are short peptides that often comprise major fractions of rattlesnake venoms and have been extensively studied for their bioactive properties. These toxins are thought to be important for rapidly immobilizing mammalian prey and are implicated in serious, and sometimes fatal, responses to envenomation in humans. While high quality reference genomes for multiple venomous snakes are available, the loci that encode myotoxins have not been successfully assembled in any existing genome assembly. Here, we integrate new and existing genomic and transcriptomic data from the Prairie Rattlesnake (Crotalus viridis viridis) to reconstruct, characterize, and infer the chromosomal locations of myotoxin-encoding loci. We integrate long-read transcriptomics (Pacific Bioscience's Iso-Seq) and short-read RNA-seq to infer gene sequence diversity and characterize patterns of myotoxin and paralogous β-defensin expression across multiple tissues. We also identify two long non-coding RNA sequences which both encode functional myotoxins, demonstrating a newly discovered source of venom coding sequence diversity. We also integrate long-range mate-pair chromatin contact data and linked-read sequencing to infer the structure and chromosomal locations of the three myotoxin-like loci. Further, we conclude that the venom-associated myotoxin is located on chromosome 1 and is adjacent to non-venom paralogs. Consistent with this locus contributing to venom composition, we find evidence that the promoter of this gene is selectively open in venom gland tissue and contains transcription factor binding sites implicated in broad trans-regulatory pathways that regulate snake venoms. This study provides the best genomic reconstruction of myotoxin loci to date and raises questions about the physiological roles and interplay between myotoxin and related genes, as well as the genomic origins of snake venom variation.

2019 Novel Human Coronavirus Sars-Cov-2 And Covid-19: A Brief Review

The emergence of the 2019 novel coronavirus SARS-CoV-2, causing a highly contagious disease COVID-19 poses a tremendous global public health concern. As the disease is quickly developing into a worldwide pandemic, a careful analysis of its origin, evolution, transmission, and cellular mechanism is urgently needed to combat the deadly virus. Successful isolation of the 2019-nCoV has promoted some bioinformatic studies to understand the viral origin and the feature of its infectivity. However, at this stage, much remains unclear about the origin of the novel virus and to be investigated to develop ways to control its spread. The mechanisms associated with the infectiousness of SARS-CoV-2 are not entirely understood. The current knowledge in 2019-nCoV pathogenicity and transmissibility along with several commonly known emerging viruses and information is very much important for better control of the disease. So far, the treatment is only supportive. Any antiviral agent is yet to emerge. The scientific community raced to understand the pathogenesis of the disease for developing treatment options. Extensive research is going on to understand the host response to the pandemic virus to develop the disease therapeutic. To fill the knowledge gap about the human immune response to SARS-CoV- 2 infections that may help in designing the appropriate immune intervention for treatment, diagnosis, and prophylactic/therapeutic vaccines against COVID-19, we present hereby a brief review on the genomic organization, origin, and evolution, transmission and pathogenesis of the novel virus, and clinical spectrum and possible preventive measures against COVID-19 which may be used for future references.

Abstract 1532: Short RNAs v2.0: isomiRs, tRFs, and rRFs represent numerous novel opportunities for developing potent, precision-medicine-aware diagnostics, prognostics, and therapeutics

Abstract Background: It has been twenty-eight years since the first microRNA (miRNA) was reported in the worm C. elegans. Since then, researchers have discovered miRNAs in animals, plants, and viruses. They also showed that miRNAs are critical regulators of many cellular processes in homeostasis and disease. The advent of deep sequencing helped discover many novel miRNAs, pinpointed miRNAs that are dysregulated in cancers and other diseases, and established that miRNAs have isoforms (isomiRs). Deep sequencing also proved instrumental in discovering two more populous classes of short RNAs, the tRNA-derived fragments (tRFs) and the rRNA-derived fragments (rRFs). In the last decade, others and we generated strong evidence that isomiRs, tRFs, and rRFs regulate messenger RNA (mRNA) and protein abundance, directly and indirectly. Methods: IsomiRs, tRFs, and rRFs arise from precursors whose genomic distribution is non-random. Many of these precursors contain sequence segments that resemble one another's or are present in unrelated genomic sequences. These factors obfuscate the molecules’ identity and genomic origin, muddle their biogenesis details, and complicate the design of experiments. Thus, to help study isomiRs, tRFs, and rRFs we had to develop specialized tools that tackle these complications. We analyzed tens of thousands of human samples using these tools, including specimens of primary human tissues, cancer cell lines, and biofluids. Results: We were the first to show that the expression of isomiRs, tRFs, and rRFs is regimented and depends on tissue type and disease. We were also first to show that their expression additionally depends on genetic ancestry, sex, and age. Our findings highlight the involvement of isomiRs, tRFs, and rRFs in previously unsuspected regulatory interactions. These interactions are currently uncharacterized and depend on cell context and, unexpectedly, on personal attributes. These dependencies are consequential. Indeed, we showed for many cancers that specific isomiRs, tRFs, and rRFs predispose patients of a given sex or ancestry to aggressive cancer biology. Conclusions: Our pan-cancer analyses show that isomiRs, tRFs, and rRFs offer unparalleled granularity when studying normal physiology, cancer, and the molecular biology of health disparities. Improving our understanding of the molecules’ dependence on tissue, disease, and personal attributes (e.g., sex, ancestry, age) will help develop better diagnostics, prognostics, and new approaches to personalized therapies. Preliminary reports provide strong evidence in this regard. Because very little is known currently about these three classes of short RNAs, data-driven approaches will be essential to all these efforts. Citation Format: Isidore Rigoutsos. Short RNAs v2.0: isomiRs, tRFs, and rRFs represent numerous novel opportunities for developing potent, precision-medicine-aware diagnostics, prognostics, and therapeutics [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2022; 2022 Apr 8-13. Philadelphia (PA): AACR; Cancer Res 2022;82(12_Suppl):Abstract nr 1532.

Algorithms for Processing Coronavirus Genomes for the Goals and Objectives of Modern Immunoinformatics, Vaccinomics, and Virology

The novel coronavirus pandemic has stimulated the scientific activity of virology and interdisciplinary sciences: medical cybernetics and bioinformatics. The article is focused on the study of algorithms for processing bioinformatic data of genomic origin predominantly for the purposes of predominantly immunoinformatics and computational vaccinology. The schemes of algorithms developed by the authors for the analysis of bioinformatic data are presented. The algorithms for processing genomic information developed by the authors based on the analysis of the available literature and many years of experience in computational and laboratory experiments can be used not only for the design and analysis of epitope vaccine components, but also for the other tasks of computational virology and microbiology. In silico experiments on the analysis of bioinformatic data are relatively low-cost and multi-informative, but they require highly qualified scientists with extensive experience, interdisciplinary training, and, accordingly, a wide range of knowledge and skills. However, for the complete analysis and implementation of, for example, the epitope vaccines, subsequent validation by the laboratory and in vivo experiments are required.

Something Fishy about Siamese Fighting Fish (Betta splendens) Sex: Polygenic Sex Determination or a Newly Emerged Sex-Determining Region?

Fishes provide a unique and intriguing model system for studying the genomic origin and evolutionary mechanisms underlying sex determination and high sex-chromosome turnover. In this study, the mode of sex determination was investigated in Siamese fighting fish, a species of commercial importance. Genome-wide SNP analyses were performed on 75 individuals (40 males and 35 females) across commercial populations to determine candidate sex-specific/sex-linked loci. In total, 73 male-specific loci were identified and mapped to a 5.6 kb region on chromosome 9, suggesting a putative male-determining region (pMDR) containing localized dmrt1 and znrf3 functional sex developmental genes. Repeat annotations of the pMDR revealed an abundance of transposable elements, particularly Ty3/Gypsy and novel repeats. Remarkably, two out of the 73 male-specific loci were located on chromosomes 7 and 19, implying the existence of polygenic sex determination. Besides male-specific loci, five female-specific loci on chromosome 9 were also observed in certain populations, indicating the possibility of a female-determining region and the polygenic nature of sex determination. An alternative explanation is that male-specific loci derived from other chromosomes or female-specific loci in Siamese fighting fish recently emerged as new sex-determining loci during domestication and repeated hybridization.

Coffea cytogenetics: from the first karyotypes to the meeting with genomics.

Coffea karyotype organization and evolution has been uncovered by classical cytogenetics and cytogenomics. We revisit these discoveries and present new karyotype data. Coffea possesses ~ 124 species, including C. arabica and C. canephora responsible for commercial coffee production. We reviewed the Coffea cytogenetics, from the first chromosome counting, encompassing the karyotype characterization, chromosome DNA content, and mapping of chromosome portions and DNA sequences, until the integration with genomics. We also showed new data about Coffea karyotype. The 2n chromosome number evidenced the diploidy of almost all Coffea, and the C. arabica tetraploidy, as well as the polyploidy of other hybrids. Since then, other genomic similarities and divergences among the Coffea have been shown by karyotype morphology, nuclear and chromosomal C-value, AT and GC rich chromosome portions, and repetitive sequence and gene mapping. These cytogenomic data allowed us to know and understand the phylogenetic relations in Coffea, as well as their ploidy level and genomic origin, highlighting the relatively recent allopolyploidy. In addition to the euploidy, the role of the mobile elements in Coffea diversification is increasingly more evident, and the comparative analysis of their structure and distribution on the genome of different species is in the spotlight for future research. An integrative look at all these data is fundamental for a deeper understanding of Coffea karyotype evolution, including the key role of polyploidy in C. arabica origin. The 'Híbrido de Timor', a recent natural allotriploid, is also in the spotlight for its potential as a source of resistance genes and model for plant polyploidy research. Considering this, we also present some unprecedented results about the exciting evolutionary history of these polyploid Coffea.

The genomic origins of the world’s first farmers

SummaryThe precise genetic origins of the first Neolithic farming populations in Europe and Southwest Asia, as well as the processes and the timing of their differentiation, remain largely unknown. Demogenomic modeling of high-quality ancient genomes reveals that the early farmers of Anatolia and Europe emerged from a multiphase mixing of a Southwest Asian population with a strongly bottlenecked western hunter-gatherer population after the last glacial maximum. Moreover, the ancestors of the first farmers of Europe and Anatolia went through a period of extreme genetic drift during their westward range expansion, contributing highly to their genetic distinctiveness. This modeling elucidates the demographic processes at the root of the Neolithic transition and leads to a spatial interpretation of the population history of Southwest Asia and Europe during the late Pleistocene and early Holocene.

Discovery and characterization of sweetpotato's closest tetraploid relative.

Summary The origin of sweetpotato, a hexaploid species, is poorly understood, partly because the identity of its tetraploid progenitor remains unknown. In this study, we identify, describe and characterize a new species of Ipomoea that is sweetpotato’s closest tetraploid relative known to date and probably a direct descendant of its tetraploid progenitor.We integrate morphological, phylogenetic, and genomic analyses of herbarium and germplasm accessions of the hexaploid sweetpotato, its closest known diploid relative Ipomoea trifida, and various tetraploid plants closely related to them from across the American continent.We identify wild autotetraploid plants from Ecuador that are morphologically distinct from Ipomoea batatas and I. trifida, but monophyletic and sister to I. batatas in phylogenetic analysis of nuclear data.We describe this new species as Ipomoea aequatoriensis T. Wells & P. Muñoz sp. nov., distinguish it from hybrid tetraploid material collected in Mexico; and show that it likely played a direct role in the origin of sweetpotato’s hexaploid genome. This discovery transforms our understanding of sweetpotato’s origin.

Long-read sequencing reveals atypical mitochondrial genome structure in a New Zealand marine isopod.

Most animal mitochondrial genomes are small, circular and structurally conserved. However, recent work indicates that diverse taxa possess unusual mitochondrial genomes. In Isopoda, species in multiple lineages have atypical and rearranged mitochondrial genomes. However, more species of this speciose taxon need to be evaluated to understand the evolutionary origins of atypical mitochondrial genomes in this group. In this study, we report the presence of an atypical mitochondrial structure in the New Zealand endemic marine isopod, Isocladus armatus. Data from long- and short-read DNA sequencing suggest that I. armatus has two mitochondrial chromosomes. The first chromosome consists of two mitochondrial genomes that have been inverted and fused together in a circular form, and the second chromosome consists of a single mitochondrial genome in a linearized form. This atypical mitochondrial structure has been detected in other isopod lineages, and our data from an additional divergent isopod lineage (Sphaeromatidae) lends support to the hypothesis that atypical structure evolved early in the evolution of Isopoda. Additionally, we find that an asymmetrical site previously observed across many species within Isopoda is absent in I. armatus, but confirm the presence of two asymmetrical sites recently reported in two other isopod species.

Assembled and annotated 26.5 Gbp coast redwood genome: a resource for estimating evolutionary adaptive potential and investigating hexaploid origin

Sequencing, assembly, and annotation of the 26.5 Gbp hexaploid genome of coast redwood (Sequoia sempervirens) was completed leading toward discovery of genes related to climate adaptation and investigation of the origin of the hexaploid genome. Deep-coverage short-read Illumina sequencing data from haploid tissue from a single seed were combined with long-read Oxford Nanopore Technologies sequencing data from diploid needle tissue to create an initial assembly, which was then scaffolded using proximity ligation data to produce a highly contiguous final assembly, SESE 2.1, with a scaffold N50 size of 44.9 Mbp. The assembly included several scaffolds that span entire chromosome arms, confirmed by the presence of telomere and centromere sequences on the ends of the scaffolds. The structural annotation produced 118,906 genes with 113 containing introns that exceed 500 Kbp in length and one reaching 2 Mb. Nearly 19 Gbp of the genome represented repetitive content with the vast majority characterized as long terminal repeats, with a 2.9:1 ratio of Copia to Gypsy elements that may aid in gene expression control. Comparison of coast redwood to other conifers revealed species-specific expansions for a plethora of abiotic and biotic stress response genes, including those involved in fungal disease resistance, detoxification, and physical injury/structural remodeling and others supporting flavonoid biosynthesis. Analysis of multiple genes that exist in triplicate in coast redwood but only once in its diploid relative, giant sequoia, supports a previous hypothesis that the hexaploidy is the result of autopolyploidy rather than any hybridizations with separate but closely related conifer species.

Quantitative assessment reveals the dominance of duplicated sequences in germline-derived extrachromosomal circular DNA

Extrachromosomal circular DNA (eccDNA) originates from linear chromosomal DNA in various human tissues under physiological and disease conditions. The genomic origins of eccDNA have largely been investigated using invitro-amplified DNA. However, invitro amplification obscures quantitative information by skewing the total population stoichiometry. In addition, the analyses have focused on eccDNA stemming from single-copy genomic regions, leaving eccDNA from multicopy regions unexamined. To address these issues, we isolated eccDNA without invitro amplification (naïve small circular DNA, nscDNA) and assessed the populations quantitatively by integrated genomic, molecular, and cytogenetic approaches. nscDNA of up to tens of kilobases were successfully enriched by our approach and were predominantly derived from multicopy genomic regions including segmental duplications (SDs). SDs, which account for 5% of the human genome and are hotspots for copy number variations, were significantly overrepresented in sperm nscDNA, with three times more sequencing reads derived from SDs than from the entire single-copy regions. SDs were also overrepresented in mouse sperm nscDNA, which we estimated to comprise 0.2% of nuclear DNA. Considering that eccDNA can be integrated into chromosomes, germline-derived nscDNA may be a mediator of genome diversity.

Propagation of a De Novo Gene under Natural Selection: Antifreeze Glycoprotein Genes and Their Evolutionary History in Codfishes.

The de novo birth of functional genes from non-coding DNA as an important contributor to new gene formation is increasingly supported by evidence from diverse eukaryotic lineages. However, many uncertainties remain, including how the incipient de novo genes would continue to evolve and the molecular mechanisms underlying their evolutionary trajectory. Here we address these questions by investigating evolutionary history of the de novo antifreeze glycoprotein (AFGP) gene and gene family in gadid (codfish) lineages. We examined AFGP phenotype on a phylogenetic framework encompassing a broad sampling of gadids from freezing and non-freezing habitats. In three select species representing different AFGP-bearing clades, we analyzed all AFGP gene family members and the broader scale AFGP genomic regions in detail. Codon usage analyses suggest that motif duplication produced the intragenic AFGP tripeptide coding repeats, and rapid sequence divergence post-duplication stabilized the recombination-prone long repetitive coding region. Genomic loci analyses support AFGP originated once from a single ancestral genomic origin, and shed light on how the de novo gene proliferated into a gene family. Results also show the processes of gene duplication and gene loss are distinctive in separate clades, and both genotype and phenotype are commensurate with differential local selective pressures.

Heterogeneity among enhancer RNAs: origins, consequences and perspectives.

Enhancer RNAs (eRNAs) are non-coding RNAs transcribed from distal cis-regulatory elements (i.e. enhancers), which are stereotyped as short, rarely spliced and unstable. In fact, a non-negligible fraction of eRNAs seems to be longer, spliced and more stable, and their cognate enhancers are epigenomically and functionally distinguishable from typical enhancers. In this review, we first summarized the genomic and molecular origins underlying the observed heterogeneity among eRNAs. Then, we discussed how their heterogeneous properties (e.g. stability) affect the modes of interaction with their regulatory partners, from promiscuous cis-interactions to specific trans-interactions. Finally, we highlighted the existence of a seemingly continuous spectrum of eRNA properties and its implications in the genomic origins of non-coding RNA genes from an evolutionary perspective.

Genome variability, species diversity, phylogenetic relationships, origin, and geographical distribution of badnaviruses involved in the cacao swollen shoot disease: the case of West Africa

Genome variability, species diversity, phylogenetic relationships, origin, and geographical distribution of badnaviruses involved in the cacao swollen shoot disease: the case of West Africa

DNA Repair Inhibition Leads to Active Export of Repetitive Sequences to the Cytoplasm Triggering an Inflammatory Response.

Adult-onset neurodegenerative diseases are often accompanied by evidence of a chronic inflammation that includes activation of microglial cells and altered levels of brain cytokines. Aspects of this response are likely secondary reactions to neurodegeneration, but for many illnesses the inflammation may itself be an early and even causative disease event. In such cases, the inflammation is referred to as "sterile" as it occurs in the absence of an actual bacterial or viral pathogen. A potent trigger of sterile inflammation in CNS microglia has been shown to be the presence of DNA in the cytoplasm (cytoDNA) induced either by direct DNA damage or by inhibited DNA repair. We have shown that cytoDNA comes from the cell nucleus as a result of insufficient DNA damage repair. Using wild-type and Atm-/- mouse microglia, we extend these observations here by showing that its genomic origins are not random, but rather are heavily biased toward transcriptionally inactive, intergenic regions, in particular repetitive elements and AT-rich sequences. Once released from the genome, in both males and females, we show that cytoDNA is actively exported to the cytoplasm by a CRM1-dependent mechanism. In the cytoplasm, it is degraded either by a cytosolic exonuclease, Trex1, or an autophagy pathway that ends with degradation in the lysosome. Blocking the accumulation of cytoDNA prevents the emergence of the sterile inflammation reaction. These findings offer new insights into the emergence of sterile inflammation and offer novel approaches that may be of use in combatting a wide range of neurodegenerative conditions.SIGNIFICANCE STATEMENT Sterile inflammation describes a state where the defenses of the immune system are activated in the absence of a true pathogen. A potent trigger of this unorthodox response is the presence of DNA in the cytoplasm, which immune cells interpret as an invading virus or pathogen. We show that when DNA damage increases, fragments of the cell's own genome are actively exported to the cytoplasm where they are normally degraded. If this degradation is incomplete an immune reaction is triggered. Both age and stress increase DNA damage, and as age-related neurodegenerative diseases are frequently accompanied by a chronic low-level inflammation, strategies that reduce the induction of cytoplasmic DNA or speed its clearance become attractive therapeutic targets.

Pseudomonas Phage MD8: Genetic Mosaicism and Challenges of Taxonomic Classification of Lambdoid Bacteriophages.

Pseudomonas phage MD8 is a temperate phage isolated from the freshwater lake Baikal. The organisation of the MD8 genome resembles the genomes of lambdoid bacteriophages. However, MD8 gene and protein sequences have little in common with classified representatives of lambda-like phages. Analysis of phage genomes revealed a group of other Pseudomonas phages related to phage MD8 and the genomic layout of MD8-like phages indicated extensive gene exchange involving even the most conservative proteins and leading to a high degree of genomic mosaicism. Multiple horizontal transfers and mosaicism of the genome of MD8, related phages and other λ-like phages raise questions about the principles of taxonomic classification of the representatives of this voluminous phage group. Comparison and analysis of various bioinformatic approaches applied to λ-like phage genomes demonstrated different efficiency and contradictory results in the estimation of genomic similarity and relatedness. However, we were able to make suggestions for the possible origin of the MD8 genome and the basic principles for the taxonomic classification of lambdoid phages. The group comprising 26 MD8-related phages was proposed to classify as two close genera belonging to a big family of λ-like phages.

Ancient genomics reveals tripartite origins of Japanese populations.

Prehistoric Japan underwent rapid transformations in the past 3000 years, first from foraging to wet rice farming and then to state formation. A long-standing hypothesis posits that mainland Japanese populations derive dual ancestry from indigenous Jomon hunter-gatherer-fishers and succeeding Yayoi farmers. However, the genomic impact of agricultural migration and subsequent sociocultural changes remains unclear. We report 12 ancient Japanese genomes from pre- and postfarming periods. Our analysis finds that the Jomon maintained a small effective population size of ~1000 over several millennia, with a deep divergence from continental populations dated to 20,000 to 15,000 years ago, a period that saw the insularization of Japan through rising sea levels. Rice cultivation was introduced by people with Northeast Asian ancestry. Unexpectedly, we identify a later influx of East Asian ancestry during the imperial Kofun period. These three ancestral components continue to characterize present-day populations, supporting a tripartite model of Japanese genomic origins.