Comprehensive multiomics analysis reveals distinct differences between pediatric choroid plexus papilloma and carcinoma.

Choroid plexus tumors (CPTs) are intraventricular tumors derived from the choroid plexus epithelium and occur frequently in children. The aim of this study was to investigate the genomic and epigenomic characteristics of CPT and identify the differences between choroid plexus papilloma (CPP) and choroid plexus carcinoma (CPC). We conducted multiomics analyses of 20 CPT patients including CPP and CPC. Multiomics analysis included whole-genome sequencing, whole-transcriptome sequencing, and methylation sequencing. Mutually exclusive TP53 and EPHA7 point mutations, coupled with the amplification of chromosome 1, were exclusively identified in CPC. In contrast, amplification of chromosome 9 was specific to CPP. Differential gene expression analysis uncovered a significant overexpression of genes related to cell cycle regulation and epithelial-mesenchymal transition pathways in CPC compared to CPP. Overexpression of genes associated with tumor metastasis and progression was observed in the CPC subgroup with leptomeningeal dissemination. Furthermore, methylation profiling unveiled hypomethylation in major repeat regions, including long interspersed nuclear elements, short interspersed nuclear elements, long terminal repeats, and retrotransposons in CPC compared to CPP, implying that the loss of epigenetic silencing of transposable elements may play a role in tumorigenesis of CPC. Finally, the differential expression of AK1, regulated by both genomic and epigenomic factors, emerged as a potential contributing factor to the histological difference of CPP against CPC. Our results suggest pronounced genomic and epigenomic disparities between CPP and CPC, providing insights into the pathogenesis of CPT at the molecular level.

Korea4K: whole genome sequences of 4,157 Koreans with 107 phenotypes derived from extensive health check-ups.

BACKGROUND: Phenome-wide association studies (PheWASs) have been conducted on Asian populations, including Koreans, but many were based on chip or exome genotyping data. Such studies have limitations regarding whole genome-wide association analysis, making it crucial to have genome-to-phenome association information with the largest possible whole genome and matched phenome data to conduct further population-genome studies and develop health care services based on population genomics. RESULTS: Here, we present 4,157 whole genome sequences (Korea4K) coupled with 107 health check-up parameters as the largest genomic resource of the Korean Genome Project. It encompasses most of the variants with allele frequency >0.001 in Koreans, indicating that it sufficiently covered most of the common and rare genetic variants with commonly measured phenotypes for Koreans. Korea4K provides 45,537,252 variants, and half of them were not present in Korea1K (1,094 samples). We also identified 1,356 new genotype-phenotype associations that were not found by the Korea1K dataset. Phenomics analyses further revealed 24 significant genetic correlations, 14 pleiotropic associations, and 127 causal relationships based on Mendelian randomization among 37 traits. In addition, the Korea4K imputation reference panel, the largest Korean variants reference to date, showed a superior imputation performance to Korea1K across all allele frequency categories. CONCLUSIONS: Collectively, Korea4K provides not only the largest Korean genome data but also corresponding health check-up parameters and novel genome-phenome associations. The large-scale pathological whole genome-wide omics data will become a powerful set for genome-phenome level association studies to discover causal markers for the prediction and diagnosis of health conditions in future studies.

Conservation implications of elucidating the Korean wolf taxonomic ambiguity through whole-genome sequencing.

The taxonomic status of the now likely extirpated Korean Peninsula wolf has been extensively debated, with some arguing it represents an independent wolf lineage, Canis coreanus. To investigate the Korean wolf's genetic affiliations and taxonomic status, we sequenced and analysed the genomes of a Korean wolf dated to the beginning of the 20th century, and a captive wolf originally from the Pyongyang Central Zoo. Our results indicated that the Korean wolf bears similar genetic ancestry to other regional East Asian populations, therefore suggesting it is not a distinct taxonomic lineage. We identified regional patterns of wolf population structure and admixture in East Asia with potential conservation consequences in the Korean Peninsula and on a regional scale. We find that the Korean wolf has similar genomic diversity and inbreeding to other East Asian wolves. Finally, we show that, in contrast to the historical sample, the captive wolf is genetically more similar to wolves from the Tibetan Plateau; hence, Korean wolf conservation programmes might not benefit from the inclusion of this specimen.

The complete mitochondrial genome of Rana coreana (Anura: Ranidae).

