The ChinaMAP analytics of deep whole genome sequences in 10,588 individuals

Understanding the genetic epidemiology of hereditary breast cancer in India using whole genome data from 1029 healthy individuals

Breast cancer is the most highly reported cancer in India. Genetic testing could help tackle the increasing cancer burden by enabling carriers obtain early diagnosis through increased surveillance, and help guide treatment. However, accurate interpretation of variant pathogenicity must be established in a population-specific manner to ensure effective use of genetic testing. Here we query IndiGen data obtained from sequencing 1029 Indian individuals, and perform variant classification of all reported BRCA variants using gold-standard ACMG & AMP guidelines to establish disease epidemiology. To address the high number of VUS variants thus obtained, we further utilize the brca-NOVUS ML tool to obtain pathogenicity predictions in a manner close to ACMG guidelines at scale. Through the manual application of ACMG & AMP guidelines, we determined the genetic prevalence to be the following: 1 in 342 carriers of BRCA1, and 1 in 256 carriers of BRCA2 pathogenic/likely pathogenic variants bear a significant lifetime risk of developing breast / ovarian cancer in India. The high population prevalence and unique variant landscape emphasizes a need for population-scale studies of causative variants to enable effective screening. We advise cautious clinical interpretation, given incomplete penetrance and other complex factors that result in cancer disease. To the best of our knowledge, this is the first and most comprehensive population-scale genetic epidemiological study of BRCA-linked breast cancer variants reported from India.

Nuciferine activates intestinal TAS2R46 to attenuate metabolic disorders and hyperlipidemia via hepatic VLDL regulation

BACKGROUND

Dysregulated blood lipid metabolism, a primary driver of hyperlipidemia, is closely associated with excessive very low-density lipoprotein (VLDL) synthesis and secretion. Nuciferine, a bioactive compound isolated from lotus leaves, demonstrates remarkable hypolipidemic efficacy; however, its limited bioavailability challenges existing mechanistic explanations for this pronounced therapeutic effect.

PURPOSE

This research aims to investigate the disease-modifying effect and underlying mechanism of nuciferine against hyperlipidemia from a novel perspective by modulating the synthesis and secretion of VLDL.

METHODS

High-fat diet-induced hyperlipidemic rats were assessed by biochemical assays and histopathological examinations to assess the therapeutic effect of nuciferine. Untargeted metabo-lipidomics were launched to obtain the metabolic and lipid profiles, and explainable machine learning algorithms were innovatively utilized in screening differentially expressed metabolites for pathway analysis. A hyperlipidemic two-layer cell co-culture model was analyzed using quantitative polymerase chain reaction, molecular docking, immunofluorescence, cellular thermal shift assay, western blotting, and flow cytometry to delineate VLDL regulatory mechanisms.

RESULTS

Nuciferine significantly attenuated lipid accumulation and metabolic dysfunction in a hyperlipidemic rat model, with the biosynthesis and metabolism of phenylalanine, tyrosine, and tryptophan as pivotal metabolic pathways. Mechanistically, nuciferine activated intestinal type 2 taste receptor 46 (TAS2R46), promoting Ca²⁺-dependent secretion of glucagon-like peptide-1 (GLP-1). Subsequently, the hepatic GLP-1 receptor was cascaded to upregulate expression of liver X receptor alpha and ATP-binding cassette transporter A1, thereby reducing pathological VLDL overproduction.

CONCLUSION

This investigation establishes nuciferine's therapeutic potential in metabolic disorders and hyperlipidemia by activating intestinal TAS2R46 to regulate hepatic VLDL synthesis and secretion.

Analysis of germline-somatic mutational connections in colorectal cancer reveals differential tumorigenic patterns and a novel predictive marker for germline mutation carriers

Colorectal cancer (CRC) genetic testing of regions beyond clinical guidelines has revealed a substantial number of likely pathogenic germline mutations (GMs). It remains largely undetermined whether and how these GMs, typically located in non-mismatch repair (non-MMR) genes, are associated with the tumorigenesis of CRC. This study aimed to identify CRC-predisposing GMs among 93 cancer susceptibility genes and investigate their potential influences on CRC somatic mutational features. We secondarily aimed to investigate whether somatic ERBB2 amplification contributes to identifying GM carriers. This study incorporated a total of 3,240 Chinese CRC patients and 10,588 control individuals. CRC patients were subjected to paired tumor-normal sequencing with a 1,021-gene panel. A case-control analysis was conducted to profile the GM-associated CRC risk. A comprehensive germline-somatic association analysis was performed among 2,405 patients, with key findings subsequently validated in an independent 835-patient cohort and the TCGA CRC cohort. The case-control results supported CRC-predisposing effects of GMs in certain homologous recombination repair (HRR) and DNA damage checkpoint factor (CPF) genes, such as BRCA1/2, RecQ helicase genes, ATM, and CHEK2. HRR GMs were associated with an increased copy number alteration burden, more TP53 clonal mutations, and a higher probability of carrying somatic ERBB2 amplification. CPF GMs were inferred to have synergistic effects with ARID1A and KDM6A somatic mutations in CRC tumorigenesis. Among patients with onset age ≥55 years, stable microsatellites, and no cancer family history, ERBB2 amplification was significantly predictive of GM carriers. Our findings elucidate different germline tumorigenic patterns not driven by deficient MMR. Somatic ERBB2 amplification in CRC can serve as an indicator for germline genetic testing when traditional risk features are absent.

Rare damaging variants in the sex differences of congenital heart disease: an exome sequencing study

BACKGROUND

Congenital heart disease (CHD) exhibits a marked male predominance in birth prevalence, yet the genetic mechanisms underlying this sex disparity remain poorly understood. This study investigates the contribution of rare damaging variants on autosomes and the X chromosome to sex differences in foetal CHD.

METHODS

Parents of foetuses with CHD were recruited for the study. Rare damaging variants were identified by analysing whole-exome sequencing data from foetus-parental trios, and their contributions to sex differences were estimated through case-control studies. Functional enrichment analysis was conducted to assess functional differences in genetic variants between sexes.

FINDINGS

820 foetal probands with CHD were recruited, including 487 males and 333 females. We identified a significant enrichment of X-linked rare damaging variants, primarily driven by maternally inherited hemizygous variants (MIHVs) in male probands (OR = 1·84, P < 0·05), which accounted for 7·2% of male cases in our cohort. These variants were not found to be enriched in female probands. Additionally, X-linked rare damaging de novo variants (DNVs) were not enriched in either male or female probands (female probands: 1·8% versus female controls: 0·7%, P = 0·6789; no DNVs observed in males). Gene-level variant burden analysis revealed three X-linked CHD candidate genes: DCX, CACNA1F, and MAP3K15. Autosomal rare variants did not show significant differences in variant burdens between sexes. Notably, male probands showed specific functional enrichments in cilia-related pathways for autosomal recessive variants, as well as in chromatin remodelling and nervous system development pathways for autosomal DNVs.

INTERPRETATION

Male and female foetal CHD have significantly different genetic landscapes. The enrichment of X-linked rare damaging MIHVs in males provides a genetic explanation for the higher prevalence of CHD in males. This finding highlights the importance of incorporating sex-stratified approaches in clinical diagnostics and research.

FUNDING

Natural Science Foundation of China and Key research and development project of Hebei Province.

Innovative integration of biometric data and blockchain to enhance ownership and trust with NFTs

Biometric data has the potential to revolutionize health analytics and pharmacology by providing personalized insights into drug efficacy and health trajectories. However, its governance presents significant ethical challenges, particularly around individual ownership and privacy. This study addresses these challenges by proposing a sustainable and ethical framework that integrates biometric data with non-fungible tokens (NFTs). We developed a customized NFT framework with advanced smart contract functionalities that enhance privacy protection and decentralized authentication of biometric data ownership. This approach ensures the secure and ethical management of digital health data while reinforcing individuals' control over their biometric information. By leveraging cryptographic techniques for privacy protection, this framework enhances both the security and efficiency of personal health data management, offering a new perspective on ownership in digital health. Furthermore, a sustainable economic model is proposed to facilitate ethical transactions of tokenized biometric data within the NFT marketplace. The implications of this study extend beyond technology and commerce, offering valuable insights into human behavior in emerging digital economies and contributing to the creation of a more sustainable and equitable digital health ecosystem.

Loss of miRNA-Mediated VEGFA Regulation by SNP-Induced Impairment: A Bioinformatic Analysis in Diabetic Complications

