Whole-Genome Sequencing of a Healthy Aging Cohort

The dual role of the DREAM/G2M pathway in non-tumorigenic immortalization of senescent cells

Anti-aging and tumorigenesis share common genes and pathways, and thus targeting these genes as part of anti-aging interventions carries the risk of tumorigenesis. It is essential to understand the gene signatures that balance tumorigenesis and aging. To achieve this goal, we analyzed RNA-sequencing data from three non-tumorigenic immortalized cell lines that spontaneously escaped from senescence. By single sample gene set enrichment assay (ssGSEA) and GSEA analysis, we found that both cell growth signaling (E2F targets, MYC targets) and tumor surveillance mechanisms (DNA repair, G2M checkpoint, mitotic spindle) were up-regulated in all three cell lines, suggesting that these genes are potential signatures for non-tumorigenic immortalization. Further analysis revealed that the 182 commonly up-regulated genes in these three cell lines overlapped with the DREAM/G2M pathway, which is known to be the upstream regulator of E2F, Myc targets, DNA repair, G2M checkpoint and mitotic spindle pathways in its cell cycle activation or inhibitory form. By western blotting, quantitative PCR and co-immunoprecipitation, we verified that both forms of the DREAM pathway are up-regulated in all three cell lines; this pathway facilitates control of cell cycle progression, supporting a new mechanism for non-tumorigenic immortalization. Thus, we propose that the DREAM/G2M pathway plays important dual roles with respect to preventing tumorigenesis in the process of immortalization. Our data might serve as the basis for the identification of new signature pathways or gene biomarkers for non-tumorigenic immortalization, and may aid in the discovery of new targets for tumor-free anti-aging drug screening.

Polygenic prediction of human longevity on the supposition of pervasive pleiotropy

The highly polygenic nature of human longevity renders cross-trait pleiotropy an indispensable feature of its genetic architecture. Leveraging the genetic correlation between the aging-related traits (ARTs), we sought to model the additive variance in lifespan as a function of cumulative liability from pleiotropic segregating variants. We tracked allele frequency changes as a function of viability across different age bins and prioritized 34 variants with an immediate implication on lipid metabolism, body mass index (BMI), and cognitive performance, among other traits, revealed by PheWAS analysis in the UK Biobank. Given the highly complex and non-linear interactions between the genetic determinants of longevity, we reasoned that a composite polygenic score would approximate a substantial portion of the variance in lifespan and developed the integrated longevity genetic scores (iLGSs) for distinguishing exceptional survival. We showed that coefficients derived from our ensemble model could potentially reveal an interesting pattern of genomic pleiotropy specific to lifespan. We assessed the predictive performance of our model for distinguishing the enrichment of exceptional longevity among long-lived individuals in two replication cohorts and showed that the median lifespan in the highest decile of our composite prognostic index is up to 4.8 years longer. Finally, using the proteomic correlates of i L G S , we identified protein markers associated with exceptional longevity irrespective of chronological age and prioritized drugs with repurposing potentials for gerotherapeutics. Together, our approach demonstrates a promising framework for polygenic modeling of additive liability conferred by ARTs in defining exceptional longevity and assisting the identification of individuals at higher risk of mortality for targeted lifestyle modifications earlier in life. Furthermore, the proteomic signature associated with i L G S highlights the functional pathway upstream of the PI3K-Akt that can be effectively targeted to slow down aging and extend lifespan.

Quartet DNA reference materials and datasets for comprehensively evaluating germline variant calling performance

BACKGROUND

Genomic DNA reference materials are widely recognized as essential for ensuring data quality in omics research. However, relying solely on reference datasets to evaluate the accuracy of variant calling results is incomplete, as they are limited to benchmark regions. Therefore, it is important to develop DNA reference materials that enable the assessment of variant detection performance across the entire genome.

RESULTS

We established a DNA reference material suite from four immortalized cell lines derived from a family of parents and monozygotic twins. Comprehensive reference datasets of 4.2 million small variants and 15,000 structural variants were integrated and certified for evaluating the reliability of germline variant calls inside the benchmark regions. Importantly, the genetic built-in-truth of the Quartet family design enables estimation of the precision of variant calls outside the benchmark regions. Using the Quartet reference materials along with study samples, batch effects are objectively monitored and alleviated by training a machine learning model with the Quartet reference datasets to remove potential artifact calls. Moreover, the matched RNA and protein reference materials and datasets from the Quartet project enables cross-omics validation of variant calls from multiomics data.

CONCLUSIONS

The Quartet DNA reference materials and reference datasets provide a unique resource for objectively assessing the quality of germline variant calls throughout the whole-genome regions and improving the reliability of large-scale genomic profiling.

Large-scale rare variant burden testing in Parkinson's disease

Parkinson's disease has a large heritable component and genome-wide association studies have identified over 90 variants with disease-associated common variants, providing deeper insights into the disease biology. However, there have not been large-scale rare variant analyses for Parkinson's disease. To address this gap, we investigated the rare genetic component of Parkinson's disease at minor allele frequencies <1%, using whole genome and whole exome sequencing data from 7184 Parkinson's disease cases, 6701 proxy cases and 51 650 healthy controls from the Accelerating Medicines Partnership Parkinson's disease (AMP-PD) initiative, the National Institutes of Health, the UK Biobank and Genentech. We performed burden tests meta-analyses on small indels and single nucleotide protein-altering variants, prioritized based on their predicted functional impact. Our work identified several genes reaching exome-wide significance. Two of these genes, GBA1 and LRRK2, have variants that have been previously implicated as risk factors for Parkinson's disease, with some variants in LRRK2 resulting in monogenic forms of the disease. We identify potential novel risk associations for variants in B3GNT3, AUNIP, ADH5, TUBA1B, OR1G1, CAPN10 and TREML1 but were unable to replicate the observed associations across independent datasets. Of these, B3GNT3 and TREML1 could provide new evidence for the role of neuroinflammation in Parkinson's disease. To date, this is the largest analysis of rare genetic variants in Parkinson's disease.

Whole exome sequencing of patients with varicella-zoster virus and herpes simplex virus induced acute retinal necrosis reveals rare disease-associated genetic variants

Purpose

Herpes simplex virus (HSV) and varicella-zoster virus (VZV) are neurotropic human alphaherpesviruses endemic worldwide. Upon primary infection, both viruses establish lifelong latency in neurons and reactivate intermittently to cause a variety of mild to severe diseases. Acute retinal necrosis (ARN) is a rare, sight-threatening eye disease induced by ocular VZV or HSV infection. The virus and host factors involved in ARN pathogenesis remain incompletely described. We hypothesize an underlying genetic defect in at least part of ARN cases.

Methods

We collected blood from 17 patients with HSV-or VZV-induced ARN, isolated DNA and performed Whole Exome Sequencing by Illumina followed by analysis in Varseq with criteria of CADD score > 15 and frequency in GnomAD < 0.1% combined with biological filters. Gene modifications relative to healthy control genomes were filtered according to high quality and read-depth, low frequency, high deleteriousness predictions and biological relevance.

Results

We identified a total of 50 potentially disease-causing genetic variants, including missense, frameshift and splice site variants and on in-frame deletion in 16 of the 17 patients. The vast majority of these genes are involved in innate immunity, followed by adaptive immunity, autophagy, and apoptosis; in several instances variants within a given gene or pathway was identified in several patients.

Discussion

We propose that the identified variants may contribute to insufficient viral control and increased necrosis ocular disease presentation in the patients and serve as a knowledge base and starting point for the development of improved diagnostic, prophylactic, and therapeutic applications.

Targeted sequencing of the 9p21.3 region reveals association with reduced disease risks in Ashkenazi Jewish centenarians

Genome-wide association studies (GWAS) have pinpointed the chromosomal locus 9p21.3 as a genetic hotspot for various age-related disorders. Common genetic variants in this locus are linked to multiple traits, including coronary artery diseases, cancers, and diabetes. Centenarians are known for their reduced risk and delayed onset of these conditions. To investigate whether this evasion of disease risks involves diminished genetic risks in the 9p21.3 locus, we sequenced this region in an Ashkenazi Jewish centenarian cohort (centenarians: n = 450, healthy controls: n = 500). Risk alleles associated with cancers, glaucoma, CAD, and T2D showed a significant depletion in centenarians. Furthermore, the risk and non-risk genotypes are linked to two distinct low-frequency variant profiles, enriched in controls and centenarians, respectively. Our findings provide evidence that the extreme longevity cohort is associated with collectively lower risks of multiple age-related diseases in the 9p21.3 locus.