Rana coreana is a brown frog species native to the Korean Peninsula. We characterized the complete mitochondrial genome of the species. The mitochondrial genome sequence of R. coreana is 22,262 bp and comprises 13 protein-coding genes, two ribosomal RNA (rRNA) genes, 22 transfer RNA (tRNA) genes, and two control regions (CRs). The CR duplication and gene organization were identical to those observed in Rana kunyuensis and Rana amurensis. A total of 13 protein-coding genes were used to examine the phylogenetic relationships between this species and the genus Rana. R. coreana living on the Korean Peninsula, formed a cluster with R. kunyuensis and R. amurensis, with R. coreana showing the closest phylogenetic affinity for R. kunyuensis.

Identification and validation of six acute myocardial infarction-associated variants, including a novel prognostic marker for cardiac mortality.

Background: Acute myocardial infarction (AMI) is one of the leading causes of death worldwide, and approximately half of AMI-related deaths occur before the affected individual reaches the hospital. The present study aimed to identify and validate genetic variants associated with AMI and their role as prognostic markers. Materials and methods: We conducted a replication study of 29 previously identified novel loci containing 85 genetic variants associated with early-onset AMI using a new independent set of 2,920 Koreans [88 patients with early- and 1,085 patients with late-onset AMI, who underwent percutaneous coronary intervention (PCI), and 1,747 healthy controls]. Results: Of the 85 previously reported early-onset variants, six were confirmed in our genome-wide association study with a false discovery rate of less than 0.05. Notably, rs12639023, a cis-eQTL located in the intergenic region between LINC02005 and CNTN3, significantly increased longitudinal cardiac mortality and recurrent AMI. CNTN3 is known to play a role in altering vascular permeability. Another variant, rs78631167, located upstream of PLAUR and known to function in fibrinolysis, was moderately replicated in this study. By surveying the nearby genomic region around rs78631167, we identified a significant novel locus (rs8109584) located 13 bp downstream of rs78631167. The present study showed that six of the early-onset variants of AMI are applicable to both early- and late-onset cases. Conclusion: Our results confirm markers that can potentially be utilized to predict, screen, prevent, and treat candidate patients with AMI and highlight the potential of rs12639023 as a prognostic marker for cardiac mortality in AMI.

Germline DNA-Repair Genes and HOXB13 Mutations in Korean Men with Metastatic Prostate Cancer: Data from a Large Korean Cohort.