Background/Objectives: MicroRNAs (miRNAs) are molecules involved in biological regulation processes, including type 2 diabetes and its complications development. Single nucleotide polymorphisms (SNPs) can alter miRNA mechanisms, resulting in loss or gain effects. VEGFA is recognized for its role in angiogenesis. However, its overexpression can lead to deleterious effects, such as disorganized and inefficient vasculature. Under hyperglycemic conditions, VEGFA expression seems to increase, which may contribute to the development of microvascular and macrovascular diabetic complications. Several miRNAs are associated with VEGFA regulation and seem to act in the prevention of dysregulated expression. This study aimed to investigate SNPs in miRNA regions related to the loss effect in VEGFA regulation, examining their frequency and potential physiological effects in the development of diabetic complications. Methods: VEGFA-targeting miRNAs were identified using the R package multimiR, with validated and predicted results. Tissue expression analysis and SNP search were data-mined with Python 3 for miRNASNP-v3 SNP raw databases. Allele frequencies were obtained from dbSNP. The miRNA-mRNA interaction comparison was obtained in the miRmap tool through Python 3. MalaCards were used to infer physiological disease association. Results: The variant rs371699284 was selected in hsa-miR-654-3p among 103 potential VEGFA-targeting miRNAs. This selected SNP demonstrated promising results in bioinformatics predictions, tissue-specific expression, and population frequency, highlighting its potential role in miRNA regulation and the resulting loss in VEGFA-silencing efficiency. Conclusions: Our findings suggest that carriers of rs1238947970 may increase susceptibility to diabetic microvascular and macrovascular complications. Furthermore, in vitro and in silico studies are necessary to better understand these processes.

Contrastive sequence signatures between the both sides of a recombination spot reveal an adaptation at PPARD locus from standing variation for pleiotropy since out-of-Africa dispersal

BACKGROUND

Drug metabolism and transporter genes are a specialized class of genes involved in absorption, distribution, metabolism and excretion. They easily present distinct genetic population differentiation and are vulnerable to natural selection.

RESULTS

We initiated a study using a special panel of informative genetic markers in such genes and dissected the genetic structure in representative Chinese and worldwide populations. A distinctive sub-population stratification was discovered in extensive Eurasians and resulted from divergence at the PPARD locus. The contrastive sequence signatures between the both sides of a recombination spot prove a selective sweep on this locus for genetic hitchhiking effect. A genealogy-based framework demonstrates the positive selection acting from standing variation exerted a moderate pressure in Eurasians, and drove the adaptive allele up to a high frequency. The timing and tempo estimations for the genetic adaptation indicate its onset coincided with the early out-of-Africa migration of modern humans and it lasted over a prolonged evolutionary history. A phenome-wide association analysis reveals an extended cis-regulation on the local gene expression and the pleiotropy implicated in a variety of complex traits. The colocalization analyses between the genetic associations from cis-acting gene expression and complex traits signify the most likely selective pressure from physical capacity, energy metabolism, and immune-related involvement, and provide prioritization for the effective genes and casual variants.

CONCLUSIONS

This work has laid a foundation for following efforts to make full sense of the biological mechanisms underlying the genetic adaptation.

Population genomics advances in frontier ethnic minorities in China

China, with its large geographic span, possesses rich genetic diversity across vast frontier regions in addition to the Han Chinese majority. Importantly, demographic events and various natural and cultural environments in Chinese frontier regions have shaped the genomic diversity of ethnic minorities via local adaptations. Thus, insights into the genetic diversity and adaptive evolution of these under-represented ethnic groups are crucial for understanding evolutionary scenarios and biomedical implications in East Asian populations. Here, we focus on ethnic minorities in Chinese frontier regions and review research advances regarding genomic diversity, genetic structure, population history, genetic admixture, and local adaptation. We first provide an overview of the extensive genetic diversity across populations in different Chinese frontier regions. Next, we summarize research progress regarding genetic ancestry, demographic history, the adaptive process, and the archaic identification of multiple ethnic minorities in different Chinese frontier regions. Finally, we discuss the gaps and opportunities in genomic studies of Chinese populations and the need for a more comprehensive understanding of genomic diversity and the evolution of populations of East Asian ancestry in the post-genomic era.

Mutation Screening of ATXN1, ATXN2, and ATXN3 in Amyotrophic Lateral Sclerosis

Emerging evidence suggests potential disease modifying roles of ATXN1, ATXN2, and ATXN3 in amyotrophic lateral sclerosis (ALS). We aimed to provide a comprehensive variants profile of the ATXN1, ATXN2, and ATXN3 genes and examine the association of these variants with the risk and clinical characteristics of ALS. We screened and analyzed the rare variants in a cohort of 2220 ALS patients from Southwest China, using controls from the Genome Aggregation Database (gnomAD) and the China Metabolic Analytics Project (ChinaMAP). The over-representation of rare variants and their association with disease risk in ALS patients were assessed using Fisher's exact test with Bonferroni correction at both allele and gene levels. Kaplan-Meier analysis was employed to explore the relationship between the distribution of variants and survival. A total of 62 eligible rare missense variants were identified, comprising 32 from ATXN1, 21 from ATXN2, and 9 from ATXN3. Allelic association testing revealed a significant enrichment of the ATXN1 (c.2122C > G, p.Leu708Val) variant and the ATXN2 (c.3778C > G, p.Pro1260Ala) variant in ALS. Gene burden analysis indicated that variants in the ATXN1 and ATXN3 genes had a higher burden in ALS. Substantial heterogeneity in survival time was observed among patients carrying different variants within the same gene. However, there were no significant differences in survival between ALS patients grouped by N-terminal or C-terminal distribution. Our results provided a genetic variation profile of ATXN1, ATXN2, and ATXN3 in ALS patients, along with the clinical characteristics of individuals carrying these variations. This information might offer valuable insights for the ongoing ALS disease-modifying treatments.

Pilot work of the 10K Chinese People Genomic Diversity Project along the Silk Road suggests a complex east-west admixture landscape and biological adaptations

Genomic sources from China are underrepresented in the population-specific reference database. We performed whole-genome sequencing or genome-wide genotyping on 1,207 individuals from four linguistically diverse groups (1,081 Sinitic, 56 Mongolic, 40 Turkic, and 30 Tibeto-Burman people) living in North China included in the 10K Chinese People Genomic Diversity Project (10K_CPGDP) to characterize the genetic architecture and adaptative history of ethnic groups in the Silk Road Region of China. We observed a population split between Northwest Chinese minorities (NWCMs) and Han Chinese since the Upper Paleolithic and later Neolithic genetic differentiation within NWCMs. The observed population substructures among ethnically/linguistically diverse NWCMs suggested that differentiated admixture events contributed to the differences in their genomic and phenotypic diversity. We estimated that the Dongxiang, Tibetan, and Yugur people inherited more than 10% of the Western Eurasian ancestry, which is much greater than that of the Salar and Tu people (<7%), while Han neighbors showed less West Eurasian ancestry (∼1%-3%). Male-biased admixture introduced Western Eurasian ancestry in the Dongxiang, Tibetan, and Yugur populations. We found that the eastern-western admixture in NWCMs occurred ∼800-1,100 years ago, coinciding with intensive economic and cultural exchanges during the Tang and Song dynasties. Additionally, we identified the signatures of natural selection associated with cardiovascular system diseases or lipid metabolism and developmental/neurogenetic disorders. Moreover, the EPAS1 gene showed relatively high population branch statistic values in NWCMs. The well-fitted demographical models presented the vast landscape of complex admixture processes of the Silk Road people, and the newly reported functionally important variations suggested the importance of including ethnolinguistically diverse populations in Chinese genetic studies for uncovering the genetic basis of complex traits/diseases.

Tailoring monogenic disease carrier screening panels for Chinese populations: The importance of considering regional differences

Carrier screening for monogenic diseases is becoming increasingly important in preventive medicine, yet selecting appropriate target genes remains a complex task, especially in countries with significant ethnic and geographic diversity such as China. This study aimed to develop a strategy to screen for carrier screening target genes suitable for the Chinese population, considering regional variations in carrier frequencies (CFs). We analyzed a dataset from a large-scale, multicenter carrier screening study, encompassing 33,104 individuals from different regions of China and carrier status for 223 genes. We focused on the CFs of these genes across regions. The study first stratified the population based on participants' self-reported ancestral places and then applied consensus k-means clustering analysis to the CF characteristics of these regions. This approach enabled us to identify distinct regional subpopulations with shared genetic backgrounds. The results showed that the regions clustered into three subpopulations (North, South, and Far South) based on CF characteristics, and 44 genes exhibited significant CF differences across these subpopulations (α = 0.05). Applying an overall CF threshold without considering regional diversity would have excluded 11 regionally prevalent genes from the screening panel. By incorporating regional variations, we accurately identified 58 genes that met the recommended CF criteria (autosomal gene CF > 1/200, X-linked gene CF > 1/40,000) in at least one subpopulation. This study emphasizes the importance of considering regional diversity when designing carrier screening panels for monogenic diseases in China. Our proposed strategy, combining regional stratification and clustering analysis, provides a more precise method for selecting target genes, thereby enhancing the effectiveness and relevance of screening programs across different Chinese populations.

