Whole Genome Analyses of Chinese Population and De Novo Assembly of A Northern Han Genome

KAT8 beyond Acetylation: A Survey of Its Epigenetic Regulation, Genetic Variability, and Implications for Human Health

Lysine acetyltransferase 8, also known as KAT8, is an enzyme involved in epigenetic regulation, primarily recognized for its ability to modulate histone acetylation. This review presents an overview of KAT8, emphasizing its biological functions, which impact many cellular processes and range from chromatin remodeling to genetic and epigenetic regulation. In many model systems, KAT8's acetylation of histone H4 lysine 16 (H4K16) is critical for chromatin structure modification, which influences gene expression, cell proliferation, differentiation, and apoptosis. Furthermore, this review summarizes the observed genetic variability within the KAT8 gene, underscoring the implications of various single nucleotide polymorphisms (SNPs) that affect its functional efficacy and are linked to diverse phenotypic outcomes, ranging from metabolic traits to neurological disorders. Advanced insights into the structural biology of KAT8 reveal its interaction with multiprotein assemblies, such as the male-specific lethal (MSL) and non-specific lethal (NSL) complexes, which regulate a wide range of transcriptional activities and developmental functions. Additionally, this review focuses on KAT8's roles in cellular homeostasis, stem cell identity, DNA damage repair, and immune response, highlighting its potential as a therapeutic target. The implications of KAT8 in health and disease, as evidenced by recent studies, affirm its importance in cellular physiology and human pathology.

Human pangenome analysis of sequences missing from the reference genome reveals their widespread evolutionary, phenotypic, and functional roles

Nonreference sequences (NRSs) are DNA sequences present in global populations but absent in the current human reference genome. However, the extent and functional significance of NRSs in the human genomes and populations remains unclear. Here, we de novo assembled 539 genomes from five genetically divergent human populations using long-read sequencing technology, resulting in the identification of 5.1 million NRSs. These were merged into 45284 unique NRSs, with 29.7% being novel discoveries. Among these NRSs, 38.7% were common across the five populations, and 35.6% were population specific. The use of a graph-based pangenome approach allowed for the detection of 565 transcript expression quantitative trait loci on NRSs, with 426 of these being novel findings. Moreover, 26 NRS candidates displayed evidence of adaptive selection within human populations. Genes situated in close proximity to or intersecting with these candidates may be associated with metabolism and type 2 diabetes. Genome-wide association studies revealed 14 NRSs to be significantly associated with eight phenotypes. Additionally, 154 NRSs were found to be in strong linkage disequilibrium with 258 phenotype-associated SNPs in the GWAS catalogue. Our work expands the understanding of human NRSs and provides novel insights into their functions, facilitating evolutionary and biomedical researches.

Haplotype-resolved assemblies and variant benchmark of a Chinese Quartet

BACKGROUND

Recent state-of-the-art sequencing technologies enable the investigation of challenging regions in the human genome and expand the scope of variant benchmarking datasets. Herein, we sequence a Chinese Quartet, comprising two monozygotic twin daughters and their biological parents, using four short and long sequencing platforms (Illumina, BGI, PacBio, and Oxford Nanopore Technology).

RESULTS

The long reads from the monozygotic twin daughters are phased into paternal and maternal haplotypes using the parent-child genetic map and for each haplotype. We also use long reads to generate haplotype-resolved whole-genome assemblies with completeness and continuity exceeding that of GRCh38. Using this Quartet, we comprehensively catalogue the human variant landscape, generating a dataset of 3,962,453 SNVs, 886,648 indels (< 50 bp), 9726 large deletions (≥ 50 bp), 15,600 large insertions (≥ 50 bp), 40 inversions, 31 complex structural variants, and 68 de novo mutations which are shared between the monozygotic twin daughters. Variants underrepresented in previous benchmarks owing to their complexity-including those located at long repeat regions, complex structural variants, and de novo mutations-are systematically examined in this study.

CONCLUSIONS