Removal of extracellular human amyloid beta aggregates by extracellular proteases in C. elegans

The amyloid beta (Aβ) plaques found in Alzheimer's disease (AD) patients' brains contain collagens and are embedded extracellularly. Several collagens have been proposed to influence Aβ aggregate formation, yet their role in clearance is unknown. To investigate the potential role of collagens in forming and clearance of extracellular aggregates in vivo, we created a transgenic Caenorhabditis elegans strain that expresses and secretes human Aβ1-42. This secreted Aβ forms aggregates in two distinct places within the extracellular matrix. In a screen for extracellular human Aβ aggregation regulators, we identified different collagens to ameliorate or potentiate Aβ aggregation. We show that a disintegrin and metalloprotease a disintegrin and metalloprotease 2 (ADM-2), an ortholog of ADAM9, reduces the load of extracellular Aβ aggregates. ADM-2 is required and sufficient to remove the extracellular Aβ aggregates. Thus, we provide in vivo evidence of collagens essential for aggregate formation and metalloprotease participating in extracellular Aβ aggregate removal.

Genome-wide structural variant analysis identifies risk loci for non-Alzheimer's dementias

We characterized the role of structural variants, a largely unexplored type of genetic variation, in two non-Alzheimer's dementias, namely Lewy body dementia (LBD) and frontotemporal dementia (FTD)/amyotrophic lateral sclerosis (ALS). To do this, we applied an advanced structural variant calling pipeline (GATK-SV) to short-read whole-genome sequence data from 5,213 European-ancestry cases and 4,132 controls. We discovered, replicated, and validated a deletion in TPCN1 as a novel risk locus for LBD and detected the known structural variants at the C9orf72 and MAPT loci as associated with FTD/ALS. We also identified rare pathogenic structural variants in both LBD and FTD/ALS. Finally, we assembled a catalog of structural variants that can be mined for new insights into the pathogenesis of these understudied forms of dementia.

Contribution of A-to-I RNA editing, M6A RNA Methylation, and Alternative Splicing to physiological brain aging and neurodegenerative diseases

Aging is a physiological and progressive phenomenon in all organisms' life cycle, characterized by the accumulation of degenerative processes triggered by several alterations within molecular pathways. These changes compromise cell fate, resulting in the loss of functions in tissues throughout the body, including the brain. Physiological brain aging has been linked to structural and functional alterations, as well as to an increased risk of neurodegenerative diseases. Post-transcriptional RNA modifications modulate mRNA coding properties, stability, translatability, expanding the coding capacity of the genome, and are involved in all cellular processes. Among mRNA post-transcriptional modifications, the A-to-I RNA editing, m6A RNA Methylation and Alternative Splicing play a critical role in all the phases of a neuronal cell life cycle and alterations in their mechanisms of action significantly contribute to aging and neurodegeneration. Here we review our current understanding of the contribution of A-to-I RNA editing, m6A RNA Methylation, and Alternative Splicing to physiological brain aging process and neurodegenerative diseases.

A germline STAT6 gain-of-function variant is associated with early-onset allergies

BACKGROUND

The signal transducer and activator of transcription 6 (STAT6) signaling pathway plays a central role in allergic inflammation. To date, however, there have been no descriptions of STAT6 gain-of-function variants leading to allergies in humans.

OBJECTIVE

We report a STAT6 gain-of-function variant associated with early-onset multiorgan allergies in a family with 3 affected members.

METHODS

Exome sequencing and immunophenotyping of T-helper cell subsets were conducted. The function of the STAT6 protein was analyzed by Western blot, immunofluorescence, electrophoretic mobility shift assays, and luciferase assays. Gastric organoids obtained from the index patient were used to study downstream effector cytokines.

RESULTS

We identified a heterozygous missense variant (c.1129G>A;p.Glu377Lys) in the DNA binding domain of STAT6 that was de novo in the index patient's father and was inherited by 2 of his 3 children. Severe atopic dermatitis and food allergy were key presentations. Clinical heterogeneity was observed among the affected individuals. Higher levels of peripheral blood T_H2 lymphocytes were detected. The mutant STAT6 displayed a strong preference for nuclear localization, increased DNA binding affinity, and spontaneous transcriptional activity. Moreover, gastric organoids showed constitutive activation of STAT6 downstream signaling molecules.

CONCLUSIONS

A germline STAT6 gain-of-function variant results in spontaneous activation of the STAT6 signaling pathway and is associated with an early-onset and severe allergic phenotype in humans. These observations enhance our knowledge of the molecular mechanisms underlying allergic diseases and will potentially contribute to novel therapeutic interventions.

SP1 impacts the primordial to primary follicle transition by regulating cholesterol metabolism in granulosa cells

The primordial to primary follicle transition (PPT) in the ovary is critical to maintain sustainable reproductive resources in female mammals. However, it is unclear how granulosa cells (GCs) of the primary follicle participate in regulating PPT. This study focused on exploring the role of transcription factor Sp1 (SP1) in regulating PPT based on the fact that SP1 is pivotal for pregranulosa cell proliferation before primordial follicle formation. The results showed that mice fertility was prolonged when Sp1 was specifically depleted from GCs (GC- Sp1 ^-/- ). Besides, the PPT in GC- Sp1 ^-/- mice was reduced, resulting in more primordial follicles being preserved. Single-cell RNA-seq also indicated that the level of cholesterol metabolism was downregulated in GC- Sp1 ^-/- mice. Additionally, the PPT was promoted by either overexpression of ferredoxin-1 (FDX1), one of the key genes in mediating cholesterol metabolism or supplementing cholesterol for cultured fetal ovaries. Collectively, SP1 in GCs participates in the metabolism of cholesterol partially by regulating the transcription of Fdx1 during the PPT.

Heritable defects in telomere and mitotic function selectively predispose to sarcomas

Cancer genetics has to date focused on epithelial malignancies, identifying multiple histotype-specific pathways underlying cancer susceptibility. Sarcomas are rare malignancies predominantly derived from embryonic mesoderm. To identify pathways specific to mesenchymal cancers, we performed whole-genome germline sequencing on 1644 sporadic cases and 3205 matched healthy elderly controls. Using an extreme phenotype design, a combined rare-variant burden and ontologic analysis identified two sarcoma-specific pathways involved in mitotic and telomere functions. Variants in centrosome genes are linked to malignant peripheral nerve sheath and gastrointestinal stromal tumors, whereas heritable defects in the shelterin complex link susceptibility to sarcoma, melanoma, and thyroid cancers. These studies indicate a specific role for heritable defects in mitotic and telomere biology in risk of sarcomas.

Hallmarks of aging: An expanding universe

Aging is driven by hallmarks fulfilling the following three premises: (1) their age-associated manifestation, (2) the acceleration of aging by experimentally accentuating them, and (3) the opportunity to decelerate, stop, or reverse aging by therapeutic interventions on them. We propose the following twelve hallmarks of aging: genomic instability, telomere attrition, epigenetic alterations, loss of proteostasis, disabled macroautophagy, deregulated nutrient-sensing, mitochondrial dysfunction, cellular senescence, stem cell exhaustion, altered intercellular communication, chronic inflammation, and dysbiosis. These hallmarks are interconnected among each other, as well as to the recently proposed hallmarks of health, which include organizational features of spatial compartmentalization, maintenance of homeostasis, and adequate responses to stress.

A review on the application of the exposome paradigm to unveil the environmental determinants of age-related diseases

Age-related diseases account for almost half of all diseases among adults worldwide, and their incidence is substantially affected by the exposome, which is the sum of all exogenous and endogenous environmental exposures and the human body’s response to these exposures throughout the entire lifespan. Herein, we perform a comprehensive review of the epidemiological literature to determine the key elements of the exposome that affect the development of age-related diseases and the roles of aging hallmarks in this process. We find that most exposure assessments in previous aging studies have used a reductionist approach, whereby the effect of only a single environmental factor or a specific class of environmental factors on the development of age-related diseases has been examined. As such, there is a lack of a holistic and unbiased understanding of the effect of multiple environmental factors on the development of age-related diseases. To address this, we propose several research strategies based on an exposomic framework that could advance our understanding—in particular, from a mechanistic perspective—of how environmental factors affect the development of age-related diseases. We discuss the statistical methods and other methods that have been used in exposome-wide association studies, with a particular focus on multiomics technologies. We also address future challenges and opportunities in the realm of multidisciplinary approaches and genome–exposome epidemiology. Furthermore, we provide perspectives on precise public health services for vulnerable populations, public communications, the integration of risk exposure information, and the bench-to-bedside translation of research on age-related diseases.