PURPOSE: Germline mutations in DNA damage repair (DDR) genes such as BRCA2 have been associated with prostate cancer (PC) risk but has not been thoroughly evaluated for metastatic prostate cancer (mPC) in Asian men. This study attempts to evaluate frequency of DDR mutations in the largest cohort of Koreans. MATERIALS AND METHODS: We recruited 340 patients with mPC unselected for family history of cancer and compared to 495 controls. Whole genome sequencing was applied to assess germline pathogenic/likely pathogenic variants (PV/LPVs) in 26 DDR genes and HOXB13, including 7 genes (ATM, BRCA1/2, CHEK2, BRIP1, PALB2, and NBN) associated with hereditary PC. Comparisons to published Caucasian and Japanese cohorts were performed. RESULTS: Total of 28 PV/LPVs were identified in 30 (8.8%) patients; mutations were found in 13 genes, including BRCA2 (15 men [4.41%]), ATM (2 men [0.59%]), NBN (2 men [0.59%], and BRIP1 (2 men [0.59%]). Only one patient had HOXB13 mutation (0.29%). A lower rate of overall germline variant frequency was observed in Korean mPC compared to Caucasians (8.8% vs. 11.8%), but individual variants notably differed from Caucasian and geographically similar Japanese cohorts. PV/LPVs in DDR genes tended to increase gradually with higher Gleason scores (GS 7, 7.1%; GS 8, 7.5%; GS 9-10, 9.9%). CONCLUSIONS: BRCA2 was the most frequently mutated gene common to different cohorts supporting its importance, but differences in variant distribution in Korean mPC underscore the need for ethnic-specific genetic models. Future ethnic-specific analyses are warranted to verify our findings.

Genome-wide analyses of early-onset acute myocardial infarction identify 29 novel loci by whole genome sequencing.

Early-onset acute myocardial infarction (AMI) may have a higher genetic predisposition than late-onset AMI. The present study aimed to identify and characterize germline variants that affect early-onset AMI using whole-genome sequencing (WGS). We performed a genome-wide association study based on the WGS of 1239 Koreans, including 596 early-onset AMI patients and 643 healthy individuals. Patients with AMI who underwent percutaneous coronary intervention (PCI) caused by atherothrombotic occlusive lesions were included in the study. A total of 29 novel loci were found to be associated with early-onset AMI. These loci are involved in thrombosis, fibrinolysis, inflammation, and lipid metabolism. One of the associated single nucleotide variants (SNVs), rs1614576, located upstream of PRKCB, is known to be associated with thrombus formation. Additionally, the results revealed a novel locus, rs78631167, located upstream of PLAUR which plays a critical role in regulating plasminogen activation and is related to fibrinolysis. The association between early-onset AMI and rs9357455, which is located upstream of PHACTR1 and regulates inflammation in AMI, was found. Moreover, we identified a lipid metabolism related genetic risk locus, rs5072, in the APOA1-AS gene. This study provides new evidence supporting the genetic association between early-onset AMI and thrombosis and fibrinolysis, as well as inflammation and lipid metabolism, by analyzing the whole-genome of 596 patients with early-onset AMI who have been treated with PCI. Our findings highlight potential genetic markers for the prediction and management of AMI, as well as for understanding the etiology of AMI.

A database of 5305 healthy Korean individuals reveals genetic and clinical implications for an East Asian population.

Despite substantial advances in disease genetics, studies to date have largely focused on individuals of European descent. This limits further discoveries of novel functional genetic variants in other ethnic groups. To alleviate the paucity of East Asian population genome resources, we established the Korean Variant Archive 2 (KOVA 2), which is composed of 1896 whole-genome sequences and 3409 whole-exome sequences from healthy individuals of Korean ethnicity. This is the largest genome database from the ethnic Korean population to date, surpassing the 1909 Korean individuals deposited in gnomAD. The variants in KOVA 2 displayed all the known genetic features of those from previous genome databases, and we compiled data from Korean-specific runs of homozygosity, positively selected intervals, and structural variants. In doing so, we found loci, such as the loci of ADH1A/1B and UHRF1BP1, that are strongly selected in the Korean population relative to other East Asian populations. Our analysis of allele ages revealed a correlation between variant functionality and evolutionary age. The data can be browsed and downloaded from a public website ( https://www.kobic.re.kr/kova/ ). We anticipate that KOVA 2 will serve as a valuable resource for genetic studies involving East Asian populations.

A method for early diagnosis of lung cancer from tumor originated DNA fragments using plasma cfDNA methylome and fragmentome profiles.

Early detection is critical for minimizing mortality from cancer. Plasma cell-free DNA (cfDNA) contains the signatures of tumor DNA, allowing us to quantify the signature and diagnose early-stage tumors. Here, we report a novel tumor fragment quantification method, TOF (Tumor Originated Fragment) for the diagnosis of lung cancer by quantifying and analyzing both the plasma cfDNA methylation patterns and fragmentomic signatures. TOF utilizes the amount of ctDNA predicted from the methylation density information of each cfDNA read mapped on 6243 lung-tumor-specific CpG markers. The 6243 tumor-specific markers were derived from lung tumor tissues by comparing them with corresponding normal tissues and healthy blood from public methylation data. TOF also utilizes two cfDNA fragmentomic signatures: 1) the short fragment ratio, and 2) the 5' end-motif profile. We used 298 plasma samples to analyze cfDNA signatures using enzymatic methyl-sequencing data from 201 lung cancer patients and 97 healthy controls. The TOF score showed 0.98 of the area under the curve in correctly classifying lung cancer from normal samples. The TOF score resolution was high enough to clearly differentiate even the early-stage non-small cell lung cancer patients from the healthy controls. The same was true for small cell lung cancer patients.

The application of circulating tumor cell and cell-free DNA liquid biopsies in ovarian cancer.

