Genomic insights into the formation of human populations in East Asia.

The deep population history of East Asia remains poorly understood owing to a lack of ancient DNA data and sparse sampling of present-day people1,2. Here we report genome-wide data from 166 East Asian individuals dating to between 6000 BC and AD 1000 and 46 present-day groups. Hunter-gatherers from Japan, the Amur River Basin, and people of Neolithic and Iron Age Taiwan and the Tibetan Plateau are linked by a deeply splitting lineage that probably reflects a coastal migration during the Late Pleistocene epoch. We also follow expansions during the subsequent Holocene epoch from four regions. First, hunter-gatherers from Mongolia and the Amur River Basin have ancestry shared by individuals who speak Mongolic and Tungusic languages, but do not carry ancestry characteristic of farmers from the West Liao River region (around 3000 BC), which contradicts theories that the expansion of these farmers spread the Mongolic and Tungusic proto-languages. Second, farmers from the Yellow River Basin (around 3000 BC) probably spread Sino-Tibetan languages, as their ancestry dispersed both to Tibet-where it forms approximately 84% of the gene pool in some groups-and to the Central Plain, where it has contributed around 59-84% to modern Han Chinese groups. Third, people from Taiwan from around 1300 BC to AD 800 derived approximately 75% of their ancestry from a lineage that is widespread in modern individuals who speak Austronesian, Tai-Kadai and Austroasiatic languages, and that we hypothesize derives from farmers of the Yangtze River Valley. Ancient people from Taiwan also derived about 25% of their ancestry from a northern lineage that is related to, but different from, farmers of the Yellow River Basin, which suggests an additional north-to-south expansion. Fourth, ancestry from Yamnaya Steppe pastoralists arrived in western Mongolia after around 3000 BC but was displaced by previously established lineages even while it persisted in western China, as would be expected if this ancestry was associated with the spread of proto-Tocharian Indo-European languages. Two later gene flows affected western Mongolia: migrants after around 2000 BC with Yamnaya and European farmer ancestry, and episodic influences of later groups with ancestry from Turan.

Enhancer variants on chromosome 2p14 regulating SPRED2 and ACTR2 act as a signal amplifier to protect against rheumatoid arthritis.

Genome-wide association studies (GWASs) have repeatedly reported multiple non-coding single-nucleotide polymorphisms (SNPs) at 2p14 associated with rheumatoid arthritis (RA), but their functional roles in the pathological mechanisms of RA remain to be explored. In this study, we integrated a series of bioinformatics and functional experiments and identified three intronic RA SNPs (rs1876518, rs268131, and rs2576923) within active enhancers that can regulate the expression of SPRED2 directly. At the same time, SPRED2 and ACTR2 influence each other as a positive feedback signal amplifier to strengthen the protective role in RA by inhibiting the migration and invasion of rheumatoid fibroblast-like synoviocytes (FLSs). In particular, the transcription factor CEBPB preferentially binds to the rs1876518-T allele to increase the expression of SPRED2 in FLSs. Our findings decipher the molecular mechanisms behind the GWAS signals at 2p14 for RA and emphasize SPRED2 as a potential candidate gene for RA, providing a potential target and direction for precise treatment of RA.

High-throughput functional dissection of noncoding SNPs with biased allelic enhancer activity for insulin resistance-relevant phenotypes.