A gain-of-function variant in RICTOR predisposes to human obesity

mTORC1/2 play central roles as signaling hubs of cell growth and metabolism and are therapeutic targets for several diseases. However, the human genetic evidence linking mutations of mTORC1/2 to obesity remains elusive. Using whole-exome sequencing of 1944 cases with severe obesity and 2161 healthy lean controls, we identify a rare RICTOR p.I116V variant enriched in 9 unrelated cases. In Rictor null mouse embryonic fibroblasts, overexpression of the RICTOR p.I116V mutant increases phosphorylation of AKT, a canonical mTORC2 substrate, compared with wild-type RICTOR, indicating a gain-of-function change. Consistent with the human obesity phenotype, the knock-in mice carrying homogenous Rictor p.I116V variants gain more body weight under a high-fat diet. Additionally, the stromal vascular fraction cells derived from inguinal white adipose tissue of knock-in mice display an enhanced capacity for adipocyte differentiation via AKT activity. These findings demonstrate that the rare gain-of-function RICTOR p.I116V mutation activates AKT signaling, promotes adipogenesis, and contributes to obesity in humans.

Paleolithic divergence and multiple Neolithic expansions of ancestral nomadic emperor-related paternal lineages

The reconstruction of demographic history using ancient and modern genomic resources reveals extensive interactions and admixture between ancient nomadic pastoralists and the social organizations of the Chinese Central Plain. However, the extent to which Y-chromosome genetic legacies from nomadic emperor-related ancestral lineages influence the Chinese paternal gene pool remains unclear. Here, we genotype 2717 ethnolinguistically diverse samples belonging to C2a lineages, perform whole-genome sequencing on 997 representative samples, and integrate these data with ancient genomic sequences. We reconstruct the evolutionary histories of Northern Zhou-, Qing emperor-, and pastoralist-related lineages to assess their genetic impact on modern Chinese populations. This reassembled fine-scale Y-chromosome phylogeny identifies deep divergence and five Neolithic expansion events contributing differently to the formation of northern Chinese populations. Phylogeographic modeling indicates that the nomadic empires of the Northern Zhou and Qing dynasties genetically originated from the Mongolian Plateau. Phylogenetic topology and shared haplotype patterns show that three upstream ancestors of Northern Zhou (C2a1a1b1a2a1b-FGC28857), Donghu tribe (C2a1a1b1-F1756), and Qing (C2a1a3a2-F10283) emperor-related lineages expanded during the middle Neolithic, contributing significantly to genetic flow between ancient northeastern Asians and modern East Asians. Notably, this study reveals limited direct contributions of Emperor Wu of Northern Zhou's lineages to modern East Asians.

Identification of novel germline mutations in FUT7 and EXT1 linked with hereditary multiple exostoses

Hereditary multiple exostoses (HME) is an autosomal dominant skeletal disorder primarily linked with mutations in Exostosin-1 (EXT1) and Exostosin-2 (EXT2) genes. However, not all HME cases can be explained by these mutations, and its pathogenic mechanisms are not fully understood. Herein, utilizing whole-exome sequencing and genetic screening with a family trio design, we identify two novel rare mutations co-segregating with HME in a Chinese family, including a nonsense mutation (c.204G>A, p.Trp68*) in EXT1 and a missense mutation (c.893T>G, p.Phe298Cys) in FUT7. Functional assays reveal that the FUT7 mutation affects the cellular localization of FUT7 protein and regulates cell proliferation. Notably, the simultaneous loss of fut7 and ext1 in a zebrafish model results in severe chondrodysplasia, indicating a functional link between FUT7 and EXT1 in chondrocyte regulation. Additionally, we unveil that FUT7 p.Phe298Cys reduces EXT1 expression through IL6/STAT3/SLUG axis at the transcription level and through ubiquitination-related proteasomal degradation at the protein level. Together, our findings not only identify novel germline mutations in FUT7 and EXT1 genes, but also highlight the critical interaction between these genes, suggesting a potential 'second-hit' mechanism over EXT1 mutations in HME pathogenesis. This insight enhances our understanding of the mechanisms underlying HME and opens new avenues for potential therapeutic interventions.

Cell-free DNA from clinical testing as a resource of population genetic analysis

As a noninvasive biomarker, cell-free DNA (cfDNA) has achieved remarkable success in clinical applications. Notably, cfDNA is essentially DNA, and conducting whole-genome sequencing (WGS) can yield a wealth of genetic information. These invaluable data should not be confined to one-time use; instead, they should be leveraged for more comprehensive population genetic analysis, including genetic variation spectrum, population structure and genetic selection, and genome-wide association studies (GWASs), among others. Such research findings can, in turn, facilitate clinical practice, enabling more advanced and accurate disease predictions. This review explores the advantages, challenges, and current research areas of cfDNA in population genetics. We hope that this review can serve as a new chapter in the repurposing of cfDNA sequence data generated from clinical testing in population genetics.

Rare and common genetic variants underlying the risk of Hirschsprung's disease

Hirschsprung's disease (HSCR) is a congenital enteric neuropathic disorder characterized by high heritability (>80%) and polygenic inheritance (>20 genes). The previous genome-wide association studies (GWAS) identified several common variants associated with HSCR and demonstrated increased predictive performance for HSCR risk in Europeans using a genetic risk score, there remains a notable gap in knowledge regarding Chinese populations. We conducted whole exome sequencing in a HSCR case cohort in Chinese. By using the common controls (505 controls from 1KG EAS and 10 588 controls from ChinaMAP), we conducted GWAS for the common variants in the exome and gene-based association for rare variants. We further validated the associated variants and genes in replicated samples and in vitro and vivo experiments. We identified one novel gene PLK5 by GWAS and suggested 45 novel putative genes based the gene-based test. By using genetic variant at RET and PLK5, we constructed a genetic risk score that could identify the individuals with very high genetic risk for HSCR. Compared with patients with zero or one risk allele from the three variants, the risk for HSCR was 36.61 times higher with six alleles. In addition, we delineated a HSCR risk gene landscape that encompasses 57 genes, which explains 88.5% and 54.5% of HSCR in Chinese and European, respectively. In summary, this study improved the understanding of genetic architecture of HSCR and provided a risk prediction approach for HSCR in the Chinese.

Associations between RetNet gene polymorphisms and the efficacy of orthokeratology for myopia control: a retrospective clinical study

BACKGROUND

This study investigated how clinical and genetic factors impact the effectiveness of orthokeratology lenses in myopia.

METHODS

A retrospective clinical study was conducted with a sample of 545 children aged 8-12 years who had myopia and have initially worn orthokeratology lenses for one year. Whole-genome sequencing (WGS) was also performed on 60 participants in two groups, one with rapid axial length (AL) progression of larger than 0.33 mm and the other with slow AL progression of less than 0.09 mm. The RetNet database was used to screen candidate genes that may contribute to the effectiveness of orthokeratology lenses in controlling myopia.

RESULTS

Children with greater baseline AL, greater spherical equivalent (SE) and greater age had better myopia control with orthokeratology lenses. A significant excess of nonsynonymous variants was observed among those with slow myopia progression, and these were prominently enriched in retinal disease-related genes. Subsequently, RIMS2 [odds ratio (OR) = 0.01, P = 0.0097] and LCA5 (OR = 9.27, P = 0.0089) were found to harbor an excess number of nonsynonymous variants in patients with slow progression of high myopia. Two intronic common variants rs36006402 in SLC7A14 and rs2285814 in CLUAP1 were strongly associated with AL growth. The identification of these novel genes associated with the effectiveness of orthokeratology lens therapy in myopic children provides insight into the genetic mechanism of orthokeratology treatment.

CONCLUSION

The effectiveness of orthokeratology lens treatment relates to interindividual variability in the control of AL growth in myopic eyes. The efficacy increased when patients carried more nonsynonymous variants in retinal disease-related gene sets. These data serve as reference for genetic counselling and the management of patients who choose orthokeratology lenses to control myopia.

A scoping review of human genetic resources management policies and databases in high- and middle-low-income countries

BACKGROUND

This review examines global human genetic resources management, focusing on genetic data policies and repositories in high- and middle-low-income countries.

METHODS

A comprehensive search strategy was employed across multiple databases, including official government websites and Google, to gather relevant literature on human genetic resources management policies and genetic resource databases. Documents were screened for relevance, focusing on high-income countries (United States, United Kingdom, Japan) and middle-low-income countries (China, India, Kenya). Data were extracted, coded, and analyzed to identify common themes and differences in genetic resource management practices.

RESULTS