In summary, this study provides high-quality haplotype-resolved assemblies and a comprehensive set of benchmarking resources for two Chinese monozygotic twin samples which, relative to existing benchmarks, offers expanded genomic coverage and insight into complex variant categories.

T2T-YAO: A Telomere-to-telomere Assembled Diploid Reference Genome for Han Chinese

Since its initial release in 2001, the human reference genome has undergone continuous improvement in quality, and the recently released telomere-to-telomere (T2T) version - T2T-CHM13 - reaches its highest level of continuity and accuracy after 20 years of effort by working on a simplified, nearly homozygous genome of a hydatidiform mole cell line. Here, to provide an authentic complete diploid human genome reference for the Han Chinese, the largest population in the world, we assembled the genome of a male Han Chinese individual, T2T-YAO, which includes T2T assemblies of all the 22 + X + M and 22 + Y chromosomes in both haploids. The quality of T2T-YAO is much better than those of all currently available diploid assemblies, and its haploid version, T2T-YAO-hp, generated by selecting the better assembly for each autosome, reaches the top quality of fewer than one error per 29.5 Mb, even higher than that of T2T-CHM13. Derived from an individual living in the aboriginal region of the Han population, T2T-YAO shows clear ancestry and potential genetic continuity from the ancient ancestors. Each haplotype of T2T-YAO possesses ∼ 330-Mb exclusive sequences, ∼ 3100 unique genes, and tens of thousands of nucleotide and structural variations as compared with CHM13, highlighting the necessity of a population-stratified reference genome. The construction of T2T-YAO, an accurate and authentic representative of the Chinese population, would enable precise delineation of genomic variations and advance our understandings in the hereditability of diseases and phenotypes, especially within the context of the unique variations of the Chinese population.

Recent positive selection signatures reveal phenotypic evolution in the Han Chinese population