Ribosomal protein S6 kinase beta-1 gene variants cause hypertrophic cardiomyopathy

BACKGROUND

Hypertrophic cardiomyopathy (HCM) is a genetic heart muscle disease with preserved or increased ejection fraction in the absence of secondary causes. Mutations in the sarcomeric protein-encoding genes predominantly cause HCM. However, relatively little is known about the genetic impact of signalling proteins on HCM.

METHODS AND RESULTS

Here, using exome and targeted sequencing methods, we analysed two independent cohorts comprising 401 Indian patients with HCM and 3521 Indian controls. We identified novel variants in ribosomal protein S6 kinase beta-1 (RPS6KB1 or S6K1) gene in two unrelated Indian families as a potential candidate gene for HCM. The two unrelated HCM families had the same heterozygous missense S6K1 variant (p.G47W). In a replication association study, we identified two S6K1 heterozygotes variants (p.Q49K and p.Y62H) in the UK Biobank cardiomyopathy cohort (n=190) compared with matched controls (n=16 479). These variants are neither detected in region-specific controls nor in the human population genome data. Additionally, we observed an S6K1 variant (p.P445S) in an Arab patient with HCM. Functional consequences were evaluated using representative S6K1 mutated proteins compared with wild type in cellular models. The mutated proteins activated the S6K1 and hyperphosphorylated the rpS6 and ERK1/2 signalling cascades, suggesting a gain-of-function effect.

CONCLUSIONS

Our study demonstrates for the first time that the variants in the S6K1 gene are associated with HCM, and early detection of the S6K1 variant carriers can help to identify family members at risk and subsequent preventive measures. Further screening in patients with HCM with different ethnic populations will establish the specificity and frequency of S6K1 gene variants.

Rapid, Reference-Free human genotype imputation with denoising autoencoders

Genotype imputation is a foundational tool for population genetics. Standard statistical imputation approaches rely on the co-location of large whole-genome sequencing-based reference panels, powerful computing environments, and potentially sensitive genetic study data. This results in computational resource and privacy-risk barriers to access to cutting-edge imputation techniques. Moreover, the accuracy of current statistical approaches is known to degrade in regions of low and complex linkage disequilibrium. Artificial neural network-based imputation approaches may overcome these limitations by encoding complex genotype relationships in easily portable inference models. Here, we demonstrate an autoencoder-based approach for genotype imputation, using a large, commonly used reference panel, and spanning the entirety of human chromosome 22. Our autoencoder-based genotype imputation strategy achieved superior imputation accuracy across the allele-frequency spectrum and across genomes of diverse ancestry, while delivering at least fourfold faster inference run time relative to standard imputation tools.

The pZRS non-coding regulatory mutation resulting in triphalangeal thumb-polysyndactyly syndrome changes the pattern of local interactions

Herein, we report on a large Polish family presenting with a classical triphalangeal thumb-polysyndactyly syndrome (TPT-PS). This rare congenital limb anomaly is generally caused by microduplications encompassing the Sonic Hedgehog (SHH) limb enhancer, termed the zone of polarizing activity (ZPA) regulatory sequence (ZRS). Recently, a pathogenic variant in the pre-ZRS (pZRS), a conserved sequence located near the ZRS, has been described in a TPT-PS Dutch family. We performed targeted ZRS sequencing, array comparative genomic hybridization, and whole-exome sequencing. Next, we sequenced the recently described pZRS region. Finally, we performed a circular chromatin conformation capture-sequencing (4C-seq) assay on skin fibroblasts of one affected family member and control samples to examine potential alterations in the SHH regulatory domain and functionally characterize the identified variant. We found that all affected individuals shared a recently identified pathogenic point mutation in the pZRS region: NC_000007.14:g.156792782C>G (GRCh38/hg38), which is the same as in the Dutch family. The results of 4C-seq experiments revealed increased interactions within the whole SHH regulatory domain (SHH-LMBR1 TAD) in the patient compared to controls. Our study expands the number of TPT-PS families carrying a pathogenic alteration of the pZRS and underlines the importance of routine pZRS sequencing in the genetic diagnostics of patients with TPT-PS or similar phenotypes. The pathogenic mutation causative for TPT-PS in our patient gave rise to increased interactions within the SHH regulatory domain in yet unknown mechanism.

Identification of ultra-rare genetic variants in pediatric acute onset neuropsychiatric syndrome (PANS) by exome and whole genome sequencing

Abrupt onset of severe neuropsychiatric symptoms including obsessive-compulsive disorder, tics, anxiety, mood swings, irritability, and restricted eating is described in children with Pediatric Acute-Onset Neuropsychiatric Syndrome (PANS). Symptom onset is often temporally associated with infections, suggesting an underlying autoimmune/autoinflammatory etiology, although direct evidence is often lacking. The pathological mechanisms are likely heterogeneous, but we hypothesize convergence on one or more biological pathways. Consequently, we conducted whole exome sequencing (WES) on a U.S. cohort of 386 cases, and whole genome sequencing (WGS) on ten cases from the European Union who were selected because of severe PANS. We focused on identifying potentially deleterious genetic variants that were de novo or ultra-rare (MAF) < 0.001. Candidate mutations were found in 11 genes (PPM1D, SGCE, PLCG2, NLRC4, CACNA1B, SHANK3, CHK2, GRIN2A, RAG1, GABRG2, and SYNGAP1) in 21 cases, which included two or more unrelated subjects with ultra-rare variants in four genes. These genes converge into two broad functional categories. One regulates peripheral immune responses and microglia (PPM1D, CHK2, NLRC4, RAG1, PLCG2). The other is expressed primarily at neuronal synapses (SHANK3, SYNGAP1, GRIN2A, GABRG2, CACNA1B, SGCE). Mutations in these neuronal genes are also described in autism spectrum disorder and myoclonus-dystonia. In fact, 12/21 cases developed PANS superimposed on a preexisting neurodevelopmental disorder. Genes in both categories are also highly expressed in the enteric nervous system and the choroid plexus. Thus, genetic variation in PANS candidate genes may function by disrupting peripheral and central immune functions, neurotransmission, and/or the blood-CSF/brain barriers following stressors such as infection.

Genetic trade-offs between complex diseases and longevity

Longevity was influenced by many complex diseases and traits. However, the relationships between human longevity and genetic risks of complex diseases were not broadly studied. Here, we constructed polygenic risk scores (PRSs) for 225 complex diseases/traits and evaluated their relationships with human longevity in a cohort with 2178 centenarians and 2299 middle‐aged individuals. Lower genetic risks of stroke and hypotension were observed in centenarians, while higher genetic risks of schizophrenia (SCZ) and type 2 diabetes (T2D) were detected in long‐lived individuals. We further stratified PRSs into cell‐type groups and significance‐level groups. The results showed that the immune component of SCZ genetic risk was positively linked to longevity, and the renal component of T2D genetic risk was the most deleterious. Additionally, SNPs with very small p‐values (p ≤ 1x10^‐5) for SCZ and T2D were negatively correlated with longevity. While for the less significant SNPs (1x10^‐5 < p ≤ 0.05), their effects on disease and longevity were positively correlated. Overall, we identified genetically informed positive and negative factors for human longevity, gained more insights on the accumulation of disease risk alleles during evolution, and provided evidence for the theory of genetic trade‐offs between complex diseases and longevity.

ATXN2 intermediate expansions in amyotrophic lateral sclerosis

