With expanding applications of next-generation sequencing in medical genetics, increasing computational methods are being developed to predict the pathogenicity of missense variants. Selecting optimal methods can accelerate the identification of candidate genes. However, the performances of different computational methods under various conditions have not been completely evaluated. Here, we compared 12 performance measures of 23 methods based on three independent benchmark datasets: (i) clinical variants from the ClinVar database related to genetic diseases, (ii) somatic variants from the IARC TP53 and ICGC databases related to human cancers and (iii) experimentally evaluated PPARG variants. Some methods showed different performances under different conditions, suggesting that they were not always applicable for different conditions. Furthermore, the specificities were lower than the sensitivities for most methods (especially, for the experimentally evaluated benchmark datasets), suggesting that more rigorous cutoff values are necessary to distinguish pathogenic variants. Furthermore, REVEL, VEST3 and the combination of both methods (i.e. ReVe) showed the best overall performances with all the benchmark data. Finally, we evaluated the performances of these methods with de novo mutations, finding that ReVe consistently showed the best performance. We have summarized the performances of different methods under various conditions, providing tentative guidance for optimal tool selection.

Autism spectrum disorder (ASD) is a complex neurodevelopmental disorder with a male-to-female prevalence of 4:1. However, the genetic mechanisms underlying this gender difference remain unclear. Mutation burden analysis, a TADA model, and co-expression and functional network analyses were performed on de novo mutations (DNMs) and corresponding candidate genes. We found that the prevalence of putative functional DNMs (loss-of-function and predicted deleterious missense mutations) in females was significantly higher than that in males, suggesting that a higher genetic load was required in females to reach the threshold for a diagnosis. We then prioritized 174 candidate genes, including 60 shared genes, 91 male-specific genes, and 23 female-specific genes. All of the three subclasses of candidate genes were significantly more frequently co-expressed in female brains than male brains, suggesting that compensation effects of the deficiency of ASD candidate genes may be more likely in females. Nevertheless, the three subclasses of candidate genes were co-expressed with each other, suggesting a convergent functional network of male and female-specific genes. Our analysis of different aspects of genetic components provides suggestive evidence supporting the female-protective effect in ASD. Moreover, further study is needed to integrate neuronal and hormonal data to elucidate the underlying gender difference in ASD.

This study aimed to determine the mutational spectrum of familial Parkinson’s disease and sporadic early-onset Parkinson’s disease (sEOPD) in a mainland Chinese population and the clinical features of mutation carriers. We performed multiplex ligation-dependent probe amplification assays and whole-exome sequencing for 1676 unrelated patients with Parkinson’s disease in a mainland Chinese population, including 192 probands from families with autosomal-recessive Parkinson’s disease, 242 probands from families with autosomal-dominant Parkinson’s disease, and 1242 sEOPD patients (age at onset ≤ 50). According to standards and guidelines from the American College of Medical Genetics and Genomics, pathogenic/likely pathogenic variants in 23 known Parkinson’s disease-associated genes occurred more frequently in the autosomal-recessive Parkinson’s disease cohort (65 of 192, 33.85%) than in the autosomal-dominant Parkinson’s disease cohort (10 of 242, 4.13%) and the sEOPD cohort (57 of 1242, 4.59%), which leads to an overall molecular diagnostic yield of 7.88% (132 of 1676). We found that PRKN was the most frequently mutated gene (n = 83, 4.95%) and present the first evidence of an SNCA duplication and LRRK2 p.N1437D variant in mainland China. In addition, several novel pathogenic/likely pathogenic variants including LRRK2 (p.V1447M and p.Y1645S), ATP13A2 (p.R735X and p.A819D), FBXO7 (p.G67E), LRP10 (c.322dupC/p.G109Rfs*51) and TMEM230 (c.429delT/p.P144Qfs*2) were identified in our cohort. Furthermore, the age at onset of the 132 probands with genetic diagnoses (median, 31.5 years) was about 14.5 years earlier than that of patients without molecular diagnoses (i.e. non-carriers, median 46.0 years). Specifically, the age at onset of Parkinson’s disease patients with pathogenic/likely pathogenic variants in ATP13A2, PLA2G6, PRKN, or PINK1 was significantly lower than that of non-carriers, while the age at onset of carriers with other gene pathogenic/likely pathogenic variants was similar to that of non-carriers. The clinical spectrum of Parkinson’s disease-associated gene carriers in this mainland Chinese population was similar to that of other populations. We also detected 61 probands with GBA possibly pathogenic variants (3.64%) and 59 probands with GBA p.L444P (3.52%). These results shed insight into the genetic spectrum and clinical manifestations of Parkinson’s disease in mainland China and expand the existing repertoire of pathogenic or likely pathogenic variants involved in known Parkinson’s disease-associated genes. Our data highlight the importance of genetic testing in Parkinson’s disease patients with age at onset < 40 years, especially in those from families with a recessive inheritance pattern, who may benefit from early diagnosis and treatment.