Characterizing natural selection signatures and relationships with phenotype spectra is important for understanding human evolution and both biological and pathological mechanisms. Here, we identified 24 genetic loci under recent selection by analyzing rare singletons in 3946 high-depth whole-genome sequencing data of Han Chinese. The loci include immune-related gene regions (MHC cluster, IGH cluster, STING1, and PSG), alcohol metabolism-related gene regions (ADH1B, ALDH2, and ALDH3B2), and the olfactory perception gene OR4C16, in which the MHC cluster, ADH1B, and ALDH2 were also identified by TOPMed and WestLake Biobank. Among the signals, the IGH cluster is particularly interesting, in which the favored allele of variant 14_105737776_C_T (rs117518546, IgG1-G396R) promotes immune response, but also increases the risk of an autoimmune disease systemic lupus erythematosus (SLE). It is also surprising that our newly discovered ALDH3B2 evolved in the opposite direction to ALDH2 for alcohol metabolism. Besides monogenic traits, we found that multiple complex traits experienced polygenic adaptation. Particularly, multi-methods consistently revealed that lower blood pressure was favored in natural selection. Finally, we built a database named RePoS (recent positive selection, http://bigdata.ibp.ac.cn/RePoS/) to integrate and display multi-population selection signals. Our study extended our understanding of natural evolution and phenotype adaptation in Han Chinese as well as other populations.

A pangenome reference of 36 Chinese populations

Human genomics is witnessing an ongoing paradigm shift from a single reference sequence to a pangenome form, but populations of Asian ancestry are underrepresented. Here we present data from the first phase of the Chinese Pangenome Consortium, including a collection of 116 high-quality and haplotype-phased de novo assemblies based on 58 core samples representing 36 minority Chinese ethnic groups. With an average 30.65× high-fidelity long-read sequence coverage, an average contiguity N50 of more than 35.63 megabases and an average total size of 3.01 gigabases, the CPC core assemblies add 189 million base pairs of euchromatic polymorphic sequences and 1,367 protein-coding gene duplications to GRCh38. We identified 15.9 million small variants and 78,072 structural variants, of which 5.9 million small variants and 34,223 structural variants were not reported in a recently released pangenome reference¹. The Chinese Pangenome Consortium data demonstrate a remarkable increase in the discovery of novel and missing sequences when individuals are included from underrepresented minority ethnic groups. The missing reference sequences were enriched with archaic-derived alleles and genes that confer essential functions related to keratinization, response to ultraviolet radiation, DNA repair, immunological responses and lifespan, implying great potential for shedding new light on human evolution and recovering missing heritability in complex disease mapping.

The first gapless, reference-quality, fully annotated genome from a Southern Han Chinese individual

We used long-read DNA sequencing to assemble the genome of a Southern Han Chinese male. We organized the sequence into chromosomes and filled in gaps using the recently completed T2T-CHM13 genome as a guide, yielding a gap-free genome, Han1, containing 3,099,707,698 bases. Using the T2T-CHM13 annotation as a reference, we mapped all genes onto the Han1 genome and identified additional gene copies, generating a total of 60,708 putative genes, of which 20,003 are protein-coding. A comprehensive comparison between the genes revealed that 235 protein-coding genes were substantially different between the individuals, with frameshifts or truncations affecting the protein-coding sequence. Most of these were heterozygous variants in which one gene copy was unaffected. This represents the first gene-level comparison between two finished, annotated individual human genomes.

Core promoter in TNBC is highly mutated with rich ethnic signature

The core promoter plays an essential role in regulating transcription initiation by controlling the interaction between transcriptional factors and sequence motifs in the core promoter. Although mutation in core promoter sequences is expected to cause abnormal gene expression leading to pathogenic consequences, limited supporting evidence showed the involvement of core promoter mutation in diseases. Our previous study showed that the core promoter is highly polymorphic in worldwide human ethnic populations in reflecting human history and adaptation. Our recent characterization of the core promoter in triple-negative breast cancer (TNBC), a subtype of breast cancer, in a Chinese TNBC cohort revealed the wide presence of core promoter mutation in TNBC. In the current study, we analyzed the core promoter in a Thai TNBC cohort. We also observed rich core promoter mutation in the Thai TNBC patients. We compared the core promoter mutations between Chinese and Thai TNBC cohorts. We observed substantial differences of core promoter mutation in TNBC between the two cohorts, as reflected by the mutation spectrum, mutation-effected gene and functional category, and altered gene expression. Our study confirmed that the core promoter in TNBC is highly mutable, and is highly ethnic-specific.

Haplotype-resolved de novo assembly of a Tujia genome suggests the necessity for high-quality population-specific genome references

Even though the human reference genome assembly is continually being improved, it remains debatable whether a population-specific reference is necessary for every ethnic group. Here, we de novo assembled an individual genome (TJ1) from the Tujia population, an ethnic minority group most closely related to the Han Chinese. TJ1 provided a high-quality haplotype-resolved assembly of chromosome-scale with a scaffold N50 size >78 Mb. Compared with GRCh38 and other de novo assemblies, TJ1 improved short-read mapping, enhanced calling precision for structural variants, and detected rare and low-frequency variants. This revealed fine-scale differences between the closely related Han and Tujia populations, such as population-stratified variants of LCT and UBXN8, and improved screening for ancestry informative markers. We demonstrated that TJ1 could reduce false positives in clinical diagnosis and analyzed the PRSS1-PRSS2 locus as a test case. Our results suggest that population-specific assemblies are necessary for genetic and medical analysis, especially when closely related populations are studied. A record of this paper's transparent peer review process is included in the supplemental information.

NyuWa Genome resource: A deep whole-genome sequencing-based variation profile and reference panel for the Chinese population

The lack of haplotype reference panels and whole-genome sequencing resources specific to the Chinese population has greatly hindered genetic studies in the world's largest population. Here, we present the NyuWa genome resource, based on deep (26.2×) sequencing of 2,999 Chinese individuals, and construct a NyuWa reference panel of 5,804 haplotypes and 19.3 million variants, which is a high-quality publicly available Chinese population-specific reference panel with thousands of samples. Compared with other panels, the NyuWa reference panel reduces the Han Chinese imputation error rate by a margin ranging from 30% to 51%. Population structure and imputation simulation tests support the applicability of one integrated reference panel for northern and southern Chinese. In addition, a total of 22,504 loss-of-function variants in coding and noncoding genes are identified, including 11,493 novel variants. These results highlight the value of the NyuWa genome resource in facilitating genetic research in Chinese and Asian populations.

Structural variants in the Chinese population and their impact on phenotypes, diseases and population adaptation

A complete characterization of genetic variation is a fundamental goal of human genome research. Long-read sequencing has improved the sensitivity of structural variant discovery. Here, we conduct the long-read sequencing-based structural variant analysis for 405 unrelated Chinese individuals, with 68 phenotypic and clinical measurements. We discover a landscape of 132,312 nonredundant structural variants, of which 45.2% are novel. The identified structural variants are of high-quality, with an estimated false discovery rate of 3.2%. The concatenated length of all the structural variants is approximately 13.2% of the human reference genome. We annotate 1,929 loss-of-function structural variants affecting the coding sequence of 1,681 genes. We discover rare deletions in HBA1/HBA2/HBB associated with anemia. Furthermore, we identify structural variants related to immunity which differentiate the northern and southern Chinese populations. Our study describes the landscape of structural variants in the Chinese population and their contribution to phenotypes and disease.

RefRGim: an intelligent reference panel reconstruction method for genotype imputation with convolutional neural networks

Genotype imputation is a statistical method for estimating missing genotypes from a denser haplotype reference panel. Existing methods usually performed well on common variants, but they may not be ideal for low-frequency and rare variants. Previous studies showed that the population similarity between study and reference panels is one of the key factors influencing the imputation accuracy. Here, we developed an imputation reference panel reconstruction method (RefRGim) using convolutional neural networks (CNNs), which can generate a study-specified reference panel for each input data based on the genetic similarity of individuals from current study and references. The CNNs were pretrained with single nucleotide polymorphism data from the 1000 Genomes Project. Our evaluations showed that genotype imputation with RefRGim can achieve higher accuracies than original reference panel, especially for low-frequency and rare variants. RefRGim will serve as an efficient reference panel reconstruction method for genotype imputation. RefRGim is freely available via GitHub: https://github.com/shishuo16/RefRGim

Study on Environmental and Lifestyle Factors for the North-South Differential of Cardiovascular Disease in China

Human death and life span are closely related to the geographical environment and regional lifestyle. These factors considerably vary among counties and regions, leading to the geographical disparity of disease. Quantitative studies on this phenomenon are insufficient. Cerebrovascular and heart diseases are the leading causes of death. The mortality rate of cerebrovascular and heart diseases is statistically higher in northern China than in southern China; the p-value of t-test for cerebrovascular and heart diseases was 0.047 and 0.000, respectively. The population attribution fraction of 12 major risk factors for cardiovascular disease (CVD) in each province was calculated based on their exposure and relative risk. The results found that residents in northern China consume high sodium-containing food, fewer vegetables, and less sea food products, and tend to be overweight. Fine particulate matter is higher in northern China than in southern China. Cold temperatures also cause a greater number of deaths than hot temperatures. All these factors have resulted in a higher CVD mortality rate in northern China. The attributive differential for sodium, vegetable, fruit, smoking, PM_2.5, omega-3, obesity, low temperature, and high temperature of heart disease between the two parts of China is 9.1, 0.7, -2.5, 0.1, 1.4, 1.3, 2.0, 4.7, and -2.1%, respectively. Furthermore, the attributive differential for the above factors of cerebrovascular disease between the two parts of China is 8.7, 0.0, -5.2, 0.1, 1.0, 0.0, 2.4, 4.7, and -2.1%. Diet high in sodium is the leading cause of the north-south differential in CVD, resulting in 0.71 less years of life expectancy in northern compared with that in southern China.

Profiling variable-number tandem repeat variation across populations using repeat-pangenome graphs

Variable number tandem repeats (VNTRs) are composed of consecutive repetitive DNA with hypervariable repeat count and composition. They include protein coding sequences and associations with clinical disorders. It has been difficult to incorporate VNTR analysis in disease studies that use short-read sequencing because the traditional approach of mapping to the human reference is less effective for repetitive and divergent sequences. In this work, we solve VNTR mapping for short reads with a repeat-pangenome graph (RPGG), a data structure that encodes both the population diversity and repeat structure of VNTR loci from multiple haplotype-resolved assemblies. We develop software to build a RPGG, and use the RPGG to estimate VNTR composition with short reads. We use this to discover VNTRs with length stratified by continental population, and expression quantitative trait loci, indicating that RPGG analysis of VNTRs will be critical for future studies of diversity and disease.

Prevalence and spectrum of DNA mismatch repair gene variation in the general Chinese population

Background

Identifying genetic disease-susceptible individuals through population screening is considered as a promising approach for disease prevention. DNA mismatch repair (MMR) genes including MLH1, MSH2, MSH6 and PMS2 play essential roles in maintaining microsatellite stability through DNA mismatch repair, and pathogenic variation in MMR genes causes microsatellite instability and is the genetic predisposition for cancer as represented by the Lynch syndrome. While the prevalence and spectrum of MMR variation has been extensively studied in cancer, it remains largely elusive in the general population. Lack of the knowledge prevents effective prevention for MMR variation–caused cancer. In the current study, we addressed the issue by using the Chinese population as a model.

Methods

We performed extensive data mining to collect MMR variant data from 18 844 ethnic Chinese individuals and comprehensive analyses for the collected MMR variants to determine its prevalence, spectrum and features of the MMR data in the Chinese population.

Results

We identified 17 687 distinct MMR variants. We observed substantial differences of MMR variation between the general Chinese population and Chinese patients with cancer, identified highly Chinese-specific MMR variation through comparing MMR data between Chinese and non-Chinese populations, predicted the enrichment of deleterious variants in the unclassified Chinese-specific MMR variants, determined MMR pathogenic prevalence of 0.18% in the general Chinese population and determined that MMR variation in the general Chinese population is evolutionarily neutral.

Conclusion

Our study provides a comprehensive view of MMR variation in the general Chinese population, a resource for biological study of human MMR variation, and a reference for MMR-related cancer applications.

Biotech in China 2021, at the beginning of the 14th five-year period ("145")

Abstract

As China assumes a more and more dominant role in global science, this mini-review attempts to provide a bird’s eye view on how the bio-digital revolution impacts China’s biosciences and bioindustry. Triggered by top-down political programs and the buildup of an impressive infrastructure in science, information technology, and education, China’s biomedical and MedTech industries prosper. Plant and animal breeding programs transform agriculture and food supply as much as the Internet of things, and synthetic biology offers new opportunities for the manufacturing of specialty chemicals within the Chinese version of a “bioeconomy.” It is already becoming apparent that the new five-year period “145” (2021–2025) will further emphasize emission control, bioenvironmental protection, and more supply of biomass-derived energy. This review identifies key drivers in China’s government, industry, and academia behind these developments and details many access points for deeper studies.

Key points

Biotechnology in China

Biomedical technology

New five-year period

Genome wide association study identifies four loci for early onset schizophrenia

Early onset schizophrenia (EOS, defined as first onset of schizophrenia before age 18) is a rare form of schizophrenia (SCZ). Though genome-wide association studies (GWASs) have identified multiple risk variants for SCZ, most of the cases included in these GWASs were not stratified according to their first age at onset. To date, the genetic architecture of EOS remains largely unknown. To identify the risk variants and to uncover the genetic basis of EOS, we conducted a two-stage GWAS of EOS in populations of Han Chinese ancestry in this study. We first performed a GWAS using 1,256 EOS cases and 2,661 healthy controls (referred as discovery stage). The genetic variants with a P < 1.0 × 10^-04 in discovery stage were replicated in an independent sample (903 EOS cases and 3,900 controls). We identified four genome-wide significant risk loci for EOS in the combined samples (2,159 EOS cases and 6,561 controls), including 1p36.22 (rs1801133, P_meta = 4.03 × 10^-15), 1p31.1 (rs1281571, P_meta = 4.14 × 10^-08), 3p21.31 (rs7626288, P_meta = 1.57 × 10^-09), and 9q33.3 (rs592927, P_meta = 4.01 × 10^-11). Polygenic risk scoring (PRS) analysis revealed substantial genetic overlap between EOS and SCZ. These discoveries shed light on the genetic basis of EOS. Further functional characterization of the identified risk variants and genes will help provide potential targets for therapeutics and diagnostics.

EigenGWAS: An online visualizing and interactive application for detecting genomic signatures of natural selection

Detecting genetic regions under selection in structured populations is of great importance in ecology, evolutionary biology and breeding programmes. We recently proposed EigenGWAS, an unsupervised genomic scanning approach that is similar to F ST but does not require grouping information of the population, for detection of genomic regions under selection. The original EigenGWAS is designed for the random mating population, and here we extend its use to inbred populations. We also show in theory and simulation that eigenvalues, the previous corrector for genetic drift in EigenGWAS, are overcorrected for genetic drift, and the genomic inflation factor is a better option for this adjustment. Applying the updated algorithm, we introduce the new EigenGWAS online platform with highly efficient core implementation. Our online computational tool accepts plink data in a standard binary format that can be easily converted from the original sequencing data, provides the users with graphical results via the R-Shiny user-friendly interface. We applied the proposed method and tool to various data sets, and biologically interpretable results as well as caveats that may lead to an unsatisfactory outcome are given. The EigenGWAS online platform is available at www.eigengwas.com, and can be localized and scaled up via R (recommended) or docker.

KAT8, lysine acetyltransferase 8, is required for adipocyte differentiation in vitro

KAT8 is a lysine acetyltransferase (KAT) that plays a role in a variety of cellular functions ranging from DNA damage repair to apoptosis. The role of KAT8 in adipocyte development and function has not been studied. Notably, a large genome-wide association study identified KAT8 as part of a novel locus that significantly contributed to body mass index and other metabolic phenotypes. Hence, we examined the expression and regulation of KAT8 during adipocyte development. KAT8 mRNA and protein levels were examined over a time course of adipocyte development, and KAT8 was found to be present in both the cytosol and nucleus of 3T3-L1 adipocytes. Although KAT8 expression was not highly regulated by adipogenesis, its expression was required for the adipogenesis of 3T3-L1 cells. Loss of KAT8 expression in preadipocytes inhibited their ability to differentiate as judged by both lipid accumulation and adipocyte marker gene expression. However, if KAT8 was knocked down after clonal expansion, its absence did not inhibit adipocyte differentiation. Also, loss of KAT8 in adipocytes did not impact lipid accumulation or the expression of adiponectin or other fat markers. Although our data demonstrate that KAT8 is required for adipocyte differentiation, further studies are necessary to determine the functions and regulation of KAT8 in adipose tissue.

Evaluation of clinical impact of pharmacogenomics knowledge involved in CPIC guidelines on Chinese pediatric patients

Aim: To evaluate the clinical benefits of implementing pharmacogenomics testing for Chinese pediatric patients. Materials & methods : Based on the drug-gene interactions involved in the Clinical Pharmacogenetics Implementation Consortium guidelines, whole-genome sequencing data from the Chinese Academy of Sciences Precision Medicine Initiative project and the medication data of pediatric patients from a children's hospital, the prevalence of the Chinese population with actionable pharmacogenomic variants was calculated, the prescribing pattern for pediatric patients was analyzed. Results: 37.0% of the drugs involved in the Clinical Pharmacogenetics Implementation Consortium guidelines were used by Chinese pediatric patients, 8.91% inpatients and 0.89% outpatients received at least one pharmacogenomics medication, 1.24% (4803) inpatients and 0.16% (2940) outpatients were estimated to be at high risk of pharmacogenomic-related adverse therapeutic outcomes. Conclusion: Implementing pharmacogenomics testing can improve therapeutic outcomes for many Chinese pediatric patients.

Germline variants of DNA repair genes in early onset mantle cell lymphoma

Although somatic mutations of DNA repair genes are frequent in mantle cell lymphoma (MCL), our understanding of their germline defects is limited. In a Chinese family with maternal Lynch syndrome and paternal B cell non-Hodgkin lymphoma, one sibling developed both Lynch syndrome and MCL. Lynch syndrome is caused by heterozygous mutations in mismatch repair (MMR) genes. To understand the genetic predispositions in the family, we performed exome sequencing and analyses of affected individuals and their tumor samples. A novel germline indel, MLH1 Gly101fsX1, was identified as the cause of Lynch syndrome, and unstable microsatellite loci and mutational signatures as evidence of defective MMR were revealed in the MCL sample. Furthermore, we included additional 15 MCL patients with early onset, and found by exome sequencing that 11 patients carried heterozygous germline variants of 20 DNA repair genes, including MSH2 in MMR. In the MCL with MSH2 Arg359fsX16, unstable microsatellite loci and defective MMR signatures were also found. In addition, five patients also had heterozygous germline variants of genes involved in B cell functions. Thus, our study found germline variants of genes in single-strand break repair, double-strand break repair, and Fanconi anemia pathway in early onset MCL; and for the first time we identified germline defects of MMR in two MCLs.

Polyunsaturated fatty acid biosynthesis pathway and genetics. implications for interindividual variability in prothrombotic, inflammatory conditions such as COVID-19✰,✰✰,★,★★

COVID-19 symptoms vary from silence to rapid death, the latter mediated by both a cytokine storm and a thrombotic storm. SARS-CoV (2003) induces Cox-2, catalyzing the synthesis, from highly unsaturated fatty acids (HUFA), of eicosanoids and docosanoids that mediate both inflammation and thrombosis. HUFA balance between arachidonic acid (AA) and other HUFA is a likely determinant of net signaling to induce a healthy or runaway physiological response. AA levels are determined by a non-protein coding regulatory polymorphisms that mostly affect the expression of FADS1, located in the FADS gene cluster on chromosome 11. Major and minor haplotypes in Europeans, and a specific functional insertion-deletion (Indel), rs66698963, consistently show major differences in circulating AA (>50%) and in the balance between AA and other HUFA (47-84%) in free living humans; the indel is evolutionarily selective, probably based on diet. The pattern of fatty acid responses is fully consistent with specific genetic modulation of desaturation at the FADS1-mediated 20:3→20:4 step. Well established principles of net tissue HUFA levels indicate that the high linoleic acid and low alpha-linoleic acid in populations drive the net balance of HUFA for any individual. We predict that fast desaturators (insertion allele at rs66698963; major haplotype in Europeans) are predisposed to higher risk and pathological responses to SARS-CoV-2 could be reduced with high dose omega-3 HUFA.

De novo genome assembly of a Han Chinese male and genome-wide detection of structural variants using Oxford Nanopore sequencing

Advances in third-generation sequencing technologies provide an opportunity to investigate the complex organizational structure of the genome and unravel the genetic mechanisms of disease and physiological traits. Here we report the sequencing and de novo assembly of a healthy male northern Han Chinese genome and detection of structural variants using only nanopore sequencing data. We performed de novo assembly after filtering the raw data. Then, we aligned the assembled contigs to the human reference genome, and visualized chromosomes plot, which illustrated the contiguity of the nanopore assembly. Additionally, genomic structural variants were detected using a structure variation detection tool with long-read sequencing data. Median coverage depth was 30-fold and the read N50 was 27,136 bp. 96.51% of reads had at least one alignment to the human reference genome. The final assembled genome was 2.85 GB in size, with an N50 contig size of 5.4 MB. We identified 20,085 structural variants. Third-generation sequencing technologies have many advantages in de novo whole-genome assembly and detection of structural variants. Our results provide reference data for disease research, and can be used as a novel population-specific dataset of structural variants to support the efficient development of personalized precision medicine.

Aristolochic acid mutational signature defines the low-risk subtype in upper tract urothelial carcinoma

Rationale: Dietary exposure to aristolochic acids and similar compounds (collectively, AA) is a significant risk factor for nephropathy and subsequent upper tract urothelial carcinoma (UTUC). East Asian populations, who have a high prevalence of UTUC, have an unusual genome-wide AA-induced mutational pattern (COSMIC signature 22). Integrating mutational signature analysis with clinicopathological information may demonstrate great potential for risk ranking this UTUC subtype.

Methods: We performed whole-genome sequencing (WGS) on 90 UTUC Chinese patients to extract mutational signatures. Genome sequencing data for urinary cell-free DNA from 26 UTUC patients were utilized to noninvasively identify the mutational signatures. Genome sequencing for primary tumors on 8 out of 26 patients was also performed. Metastasis-free survival (MFS) and cancer-specific survival (CSS) were measured using Kaplan-Meier methods.

Results: Data analysis showed that a substantial proportion of patients harbored the AA mutational signature and were associated with AA-containing herbal drug intake, female gender, poor renal function, and multifocality. Field cancerization was found to partially contribute to multifocality. Nevertheless, AA Sig subtype UTUC patients exhibited favorable outcomes of CSS and MFS compared to the No-AA Sig subtype. Additionally, AA Sig subtype patients showed a higher tumor mutation burden, higher numbers of predicted neoantigens, and infiltrating lymphocytes, suggesting the potential for immunotherapy. We also confirmed the AA signature in AA-treated human renal tubular HK-2 cells. Notably, the AA subtype could be ascertained using a clinically applicable sequencing strategy (low coverage) in both primary tumors and urinary cell-free DNA as a basis for therapy selection.

Conclusion: The AA mutational signature as a screening tool defines low-risk UTUC with therapeutic relevance. The AA mutational signature, as a molecular prognostic marker using either ureteroscopy and/or urinary cell-free DNA, is especially useful for diagnostic uncertainty when kidney-sparing treatment and/or immune checkpoint inhibitor therapy were considered.

RETRACTED: CGVD: a genomic variation database for Chinese populations

Precision medicine calls upon deeper coverage of population-based sequencing and thorough gene-content and phenotype-based analysis, which lead to a population-associated genomic variation map or database. The Chinese Genomic Variation Database (CGVD; https://bigd.big.ac.cn/cgvd/) is such a database that has combined 48.30 million (M) SNVs and 5.77 M small indels, identified from 991 Chinese individuals of the Chinese Academy of Sciences Precision Medicine Initiative Project (CASPMI) and 301 Chinese individuals of the 1000 Genomes Project (1KGP). The CASPMI project includes whole-genome sequencing data (WGS, 25–30×) from ∼1000 healthy individuals of the CASPMI cohort. To facilitate the usage of such variations for pharmacogenomics studies, star-allele frequencies of the drug-related genes in the CASPMI and 1KGP populations are calculated and provided in CGVD. As one of the important database resources in BIG Data Center, CGVD will continue to collect more genomic variations and to curate structural and functional annotations to support population-based healthcare projects and studies in China and worldwide.

Database Resources of the National Genomics Data Center in 2020

The National Genomics Data Center (NGDC) provides a suite of database resources to support worldwide research activities in both academia and industry. With the rapid advancements in higher-throughput and lower-cost sequencing technologies and accordingly the huge volume of multi-omics data generated at exponential scales and rates, NGDC is continually expanding, updating and enriching its core database resources through big data integration and value-added curation. In the past year, efforts for update have been mainly devoted to BioProject, BioSample, GSA, GWH, GVM, NONCODE, LncBook, EWAS Atlas and IC4R. Newly released resources include three human genome databases (PGG.SNV, PGG.Han and CGVD), eLMSG, EWAS Data Hub, GWAS Atlas, iSheep and PADS Arsenal. In addition, four web services, namely, eGPS Cloud, BIG Search, BIG Submission and BIG SSO, have been significantly improved and enhanced. All of these resources along with their services are publicly accessible at https://bigd.big.ac.cn.