Intermediate CAG (polyQ) expansions in the gene ataxin-2 (ATXN2) are now recognized as a risk factor for amyotrophic lateral sclerosis. The threshold for increased risk is not yet firmly established, with reports ranging from 27 to 31 repeats. We investigated the presence of ATXN2 polyQ expansions in 9268 DNA samples collected from people with amyotrophic lateral sclerosis, amyotrophic lateral sclerosis with frontotemporal dementia, frontotemporal dementia alone, Lewy body dementia and age matched controls. This analysis confirmed ATXN2 intermediate polyQ expansions of ≥31 as a risk factor for amyotrophic lateral sclerosis with an odds ratio of 6.31. Expansions were an even greater risk for amyotrophic lateral sclerosis with frontotemporal dementia (odds ratio 27.59) and a somewhat lesser risk for frontotemporal dementia alone (odds ratio 3.14). There was no increased risk for Lewy body dementia. In a subset of 1362 patients with amyotrophic lateral sclerosis with complete clinical data, we could not confirm previous reports of earlier onset of amyotrophic lateral sclerosis or shorter survival in 25 patients with expansions. These new data confirm ≥31 polyQ repeats in ATXN2 increase the risk for amyotrophic lateral sclerosis, and also for the first time show an even greater risk for amyotrophic lateral sclerosis with frontotemporal dementia. The lack of a more aggressive phenotype in amyotrophic lateral sclerosis patients with expansions has implications for ongoing gene-silencing trials for amyotrophic lateral sclerosis.

Genetic associations with healthy ageing among Chinese adults

The genetic basis of overall healthy ageing, especially among the East-Asian population is understudied. We conducted a genome-wide association study among 1618 Singapore Chinese elderly participants (65 years or older) ascertained to have aged healthily and compared their genome-wide genotypes to 6221 participants who did not age healthily, after a 20-year follow-up. Two genetic variants were identified (PMeta < 2.59 × 10-8) to be associated with healthy aging, including the LRP1B locus previously associated in long-lived individuals without cognitive decline. Our study sheds additional insights on the genetic basis of healthy ageing.

The 90 plus: longevity and COVID-19 survival

The world population is getting older and studies aiming to enhance our comprehension of the underlying mechanisms responsible for health span are of utmost interest for longevity and as a measure for health care. In this review, we summarized previous genetic association studies (GWAS) and next-generation sequencing (NGS) of elderly cohorts. We also present the updated hypothesis for the aging process, together with the factors associated with healthy aging. We discuss the relevance of studying older individuals and build databanks to characterize the presence and resistance against late-onset disorders. The identification of about 2 million novel variants in our cohort of more than 1000 elderly Brazilians illustrates the importance of studying highly admixed populations of non-European ancestry. Finally, the ascertainment of nonagenarians and particularly of centenarians who were recovered from COVID-19 or remained asymptomatic opens new avenues of research aiming to enhance our comprehension of biological mechanisms associated with resistance against pathogens.

Identification of a novel variant of the ciliopathic gene FUZZY associated with craniosynostosis

Craniosynostosis is a birth defect occurring in approximately one in 2000 live births, where premature fusion of the cranial bones inhibits growth of the skull during critical periods of brain development. The resulting changes in skull shape can lead to compression of the brain, causing severe complications. While we have some understanding of the molecular pathology of craniosynostosis, a large proportion of cases are of unknown genetic aetiology. Based on studies in mouse, we previously proposed that the ciliopathy gene Fuz should be considered a candidate craniosynostosis gene. Here, we report a novel variant of FUZ (c.851 G > C, p.(Arg284Pro)) found in monozygotic twins presenting with craniosynostosis. To investigate whether Fuz has a direct role in regulating osteogenic fate and mineralisation, we cultured primary osteoblasts and mouse embryonic fibroblasts (MEFs) from Fuz mutant mice. Loss of Fuz resulted in increased osteoblastic mineralisation. This suggests that FUZ protein normally acts as a negative regulator of osteogenesis. We then used Fuz mutant MEFs, which lose functional primary cilia, to test whether the FUZ p.(Arg284Pro) variant could restore FUZ function during ciliogenesis. We found that expression of the FUZ p.(Arg284Pro) variant was sufficient to partially restore cilia numbers, but did not mediate a comparable response to Hedgehog pathway activation. Together, this suggests the osteogenic effects of FUZ p.(Arg284Pro) do not depend upon initiation of ciliogenesis.

Different phases of aging in mouse old skeletal muscle

With a graying population and increasing longevity, it is essential to identify life transition in later years and discern heterogeneity among older people. Subclassifying the elderly population to inspect the subdivisions for pathophysiological differences is particularly important for the investigation of age-related illnesses. For this purpose, using 24- and 28-month-old mice to represent the "young-old" and "old-old", respectively, we compared their skeletal muscle transcriptomes and found each in a distinct stage: early/gradual (E-aging) and late/accelerated aging phase (L-aging). Principal component analysis showed that the old-old transcriptomes were largely disengaged from the forward transcriptomic trajectory generated in the younger-aged group, indicating a substantial change in gene expression profiles during L-aging. By calculating the transcriptomic distance, it was found that the 28-month group was closer to the two-month group than to the 24-month group. The divergence rate per month for the transcriptomes was the highest in L-aging, twice as fast as the rate in E-aging. Indeed, many of the L-aging genes were significantly altered in transcription, although the changes did not seem random but rather coordinated in a variety of functional gene sets. Of 2,707 genes transcriptionally altered during E-aging, two-thirds were also significantly changed during L-aging, to either downturning or upturning way. The downturn genes were related to mitochondrial function and translational gene sets, while the upturn genes were linked to inflammation-associated gene sets. Our results provide a transcriptomic muscle signature that distinguishes old-old mice from young-old mice. This can help to methodically examine muscle disorders in the elderly.

Co-Inference of Data Mislabelings Reveals Improved Models in Genomics and Breast Cancer Diagnostics

Mislabeling of cases as well as controls in case–control studies is a frequent source of strong bias in prognostic and diagnostic tests and algorithms. Common data processing methods available to the researchers in the biomedical community do not allow for consistent and robust treatment of labeled data in the situations where both, the case and the control groups, contain a non-negligible proportion of mislabeled data instances. This is an especially prominent issue in studies regarding late-onset conditions, where individuals who may convert to cases may populate the control group, and for screening studies that often have high false-positive/-negative rates. To address this problem, we propose a method for a simultaneous robust inference of Lasso reduced discriminative models and of latent group-specific mislabeling risks, not requiring any exactly labeled data. We apply it to a standard breast cancer imaging dataset and infer the mislabeling probabilities (being rates of false-negative and false-positive core-needle biopsies) together with a small set of simple diagnostic rules, outperforming the state-of-the-art BI-RADS diagnostics on these data. The inferred mislabeling rates for breast cancer biopsies agree with the published purely empirical studies. Applying the method to human genomic data from a healthy-ageing cohort reveals a previously unreported compact combination of single-nucleotide polymorphisms that are strongly associated with a healthy-ageing phenotype for Caucasians. It determines that 7.5% of Caucasians in the 1000 Genomes dataset (selected as a control group) carry a pattern characteristic of healthy ageing.

Incomplete Penetrance and Variable Expressivity: From Clinical Studies to Population Cohorts

The same genetic variant found in different individuals can cause a range of diverse phenotypes, from no discernible clinical phenotype to severe disease, even among related individuals. Such variants can be said to display incomplete penetrance, a binary phenomenon where the genotype either causes the expected clinical phenotype or it does not, or they can be said to display variable expressivity, in which the same genotype can cause a wide range of clinical symptoms across a spectrum. Both incomplete penetrance and variable expressivity are thought to be caused by a range of factors, including common variants, variants in regulatory regions, epigenetics, environmental factors, and lifestyle. Many thousands of genetic variants have been identified as the cause of monogenic disorders, mostly determined through small clinical studies, and thus, the penetrance and expressivity of these variants may be overestimated when compared to their effect on the general population. With the wealth of population cohort data currently available, the penetrance and expressivity of such genetic variants can be investigated across a much wider contingent, potentially helping to reclassify variants that were previously thought to be completely penetrant. Research into the penetrance and expressivity of such genetic variants is important for clinical classification, both for determining causative mechanisms of disease in the affected population and for providing accurate risk information through genetic counseling. A genotype-based definition of the causes of rare diseases incorporating information from population cohorts and clinical studies is critical for our understanding of incomplete penetrance and variable expressivity. This review examines our current knowledge of the penetrance and expressivity of genetic variants in rare disease and across populations, as well as looking into the potential causes of the variation seen, including genetic modifiers, mosaicism, and polygenic factors, among others. We also considered the challenges that come with investigating penetrance and expressivity.

Whole genome sequencing and in vitro splice assays reveal genetic causes for inherited retinal diseases