Ovarian cancer is the deadliest gynecological cancer. 70% of the cases are diagnosed at late stages with already developed metastases due to the absence of easily noticeable symptoms. Early-stage ovarian cancer has a good prognosis with a 5-year survival rate reaching 95%, hence the identification of effective biomarkers for early diagnosis is important. Advances in liquid biopsy-based methods can have a significant impact not just on the development of an efficient screening strategy, but also in clinical decision-making with additional molecular profiling and genetic alterations linked to therapy resistance. Despite the well-known advantages of liquid biopsy, there are still challenges that need to be addressed before its routine use in clinical practice. Various liquid biopsy-based biomarkers have been investigated in ovarian cancer; however, in this review, we are concentrating on the current use of cell-free DNA (cfDNA) and circulating tumor cells (CTCs) in disease management, focusing on their emerging importance in clinical practice. We also discuss the technical aspects of these workflows. The analysis of cfDNA is often chosen for the detection of mutations, copy number aberrations, and DNA methylation changes, whereas CTC analysis provides a unique opportunity to study whole cells, thus allowing DNA, RNA, and protein-based molecular profiling as well as in vivo studies. Combined solutions which merge the strengths of cfDNA and CTC approaches should be developed to maximize the potential of liquid biopsy technology.

PharmaKoVariome database for supporting genetic testing.

Pharmacogenomics (PGx) provides information about routine precision medicine, based on the patient's genotype. However, many of the available information about human allele frequencies, and about clinical drug-gene interactions, is based on American and European populations. PharmaKoVariome database was constructed to support genetic testing for safe prescription and drug development. It consolidated and stored 2507 diseases, 11 459 drugs and 61 627 drug-target or druggable genes from public databases. PharmaKoVariome precomputed ethnic-specific abundant variants for approximately 120 M single-nucleotide variants of drug-target or druggable genes. A user can search by gene symbol, drug name, disease and reference SNP ID number (rsID) to statistically analyse the frequency of ethnical variations, such as odds ratio and P-values for related genes. In an example study, we observed five Korean-enriched variants in the CYP2B6 and CYP2D6 genes, one of which (rs1065852) is known to be incapable of metabolizing drug. It is also shown that 4-6% of North and East Asians have risk factors for drugs metabolized by the CYP2D6 gene. Therefore, PharmaKoVariome is a useful database for pharmaceutical or diagnostic companies for developing diagnostic technologies that can be applied in the Asian PGx industry. Database URL: http://www.pharmakovariome.com/.

Comparative Genome and Evolution Analyses of an Endangered Stony Coral Species Dendrophyllia cribrosa Near Dokdo Islands in the East Sea.

Stony corals often harbor intracellular photosynthetic dinoflagellate algae that receive dissolved inorganic nutrients. However, Dendrophyllia cribrosa is a nonsymbiotic stony coral distributed in the western Pacific. We assembled a chromosome-level D. cribrosa genome using PacBio and Hi-C technologies. The final assembly was 625 Mb, distributed on 14 chromosomes, and contained 30,493 protein-coding genes. The Benchmarking Universal Single-Copy Orthologs analysis revealed a percentage of 96.8 of the metazoan genome. A comparative phylogenetic analysis revealed that D. cribrosa, which lacks symbionts, evolved to acquire cellular energy by expanding genes related to acyl-CoA metabolism and carbohydrate transporters. This species also has expanded immune-related genes involved in the receptor protein tyrosine kinase signaling pathway. In addition, we observed a specific expansion of calcification genes, such as coral acid-rich proteins and carbonic anhydrase, in D. cribrosa. This high-quality reference genome and comparative analysis provides insights into the ecology and evolution of nonsymbiotic stony corals.

Chromosome-Level Genome Assembly of the Butter Clam Saxidomus purpuratus.

Herein, we provide the first whole-genome sequence of the purple butter clam (Saxidomus purpuratus), an economically important bivalve shellfish. Specifically, we sequenced and de novo assembled the genome of Sa. purpuratus based on PromethION long reads and Hi-C data. The 978-Mb genome of Sa. purpuratus comprises 19 chromosomes with 36,591 predicted protein-coding genes. The N50 length of Sa. purpuratus genome is 52 Mb, showing the highest continuous assembly among bivalve genomes. The Benchmarking by Universal Single-Copy Orthologs assessment indicated that 95.07% of complete metazoan universal single-copy orthologs (n = 954) were present in the assembly. Approximately 51% of Sa. purpuratus genome comprises repetitive sequences. Based on the high-quality Sa. purpuratus genome, we resolved half of the immune-associated genes, namely, scavenger receptor (SR) proteins, which are collinear to those in the closely related Cyclina sinensis genome. This finding suggested a high degree of conservation among immune-associated genes. Twenty-two (19%) SR proteins are tandemly duplicated in Sa. purpuratus genome, suggesting putative convergence evolution. Overall, Sa. purpuratus genome provides a new resource for the discovery of economically important traits and immune-response genes.