Most of the single-nucleotide polymorphisms (SNPs) associated with insulin resistance (IR)-relevant phenotypes by genome-wide association studies (GWASs) are located in noncoding regions, complicating their functional interpretation. Here, we utilized an adapted STARR-seq to evaluate the regulatory activities of 5,987 noncoding SNPs associated with IR-relevant phenotypes. We identified 876 SNPs with biased allelic enhancer activity effects (baaSNPs) across 133 loci in three IR-relevant cell lines (HepG2, preadipocyte, and A673), which showed pervasive cell specificity and significant enrichment for cell-specific open chromatin regions or enhancer-indicative markers (H3K4me1, H3K27ac). Further functional characterization suggested several transcription factors (TFs) with preferential allelic binding to baaSNPs. We also incorporated multi-omics data to prioritize 102 candidate regulatory target genes for baaSNPs and revealed prevalent long-range regulatory effects and cell-specific IR-relevant biological functional enrichment on them. Specifically, we experimentally verified the distal regulatory mechanism at IRS1 locus, in which rs952227-A reinforces IRS1 expression by long-range chromatin interaction and preferential binding to the transcription factor HOXC6 to augment the enhancer activity. Finally, based on our STARR-seq screening data, we predicted the enhancer activity of 227,343 noncoding SNPs associated with IR-relevant phenotypes (fasting insulin adjusted for BMI, HDL cholesterol, and triglycerides) from the largest available GWAS summary statistics. We further provided an open resource (http://www.bigc.online/fnSNP-IR) for better understanding genetic regulatory mechanisms of IR-relevant phenotypes.

Transcriptome-wide association study identifies multiple genes and pathways associated with thyroid function.

Thyroid dysfunction is a common endocrine disease measured by thyroid-stimulating hormone (TSH) level. Although >70 genetic loci associated with TSH have been reported through genome-wide association studies (GWASs), the variants can only explain a small fraction of the thyroid function heritability. To identify novel candidate genes for thyroid function, we conducted the first large-scale transcriptome-wide association study (TWAS) for thyroid function using GWAS-summary data for TSH levels in up to 119 715 individuals combined with precomputed gene expression weights of six panels from four tissue types. The candidate genes identified by TWAS were further validated by TWAS replication and gene expression profiles. We identified 74 conditionally independent genes significantly associated with thyroid function, such as PDE8B (P = 1.67 × 10-282), PDE10A (P = 7.61 × 10-119), NR3C2 (P = 1.50 × 10-92) and CAPZB (P = 3.13 × 10-79). After TWAS replication using UKBB datasets, 26 genes were replicated for significant associations with thyroid-relevant diseases/traits. Among them, 16 genes were causal for their associations to thyroid-relevant diseases/traits and further validated in differential expression analyses, including two novel genes (MFSD6 and RBM47) that did not implicate in previous GWASs. Enrichment analyses detected several pathways associated with thyroid function, such as the cAMP signaling pathway (P = 7.27 × 10-4), hemostasis (P = 3.74 × 10-4), and platelet activation, signaling and aggregation (P = 9.98 × 10-4). Our study identified multiple candidate genes and pathways associated with thyroid function, providing novel clues for revealing the genetic mechanisms of thyroid function and disease.

ExsgRNA: reduce off-target efficiency by on-target mismatched sgRNA.

Clustered regularly interspaced short palindromic repeats (CRISPR)/Cas9 gene editing technology has been widely used to facilitate efficient genome editing. Current popular sgRNA design tools only consider the sgRNA perfectly matched to the target site and provide the results without any on-target mismatch. We suppose taking on-target gRNA-DNA mismatches into consideration might provide better sgRNA with similar binding activity and reduced off-target sites. Here, we trained a seq2seq-attention model with feedback-loop architecture, to automatically generate sgRNAs with on-target mismatches. Dual-luciferase reporter experiment showed that multiple sgRNAs with three mismatches could achieve the 80% of the relative activity of the perfect matched sgRNA. Meanwhile, it could reduce the number of off-target sites using sgRNAs with on-target mismatches. Finally, we provided a freely accessible web server sgRNA design tool named ExsgRNA. Users could submit their target sequence to this server and get optimal sgRNAs with less off-targets and similar on-target activity compared with the perfect-matched sgRNA.

Two-sample Mendelian randomization analysis reveals causal relationships between blood lipids and venous thromboembolism.

Venous Thromboembolism (VTE) is a complex disease that can be classified into two subtypes: Deep Vein Thrombosis (DVT) and Pulmonary Embolism (PE). Previous observational studies have shown associations between lipids and VTE, but causality remains unclear. Hence, by utilizing 241 lipid-related traits as exposures and data from the FinnGen consortium on VTE, DVT, and PE as outcomes, we conducted two-sample Mendelian randomization (MR) analysis to investigate causal relationships between lipids and VTE, DVT and PE. The MR results identified that fatty acid (FA) unsaturation traits (Ratio of bis-allylic bonds to double bonds in lipids, and Ratio of bis-allylic bonds to total fatty acids in lipids) were associated with VTE (OR [95% CI]: 1.21 [1.15-1.27]; 1.21 [1.13-1.30]), DVT (OR [95%CI]: 1.24 [1.16-1.33]; 1.26 [1.16-1.36]) and PE (OR [95%CI]: 1.18 [1.08-1.29]; 1.18 [1.09-1.27]). Phosphatidylcholines exhibit potential causal effects on VTE and PE. Phosphatidylcholine acyl-alkyl C40:4 (PC ae C40:4) was negatively associated with VTE (OR [95% CI]: 0.79 [0.73-0.86]), while phosphatidylcholine diacyl C42:6 (PC aa C42:6) and phosphatidylcholine acyl-alkyl C36:4 (PC ae C36:4) were positively associated with PE (OR [95%CI]: 1.44 [1.20-1.72]; 1.22 [1.10-1.35]). Additionally, we found that medium LDL had a protective effect on VTE. Our study indicates that higher FA unsaturation may increase the risk of VTE, DVT, and PE. Different types of phosphatidylcholine have either promotive or inhibitory effects on VTE and PE, contributing to a better understanding of the risk factors for VTE.

Multi-tissue transcriptome-wide association study reveals susceptibility genes and drug targets for insulin resistance-relevant phenotypes.

AIM: Genome-wide association studies (GWAS) have identified multiple susceptibility loci associated with insulin resistance (IR)-relevant phenotypes. However, the genes responsible for these associations remain largely unknown. We aim to identify susceptibility genes for IR-relevant phenotypes via a transcriptome-wide association study. MATERIALS AND METHODS: We conducted a large-scale multi-tissue transcriptome-wide association study for IR (Insulin Sensitivity Index, homeostasis model assessment-IR, fasting insulin) and lipid-relevant traits (high-density lipoprotein cholesterol, triglycerides, low-density lipoprotein cholesterol and total cholesterol) using the largest GWAS summary statistics and precomputed gene expression weights of 49 human tissues. Conditional and joint analyses were implemented to identify significantly independent genes. Furthermore, we estimated the causal effects of independent genes by Mendelian randomization causal inference analysis. RESULTS: We identified 1190 susceptibility genes causally associated with IR-relevant phenotypes, including 58 genes that were not implicated in the original GWAS. Among them, 11 genes were further supported in differential expression analyses or a gene knockout mice database, such as KRIT1 showed both significantly differential expression and IR-related phenotypic effects in knockout mice. Meanwhile, seven proteins encoded by susceptibility genes were targeted by clinically approved drugs, and three of these genes (H6PD, CACNB2 and DRD2) have been served as drug targets for IR-related diseases/traits. Moreover, drug repurposing analysis identified four compounds with profiles opposing the expression of genes associated with IR risk. CONCLUSIONS: Our study provided new insights into IR aetiology and avenues for therapeutic development.

Protoilludane-Type and Related Nor-Sesquiterpenes from Phellinus hartigii and Their Anti-Hypertrophic Activities in Rat Cardiomyocytes.

Three nor-sesquiterpenes, phellinharts A-C (1-3), isolated from Phellinus hartigii, exhibited unprecedented protoilludane and cerapicane-type structures. The structures of compounds 1-3 were elucidated via spectroscopic analysis, quantum chemical calculations, and X-ray diffraction. Potential biogenic pathways involving demethylation, ring cleavage, and rearrangement were proposed. Compounds 1-3 displayed potent anti-hypertrophic activities with low cytotoxicity (CC50 > 50 µM) in rat cardiomyocytes, underscoring their therapeutic potential.

Scutellarin inhibits oleic acid induced vascular smooth muscle foam cell formation via activating autophagy and inhibiting NLRP3 inflammasome activation.

Abnormalities in vascular smooth muscle cells (VSMCs) are pivotal in the pathogenesis of cardiovascular pathologies such as atherosclerosis and hypertension. Scutellarin (Scu), a flavonoid derived from marigold flowers, exhibits a spectrum of biological activities including anti-inflammatory, antioxidant, antitumor, immunomodulatory and antimicrobial effects. Notably, Scu has demonstrated the capacity to mitigate vascular endothelial damage and prevent atherosclerosis via its antioxidative properties. Nevertheless, the influence of Scu on the formation of VSMC-derived foam cells remains underexplored. In this study, Scu was evidenced to efficaciously attenuate oleic acid (OA)-induced lipid accumulation and the upregulation of adipose differentiation-associated protein Plin2 in a dose- and time-responsive manner. We elucidated that Scu effectively diminishes OA-provoked VSMC foam cell formation. Further, it was established that Scu pretreatment augments the protein expression of LC3B-II and the mRNA levels of Map1lc3b and Becn1, concurrently diminishing the protein levels of the NLRP3 inflammasome compared to the OA group. Activation of autophagy through rapamycin attenuated NLRP3 inflammasome protein expression, intracellular lipid droplet content and Plin2 mRNA levels. Scu also counteracted the OA-induced decrement of LC3B-II levels in the presence of bafilomycin-a1, facilitating the genesis of autophagosomes and autolysosomes. Complementarily, in vivo experiments revealed that Scu administration substantially reduced arterial wall thickness, vessel wall cross-sectional area, wall-to-lumen ratio and serum total cholesterol levels in comparison to the high-fat diet model group. Collectively, our findings suggest that Scu attenuates OA-induced VSMC foam cell formation through the induction of autophagy and the suppression of NLRP3 inflammasome activation.

Assessment of bidirectional relationships between brain imaging-derived phenotypes and stroke: a Mendelian randomization study.

BACKGROUND: Stroke is a major cause of mortality and long-term disability worldwide. Whether the associations between brain imaging-derived phenotypes (IDPs) and stroke are causal is uncertain. METHODS: We performed two-sample bidirectional Mendelian randomization (MR) analyses to explore the causal associations between IDPs and stroke. Summary data of 587 brain IDPs (up to 33,224 individuals) from the UK Biobank and five stroke types (sample size range from 301,663 to 446,696, case number range from 5,386 to 40,585) from the MEGASTROKE consortium were used. RESULTS: Forward MR indicated 14 IDPs belong to projection fibers or association fibers were associated with stroke. For example, higher genetically determined mean diffusivity (MD) in the right external capsule was causally associated with an increased risk of small vessel stroke (IVW OR = 2.76, 95% CI 2.07 to 3.68, P = 5.87 × 10-12). Reverse MR indicated that genetically determined higher risk of any ischemic stroke was associated with increased isotropic or free water volume fraction (ISOVF) in body of corpus callosum (IVW ß = 0.23, 95% CI 0.14 to 0.33, P = 3.22 × 10-7). This IDP is a commissural fiber and it is not included in the IDPs identified by forward MR. CONCLUSIONS: We identified 14 IDPs with statistically significant evidence of causal effects on stroke or stroke subtypes. We also identified potential causal effects of stroke on one IDP of commissural fiber. These findings might guide further work toward identifying preventative strategies at the brain imaging levels.

LDBlockShow: a fast and convenient tool for visualizing linkage disequilibrium and haplotype blocks based on variant call format files.

The triangular correlation heatmap aiming to visualize the linkage disequilibrium (LD) pattern and haplotype block structure of SNPs is ubiquitous component of population-based genetic studies. However, current tools suffered from the problem of time and memory consuming. Here, we developed LDBlockShow, an open source software, for visualizing LD and haplotype blocks from variant call format files. It is time and memory saving. In a test dataset with 100 SNPs from 60 000 subjects, it was at least 10.60 times faster and used only 0.03-13.33% of physical memory as compared with other tools. In addition, it could generate figures that simultaneously display additional statistical context (e.g. association P-values) and genomic region annotations. It can also compress the SVG files with a large number of SNPs and support subgroup analysis. This fast and convenient tool will facilitate the visualization of LD and haplotype blocks for geneticists.

Straight and twisted open nodal-line phonon states in the CaI2 family of materials.

Based on first-principles calculations and symmetry analysis, we propose that trigonal CaI2 with the space group P3̄m1 possesses straight and twisted open nodal-line phonon states with linear dispersion. The symmetry analysis indicates that joint symmetry PT and rotational symmetry C3z protect the straight nodal lines along Γ-A and K-H while PT and mirror symmetry M010 (M110) maintain the twisted nodal lines that traverse Γ-M (Γ-K) and A-L (A-H). The calculated π Berry phase suggests that all the nodal lines are nontrivial and the corresponding drumhead-like surface states are clearly visible in the observation window, which is less than 6 THz, suggesting a significant chance for them to be measured using meV-resolution inelastic X-ray scattering. The distribution of the nodal lines in the Brillouin zone is also confirmed by the phononic tight-binding model. Furthermore, the isostructural compounds MgBr2 and MgI2 show similar phonon spectra and topological nontrivial surface states. This work provides promising candidates for investigating straight and twisted open nodal-line phonon states in a single material, which will facilitate future experimental observation.

Application of Hi-C technology in three-dimensional genomics research and disease pathogenesis analysis.

3D genomics aims to investigate the spatial structure of chromatin in the nucleus on the basis of genomic sequences, gene structures and relevant regulatory elements. The spatial organization of chromosomes is fundamental for gene expression regulation. Recent advances of high-throughput chromosome conformation capture (Hi-C) technology and its derivatives, has enabled capture of chromatin architecture with high resolution. In this review, we summarize the development and applications of various technologies of 3D genomes in disease research, particularly in the elucidation of pathogenic mechanisms in cancers and other systemic disorders.

Causal Associations Between Blood Lipids and COVID-19 Risk: A Two-Sample Mendelian Randomization Study.

[Figure: see text].

An ideal Weyl nodal ring with a large drumhead surface state in the orthorhombic compound TiS2.

Topological metals or semimetals have attracted great research attention and interest in condensed matter physics and chemistry due to their exotic properties. Different from the conventional topological insulators, topological metals or semimetals are characterized by distinct topological surface states, such as a Fermi arc or a drumhead surface state, which are often used in experiments to verify the corresponding topological properties. However, the current study in this field is strongly limited in the experimental characterization because of the extreme lack of perfect material candidates with a clean band structure and clear surface states. In this work, based on theoretical calculations, we propose a new topological semimetal TiS2, which has an orthorhombic structure and exhibits excellent stability. Calculated electronic band structures reveal that there is a single Weyl nodal ring in the ky = 0 plane. A detailed symmetry analysis is provided and the corresponding surface state is calculated, which exhibits both a large energy variation of 1.5 eV and wide space distribution without and with the spin orbit coupling effect. Besides, the surface states are well separated from the bulk state. These ideal features together make TiS2 a promising nodal line semimetal for experimental investigation. In combination with the other two isostructural compounds TiSe2 and TiTe2 with similar properties, their further experimental synthesis and characterization can be highly expected and the corresponding study for the topological nodal line state can thus be greatly facilitated.

A two-dimensional PtS2/BN heterostructure as an S-scheme photocatalyst with enhanced activity for overall water splitting.

Photocatalytic hydrogen production from water is a sustainable solution to the environmental pollution and energy crises. Encouraged by the successful synthesis of PtS2 and BN nanosheets and their suitable band edges, we have designed a PtS2/BN bilayer heterojunction and investigated its electronic and optical properties for the first time based on hybrid DFT calculations. In this system, the built-in electric field and band edge bending can retain useful electrons on the conduction band of BN and holes on the valence band of PtS2, which endow this system with a stronger redox ability. Meanwhile, this electric field can efficiently separate photoinduced electron-hole pairs and improve the photocatalytic efficiency. Compared with BN and PtS2 single layers, the PtS2/BN heterojunction with its smaller bandgap can make better use of visible and infrared light. Additionally, we have studied the effect of applied strain on the electronic and optical properties. This work aims to provide a method for constructing high-efficiency BN-based photocatalysts and illuminating the electron migration mechanism in step-scheme (S-scheme) heterostructures. We have found that the PtS2/BN bilayer heterojunction is a promising S-scheme photocatalyst for overall water decomposition.

An Osteoporosis Risk SNP at 1p36.12 Acts as an Allele-Specific Enhancer to Modulate LINC00339 Expression via Long-Range Loop Formation.

Genome-wide association studies (GWASs) have reproducibly associated variants within intergenic regions of 1p36.12 locus with osteoporosis, but the functional roles underlying these noncoding variants are unknown. Through an integrative functional genomic and epigenomic analyses, we prioritized rs6426749 as a potential causal SNP for osteoporosis at 1p36.12. Dual-luciferase assay and CRISPR/Cas9 experiments demonstrate that rs6426749 acts as a distal allele-specific enhancer regulating expression of a lncRNA (LINC00339) (∼360 kb) via long-range chromatin loop formation and that this loop is mediated by CTCF occupied near rs6426749 and LINC00339 promoter region. Specifically, rs6426749-G allele can bind transcription factor TFAP2A, which efficiently elevates the enhancer activity and increases LINC00339 expression. Downregulation of LINC00339 significantly increases the expression of CDC42 in osteoblast cells, which is a pivotal regulator involved in bone metabolism. Our study provides mechanistic insight into how a noncoding SNP affects osteoporosis by long-range interaction, a finding that could indicate promising therapeutic targets for osteoporosis.

Integrating regulatory features data for prediction of functional disease-associated SNPs.

Genome-wide association studies (GWASs) are an effective strategy to identify susceptibility loci for human complex diseases. However, missing heritability is still a big problem. Most GWASs single-nucleotide polymorphisms (SNPs) are located in noncoding regions, which has been considered to be the unexplored territory of the genome. Recently, data from the Encyclopedia of DNA Elements (ENCODE) and Roadmap Epigenomics projects have shown that many GWASs SNPs in the noncoding regions fall within regulatory elements. In this study, we developed a pipeline named functional disease-associated SNPs prediction (FDSP), to identify novel susceptibility loci for complex diseases based on the interpretation of the functional features for known disease-associated variants with machine learning. We applied our pipeline to predict novel susceptibility SNPs for type 2 diabetes (T2D) and hypertension. The predicted SNPs could explain heritability beyond that explained by GWAS-associated SNPs. Functional annotation by expression quantitative trait loci analyses showed that the target genes of the predicted SNPs were significantly enriched in T2D or hypertension-related pathways in multiple tissues. Our results suggest that combining GWASs and regulatory features data could identify additional functional susceptibility SNPs for complex diseases. We hope FDSP could help to identify novel susceptibility loci for complex diseases and solve the missing heritability problem.

Modeling circRNA expression pattern with integrated sequence and epigenetic features demonstrates the potential involvement of H3K79me2 in circRNA expression.

MOTIVATION: CircRNAs are an abundant class of non-coding RNAs with widespread, cell-/tissue-specific patterns. Previous work suggested that epigenetic features might be related to circRNA expression. However, the contribution of epigenetic changes to circRNA expression has not been investigated systematically. Here, we built a machine learning framework named CIRCScan, to predict circRNA expression in various cell lines based on the sequence and epigenetic features. RESULTS: The predicted accuracy of the expression status models was high with area under the curve of receiver operating characteristic (ROC) values of 0.89-0.92 and the false-positive rates of 0.17-0.25. Predicted expressed circRNAs were further validated by RNA-seq data. The performance of expression-level prediction models was also good with normalized root-mean-square errors of 0.28-0.30 and Pearson's correlation coefficient r over 0.4 in all cell lines, along with Spearman's correlation coefficient ρ of 0.33-0.46. Noteworthy, H3K79me2 was highly ranked in modeling both circRNA expression status and levels across different cells. Further analysis in additional nine cell lines demonstrated a significant enrichment of H3K79me2 in circRNA flanking intron regions, supporting the potential involvement of H3K79me2 in circRNA expression regulation. AVAILABILITY AND IMPLEMENTATION: The CIRCScan assembler is freely available online for academic use at https://github.com/johnlcd/CIRCScan. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Organ- and Age-Specific Differences of Dioscorea polystachya Compounds Measured by UPLC-QTOF/MS.

Dioscorea polystachya, named Chinese yam, is widely cultivated as a functional food and natural medicine in China. There is currently little information about the chemical characteristics of Dioscorea polystachya in different organs (tuber cortex and tuber flesh) and at various ages. In this study, an ultra-high-performance liquid chromatography coupled to quadrupole time-of-flight mass spectrometry (UPLC-QTOF/MS) was used to profile chemical compounds in Dioscorea polystachya. As a result, thirty-eight compounds were detected in yam tuber cortex and tuber flesh. More compounds were detected in yam tuber cortex than in tuber flesh. Compounds such as dehydroepiandrosterone, allantoin and flavonoids were selected as biomarker candidates. Dehydroepiandrosterone was found more abundant in tuber flesh, while allantoin and flavonoids showed higher levels in tuber cortex. Furthermore, the levels of dioscin, malvalic acid and sucrose differed significantly among age groups and were highest in the tubers at 2 years. While the levels of allantoin, adenosine and glutamine increased with the growing years and were highest at 4 years. Thus, 2-year old Dioscorea polystachya tubers could be harvested to prepare dioscin, malvalic acid and sucrose supplements. The 4-year-old Dioscorea polystachya tubers would be the best choice for obtaining a large amount of allantoin and adenosine in industrial production.