High-income countries benefit from robust legal frameworks and advanced technological infrastructures. The United States enforces the Health Insurance Portability and Accountability Act and the Genetic Information Nondiscrimination Act to protect privacy and facilitate data sharing, while Japan relies on the Act on the Protection of Personal Information and ethical guidelines. Additionally, high-income countries host a variety of genetic databases and biobanks that support scientific research. In contrast, middle-low-income countries like China, India, and Kenya are still developing their frameworks. China has regulations such as the Biosecurity Law and the Regulations on the Management of Human Genetic Resources, but still requires more unified standards. India's policies focus on genetic research and data protection through the Biological Diversity Act, while Kenya seeks to improve data management through the 2019 Data Protection Act.

CONCLUSION

Significant disparities exist in human genetic resources management between high-income and middle-low-income countries. High-income countries have robust systems balancing privacy protection with research facilitation, supported by comprehensive and large-scale databases for scientific research. Middle-low-income countries need to enhance legal frameworks and build population-specific databases. Promoting equitable data sharing and adopting best practices from high-income countries are essential for advancing global scientific discovery and ensuring fair management of genetic resources.

Identification of mutations associated with congenital cataracts in nineteen Chinese families

BACKGROUND

Congenital cataracts (CC) are one of the leading causes of impaired vision or blindness in children, with approximately 8.3-25% being inherited. The aim of this study is to investigate the mutation spectrum and frequency of 9 cataract-associated genes in 19 Chinese families with congenital cataracts.

PURPOSE

To identify the gene variants associated with congenital cataracts.

METHODS

This study included a total of 58 patients from 19 pedigrees with congenital cataracts. All probands were initially screened by whole-exome sequencing(WES), and then validated by co-segregation analysis using Sanger sequencing.

RESULTS

Likely pathogenic variants were detected in 8 families, with a positivity rate of 42.1%. Variants in various genes were identified, including GJA3, CRYGD, CRYBA4, BFSP2, IARS2, CRYAA, CRYBA1, ARL2 and CRYBB3. Importantly, this study identified compound heterozygous variants of IARS2 in one family.

CONCLUSIONS

Our research findings have revealed multiple gene variants associated with cataracts, providing clinical guidance for improved molecular diagnosis of congenital cataracts in the era of precision medicine.

Identification of a novel single nucleotide deletion in the NHS causing Nance-Horan syndrome

BACKGROUND

Nance-Horan syndrome (NHS) is a rare X-linked dominant disorder caused by pathogenic variants in the NHS gene on chromosome Xp22.2-Xp22.13. Clinical manifestations consist of congenital cataracts, along with dysmorphic facial features and dental anomalies and, in certain instances, intellectual disability. This study aimed to identify the genetic cause responsible for NHS in a Chinese family with four individuals primarily presenting with congenital cataracts.

METHODS

Genomic DNA was collected from six family members, including four affected individuals (three females and one male) from a two-generation family. The family history and clinical data were documented. Whole-exome sequencing was performed on the proband, and candidate pathogenic variants were filtered through a series of screening steps and validated by Sanger sequencing. Co-segregation analysis was conducted to confirm the pathogenicity of the identified variant.

RESULTS

Genetic analysis revealed a novel frameshift pathogenic variant in NHS gene (c.1735delA: p.R579Gfs*91) present in all four affected members. All affected members exhibited congenital cataracts, congenital ptosis, strabismus, high myopia as well as dental and facial anomalies, and more severe characteristic features observed in the male patient. These clinical manifestations were consistent with the phenotype of NHS.

CONCLUSION

This study identified a novel NHS pathogenic variant in a Chinese family, expanding the mutational spectrum of NHS. Contrary to previous reports of female carriers exhibiting mild symptoms, we demonstrated severe ocular phenotypes in three affected females. These findings will assist in providing genetic counseling for NHS patients.

Genetic diversity and dietary adaptations of the Central Plains Han Chinese population in East Asia

The Central Plains Han Chinese (CPHC) is the typical agricultural population of East Asia. Investigating the genome of the CPHC is crucial to understanding the genetic structure and adaptation of the modern humans in East Asia. Here, we perform whole genome sequencing of 492 CPHC individuals and obtained 22.65 million SNPs, 4.26 million INDELs and 41,959 SVs. We found the CPHC has a higher level of genetic diversity and the glycolipid metabolic genes show strong selection signals, e.g. LONP2, FADS2, FGF21 and SLC19A2. Ancient DNA analyses suggest that the domestication of crops, which drove the emergence of the candidate mutations. Notably, East Asian-specific SVs, e.g., DEL_21699 (LINC01749) and DEL_38406 (FAM102A) may be associated with the high prevalence of esophageal squamous carcinoma and primary angle-closure glaucoma. Our results provide an important genetic resource and show that dietary adaptations play an important role in phenotypic evolution in East Asian populations.

Genetic mapping of serum metabolome to chronic diseases among Han Chinese

Serum metabolites are potential regulators for chronic diseases. To explore the genetic regulation of metabolites and their roles in chronic diseases, we quantified 2,759 serum metabolites and performed genome-wide association studies (GWASs) among Han Chinese individuals. We identified 184 study-wide significant (p < 1.81 × 10-11) metabolite quantitative trait loci (metaboQTLs), 88.59% (163) of which were novel. Notably, we identified Asian-ancestry-specific metaboQTLs, including the SNP rs2296651 for taurocholic acid and taurochenodesoxycholic acid. Leveraging the GWAS for 37 clinical traits from East Asians, Mendelian randomization analyses identified 906 potential causal relationships between metabolites and clinical traits, including 27 for type 2 diabetes and 38 for coronary artery disease. Integrating genetic regulation of the transcriptome and proteome revealed putative regulators of several metabolites. In summary, we depict a landscape of the genetic architecture of the serum metabolome among Han Chinese and provide insights into the role of serum metabolites in chronic diseases.

The pharmacogenomic landscape in the Chinese: An analytics of pharmacogenetic variants in 206,640 individuals