Inherited retinal diseases (IRDs) are a major cause of visual impairment. These clinically heterogeneous disorders are caused by pathogenic variants in more than 270 genes. As 30-40% of cases remain genetically unexplained following conventional genetic testing, we aimed to obtain a genetic diagnosis in an IRD cohort in which the genetic cause was not found using whole-exome sequencing or targeted capture sequencing. We performed whole-genome sequencing (WGS) to identify causative variants in 100 unresolved cases. After initial prioritization, we performed an in-depth interrogation of all noncoding and structural variants in genes when one candidate variant was detected. In addition, functional analysis of putative splice-altering variants was performed using in vitro splice assays. We identified the genetic cause of the disease in 24 patients. Causative coding variants were observed in genes such as ATXN7, CEP78, EYS, FAM161A, and HGSNAT. Gene disrupting structural variants were also detected in ATXN7, PRPF31, and RPGRIP1. In 14 monoallelic cases, we prioritized candidate noncanonical splice sites or deep-intronic variants that were predicted to disrupt the splicing process based on in silico analyses. Of these, seven cases were resolved as they carried pathogenic splice defects. WGS is a powerful tool to identify causative variants residing outside coding regions or heterozygous structural variants. This approach was most efficient in cases with a distinct clinical diagnosis. In addition, in vitro splice assays provide important evidence of the pathogenicity of rare variants.

Long-read technologies identify a hidden inverted duplication in a family with choroideremia

The lack of molecular diagnoses in rare genetic diseases can be explained by limitations of current standard genomic technologies. Upcoming long-read techniques have complementary strengths to overcome these limitations, with a particular strength in identifying structural variants. By using optical genome mapping and long-read sequencing, we aimed to identify the pathogenic variant in a large family with X-linked choroideremia. In this family, aberrant splicing of exon 12 of the choroideremia gene CHM was detected in 2003, but the underlying genomic defect remained elusive. Optical genome mapping and long-read sequencing approaches now revealed an intragenic 1,752 bp inverted duplication including exon 12 and surrounding regions, located downstream of the wild-type copy of exon 12. Both breakpoint junctions were confirmed with Sanger sequencing and segregate with the X-linked inheritance in the family. The breakpoint junctions displayed sequence microhomology suggestive for an erroneous replication mechanism as the origin of the structural variant. The inverted duplication is predicted to result in a hairpin formation of the pre-mRNA with the wild-type exon 12, leading to exon skipping in the mature mRNA. The identified inverted duplication is deemed the hidden pathogenic cause of disease in this family. Our study shows that optical genome mapping and long-read sequencing have significant potential for the identification of (hidden) structural variants in rare genetic diseases.

Rare genetic coding variants associated with human longevity and protection against age-related diseases

Extreme longevity in humans has a strong genetic component, but whether this involves genetic variation in the same longevity pathways as found in model organisms is unclear. Using whole-exome sequences of a large cohort of Ashkenazi Jewish centenarians to examine enrichment for rare coding variants, we found most longevity-associated rare coding variants converge upon conserved insulin/insulin-like growth factor 1 signaling and AMP-activating protein kinase signaling pathways. Centenarians have a number of pathogenic rare coding variants similar to control individuals, suggesting that rare variants detected in the conserved longevity pathways are protective against age-related pathology. Indeed, we detected a pro-longevity effect of rare coding variants in the Wnt signaling pathway on individuals harboring the known common risk allele APOE4. The genetic component of extreme human longevity constitutes, at least in part, rare coding variants in pathways that protect against aging, including those that control longevity in model organisms.

The Ageing Brain: Molecular and Cellular Basis of Neurodegeneration

Ageing is an inevitable event in the lifecycle of all organisms, characterized by progressive physiological deterioration and increased vulnerability to death. Ageing has also been described as the primary risk factor of most neurodegenerative diseases, including Alzheimer's disease (AD), Parkinson's disease (PD), Huntington's disease (HD), and frontotemporal lobar dementia (FTD). These neurodegenerative diseases occur more prevalently in the aged populations. Few effective treatments have been identified to treat these epidemic neurological crises. Neurodegenerative diseases are associated with enormous socioeconomic and personal costs. Here, the pathogenesis of AD, PD, and other neurodegenerative diseases has been presented, including a summary of their known associations with the biological hallmarks of ageing: genomic instability, telomere attrition, epigenetic alterations, loss of proteostasis, mitochondrial dysfunction, cellular senescence, deregulated nutrient sensing, stem cell exhaustion, and altered intercellular communications. Understanding the central biological mechanisms that underlie ageing is important for identifying novel therapeutic targets for neurodegenerative diseases. Potential therapeutic strategies, including the use of NAD+ precursors, mitophagy inducers, and inhibitors of cellular senescence, has also been discussed.

Analysis of high-risk pedigrees identifies 11 candidate variants for Alzheimer's disease

Introduction

Analysis of sequence data in high‐risk pedigrees is a powerful approach to detect rare predisposition variants.

Methods

Rare, shared candidate predisposition variants were identified from exome sequencing 19 Alzheimer's disease (AD)‐affected cousin pairs selected from high‐risk pedigrees. Variants were further prioritized by risk association in various external datasets. Candidate variants emerging from these analyses were tested for co‐segregation to additional affected relatives of the original sequenced pedigree members.

Results

AD‐affected high‐risk cousin pairs contained 564 shared rare variants. Eleven variants spanning 10 genes were prioritized in external datasets: rs201665195 (ABCA7), and rs28933981 (TTR) were previously implicated in AD pathology; rs141402160 (NOTCH3) and rs140914494 (NOTCH3) were previously reported; rs200290640 (PIDD1) and rs199752248 (PIDD1) were present in more than one cousin pair; rs61729902 (SNAP91), rs140129800 (COX6A2, AC026471), and rs191804178 (MUC16) were not present in a longevity cohort; and rs148294193 (PELI3) and rs147599881 (FCHO1) approached significance from analysis of AD‐related phenotypes. Three variants were validated via evidence of co‐segregation to additional relatives (PELI3, ABCA7, and SNAP91).

Discussion

These analyses support ABCA7 and TTR as AD risk genes, expand on previously reported NOTCH3 variant identification, and prioritize seven additional candidate variants.

Footprints in the Sand: Deep Taxonomic Comparisons in Vertebrate Genomics to Unveil the Genetic Programs of Human Longevity

With the modern quality, quantity, and availability of genomic sequencing across species, as well as across the expanse of human populations, we can screen for shared signatures underlying longevity and lifespan. Knowledge of these mechanisms would be medically invaluable in combating aging and age-related diseases. The diversity of longevities across vertebrates is an opportunity to look for patterns of genetic variation that may signal how this life history property is regulated, and ultimately how it can be modulated. Variation in human longevity provides a unique window to look for cases of extreme lifespan within a population, as well as associations across populations for factors that influence capacity to live longer. Current large cohort studies support the use of population level analyses to identify key factors associating with human lifespan. These studies are powerful in concept, but have demonstrated limited ability to resolve signals from background variation. In parallel, the expanding catalog of sequencing and annotation from diverse species, some of which have evolved longevities well past a human lifespan, provides independent cases to look at the genomic signatures of longevity. Recent comparative genomic work has shown promise in finding shared mechanisms associating with longevity among distantly related vertebrate groups. Given the genetic constraints between vertebrates, we posit that a combination of approaches, of parallel meta-analysis of human longevity along with refined analysis of other vertebrate clades having exceptional longevity, will aid in resolving key regulators of enhanced lifespan that have proven to be elusive when analyzed in isolation.

Whole-genome sequencing analysis of semi-supercentenarians

Extreme longevity is the paradigm of healthy aging as individuals who reached the extreme decades of human life avoided or largely postponed all major age-related diseases. In this study, we sequenced at high coverage (90X) the whole genome of 81 semi-supercentenarians and supercentenarians [105+/110+] (mean age: 106.6 ± 1.6) and of 36 healthy unrelated geographically matched controls (mean age 68.0 ± 5.9) recruited in Italy. The results showed that 105+/110+ are characterized by a peculiar genetic background associated with efficient DNA repair mechanisms, as evidenced by both germline data (common and rare variants) and somatic mutations patterns (lower mutation load if compared to younger healthy controls). Results were replicated in a second independent cohort of 333 Italian centenarians and 358 geographically matched controls. The genetics of 105+/110+ identified DNA repair and clonal haematopoiesis as crucial players for healthy aging and for the protection from cardiovascular events.

Analysis of Brugada syndrome loci reveals that fine-mapping clustered GWAS hits enhances the annotation of disease-relevant variants