We initiate the Westlake BioBank for Chinese (WBBC) pilot project with 4,535 whole-genome sequencing (WGS) individuals and 5,841 high-density genotyping individuals, and identify 81.5 million SNPs and INDELs, of which 38.5% are absent in dbSNP Build 151. We provide a population-specific reference panel and an online imputation server (https://wbbc.westlake.edu.cn/) which could yield substantial improvement of imputation performance in Chinese population, especially for low-frequency and rare variants. By analyzing the singleton density of the WGS data, we find selection signatures in SNX29, DNAH1 and WDR1 genes, and the derived alleles of the alcohol metabolism genes (ADH1A and ADH1B) emerge around 7,000 years ago and tend to be more common from 4,000 years ago in East Asia. Genetic evidence supports the corresponding geographical boundaries of the Qinling-Huaihe Line and Nanling Mountains, which separate the Han Chinese into subgroups, and we reveal that North Han was more homogeneous than South Han.

The metabolomic profile of aging is complex. Here, we analyse 325 nuclear magnetic resonance (NMR) biomarkers from 250,341 UK Biobank participants, identifying 54 representative aging-related biomarkers associated with all-cause mortality. We conduct genome-wide association studies (GWAS) for these 325 biomarkers using whole-genome sequencing (WGS) data from 95,372 individuals and perform multivariable Mendelian randomization (MVMR) analyses, discovering 439 candidate “biomarker - disease” causal pairs at the nominal significance level. We develop a metabolomic aging score that outperforms other aging metrics in predicting short-term mortality risk and exhibits strong potential for discriminating aging-accelerated populations and improving disease risk prediction. A longitudinal analysis of 13,263 individuals enables us to calculate a metabolomic aging rate which provides more refined aging assessments and to identify candidate anti-aging and pro-aging NMR biomarkers. Taken together, our study has presented a comprehensive aging-related metabolomic profile and highlighted its potential for personalized aging monitoring and early disease intervention.

Long single-molecular sequencing technologies, such as PacBio circular consensus sequencing (CCS) and nanopore sequencing, are advantageous in detecting DNA 5-methylcytosine in CpGs (5mCpGs), especially in repetitive genomic regions. However, existing methods for detecting 5mCpGs using PacBio CCS are less accurate and robust. Here, we present ccsmeth, a deep-learning method to detect DNA 5mCpGs using CCS reads. We sequence polymerase-chain-reaction treated and M.SssI-methyltransferase treated DNA of one human sample using PacBio CCS for training ccsmeth. Using long (≥10 Kb) CCS reads, ccsmeth achieves 0.90 accuracy and 0.97 Area Under the Curve on 5mCpG detection at single-molecule resolution. At the genome-wide site level, ccsmeth achieves >0.90 correlations with bisulfite sequencing and nanopore sequencing using only 10× reads. Furthermore, we develop a Nextflow pipeline, ccsmethphase, to detect haplotype-aware methylation using CCS reads, and then sequence a Chinese family trio to validate it. ccsmeth and ccsmethphase can be robust and accurate tools for detecting DNA 5-methylcytosines.