Northeastern Asian and Jomon-related genetic structure in the Three Kingdoms period of Gimhae, Korea.

The genetic history of prehistoric and protohistoric Korean populations is not well understood because only a small number of ancient genomes are available. Here, we report the first paleogenomic data from the Korean Three Kingdoms period, a crucial point in the cultural and historic formation of Korea. These data comprise eight shotgun-sequenced genomes from ancient Korea (0.7×-6.1× coverage). They were derived from two archeological sites in Gimhae: the Yuha-ri shell mound and the Daesung-dong tumuli, the latter being the most important funerary complex of the Gaya confederacy. All individuals are from between the 4th and 5th century CE and are best modeled as an admixture between a northern China Bronze Age genetic source and a source of Jomon-related ancestry that shares similarities with the present-day genomes from Japan. The observed substructure and proportion of Jomon-related ancestry suggest the presence of two genetic groups within the population and diversity among the Gaya population. We could not correlate the genomic differences between these two groups with either social status or sex. All the ancient individuals' genomic profiles, including phenotypically relevant SNPs associated with hair and eye color, facial morphology, and myopia, imply strong genetic and phenotypic continuity with modern Koreans for the last 1,700 years.

The Kiss of Death: Serratia marcescens Antibacterial Activities against Staphylococcus aureus Requires Both de novo Prodigiosin Synthesis and Direct Contact.

Prodigiosin possesses antibacterial activities, but as a highly hydrophobic compound, it raised the question about how Serratia marcescens introduce this compound to other microbes. Here, we demonstrate that the production of prodigiosin by newly isolated S. marcescens RH10 correlates with its antibacterial activity against a multidrug-resistant strain of S. aureus, with this pathogen's viability decreasing 6-log over 24 h. While S. marcescens RH10 does secrete membrane vesicles that carry prodigiosin, this antibiotic was not active in this form, with 5 mg/L prodigiosin leading to only a 1.22-fold reduction in the S. aureus viability while the same quantity of purified prodigiosin led to a 2800-fold reduction. Contact assays, however, showed increased activity, with a 3-log loss in the S. aureus viabilities in only 6 h as long as de novo production of prodigiosin occurred. The role of prodigiosin was confirmed further by generating an isogenic ΔpigA mutant in S. marcescens RH10, based on the draft genome sequence reported here, to inhibit the synthesis of prodigiosin. In all experiments performed, this mutant was unable to kill S. aureus. Finally, the possibility that the type VI secretion system present in S. marcescens may also be important was also explored as it is known to be used by this strain to kill other microbes. The results here, however, found no obvious activity against S. aureus. In conclusion, the results presented here show prodigiosin requires both cell-to-cell contact and de novo synthesis for it to be effective as an antibiotic for its native host. IMPORTANCE The antibacterial activities of prodigiosin are well-established but, as a hydrophobic molecule, the mechanisms used to introduce it to susceptible microbes has never been studied. We found here, in contrast to violacein, another hydrophobic antibiotic that can be transferred using membrane vesicles (MVs), prodigiosin is also carried from Serratia marcescens in MVs released but its resulting activities were severely mitigated compared to the freely added compound, suggesting it is more tightly bound to the MVs than violacein. This led us to hypothesize that cell-to-cell contact is needed, which we demonstrate here. As well, we show de novo synthesis of prodigiosin is needed for it to be effective. As violacein- and prodigiosin-producing bacterial strains are both beneficial to amphibians, where they help protect the skin against pathogens, the findings presented here provide an important ecological perspective as they show the mechanisms used differ according to the antibacterial produced.

The Analyses of Cetacean Virus-Responsive Genes Reveal Evolutionary Marks in Mucosal Immunity-Associated Genes.