Genome-wide association studies (GWASs) are instrumental in identifying loci harboring common single-nucleotide variants (SNVs) that affect human traits and diseases. GWAS hits emerge in clusters, but the focus is often on the most significant hit in each trait- or disease-associated locus. The remaining hits represent SNVs in linkage disequilibrium (LD) and are considered redundant and thus frequently marginally reported or exploited. Here, we interrogate the value of integrating the full set of GWAS hits in a locus repeatedly associated with cardiac conduction traits and arrhythmia, SCN5A-SCN10A. Our analysis reveals 5 common 7-SNV haplotypes (Hap1-5) with 2 combinations associated with life-threatening arrhythmia-Brugada syndrome (the risk Hap1/1 and protective Hap2/3 genotypes). Hap1 and Hap2 share 3 SNVs; thus, this analysis suggests that assuming redundancy among clustered GWAS hits can lead to confounding disease-risk associations and supports the need to deconstruct GWAS data in the context of haplotype composition.

The Impact of Modern Technologies on Molecular Diagnostic Success Rates, with a Focus on Inherited Retinal Dystrophy and Hearing Loss

The identification of pathogenic variants in monogenic diseases has been of interest to researchers and clinicians for several decades. However, for inherited diseases with extremely high genetic heterogeneity, such as hearing loss and retinal dystrophies, establishing a molecular diagnosis requires an enormous effort. In this review, we use these two genetic conditions as examples to describe the initial molecular genetic identification approaches, as performed since the early 90s, and subsequent improvements and refinements introduced over the years. Next, the history of DNA sequencing from conventional Sanger sequencing to high-throughput massive parallel sequencing, a.k.a. next-generation sequencing, is outlined, including their advantages and limitations and their impact on identifying the remaining genetic defects. Moreover, the development of recent technologies, also coined "third-generation" sequencing, is reviewed, which holds the promise to overcome these limitations. Furthermore, we outline the importance and complexity of variant interpretation in clinical diagnostic settings concerning the massive number of different variants identified by these methods. Finally, we briefly mention the development of novel approaches such as optical mapping and multiomics, which can help to further identify genetic defects in the near future.

Integrated genetic analyses revealed novel human longevity loci and reduced risks of multiple diseases in a cohort study of 15,651 Chinese individuals

There is growing interest in studying the genetic contributions to longevity, but limited relevant genes have been identified. In this study, we performed a genetic association study of longevity in a total of 15,651 Chinese individuals. Novel longevity loci, BMPER (rs17169634; p = 7.91 × 10^-15 ) and TMEM43/XPC (rs1043943; p = 3.59 × 10^-8 ), were identified in a case-control analysis of 11,045 individuals. BRAF (rs1267601; p = 8.33 × 10^-15 ) and BMPER (rs17169634; p = 1.45 × 10^-10 ) were significantly associated with life expectancy in 12,664 individuals who had survival status records. Additional sex-stratified analyses identified sex-specific longevity genes. Notably, sex-differential associations were identified in two linkage disequilibrium blocks in the TOMM40/APOE region, indicating potential differences during meiosis between males and females. Moreover, polygenic risk scores and Mendelian randomization analyses revealed that longevity was genetically causally correlated with reduced risks of multiple diseases, such as type 2 diabetes, cardiovascular diseases, and arthritis. Finally, we incorporated genetic markers, disease status, and lifestyles to classify longevity or not-longevity groups and predict life span. Our predictive models showed good performance (AUC = 0.86 for longevity classification and explained 19.8% variance of life span) and presented a greater predictive efficiency in females than in males. Taken together, our findings not only shed light on the genetic contributions to longevity but also elucidate correlations between diseases and longevity.

Genome sequencing analysis identifies new loci associated with Lewy body dementia and provides insights into its genetic architecture

The genetic basis of Lewy body dementia (LBD) is not well understood. Here, we performed whole-genome sequencing in large cohorts of LBD cases and neurologically healthy controls to study the genetic architecture of this understudied form of dementia, and to generate a resource for the scientific community. Genome-wide association analysis identified five independent risk loci, whereas genome-wide gene-aggregation tests implicated mutations in the gene GBA. Genetic risk scores demonstrate that LBD shares risk profiles and pathways with Alzheimer's disease and Parkinson's disease, providing a deeper molecular understanding of the complex genetic architecture of this age-related neurodegenerative condition.

Reliability of genomic variants across different next-generation sequencing platforms and bioinformatic processing pipelines

Background

Next Generation Sequencing (NGS) is the fundament of various studies, providing insights into questions from biology and medicine. Nevertheless, integrating data from different experimental backgrounds can introduce strong biases. In order to methodically investigate the magnitude of systematic errors in single nucleotide variant calls, we performed a cross-sectional observational study on a genomic cohort of 99 subjects each sequenced via (i) Illumina HiSeq X, (ii) Illumina HiSeq, and (iii) Complete Genomics and processed with the respective bioinformatic pipeline. We also repeated variant calling for the Illumina cohorts with GATK, which allowed us to investigate the effect of the bioinformatics analysis strategy separately from the sequencing platform’s impact.

Results

The number of detected variants/variant classes per individual was highly dependent on the experimental setup. We observed a statistically significant overrepresentation of variants uniquely called by a single setup, indicating potential systematic biases. Insertion/deletion polymorphisms (indels) were associated with decreased concordance compared to single nucleotide polymorphisms (SNPs). The discrepancies in indel absolute numbers were particularly prominent in introns, Alu elements, simple repeats, and regions with medium GC content. Notably, reprocessing sequencing data following the best practice recommendations of GATK considerably improved concordance between the respective setups.

Conclusion

We provide empirical evidence of systematic heterogeneity in variant calls between alternative experimental and data analysis setups. Furthermore, our results demonstrate the benefit of reprocessing genomic data with harmonized pipelines when integrating data from different studies.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-020-07362-8.

Precision medicine for mood disorders: objective assessment, risk prediction, pharmacogenomics, and repurposed drugs

Mood disorders (depression, bipolar disorders) are prevalent and disabling. They are also highly co-morbid with other psychiatric disorders. Currently there are no objective measures, such as blood tests, used in clinical practice, and available treatments do not work in everybody. The development of blood tests, as well as matching of patients with existing and new treatments, in a precise, personalized and preventive fashion, would make a significant difference at an individual and societal level. Early pilot studies by us to discover blood biomarkers for mood state were promising [1], and validated by others [2]. Recent work by us has identified blood gene expression biomarkers that track suicidality, a tragic behavioral outcome of mood disorders, using powerful longitudinal within-subject designs, validated them in suicide completers, and tested them in independent cohorts for ability to assess state (suicidal ideation), and ability to predict trait (future hospitalizations for suicidality) [3-6]. These studies showed good reproducibility with subsequent independent genetic studies [7]. More recently, we have conducted such studies also for pain [8], for stress disorders [9], and for memory/Alzheimer's Disease [10]. We endeavored to use a similar comprehensive approach to identify more definitive biomarkers for mood disorders, that are transdiagnostic, by studying mood in psychiatric disorders patients. First, we used a longitudinal within-subject design and whole-genome gene expression approach to discover biomarkers which track mood state in subjects who had diametric changes in mood state from low to high, from visit to visit, as measured by a simple visual analog scale that we had previously developed (SMS-7). Second, we prioritized these biomarkers using a convergent functional genomics (CFG) approach encompassing in a comprehensive fashion prior published evidence in the field. Third, we validated the biomarkers in an independent cohort of subjects with clinically severe depression (as measured by Hamilton Depression Scale, (HAMD)) and with clinically severe mania (as measured by the Young Mania Rating Scale (YMRS)). Adding the scores from the first three steps into an overall convergent functional evidence (CFE) score, we ended up with 26 top candidate blood gene expression biomarkers that had a CFE score as good as or better than SLC6A4, an empirical finding which we used as a de facto positive control and cutoff. Notably, there was among them an enrichment in genes involved in circadian mechanisms. We further analyzed the biological pathways and networks for the top candidate biomarkers, showing that circadian, neurotrophic, and cell differentiation functions are involved, along with serotonergic and glutamatergic signaling, supporting a view of mood as reflecting energy, activity and growth. Fourth, we tested in independent cohorts of psychiatric patients the ability of each of these 26 top candidate biomarkers to assess state (mood (SMS-7), depression (HAMD), mania (YMRS)), and to predict clinical course (future hospitalizations for depression, future hospitalizations for mania). We conducted our analyses across all patients, as well as personalized by gender and diagnosis, showing increased accuracy with the personalized approach, particularly in women. Again, using SLC6A4 as the cutoff, twelve top biomarkers had the strongest overall evidence for tracking and predicting depression after all four steps: NRG1, DOCK10, GLS, PRPS1, TMEM161B, GLO1, FANCF, HNRNPDL, CD47, OLFM1, SMAD7, and SLC6A4. Of them, six had the strongest overall evidence for tracking and predicting both depression and mania, hence bipolar mood disorders. There were also two biomarkers (RLP3 and SLC6A4) with the strongest overall evidence for mania. These panels of biomarkers have practical implications for distinguishing between depression and bipolar disorder. Next, we evaluated the evidence for our top biomarkers being targets of existing psychiatric drugs, which permits matching patients to medications in a targeted fashion, and the measuring of response to treatment. We also used the biomarker signatures to bioinformatically identify new/repurposed candidate drugs. Top drugs of interest as potential new antidepressants were pindolol, ciprofibrate, pioglitazone and adiphenine, as well as the natural compounds asiaticoside and chlorogenic acid. The last 3 had also been identified by our previous suicidality studies. Finally, we provide an example of how a report to doctors would look for a patient with depression, based on the panel of top biomarkers (12 for depression and bipolar, one for mania), with an objective depression score, risk for future depression, and risk for bipolar switching, as well as personalized lists of targeted prioritized existing psychiatric medications and new potential medications. Overall, our studies provide objective assessments, targeted therapeutics, and monitoring of response to treatment, that enable precision medicine for mood disorders.