Pharmacogenomic landscapes and related databases are important for identifying the biomarkers of drug response and toxicity. However, these data are still lacking for the Chinese population. In this study, we constructed a pharmacogenomic landscape and an associated database using whole-genome sequencing data generated by non-invasive prenatal testing in 206,640 Chinese individuals. In total, 1,577,513 variants (including 331,610 novel variants) were identified among 3,538 pharmacogenes related to 2,086 drugs. We found that the variant spectrum in the Chinese population differed among the seven major regions. Regional differences also exist among provinces in China. The average numbers of drug enzyme, transporter, and receptor variants were 258, 557, and 632, respectively. Subsequent correlation analysis indicated that the pharmacogenes affecting multiple drugs had fewer variants. Among the 16 categories of drugs, we found that nervous system, cardiovascular system, and genitourinary system/sex hormone drugs were more likely to be affected by variants of pharmacogenes. Characteristics of the variants in the enzyme, transporter, and receptor subfamilies showed specificity. To explore the clinical utility of these data, a genetic association study was conducted on 1,019 lung cancer patients. Two novel variants, AKT2 chr19:40770621 C>G and SLC19A1 chr21:46934171 A>C, were identified as novel platinum response biomarkers. Finally, a pharmacogenomic database, named the Chinese Pharmacogenomic Knowledge Base (CNPKB: http://www.cnpkb.com.cn/), was constructed to collect all the data. In summary, a pharmacogenomic landscape and database for the Chinese population were constructed in this study, which could support personalized Chinese medicine in the future.

Associations between ABO, FUT2 and chronic pancreatitis: A comprehensive meta-analysis of multiple cohorts and public biobanks

OBJECTIVES

Associations of ABO blood group specifying transferases A/B (ABO) and fucosyltransferase 2 (FUT2) with CP remain inconclusive. We aimed to comprehensively investigate the associations by Chinese sequencing cohorts and external cohorts.

METHODS

First, we analyzed the distributions of ABO blood groups and FUT2 status, along with lead single nucleotide polymorphisms (SNPs) at ABO (rs8176693 C/T) and FUT2 (rs632111 A/G) gene loci in Chinese low-coverage whole-genome sequencing discovery cohort. Subsequently, we investigated the associations of CP with ABO and FUT2 SNPs in Chinese whole-exome sequencing validation cohort and three public biobanks (FinnGen, UK Biobank, and BioBank Japan). Finally, comprehensive meta-analysis was performed by integrating data from two Chinese cohorts, reported cohorts and public biobanks.

RESULTS

Firstly, in Chinese discovery cohort, the distribution of blood types in CP patients showed no significant difference compared to healthy controls, and similar results were observed in subgroup analyses and in meta-analysis with the reported cohorts. Secondly, results indicated no association between rs8176693 or rs632111 and CP in Chinese cohorts and meta-analysis of three biobanks, though ABO SNP was found to be significantly associated with CP in UK Biobank (odds ratio [OR] = 1.27, P = 0.04). Finally, no association was observed between rs8176693 (OR = 1.03, P = 0.29) or rs632111 (OR = 1.04, P = 0.10) and CP in comprehensive meta-analysis.

CONCLUSIONS

No association was found between lead SNPs of ABO or FUT2 and CP in meta-analysis, nor was there an association between ABO blood group or FUT2 secretor status and CP in Chinese cohort. ABO and FUT2 might play limited role in CP development.

Perspectives and opportunities in forensic human, animal, and plant integrative genomics in the Pangenome era

The Human Pangenome Reference Consortium, the Chinese Pangenome Consortium, and other plant and animal pangenome projects have announced the completion of pilot work aimed at constructing high-quality, haplotype-resolved reference graph genomes representative of global ethno-linguistically different populations or different plant and animal species. These graph-based, gapless pangenome references, which are enriched in terms of genomic diversity, completeness, and contiguity, have the potential for enhancing long-read sequencing (LRS)-based genomic research, as well as improving mappability and variant genotyping on traditional short-read sequencing platforms. We comprehensively discuss the advancements in pangenome-based genomic integrative genomic discoveries across forensic-related species (humans, animals, and plants) and summarize their applications in variant identification and forensic genomics, epigenetics, transcriptomics, and microbiome research. Recent developments in multiplexed array sequencing have introduced a highly efficient and programmable technique to overcome the limitations of short forensic marker lengths in LRS platforms. This technique enables the concatenation of short RNA transcripts and DNA fragments into LRS-optimal molecules for sequencing, assembly, and genotyping. The integration of new pangenome reference coordinates and corresponding computational algorithms will benefit forensic integrative genomics by facilitating new marker identification, accurate genotyping, high-resolution panel development, and the updating of statistical algorithms. This review highlights the necessity of integrating LRS-based platforms, pangenome-based study designs, and graph-based pangenome references in short-read mapping and LRS-based innovations to achieve precision forensic science.

YHSeqY3000 panel captures all founding lineages in the Chinese paternal genomic diversity database

BACKGROUND

The advancements in second-/third-generation sequencing technologies, alongside computational innovations, have significantly enhanced our understanding of the genomic structure of Y-chromosomes and their unique phylogenetic characteristics. These researches, despite the challenges posed by the lack of population-scale genomic databases, have the potential to revolutionize our approach to high-resolution, population-specific Y-chromosome panels and databases for anthropological and forensic applications.

OBJECTIVES

This study aimed to develop the highest-resolution Y-targeted sequencing panel, utilizing time-stamped, core phylogenetic informative mutations identified from high-coverage sequences in the YanHuang cohort. This panel is intended to provide a new tool for forensic complex pedigree search and paternal biogeographical ancestry inference, as well as explore the general patterns of the fine-scale paternal evolutionary history of ethnolinguistically diverse Chinese populations.

RESULTS

The sequencing performance of the East Asian-specific Y-chromosomal panel, including 2999-core SNP variants, was found to be robust and reliable. The YHSeqY3000 panel was designed to capture the genetic diversity of Chinese paternal lineages from 3500 years ago, identifying 408 terminal lineages in 2097 individuals across 41 genetically and geographically distinct populations. We identified a fine-scale paternal substructure that was correlating with ancient population migrations and expansions. New evidence was provided for extensive gene flow events between minority ethnic groups and Han Chinese people, based on the integrative Chinese Paternal Genomic Diversity Database.

CONCLUSIONS

This work successfully integrated Y-chromosome-related basic genomic science with forensic and anthropological translational applications, emphasizing the necessity of comprehensively characterizing Y-chromosome genomic diversity from genomically under-representative populations. This is particularly important in the second phase of our population-specific medical or anthropological genomic cohorts, where dense sampling strategies are employed.

GoFCards: an integrated database and analytic platform for gain of function variants in humans

Gain-of-function (GOF) variants, which introduce new or amplify protein functions, are essential for understanding disease mechanisms. Despite advances in genomics and functional research, identifying and analyzing pathogenic GOF variants remains challenging owing to fragmented data and database limitations, underscoring the difficulty in accessing critical genetic information. To address this challenge, we manually reviewed the literature, pinpointing 3089 single-nucleotide variants and 72 insertions and deletions in 579 genes associated with 1299 diseases from 2069 studies, and integrated these with the 3.5 million predicted GOF variants. Our approach is complemented by a proprietary scoring system that prioritizes GOF variants on the basis of the evidence supporting their GOF effects and provides predictive scores for variants that lack existing documentation. We then developed a database named GoFCards for general geneticists and clinicians to easily obtain GOF variants in humans (http://www.genemed.tech/gofcards). This database also contains data from >150 sources and offers comprehensive variant-level and gene-level annotations, with the aim of providing users with convenient access to detailed and relevant genetic information. Furthermore, GoFCards empowers users with limited bioinformatic skills to analyze and annotate genetic data, and prioritize GOF variants. GoFCards offers an efficient platform for interpreting GOF variants and thereby advancing genetic research.

Development and evaluation of a panel of newly screened Y chromosome InDels for inferring paternal ancestry information in Southwest China

Y-InDels (insertions/deletions) are genetic markers which are extremely understudied. It is unknown whether this type of markers can be utilized for genetic ancestry inference. We have developed an innovative Y chromosome ancestry inference system tailored for forensic applications. This panel amplifies 21 Y chromosome loci, encompassing Y-InDels and Y-SNPs (Single Nucleotide Polymorphism), utilizing the capillary electrophoresis (CE) platform. The system performed well at DNA concentrations greater than 0.125 ng/ul and produced accurate results at a 1:100 mixing ratio of male and female DNA. The Cumulative probability of matching (CPM) was between 0.95 and 0.97 in the experimental population. The system's efficacy in inferring ancestral origins was demonstrated through intercontinental population discrimination, revealing high discrimination power between African and East Asian populations. Population genetic analyses conducted on Han, Qiang and Hui populations in Southwest China, where the smallest FST value was 0.0002 between Han Chinese in Beijing (from 1000 Genomes Project) and Qiang Chinese from Sichuan (CQSC). Phylogenetic tree construction further illuminated distinct haplotypes among populations, with ethnically unique haplotypes observed in 34.6% of Hui and 7.1% of Qiang populations. K-fold cross-validation show the system's inference abilities at the intercontinental level. In addition, our investigations identified potential associations between the Y-InDel locus Y: 15,385,547 (GRCh37) and haplogroup R1a1a1b2a2- Z2124, as well as locus Y: 13,990,180 (GRCh37) and haplogroup F-M89. In conclusion, we have established a Y-chromosome inference system tailored for grassroots-level application, underscoring the value of incorporating Y-InDel markers in forensic analyses.

SEAD reference panel with 22,134 haplotypes boosts rare variant imputation and genome-wide association analysis in Asian populations

Limited whole genome sequencing (WGS) studies in Asian populations result in a lack of representative reference panels, thus hindering the discovery of ancestry-specific variants. Here, we present the South and East Asian reference Database (SEAD) panel ( https://imputationserver.westlake.edu.cn/ ), which integrates WGS data for 11,067 individuals from various sources across 17 Asian countries. The SEAD panel, comprising 22,134 haplotypes and 88,294,957 variants, demonstrates improved imputation accuracy for South Asian populations compared to 1000 Genomes Project, TOPMed, and ChinaMAP panels, with a higher proportion of well-imputed rare variants. For East Asian populations, SEAD shows concordance comparable to ChinaMAP, but outperforming TOPMed. Additionally, we apply the SEAD panel to conduct a genome-wide association study for total hip (Hip) and femoral neck (FN) bone mineral density (BMD) traits in 5369 genotyped Chinese samples. The single-variant test suggests that rare variants near SNTG1 are associated with Hip BMD (rs60103302, MAF = 0.0092, P = 1.67 × 10-7), and variant-set analysis further supports the association (Pslide_window = 9.08 × 10-9, Pgene_centric = 5.27 × 10-8). This association was not reported previously and can only be detected by using Asian reference panels. Preliminary in vitro experiments for one of the rare variants identified provide evidence that it upregulates SNTG1 expression, which could in turn inhibit the proliferation and differentiation of preosteoblasts.

Genomic insights into Neolithic founding paternal lineages around the Qinghai-Xizang Plateau using integrated YanHuang resource

Indigenous populations of the Qinghai-Xizang Plateau exhibit unique high-altitude adaptations, especially within Tibeto-Burman (TB) groups. However, the paternal genetic heritage of eastern Plateau regions remains less explored. We present one integrative Y chromosome dataset of 9,901 modern and ancient individuals, including whole Y chromosome sequences from 1,297 individuals and extensive Y-SNP/STR genotype data. We reveal the Paleolithic common origin and following divergence of Qinghai-Xizang Plateau ancestors from East Asian lowlands, marked by subsequent isolation and Holocene expansion involving local hunter-gatherers and millet-farming communities. We identified two key TB-related founding lineages, D-Z31591 and O-CTS4658, which underwent significant expansions around 5,000 years ago on the Qinghai-Xizang Plateau and its eastern Tibetan-Yi Corridor. The genetic legacy of these TB lineages highlights crucial migration pathways linking the Plateau and lowland southwestern China. Our findings align paternal genetic structures with East Asian geography and linguistic groups, underscoring the utility of Y chromosome analyses in unraveling complex paternal histories.

Distinct patterns of genetic variation at low-recombining genomic regions represent haplotype structure

Genomic regions sometimes show patterns of genetic variation distinct from the genome-wide population structure. Such deviations have often been interpreted to represent effects of selection. However, systematic investigation of whether and how non-selective factors, such as recombination rates, can affect distinct patterns has been limited. Here, we associate distinct patterns of genetic variation with reduced recombination rates in a songbird, the Eurasian blackcap (Sylvia atricapilla), using a new reference genome assembly, whole-genome resequencing data and recombination maps. We find that distinct patterns of genetic variation reflect haplotype structure at genomic regions with different prevalence of reduced recombination rate across populations. At low-recombining regions shared in most populations, distinct patterns reflect conspicuous haplotypes segregating in multiple populations. At low-recombining regions found only in a few populations, distinct patterns represent variance among cryptic haplotypes within the low-recombining populations. With simulations, we confirm that these distinct patterns evolve neutrally by reduced recombination rate, on which the effects of selection can be overlaid. Our results highlight that distinct patterns of genetic variation can emerge through evolutionary reduction of local recombination rate. The recombination landscape as an evolvable trait therefore plays an important role determining the heterogeneous distribution of genetic variation along the genome.

Systematic identification of pathogenic variants of non-small cell lung cancer in the promoters of DNA-damage repair genes

BACKGROUND

Deficiency in DNA-damage repair (DDR) genes, often due to disruptive coding variants, is linked to higher cancer risk. Our previous study has revealed the association between rare loss-of-function variants in DDR genes and the risk of lung cancer. However, it is still challenging to study the predisposing role of rare regulatory variants of these genes.

METHODS

Based on whole-genome sequencing data from 2984 patients with non-small cell lung cancer (NSCLC) and 3020 controls, we performed massively parallel reporter assays on 1818 rare variants located in the promoters of DDR genes. Pathway- or gene-level burden analyses were performed using Firth's logistic regression or generalized linear model.

FINDINGS

We identified 750 rare functional regulatory variants (frVars) that showed allelic differences in transcriptional activity within the promoter regions of DDR genes. Interestingly, the burden of frVars was significantly elevated in cases (odds ratio [OR] = 1.17, p = 0.026), whereas the burden of variants prioritized solely based on bioinformatics annotation was comparable between cases and controls (OR = 1.04, p = 0.549). Among the frVars, 297 were down-regulated transcriptional activity (dr-frVars) and 453 were up-regulated transcriptional activity (ur-frVars); especially, dr-frVars (OR = 1.30, p = 0.008) rather than ur-frVars (OR = 1.06, p = 0.495) were significantly associated with risk of NSCLC. Individuals with NSCLC carried more dr-frVars from Fanconi anemia, homologous recombination, and nucleotide excision repair pathways. In addition, we identified seven genes (i.e., BRCA2, GTF2H1, DDB2, BLM, ALKBH2, APEX1, and RAD51B) with promoter dr-frVars that were associated with lung cancer susceptibility.

INTERPRETATION

Our findings indicate that functional promoter variants in DDR genes, in addition to protein-truncating variants, can be pathogenic and contribute to lung cancer susceptibility.

FUNDING

National Natural Science Foundation of China, Youth Foundation of Jiangsu Province, Research Unit of Prospective Cohort of Cardiovascular Diseases and Cancer of Chinese Academy of Medical Sciences, and Natural Science Foundation of Jiangsu Province.

Complete genomic profiles of 1496 Taiwanese reveal curated medical insights

INTRODUCTION

The population of Taiwan has a long history of ethno-cultural evolution. The Taiwanese population was isolated from other large populations such as the European, Han Chinese, and Japanese population. The Taiwan Biobank (TWB) project has built a nationwide database, particularly for personal whole-genome sequence (WGS) to facilitate basic and clinical collaboration nationally and internationally, making it one of the most valuable public datasets of the East Asian population.

OBJECTIVES

This study provides comprehensive medical genomic findings from TWB WGS data, for better characterization of disease susceptibility and the choice of ideal treatment regimens in Taiwanese population.

METHODS

We reanalyzed 1496 WGS using a PrecisionFDA Truth challenge winner method Sentieon DNAscope. Single nucleotide variants (SNV) and small insertions/deletions (INDEL) were benchmarked. We also analyzed pharmacogenomic (PGx) drug-associated alleles, and copy number variants (CNV). Multiple practicing clinicians reviewed and curated the clinically significant variants. Variant annotations can be browsed at TaiwanGenomes (https://genomes.tw).

RESULTS

We found that each participant had an average of 6,870.7 globally novel variants and 75.3% (831/1103) of the participants harbored at least one PharmGKB-selected high evidence level human leukocyte antigen (HLA) risk allele. 54 PharmGKB-reported high-level instances of evidence of Cytochrome P450 variant-drug pairs, with a population frequency of over 13.2%. We also identified 23 variants in the ACMG secondary finding V3 gene list from 25 participants, suggesting that 1.67% (25/1496) of the population is harboring at least one medical actionable variant. Our carrier status analyses suggest that one in 25 couples (3.94%) would risk having offspring with at least one pathogenic variant, which is in line with rates found in Japan and Singapore. For pathogenic CNV, we detected 6.88% and 2.02% carrier rates for alpha thalassemia and spinal muscular atrophy, respectively.

CONCLUSION

Our study highlights the overall medical insights of a complete Taiwanese genomic profile.

HiOmics: A cloud-based one-stop platform for the comprehensive analysis of large-scale omics data

Analyzing the vast amount of omics data generated comprehensively by high-throughput sequencing technology is of utmost importance for scientists. In this context, we propose HiOmics, a cloud-based platform equipped with nearly 300 plugins designed for the comprehensive analysis and visualization of omics data. HiOmics utilizes the Element Plus framework to craft a user-friendly interface and harnesses Docker container technology to ensure the reliability and reproducibility of data analysis results. Furthermore, HiOmics employs the Workflow Description Language and Cromwell engine to construct workflows, ensuring the portability of data analysis and simplifying the examination of intricate data. Additionally, HiOmics has developed DataCheck, a tool based on Golang, which verifies and converts data formats. Finally, by leveraging the object storage technology and batch computing capabilities of public cloud platforms, HiOmics enables the storage and processing of large-scale data while maintaining resource independence among users.

Understanding the variant landscape, and genetic epidemiology of Multiple Endocrine Neoplasia in India

PURPOSE

Multiple Endocrine Neoplasia (MEN) is a group of familial cancer syndromes that encompasses several types of endocrine tumors differentiated by genetic mutations in RET, MEN1 and CDKN1B genes. Accurate diagnosis of MEN subtypes can thus be performed through genetic testing. However, MEN variants remain largely understudied in Indian populations. Additionally, few dedicated resources to understand these disorders currently exist.

METHODS

Using the gold-standard ACMG/AMP guidelines, we systematically classified variants reported across the three genes in the IndiGen dataset, and established the genetic epidemiology of MEN in the Indian population. We further classified ClinVar and Mastermind variants and compiled all into a database. Finally, we designed a multiplex primer panel for rapid variant identification.

RESULTS

We have established the genetic prevalence of MEN as the following: 1 in 1026 individuals is likely to be afflicted with MEN linked with pathogenic RET mutations. We have further created the MAPVar database containing 3280 ACMG-classified variants freely accessible at: https://clingen.igib.res.in/MAPVar/ . Finally, our NGS primer panel covers 33 exonic regions across two pools through 38 amplicons with a total amplified region of 65 kb.

CONCLUSION

Our work establishes that MEN is a prevalent disorder in India. The rare nature of Indian variants underscores the need of genomic and functional studies to establish a more comprehensive variant landscape. Additionally, our panel offers a means of cost-effective genetic testing, and the MAPVar database a ready reference to aid in a better understanding of variant pathogenicity in clinical as well as research settings.

The 1000 Chinese Indigenous Pig Genomes Project provides insights into the genomic architecture of pigs

Pigs play a central role in human livelihoods in China, but a lack of systematic large-scale whole-genome sequencing of Chinese domestic pigs has hindered genetic studies. Here, we present the 1000 Chinese Indigenous Pig Genomes Project sequencing dataset, comprising 1011 indigenous individuals from 50 pig populations covering approximately two-thirds of China's administrative divisions. Based on the deep sequencing (~25.95×) of these pigs, we identify 63.62 million genomic variants, and provide a population-specific reference panel to improve the imputation performance of Chinese domestic pig populations. Using a combination of methods, we detect an ancient admixture event related to a human immigration climax in the 13th century, which may have contributed to the formation of southeast-central Chinese pig populations. Analyzing the haplotypes of the Y chromosome shows that the indigenous populations residing around the Taihu Lake Basin exhibit a unique haplotype. Furthermore, we find a 13 kb region in the THSD7A gene that may relate to high-altitude adaptation, and a 0.47 Mb region on chromosome 7 that is significantly associated with body size traits. These results highlight the value of our genomic resource in facilitating genomic architecture and complex traits studies in pigs.

PPARA variant rs1800234 had a dose dependent pharmacogenetics impact on the therapeutic response to chiglitazar

BACKGROUND

Our objective was to explore the pharmacogenetic impact of three known functional variants in drug target genes and determine whether they can explain the inter-individual variation in therapeutic response.

METHODS

In a post hoc analysis of data from randomized controlled clinical trials of chiglitazar, we genotyped 481 Chinese patients with T2DM and investigated the association of variants in PPAR genes with the therapeutic outcome separated by dose using linear regression.

RESULTS

rs1800234, a gain-of-function variant of PPARA, had a dose-dependent pharmacogenetic impact on the therapeutic response to chiglitazar. The C allele was significantly associated with reduced therapeutic response in the 48 mg group, while no significant association was observed in the 32 mg group. In addition, in patients without the C allele, patients treated with 48 mg chiglitazar had a better therapeutic response than those treated with 32 mg chiglitazar. To the contrary, in patients with the C allele, patients treated with 48 mg chiglitazar had a worse therapeutic response than those treated with 32 mg of chiglitazar.

CONCLUSION

The PPARA variant rs1800234 had a dose-dependent pharmacogenetic impact on the therapeutic response to chiglitazar. It could help explain the absence of a dose effect of chiglitazar and serve as a potential biomarker for the chosen dose of chiglitazar in the future. In addition, our study provided important reference for the design and clinical application of multi-target drugs.

MGA loss-of-function variants cause premature ovarian insufficiency

Although premature ovarian insufficiency (POI), a common cause of female infertility and subfertility, has a well-established hereditary component, the genetic factors currently implicated in POI account for only a limited proportion of cases. Here, using an exome-wide, gene-based case-control analysis in a discovery cohort comprising 1,027 POI cases and 2,733 ethnically matched women controls from China, we found that heterozygous loss-of-function (LoF) variants of MAX dimerization protein (MGA) were significantly enriched in the discovery cohort, accounting for 2.6% of POI cases, while no MGA LoF variants were found in the matched control females. Further exome screening was conducted in 4 additional POI cohorts (2 from China and 2 from the United States) for replication studies, and we identified heterozygous MGA LoF variants in 1.0%, 1.4%, 1.0%, and 1.0% of POI cases, respectively. Overall, a total of 37 distinct heterozygous MGA LoF variants were discovered in 38 POI cases, accounting for approximately 2.0% of the total 1,910 POI cases analyzed in this study. Accordingly, Mga+/- female mice were subfertile, exhibiting shorter reproductive lifespan and decreased follicle number compared with WT, mimicking the observed phenotype in humans. Our findings highlight the essential role of MGA deficiency for impaired female reproductive ability.

Novel variants and genotype-phenotype correlation in a multicentre cohort of GNE myopathy in China

BACKGROUND

GlcNAc2-epimerase (GNE) myopathy is a rare autosomal recessive disorder caused by pathogenic variants in the GNE gene, which is essential for the sialic acid biosynthesis pathway.

OBJECTIVE

This multi-centre study aimed to delineate the clinical phenotype and GNE variant spectrum in Chinese patients, enhancing our understanding of the genetic diversity and clinical manifestation across different populations.

METHODS

We retrospectively analysed GNE variants from 113 patients, integrating these data with external GNE variants from online databases for a global perspective, examining their consequences, distribution, ethnicity and severity.

RESULTS

This study revealed 97 distinct GNE variants, including 35 (36.08%) novel variants. Two more patients with deep intronic variant c.862+870C>T were identified, while whole genome sequencing (WGS) uncovered another two novel intronic variants: c.52-8924G>T and c.1505-12G>A. Nanopore long reads sequencing (LRS) and further PCR analysis verified a 639 bp insertion at chr9:36249241. Missense variants predominantly located in the epimerase/kinase domain coding region, indicating the impairment of catalytic function as a key pathogenic consequence. Comparative studies with Japanese, Korean and Jewish, our cohorts showed later onset ages by 2 years. The high allele frequency of the non-catalytic GNE variant, c.620A>T, might underlie the milder phenotype of Chinese patients.

CONCLUSIONS

Comprehensive techniques such as WGS and Nanopore LRS warrants the identifying of GNE variants. Patients with the non-catalytic GNE variant, c.620A>T, had a milder disease progression and later wheelchair use.

A genome-wide association study of neonatal metabolites

Genetic factors significantly influence the concentration of metabolites in adults. Nevertheless, the genetic influence on neonatal metabolites remains uncertain. To bridge this gap, we employed genotype imputation techniques on large-scale low-pass genome data obtained from non-invasive prenatal testing. Subsequently, we conducted association studies on a total of 75 metabolic components in neonates. The study identified 19 previously reported associations and 11 novel associations between single-nucleotide polymorphisms and metabolic components. These associations were initially found in the discovery cohort (8,744 participants) and subsequently confirmed in a replication cohort (19,041 participants). The average heritability of metabolic components was estimated to be 76.2%, with a range of 69%-78.8%. These findings offer valuable insights into the genetic architecture of neonatal metabolism.

Phenome-wide association study in 25,639 pregnant Chinese women reveals loci associated with maternal comorbidities and child health

Phenome-wide association studies (PheWAS) have been less focused on maternal diseases and maternal-newborn comorbidities, especially in the Chinese population. To enhance our understanding of the genetic basis of these related diseases, we conducted a PheWAS on 25,639 pregnant women and 14,151 newborns in the Chinese Han population using ultra-low-coverage whole-genome sequence (ulcWGS). We identified 2,883 maternal trait-associated SNPs associated with 26 phenotypes, among which 99.5% were near established genome-wide association study (GWAS) loci. Further refinement delineated these SNPs to 442 unique trait-associated loci (TALs) predicated on linkage disequilibrium R2 > 0.8, revealing that 75.6% demonstrated pleiotropy and 50.9% were located in genes implicated in analogous phenotypes. Notably, we discovered 21 maternal SNPs associated with 35 neonatal phenotypes, including two SNPs associated with identical complications in both mothers and children. These findings underscore the importance of integrating ulcWGS data to enrich the discoveries derived from traditional PheWAS approaches.

A Cysteinyl-tRNA Synthetase Mutation Causes Novel Autosomal-Dominant Inheritance of a Parkinsonism/Spinocerebellar-Ataxia Complex

This study aimed to identify possible pathogenic genes in a 90-member family with a rare combination of multiple neurodegenerative disease phenotypes, which has not been depicted by the known neurodegenerative disease. We performed physical and neurological examinations with International Rating Scales to assess signs of ataxia, Parkinsonism, and cognitive function, as well as brain magnetic resonance imaging scans with seven sequences. We searched for co-segregations of abnormal repeat-expansion loci, pathogenic variants in known spinocerebellar ataxia-related genes, and novel rare mutations via whole-genome sequencing and linkage analysis. A rare co-segregating missense mutation in the CARS gene was validated by Sanger sequencing and the aminoacylation activity of mutant CARS was measured by spectrophotometric assay. This pedigree presented novel late-onset core characteristics including cerebellar ataxia, Parkinsonism, and pyramidal signs in all nine affected members. Brain magnetic resonance imaging showed cerebellar/pons atrophy, pontine-midline linear hyperintensity, decreased rCBF in the bilateral basal ganglia and cerebellar dentate nucleus, and hypo-intensities of the cerebellar dentate nuclei, basal ganglia, mesencephalic red nuclei, and substantia nigra, all of which suggested neurodegeneration. Whole-genome sequencing identified a novel pathogenic heterozygous mutation (E795V) in the CARS gene, meanwhile, exhibited none of the known repeat-expansions or point mutations in pathogenic genes. Remarkably, this CARS mutation causes a 20% decrease in aminoacylation activity to charge tRNACys with L-cysteine in protein synthesis compared with that of the wild type. All family members carrying a heterozygous mutation CARS (E795V) had the same clinical manifestations and neuropathological changes of Parkinsonism and spinocerebellar-ataxia. These findings identify novel pathogenesis of Parkinsonism-spinocerebellar ataxia and provide insights into its genetic architecture.

Association of single nucleotide polymorphisms and gene-environment interactions with major depressive disorder in Chinese

We conducted a case-control study to investigate the effects of genetics and gene-environment interactions on major depressive disorders (MDD) in the Chinese population. Using targeted-exome sequencing, we included 984 patients with MDD and 508 healthy controls in our study. A logistic regression model was employed to analyze the association between single nucleotide polymorphisms (SNPs) and MDD. Additionally, a linear regression model was utilized to examine the associations between (1) gene-environment interaction and the 17-item Hamilton Depression Rating Scale, (2) SNPs and the Beck Scale for Suicide Ideation-Chinese version, and gene-environment interaction and the Beck scale for suicide ideation-Chinese version. The association analysis between SNPs and MDD revealed that the following loci reached genome-wide significance: rs2305554 of the cholinergic receptor nicotinic alpha 7 subunit, rs9459173 of synaptojanin 2, rs372369000 of beta-1,4-galactosyltransferase 6, rs866666526 of dopa decarboxylase, rs1254882194 of calcium/calmodulin dependent protein kinase ID, rs199880487 of reelin, rs1167948188 of reelin, rs1390140186 of QKI, KH domain containing RNA binding, and rs1776342 of period circadian regulator 3. The association analysis between SNPs and the Beck Scale for Suicide Ideation-Chinese version indicated that rs264272 and rs1774784888 of piezo type mechanosensitive ion channel component 2 reached genome-wide significance. These findings may enhance our understanding of MDD and contribute to the development of new potential targets for its diagnosis and treatment.

Prevalence Estimation of the PALB2 Germline Variant in East Asians and Koreans through Population Database Analysis

PALB2 is a tumor suppressor gene. Heterozygous germline pathogenic variants of PALB2 significantly increase the lifetime risk of breast cancer and moderately increase the risk of ovarian and pancreatic cancers. This study analyzed the estimated prevalence of PALB2 variants globally, focusing on East Asian and Korean populations, where limited data were previously available. We examined 125,748 exomes from the Genome Aggregation Database (gnomAD), including 9197 East Asians, and additional data from 5305 individuals in the Korean Variant Archive and 1722 in the Korean Reference Genome Database. All PALB2 variants were interpreted according to guidelines from the American College of Medical Genetics and Genomics and the Clinical Genome Resource. The global prevalence of PALB2 variants was 0.18%, with the highest prevalence in Finnish populations (0.41%) and the lowest in Ashkenazi Jewish populations (0.04%). East Asians had a prevalence of 0.09%. By combining data from Korean genome databases and gnomAD totaling 8936 individuals, the overall prevalence of PALB2 variants in the Korean population was determined to be 0.13%. This study is the first comprehensive investigation of PALB2 variant prevalence in East Asians and Koreans using gnomAD and Korean genome databases. These findings provide essential reference data for future research and highlight the importance of region-specific genetic studies that will inform genetic counseling and hereditary cancer risk management.

A synonymous mutation of rs1137070 cause the mice Maoa gene transcription and translation to decrease

The Monoamine Oxidase-A (MAOA) EcoRV polymorphism (rs1137070) is a unique synonymous mutation (c.1409 T > C) within the MAOA gene, which plays a crucial role in Maoa gene expression and function. This study aimed to explore the relationship between the mouse Maoa rs1137070 genotype and differences in MAOA gene expression. Mice carrying the CC genotype of rs1137070 exhibited a significantly lower Maoa expression level, with an odds ratio of 2.44 compared to the T carriers. Moreover, the wild-type TT genotype of MAOA demonstrated elevated mRNA expression and a longer half-life. We also delved into the significant expression and structural disparities among genotypes. Furthermore, it was evident that different aspartic acid synonymous codons within Maoa influenced both MAOA expression and enzyme activity, highlighting the association between rs1137070 and MAOA. To substantiate these findings, a dual-luciferase reporter assay confirmed that GAC was more efficient than GAT binding. Conversely, the synonymous mutation altered Maoa gene expression in individual mice. An RNA pull-down assay suggested that this alteration could impact the interaction with RNA-binding proteins. In summary, our results illustrate that synonymous mutations can indeed regulate the downregulation of gene expression, leading to changes in MAOA function and their potential association with neurological-related diseases.

Estimating carrier rates and prevalence of porphyria-associated gene variants in the Chinese population based on genetic databases

Porphyria is a group of rare metabolic disorders caused by mutations in the genes encoding crucial enzymes in the heme biosynthetic pathway. However, the lack of comprehensive genetic analysis of porphyria patients in the Chinese population makes identifying and diagnosing carriers of the condition challenging. Using the ChinaMAP database, we determined the frequencies of P/LP porphyria-associated gene variants according to the ACMG guidelines. We also calculated the carrier rates and prevalence of each type of porphyria in the Chinese population under Hardy-Weinberg equilibrium. Compared with the variants in the gnomAD database, the genetic spectrum of porphyria-related P/LP variants in the Chinese population is distinct. In the ChinaMAP database, we identified 23 variants. We estimated the carrier rates for autosomal dominant porphyrias (AIP, HCP, VP, PCT) in the Chinese population to be 1/1059, 1/1513, 1/10588, and 1/1765, respectively. For autosomal recessive porphyrias (ADP, EPP, HEP, CEP), the estimated carrier rates were 1/5294, 1/2117, 1/1765, and 1/2647, respectively, with predicted prevalence rates of 8.92 × 10-9, 7.51 × 10-5, 8.02 × 10-8, and 3.57 × 10-8, respectively. Notably, 12 of the variants we identified were unique to the Chinese population. The predicted prevalence rate of EPP was the highest among the various types of porphyria in the Chinese population, while the others were moderate to low. This is the first comprehensive genetic study on porphyria in the Chinese population. Clarifying the genetic characteristics of various porphyria types among the Chinese population provides scientifically sound reference data for both research and genetic screening to identify porphyria carriers.

A genetic investigation in five Chinese families with keratoconus

Background

This study investigated the genetic characteristics of five Chinese families with keratoconus (KC).

Methods

In the five families affected by KC, medical records, clinical observations, and blood samples were collected from all individuals. All KC family members (n = 20) underwent both whole exome sequencing of genomic DNA and Sanger sequencing to confirm the variants. Online software was utilized to analyze all variants, and the online server I-TASSER was employed for in silico predictions of the three-dimensional protein structures of the variants. The newly discovered variants and single nucleotide polymorphisms were further examined in 322 sporadic KC patients.

Results

The Pentacam tomographic composite index in those affected first-degree family members of the probands showed a pathological change. Five new variants were detected in the five probands and other affected members in their families: a heterozygous missense variant g.19043832C>T (p.Ser145Asn) in the homer scaffolding protein 3 (HOMER3) gene; a heterozygous missense variant g.99452113G>A (p.Gly483Arg) in the insulin-like growth factor 1 receptor (IGF1R) gene; a heterozygous missense variant g.55118280G>T (p.Trp843Leu) in the echinoderm microtubule-associated protein like 6 (EML6) gene; a heterozygous frameshift variant c. 1226_1227del (p.Gln410Glufs*17) in the DOP1 leucine zipper-like protein B (DOP1B) gene; and a heterozygous splice-site variant c.7776+2T>A in the neurobeachin-like protein 2 (NBEAL2) gene. These variations were predicted to be potentially pathogenic and associated with KC.

Conclusion

Five novel variants in HOMER3, IGF1R, EML6, DOP1B, and NBEAL2 genes were identified in this study and may be associated with the pathogenesis of KC. This study provides new information about the gene variants and their protein changes in KC patients. The findings should be explored further and could potentially be applied to the early diagnosis of KC before clinical onset.

Obesity-enriched gut microbe degrades myo-inositol and promotes lipid absorption

Numerous studies have reported critical roles for the gut microbiota in obesity. However, the specific microbes that causally contribute to obesity and the underlying mechanisms remain undetermined. Here, we conducted shotgun metagenomic sequencing in a Chinese cohort of 631 obese subjects and 374 normal-weight controls and identified a Megamonas-dominated, enterotype-like cluster enriched in obese subjects. Among this cohort, the presence of Megamonas and polygenic risk exhibited an additive impact on obesity. Megamonas rupellensis possessed genes for myo-inositol degradation, as demonstrated in vitro and in vivo, and the addition of myo-inositol effectively inhibited fatty acid absorption in intestinal organoids. Furthermore, mice colonized with M. rupellensis or E. coli heterologously expressing the myo-inositol-degrading iolG gene exhibited enhanced intestinal lipid absorption, thereby leading to obesity. Altogether, our findings uncover roles for M. rupellensis as a myo-inositol degrader that enhances lipid absorption and obesity, suggesting potential strategies for future obesity management.