Viruses are the most common and abundant organisms in the marine environment. To better understand how cetaceans have adapted to this virus-rich environment, we compared cetacean virus-responsive genes to those from terrestrial mammals. We identified virus-responsive gene sequences in seven species of cetaceans, which we compared with orthologous sequences in seven terrestrial mammals. As a result of evolution analysis using the branch model and the branch-site model, 21 genes were selected using at least one model. IFN-ε, an antiviral cytokine expressed at mucous membranes, and its receptor IFNAR1 contain cetacean-specific amino acid substitutions that might change the interaction between the two proteins and lead to regulation of the immune system against viruses. Cetacean-specific amino acid substitutions in IL-6, IL-27, and the signal transducer and activator of transcription (STAT)1 are also predicted to alter the mucosal immune response of cetaceans. Since mucosal membranes are the first line of defense against the external environment and are involved in immune tolerance, our analysis of cetacean virus-responsive genes suggests that genes with cetacean-specific mutations in mucosal immunity-related genes play an important role in the protection and/or regulation of immune responses against viruses.

KOREF_S1: phased, parental trio-binned Korean reference genome using long reads and Hi-C sequencing methods.

BACKGROUND: KOREF is the Korean reference genome, which was constructed with various sequencing technologies including long reads, short reads, and optical mapping methods. It is also the first East Asian multiomic reference genome accompanied by extensive clinical information, time-series and multiomic data, and parental sequencing data. However, it was still not a chromosome-scale reference. Here, we updated the previous KOREF assembly to a new chromosome-level haploid assembly of KOREF, KOREF_S1v2.1. Oxford Nanopore Technologies (ONT) PromethION, Pacific Biosciences HiFi-CCS, and Hi-C technology were used to build the most accurate East Asian reference assembled so far. RESULTS: We produced 705 Gb ONT reads and 114 Gb Pacific Biosciences HiFi reads, and corrected ONT reads by Pacific Biosciences reads. The corrected ultra-long reads reached higher accuracy of 1.4% base errors than the previous KOREF_S1v1.0, which was mainly built with short reads. KOREF has parental genome information, and we successfully phased it using a trio-binning method, acquiring a near-complete haploid-assembly. The final assembly resulted in total length of 2.9 Gb with an N50 of 150 Mb, and the longest scaffold covered 97.3% of GRCh38's chromosome 2. In addition, the final assembly showed high base accuracy, with <0.01% base errors. CONCLUSIONS: KOREF_S1v2.1 is the first chromosome-scale haploid assembly of the Korean reference genome with high contiguity and accuracy. Our study provides useful resources of the Korean reference genome and demonstrates a new strategy of hybrid assembly that combines ONT's PromethION and PacBio's HiFi-CCS.

Sex-specific difference of in-hospital mortality from COVID-19 in South Korea.

We sought to assess the impact of sex on in-hospital mortality of patients with COVID-19 infection in South Korea. The study recruited 5,628 prospective consecutive patients who were hospitalized in South Korea with COVID-19 infection, and enrolled in the Korea Centers for Disease Control and Prevention (KCDC) dataset between January 20, 2020, and April 30, 2020. The primary endpoint was in-hospital death from COVID-19. The cohort comprised of 3,308 women (59%) and 2,320 men (41%). In-hospital death was significantly lower in women than men (3.5% vs. 5.5%, hazard ratio (HR): 0.61; 95% confidence interval (CI): 0.47 to 0.79, p <0.001). Results were consistent after multivariable regression (HR: 0.59; 95% CI: 0.41 to 0.85, p = 0.023) and propensity score matching (HR: 0.51; 95% CI: 0.30 to 0.86, p = 0.012). In South Korea, women had a significantly lower risk of in-hospital death amongst those patients hospitalized with COVID-19 infection.

LT1, an ONT long-read-based assembly scaffolded with Hi-C data and polished with short reads.

We present LT1, the first high-quality human reference genome from the Baltic States. LT1 is a female de novo human reference genome assembly, constructed using 57× nanopore long reads and polished using 47× short paired-end reads. We utilized 72 GB of Hi-C chromosomal mapping data for scaffolding, to maximize assembly contiguity and accuracy. The contig assembly of LT1 was 2.73 Gbp in length, comprising 4490 contigs with an NG50 value of 12.0 Mbp. After scaffolding with Hi-C data and manual curation, the final assembly has an NG50 value of 137 Mbp and 4699 scaffolds. Assessment of gene prediction quality using Benchmarking Universal Single-Copy Orthologs (BUSCO) identified 89.3% of the single-copy orthologous genes included in the benchmark. Detailed characterization of LT1 suggests it has 73,744 predicted transcripts, 4.2 million autosomal SNPs, 974,616 short indels, and 12,079 large structural variants. These data may be used as a benchmark for further in-depth genomic analyses of Baltic populations.