Collagenous Alzheimer amyloid plaque component impacts on the compaction of amyloid-β plaques

Massive deposition of amyloid β peptides (Aβ) as senile plaques (SP) characterizes the brain pathology of Alzheimer’s disease (AD). SPs exhibit a variety of morphologies, although little is known about the SP components that determine their morphology. Collagenous Alzheimer amyloid plaque component (CLAC) is one of the major non-Aβ proteinaceous components of SP amyloid in AD brains. Here we show that overexpression of CLAC precursor (CLAC-P) in the brains of APP transgenic mice results in a significant remodeling of amyloid pathology, i.e., reduction in diffuse-type amyloid plaques and an increase in compact plaques laden with thioflavin S-positive amyloid cores. In vivo microdialysis revealed a significant decrease in Aβ in the brain interstitial fluid of CLAC-P/APP double transgenic mice compared with APP transgenic mice. These findings implicate CLAC in the compaction of Aβ in amyloid plaques and the brain dynamics of Aβ.

Similar burden of pathogenic coding variants in exceptionally long-lived individuals and individuals without exceptional longevity

Centenarians (exceptionally long‐lived individuals—ELLI) are a unique segment of the population, exhibiting long human lifespan and healthspan, despite generally practicing similar lifestyle habits as their peers. We tested disease‐associated mutation burden in ELLI genomes by determining the burden of pathogenic variants reported in the ClinVar and HGMD databases using data from whole exome sequencing (WES) conducted in a cohort of ELLI, their offspring, and control individuals without antecedents of familial longevity (n = 1879), all descendent from the founder population of Ashkenazi Jews. The burden of pathogenic variants did not differ between the three groups. Additional analyses of variants subtypes and variant effect predictor (VEP) biotype frequencies did not reveal a decrease of pathogenic or loss‐of‐function (LoF) variants in ELLI and offspring compared to the control group. Case–control pathogenic variants enrichment analyses conducted in ELLI and controls also did not identify significant differences in any of the variants between the groups and polygenic risk scores failed to provide a predictive model. Interestingly, cancer and Alzheimer's disease‐associated variants were significantly depleted in ELLI compared to controls, suggesting slower accumulation of mutation. That said, polygenic risk score analysis failed to find any predictive variants among the functional variants tested. The high similarity in the burden of pathogenic variation between ELLI and individuals without familial longevity supports the notion that extension of lifespan and healthspan in ELLI is not a consequence of pathogenic variant depletion but rather a result of other genomic, epigenomic, or potentially nongenomic properties.

Identification and genomic analysis of pedigrees with exceptional longevity identifies candidate rare variants

BACKGROUND

Longevity as a phenotype entails living longer than average and typically includes living without chronic age-related diseases. Recently, several common genetic components to longevity have been identified. This study aims to identify additional genetic variants associated with longevity using unique and powerful analyses of pedigrees with a statistical excess of healthy elderly individuals identified in the Utah Population Database (UPDB).

METHODS

From an existing biorepository of Utah pedigrees, six independent cousin pairs were selected from four extended pedigrees that exhibited an excess of healthy elderly individuals; whole exome sequencing (WES) was performed on two elderly individuals from each pedigree who were either first cousins or first cousins once removed. Rare (<.01 population frequency) variants shared by at least one elderly cousin pair in a region likely to be identical by descent were identified as candidates. Ingenuity Variant Analysis was used to prioritize putative causal variants based on quality control, frequency, and gain or loss of function. The variant frequency was compared in healthy cohorts and in an Alzheimer's disease cohort. Remaining variants were filtered based on their presence in genes reported to have an effect on the aging process, aging of cells, or the longevity process. Validation of these candidate variants included tests of segregation on other elderly relatives.

RESULTS

Fifteen rare candidate genetic variants spanning 17 genes shared within cousins were identified as having passed prioritization criteria. Of those variants, six were present in genes that are known or predicted to affect the aging process: rs78408340 (PAM), rs112892337 (ZFAT), rs61737629 (ESPL1), rs141903485 (CEBPE), rs144369314 (UTP4), and rs61753103 (NUP88 and RABEP1). ESPL1 rs61737629 and CEBPE rs141903485 show additional evidence of segregation with longevity in expanded pedigree analyses (p-values = .001 and .0001, respectively).

DISCUSSION

This unique pedigree analysis efficiently identified several novel rare candidate variants that may affect the aging process and added support to seven genes that likely contribute to longevity. Further analyses showed evidence for segregation for two rare variants, ESPL1 rs61737629 and CEBPE rs141903485, in the original longevity pedigrees in which they were initially observed. These candidate genes and variants warrant further investigation.

Transfer of a human gene variant associated with exceptional longevity improves cardiac function in obese type 2 diabetic mice through induction of the SDF-1/CXCR4 signalling pathway

AIMS

Homozygosity for a four-missense single-nucleotide polymorphism haplotype of the human BPIFB4 gene is enriched in long-living individuals. Delivery of this longevity-associated variant (LAV) improved revascularisation and reduced endothelial dysfunction and atherosclerosis in mice through a mechanism involving the stromal cell-derived factor-1 (SDF-1). Here, we investigated if delivery of the LAV-BPIFB4 gene may attenuate the progression of diabetic cardiomyopathy.

METHODS AND RESULTS

Compared with age-matched lean controls, diabetic db/db mice showed altered echocardiographic indices of diastolic and systolic function and histological evidence of microvascular rarefaction, lipid accumulation, and fibrosis in the myocardium. All these alterations, as well as endothelial dysfunction, were prevented by systemic LAV-BPIFB4 gene therapy using an adeno-associated viral vector serotype 9 (AAV9). In contrast, AAV9 wild-type-BPIFB4 exerted no benefit. Interestingly, LAV-BPIFB4-treated mice showed increased SDF-1 levels in peripheral blood and myocardium and up-regulation of the cardiac myosin heavy chain isoform alpha, a contractile protein that was reduced in diabetic hearts. SDF-1 up-regulation was instrumental to LAV-BPIFB4-induced benefit as both haemodynamic and structural improvements were inhibited by an orally active antagonist of the SDF-1 CXCR4 receptor.

CONCLUSIONS

In mice with type-2 diabetes, LAV-BPIFB4 gene therapy promotes an advantageous remodelling of the heart, allowing it to better withstand diabetes-induced stress. These results support the viability of transferring healthy characteristics of longevity to attenuate diabetic cardiac disease.

Whole-Exome Sequencing of an Exceptional Longevity Cohort

Centenarians represent a unique cohort to study the genetic basis for longevity and factors determining the risk of neurodegenerative disorders, including Alzheimer's disease (AD). The estimated genetic contribution to longevity is highest in centenarians and super-cententenarians, but few genetic variants have been shown to clearly impact this phenotype. While the genetic risk for AD and other dementias is now well understood, the frequency of known dementia risk variants in centenarians is not fully characterized. To address these questions, we performed whole-exome sequencing on 100 individuals of 98-108 years age in search of genes with large effect sizes towards the exceptional aging phenotype. Overall, we were unable to identify a rare protein-altering variant or individual genes with an increased burden of rare variants associated with exceptional longevity. Gene burden analysis revealed three genes of nominal statistical significance associated with extreme aging, including LYST, MDN1, and RBMXL1. Several genes with variants conferring an increased risk for AD and other dementias were identified, including TREM2, EPHA1, ABCA7, PLD3, MAPT, and NOTCH3. Larger centenarian studies will be required to further elucidate the genetic basis for longevity, and factors conferring protection against age-dependent neurodegenerative syndromes.

TraPS-VarI: Identifying genetic variants altering phosphotyrosine based signalling motifs

Patient stratification and individualized therapeutic strategies rely on the established knowledge of genotype-specific molecular and cellular alterations of biological and therapeutic significance. Whilst almost all approved drugs have been developed based on the Reference Sequence protein database (RefSeq), the latest genome sequencing studies establish the substantial prevalence of non-synonymous genetic mutations in the general population, including stop-insertion and frame shift mutations within the coding regions of membrane proteins. While the availability of individual genotypes are becoming increasingly common, the biological and clinical interpretations of mutations among individual genomes is largely lagging behind. Lately, transmembrane proteins of haematopoietic (myeloid and lymphoid) derived immune cells have attracted much attention as important targets for cancer immunotherapies. As such, the signalling properties of haematological transmembrane receptors rely on the membrane-proximal phosphotyrosine based sequence motifs (TBSMs) such as ITAM (immunoreceptor tyrosine-based activation motif), ITIM (immunoreceptor tyrosine-based inhibition motif) and signal transducer and activator of transcription 3 (STAT3)-recruiting YxxQ motifs. However, mutations that alter the coding regions of transmembrane proteins, resulting in either insertion or deletion of crucial signal modulating TBSMs, remains unknown. To conveniently identify individual cell line-specific or patient-specific membrane protein altering mutations, we present the Transmembrane Protein Sequence Variant Identifier (TraPS-VarI). TraPS-VarI is an annotation tool for accurate mapping of the effect of an individual’s mutation in the transmembrane protein sequence, and to identify the prevalence of TBSMs. TraPS-VarI is a biologist and clinician-friendly algorithm with a web interface and an associated database browser (https://www.traps-vari.org/).

ExpansionHunter Denovo: a computational method for locating known and novel repeat expansions in short-read sequencing data

Repeat expansions are responsible for over 40 monogenic disorders, and undoubtedly more pathogenic repeat expansions remain to be discovered. Existing methods for detecting repeat expansions in short-read sequencing data require predefined repeat catalogs. Recent discoveries emphasize the need for methods that do not require pre-specified candidate repeats. To address this need, we introduce ExpansionHunter Denovo, an efficient catalog-free method for genome-wide repeat expansion detection. Analysis of real and simulated data shows that our method can identify large expansions of 41 out of 44 pathogenic repeats, including nine recently reported non-reference repeat expansions not discoverable via existing methods.

Genome-wide Profiling Identifies DNA Methylation Signatures of Aging in Rod Photoreceptors Associated with Alterations in Energy Metabolism

Aging-associated functional decline is accompanied by alterations in the epigenome. To explore DNA modifications that could influence visual function with age, we perform whole-genome bisulfite sequencing of purified mouse rod photoreceptors at four ages and identify 2,054 differentially methylated regions (DMRs). We detect many DMRs during early stages of aging and in rod regulatory regions, and some of these cluster at chromosomal hotspots, especially on chromosome 10, which includes a longevity interactome. Integration of methylome to age-related transcriptome changes, chromatin signatures, and first-order protein-protein interactions uncover an enrichment of DMRs in altered pathways that are associated with rod function, aging, and energy metabolism. In concordance, we detect reduced basal mitochondrial respiration and increased fatty acid dependency with retinal age in ex vivo assays. Our study reveals age-dependent genomic and chromatin features susceptible to DNA methylation changes in rod photoreceptors and identifies a link between DNA methylation and energy metabolism in aging.

Essentiality-specific pathogenicity prioritization gene score to improve filtering of disease sequence data

The causal genetic variants underlying more than 50% of single gene (monogenic) disorders are yet to be discovered. Many patients with conditions likely to have a monogenic basis do not receive a confirmed molecular diagnosis which has potential impacts on clinical management. We have developed a gene-specific score, essentiality-specific pathogenicity prioritization (ESPP), to guide the recognition of genes likely to underlie monogenic disease variation to assist in filtering of genome sequence data. When a patient genome is sequenced, there are frequently several plausibly pathogenic variants identified in different genes. Recognition of the single gene most likely to include pathogenic variation can guide the identification of a causal variant. The ESPP score integrates gene-level scores which are broadly related to gene essentiality. Previous work towards the recognition of monogenic disease genes proposed a model with increasing gene essentiality from 'non-essential' to 'essential' genes (for which pathogenic variation may be incompatible with survival) with genes liable to contain disease variation positioned between these two extremes. We demonstrate that the ESPP score is useful for recognizing genes with high potential for pathogenic disease-related variation. Genes classed as essential have particularly high scores, as do genes recently recognized as strong candidates for developmental disorders. Through the integration of individual gene-specific scores, which have different properties and assumptions, we demonstrate the utility of an essentiality-based gene score to improve sequence genome filtering.

Variable expression of subclinical phenotypes instead of reduced penetrance in families with mild triphalangeal thumb phenotypes

Background

The of zone of polarizing activity regulatory sequence (ZRS) is a regulatory element residing in intron 5 of LMBR1 and regulates Sonic Hedgehog expression in the limb bud. Variants in the ZRS are generally fully penetrant and can cause triphalangeal thumb (TPT) and polydactyly in affected families.

Objective

In this report, we describe two families with mild phenotypical presentation.

Methods

We performed a field study for clinical evaluation and sequenced the ZRS for variantsusing Sanger sequencing.

Results

In family I, a novel 165A>G variant in the ZRS (g.156584405A>G, GRCh37/Hg19) was found. In family II, we identified a 295T>C variant in the ZRS (g.156584535T>C, GRCh37/Hg19). Family members of both families who were presumed to be unaffected shared the variant in the ZRS with affected family members, suggesting reduced penetrance of the genotype. However, clinical examination of these unaffected family members revealed minor anomalies like broad thumbs and lack of thumb opposition. As the phenotype in affected patients is remarkably mild, we suggest that these ZRS variants are minimally disruptive for Sonic Hedgehog expression and therefore can result in subclinical phenotypes.

Conclusion

Our study underlines the importance of accurate clinical examination and appropriate genetic counselling in families with mild cases of TPT.

A novel nicastrin mutation in a three-generation Dutch family with hidradenitis suppurativa: a search for functional significance

Background

Mutations in the γ‐secretase enzyme subunits have been described in multiple kindreds with familial hidradenitis suppurativa (HS).

Objective

In this study, we report a novel nicastrin (NCSTN) mutation causing HS in a Dutch family. We sought to explore the immunobiological function of NCSTN mutations using data of the Immunological Genome Project.

Methods

Blood samples of three affected and two unaffected family members were collected. Whole‐genome sequencing was performed using genomic DNA isolated from peripheral blood leucocytes. Sanger sequencing was done to confirm the causative NCSTN variant and the familial segregation. The microarray data set of the Immunological Genome Project was used for thorough dissection of the expression and function of wildtype NCSTN in the immune system.

Results

In a family consisting of 23 members, we found an autosomal dominant inheritance pattern of HS and detected a novel splice site mutation (c.1912_1915delCAGT) in the NCSTN gene resulting in a frameshift and subsequent premature stop. All affected individuals had HS lesions on non‐flexural and atypical locations. Wildtype NCSTN appears to be upregulated in myeloid cells like monocytes and macrophages, and in mesenchymal cells such as fibroblastic reticular cells and fibroblasts. In addition, within the 25 highest co‐expressed genes with NCSTN we identified CAPNS1,ARNT and PPARD.

Conclusion

This study reports the identification a novel NCSTN gene splice site mutation which causes familial HS. The associated immunobiological functions of NCSTN and its co‐expressed genes ARNT and PPARD link genetics to the most common environmental and metabolic HS risk factors which are smoking and obesity.