Scalable and efficient DNA sequencing analysis on different compute infrastructures aiding variant discovery

DNA variation analysis has become indispensable in many aspects of modern biomedicine, most prominently in the comparison of normal and tumor samples. Thousands of samples are collected in local sequencing efforts and public databases requiring highly scalable, portable, and automated workflows for streamlined processing. Here, we present nf-core/sarek 3, a well-established, comprehensive variant calling and annotation pipeline for germline and somatic samples. It is suitable for any genome with a known reference. We present a full rewrite of the original pipeline showing a significant reduction of storage requirements by using the CRAM format and runtime by increasing intra-sample parallelization. Both are leading to a 70% cost reduction in commercial clouds enabling users to do large-scale and cross-platform data analysis while keeping costs and CO2 emissions low. The code is available at https://nf-co.re/sarek.

The promise of explainable deep learning for omics data analysis: Adding new discovery tools to AI

Deep learning has already revolutionised the way a wide range of data is processed in many areas of daily life. The ability to learn abstractions and relationships from heterogeneous data has provided impressively accurate prediction and classification tools to handle increasingly big datasets. This has a significant impact on the growing wealth of omics datasets, with the unprecedented opportunity for a better understanding of the complexity of living organisms. While this revolution is transforming the way these data are analyzed, explainable deep learning is emerging as an additional tool with the potential to change the way biological data is interpreted. Explainability addresses critical issues such as transparency, so important when computational tools are introduced especially in clinical environments. Moreover, it empowers artificial intelligence with the capability to provide new insights into the input data, thus adding an element of discovery to these already powerful resources. In this review, we provide an overview of the transformative effects explainable deep learning is having on multiple sectors, ranging from genome engineering and genomics, from radiomics to drug design and clinical trials. We offer a perspective to life scientists, to better understand the potential of these tools, and a motivation to implement them in their research, by suggesting learning resources they can use to move their first steps in this field.

Robertsonian fusion and centromere repositioning contributed to the formation of satellite-free centromeres during the evolution of zebras

Centromeres are epigenetically specified by the histone H3 variant CENP-A and typically associated to highly repetitive satellite DNA. We previously discovered natural satellite-free neocentromeres in Equus caballus and E. asinus. Here, through ChIP-seq with an anti-CENP-A antibody, we found an extraordinarily high number of centromeres lacking satellite DNA in the zebras E. burchelli (15 of 22) and E. grevyi (13 of 23), demonstrating that the absence of satellite DNA at the majority of centromeres is compatible with genome stability and species survival and challenging the role of satellite DNA in centromere function. Nine satellite-free centromeres are shared between the two species in agreement with their recent separation. We assembled all centromeric regions and improved the reference genome of E. burchelli. Sequence analysis of the CENP-A binding domains revealed that they are LINE-1 and AT-rich with four of them showing DNA amplification. In the two zebras, satellite-free centromeres emerged from centromere repositioning or following Robertsonian fusion. In five chromosomes, the centromeric function arose near the fusion points, which are located within regions marked by traces of ancestral pericentromeric sequences. Therefore, besides centromere repositioning, Robertsonian fusions are an important source of satellite-free centromeres during evolution. Finally, in one case, a satellite-free centromere was seeded on an inversion breakpoint. At eleven chromosomes, whose primary constrictions seemed to be associated to satellite repeats by cytogenetic analysis, satellite-free neocentromeres were instead located near the ancestral inactivated satellite-based centromeres, therefore, the centromeric function has shifted away from a satellite repeat containing locus to a satellite-free new position.

Rare coding variants in ten genes confer substantial risk for schizophrenia

Rare coding variation has historically provided the most direct connections between gene function and disease pathogenesis. By meta-analysing the whole exomes of 24,248 schizophrenia cases and 97,322 controls, we implicate ultra-rare coding variants (URVs) in 10 genes as conferring substantial risk for schizophrenia (odds ratios of 3-50, P < 2.14 × 10-6) and 32 genes at a false discovery rate of <5%. These genes have the greatest expression in central nervous system neurons and have diverse molecular functions that include the formation, structure and function of the synapse. The associations of the NMDA (N-methyl-D-aspartate) receptor subunit GRIN2A and AMPA (α-amino-3-hydroxy-5-methyl-4-isoxazole propionic acid) receptor subunit GRIA3 provide support for dysfunction of the glutamatergic system as a mechanistic hypothesis in the pathogenesis of schizophrenia. We observe an overlap of rare variant risk among schizophrenia, autism spectrum disorders1, epilepsy and severe neurodevelopmental disorders2, although different mutation types are implicated in some shared genes. Most genes described here, however, are not implicated in neurodevelopment. We demonstrate that genes prioritized from common variant analyses of schizophrenia are enriched in rare variant risk3, suggesting that common and rare genetic risk factors converge at least partially on the same underlying pathogenic biological processes. Even after excluding significantly associated genes, schizophrenia cases still carry a substantial excess of URVs, which indicates that more risk genes await discovery using this approach.

Omics in a Digital World: The Role of Bioinformatics in Providing New Insights Into Human Aging

Background

Aging is a complex phenotype influenced by a combination of genetic and environmental factors. Although many studies addressed its cellular and physiological age-related changes, the molecular causes of aging remain undetermined. Considering the biological complexity and heterogeneity of the aging process, it is now clear that full understanding of mechanisms underlying aging can only be achieved through the integration of different data types and sources, and with new computational methods capable to achieve such integration.

Recent Advances

In this review, we show that an omics vision of the age-dependent changes occurring as the individual ages can provide researchers with new opportunities to understand the mechanisms of aging. Combining results from single-cell analysis with systems biology tools would allow building interaction networks and investigate how these networks are perturbed during aging and disease. The development of high-throughput technologies such as next-generation sequencing, proteomics, metabolomics, able to investigate different biological markers and to monitor them simultaneously during the aging process with high accuracy and specificity, represents a unique opportunity offered to biogerontologists today.

Critical Issues

Although the capacity to produce big data drastically increased over the years, integration, interpretation and sharing of high-throughput data remain major challenges. In this paper we present a survey of the emerging omics approaches in aging research and provide a large collection of datasets and databases as a useful resource for the scientific community to identify causes of aging. We discuss their peculiarities, emphasizing the need for the development of methods focused on the integration of different data types.

Future Directions

We critically review the contribution of bioinformatics into the omics of aging research, and we propose a few recommendations to boost collaborations and produce new insights. We believe that significant advancements can be achieved by following major developments in bioinformatics, investing in diversity, data sharing and community-driven portable bioinformatics methods. We also argue in favor of more engagement and participation, and we highlight the benefits of new collaborations along these lines. This review aims at being a useful resource for many researchers in the field, and a call for new partnerships in aging research.

Autism spectrum disorder and attention deficit hyperactivity disorder have a similar burden of rare protein-truncating variants

We analyze the exome sequences of approximately 8,000 children with autism spectrum disorder (ASD) and/or attention-deficit/hyperactivity disorder (ADHD) and 5,000 controls, and we find that ASD and ADHD have a similar burden of rare protein-truncating variants in evolutionarily constrained genes, both significantly higher than controls. This motivates a combined analysis across ASD and ADHD, which identifies microtubule-associated protein 1A (MAP1A) as a novel exome-wide significant gene conferring risk for childhood psychiatric disorders.

Quantifying the Impact of Rare and Ultra-rare Coding Variation across the Phenotypic Spectrum

There is a limited understanding about the impact of rare protein-truncating variants across multiple phenotypes. We explore the impact of this class of variants on 13 quantitative traits and 10 diseases using whole-exome sequencing data from 100,296 individuals. Protein-truncating variants in genes intolerant to this class of mutations increased risk of autism, schizophrenia, bipolar disorder, intellectual disability, and ADHD. In individuals without these disorders, there was an association with shorter height, lower education, increased hospitalization, and reduced age at enrollment. Gene sets implicated from GWASs did not show a significant protein-truncating variants burden beyond what was captured by established Mendelian genes. In conclusion, we provide a thorough investigation of the impact of rare deleterious coding variants on complex traits, suggesting widespread pleiotropic risk.

Sequencing and de novo assembly of 150 genomes from Denmark as a population reference

Hundreds of thousands of human genomes are now being sequenced to characterize genetic variation and use this information to augment association mapping studies of complex disorders and other phenotypic traits. Genetic variation is identified mainly by mapping short reads to the reference genome or by performing local assembly. However, these approaches are biased against discovery of structural variants and variation in the more complex parts of the genome. Hence, large-scale de novo assembly is needed. Here we show that it is possible to construct excellent de novo assemblies from high-coverage sequencing with mate-pair libraries extending up to 20 kilobases. We report de novo assemblies of 150 individuals (50 trios) from the GenomeDenmark project. The quality of these assemblies is similar to those obtained using the more expensive long-read technology. We use the assemblies to identify a rich set of structural variants including many novel insertions and demonstrate how this variant catalogue enables further deciphering of known association mapping signals. We leverage the assemblies to provide 100 completely resolved major histocompatibility complex haplotypes and to resolve major parts of the Y chromosome. Our study provides a regional reference genome that we expect will improve the power of future association mapping studies and hence pave the way for precision medicine initiatives, which now are being launched in many countries including Denmark.

Whole-exome sequencing of individuals from an isolated population implicates rare risk variants in bipolar disorder

Bipolar disorder affects about 1% of the world's population, and its estimated heritability is about 75%. Only few whole genome or whole-exome sequencing studies in bipolar disorder have been reported, and no rare coding variants have yet been robustly identified. The use of isolated populations might help finding variants with a recent origin, more likely to have drifted to higher frequency by chance. Following this approach, we investigated 28 bipolar cases and 214 controls from the Faroe Islands by whole exome sequencing, and the results were followed-up in a British sample of 2025 cases and 1358 controls. Seventeen variants in 16 genes in the single-variant analysis, and 3 genes in the gene-based statistics surpassed exome-wide significance in the discovery phase. The discovery findings were supported by enrichment analysis of common variants from genome-wide association studies (GWAS) data and interrogation of protein-protein interaction networks. The replication in the British sample confirmed the association with NOS1 (missense variant rs79487279) and NCL (gene-based test). A number of variants from the discovery set were not present in the replication sample, including a novel PITPNM2 missense variant, which is located in a highly significant schizophrenia GWAS locus. Likewise, PIK3C2A identified in the gene-based analysis is located in a combined bipolar and schizophrenia GWAS locus. Our results show support both for existing findings in the literature, as well as for new risk genes, and identify rare variants that might provide additional information on the underlying biology of bipolar disorder.

Whole-exome sequencing implicates DGKH as a risk gene for panic disorder in the Faroese population

The demographic history of the isolated population of the Faroe Islands may have induced enrichment of variants rarely seen in outbred European populations, including enrichment of risk variants for panic disorder (PD). PD is a common mental disorder, characterized by recurring and unprovoked panic attacks, and genetic factors have been estimated to explain around 40% of the risk. In this study the potential enrichment of PD risk variants was explored based on whole-exome sequencing of 54 patients with PD and 211 control individuals from the Faroese population. No genome-wide significant associations were found, however several single variants and genes showed strong association with PD, where DGKH was found to be the strongest PD associated gene. Interestingly DGKH has previously demonstrated genome-wide significant association with bipolar disorder as well as evidence of association to other mental disorders. Additionally, we found an enrichment of PD risk variants in the Faroese population; variants with otherwise low frequency in more outbreed European populations. © 2016 Wiley Periodicals, Inc.

High-Quality Exome Sequencing of Whole-Genome Amplified Neonatal Dried Blood Spot DNA

Stored neonatal dried blood spot (DBS) samples from neonatal screening programmes are a valuable diagnostic and research resource. Combined with information from national health registries they can be used in population-based studies of genetic diseases. DNA extracted from neonatal DBSs can be amplified to obtain micrograms of an otherwise limited resource, referred to as whole-genome amplified DNA (wgaDNA). Here we investigate the robustness of exome sequencing of wgaDNA of neonatal DBS samples. We conducted three pilot studies of seven, eight and seven subjects, respectively. For each subject we analysed a neonatal DBS sample and corresponding adult whole-blood (WB) reference sample. Different DNA sample types were prepared for each of the subjects. Pilot 1: wgaDNA of 2x3.2mm neonatal DBSs (DBS_2x3.2) and raw DNA extract of the WB reference sample (WB_ref). Pilot 2: DBS_2x3.2, WB_ref and a WB_ref replica sharing DNA extract with the WB_ref sample. Pilot 3: DBS_2x3.2, WB_ref, wgaDNA of 2x1.6 mm neonatal DBSs and wgaDNA of the WB reference sample. Following sequencing and data analysis, we compared pairwise variant calls to obtain a measure of similarity—the concordance rate. Concordance rates were slightly lower when comparing DBS vs WB sample types than for any two WB sample types of the same subject before filtering of the variant calls. The overall concordance rates were dependent on the variant type, with SNPs performing best. Post-filtering, the comparisons of DBS vs WB and WB vs WB sample types yielded similar concordance rates, with values close to 100%. WgaDNA of neonatal DBS samples performs with great accuracy and efficiency in exome sequencing. The wgaDNA performed similarly to matched high-quality reference—whole-blood DNA—based on concordance rates calculated from variant calls. No differences were observed substituting 2x3.2 with 2x1.6 mm discs, allowing for additional reduction of sample material in future projects.

Experimental validation of methods for differential gene expression analysis and sample pooling in RNA-seq

BACKGROUND

Massively parallel cDNA sequencing (RNA-seq) experiments are gradually superseding microarrays in quantitative gene expression profiling. However, many biologists are uncertain about the choice of differentially expressed gene (DEG) analysis methods and the validity of cost-saving sample pooling strategies for their RNA-seq experiments. Hence, we performed experimental validation of DEGs identified by Cuffdiff2, edgeR, DESeq2 and Two-stage Poisson Model (TSPM) in a RNA-seq experiment involving mice amygdalae micro-punches, using high-throughput qPCR on independent biological replicate samples. Moreover, we sequenced RNA-pools and compared their results with sequencing corresponding individual RNA samples.

RESULTS

False-positivity rate of Cuffdiff2 and false-negativity rates of DESeq2 and TSPM were high. Among the four investigated DEG analysis methods, sensitivity and specificity of edgeR was relatively high. We documented the pooling bias and that the DEGs identified in pooled samples suffered low positive predictive values.

CONCLUSIONS

Our results highlighted the need for combined use of more sensitive DEG analysis methods and high-throughput validation of identified DEGs in future RNA-seq experiments. They indicated limited utility of sample pooling strategies for RNA-seq in similar setups and supported increasing the number of biological replicate samples.

STAG3 truncating variant as the cause of primary ovarian insufficiency

Primary ovarian insufficiency (POI) is a distressing cause of infertility in young women. POI is heterogeneous with only a few causative genes having been discovered so far. Our objective was to determine the genetic cause of POI in a consanguineous Lebanese family with two affected sisters presenting with primary amenorrhoea and an absence of any pubertal development. Multipoint parametric linkage analysis was performed. Whole-exome sequencing was done on the proband. Linkage analysis identified a locus on chromosome 7 where exome sequencing successfully identified a homozygous two base pair duplication (c.1947_48dupCT), leading to a truncated protein p.(Y650Sfs*22) in the STAG3 gene, confirming it as the cause of POI in this family. Exome sequencing combined with linkage analyses offers a powerful tool to efficiently find novel genetic causes of rare, heterogeneous disorders, even in small single families. This is only the second report of a STAG3 variant; the first STAG3 variant was recently described in a phenotypically similar family with extreme POI. Identification of an additional family highlights the importance of STAG3 in POI pathogenesis and suggests it should be evaluated in families affected with POI.

The use of whole-exome sequencing to disentangle complex phenotypes

The success of whole-exome sequencing to identify mutations causing single-gene disorders has been well documented. In contrast whole-exome sequencing has so far had limited success in the identification of variants causing more complex phenotypes that seem unlikely to be due to the disruption of a single gene. We describe a family where two male offspring of healthy first cousin parents present a complex phenotype consisting of peripheral neuropathy and bronchiectasis that has not been described previously in the literature. Due to the fact that both children had the same problems in the context of parental consanguinity we hypothesised illness resulted from either X-linked or autosomal recessive inheritance. Through the use of whole-exome sequencing we were able to simplify this complex phenotype and identified a causative mutation (p.R1070*) in the gene periaxin (PRX), a gene previously shown to cause peripheral neuropathy (Dejerine-Sottas syndrome) when this mutation is present. For the bronchiectasis phenotype we were unable to identify a causal single mutation or compound heterozygote, reflecting the heterogeneous nature of this phenotype. In conclusion, in this study we show that whole-exome sequencing has the power to disentangle complex phenotypes through the identification of causative genetic mutations for distinct clinical disorders that were previously masked.

Novel variation and de novo mutation rates in population-wide de novo assembled Danish trios

Building a population-specific catalogue of single nucleotide variants (SNVs), indels and structural variants (SVs) with frequencies, termed a national pan-genome, is critical for further advancing clinical and public health genetics in large cohorts. Here we report a Danish pan-genome obtained from sequencing 10 trios to high depth (50 × ). We report 536k novel SNVs and 283k novel short indels from mapping approaches and develop a population-wide de novo assembly approach to identify 132k novel indels larger than 10 nucleotides with low false discovery rates. We identify a higher proportion of indels and SVs than previous efforts showing the merits of high coverage and de novo assembly approaches. In addition, we use trio information to identify de novo mutations and use a probabilistic method to provide direct estimates of 1.27e−8 and 1.5e−9 per nucleotide per generation for SNVs and indels, respectively.

The generation of a national pan-genome, a population-specific catalogue of genetic variation, may advance the impact of clinical genetics studies. Here the Besenbacher et al. carry out deep sequencing and de novo assembly of 10 parent–child trios to generate a Danish pan-genome that provides insight into structural variation, de novo mutation rates and variant calling.

The kinetochore protein, CENPF, is mutated in human ciliopathy and microcephaly phenotypes

Background

Mutations in microtubule-regulating genes are associated with disorders of neuronal migration and microcephaly. Regulation of centriole length has been shown to underlie the pathogenesis of certain ciliopathy phenotypes. Using a next-generation sequencing approach, we identified mutations in a novel centriolar disease gene in a kindred with an embryonic lethal ciliopathy phenotype and in a patient with primary microcephaly.

Methods and results

Whole exome sequencing data from a non-consanguineous Caucasian kindred exhibiting mid-gestation lethality and ciliopathic malformations revealed two novel non-synonymous variants in CENPF, a microtubule-regulating gene. All four affected fetuses showed segregation for two mutated alleles [IVS5-2A>C, predicted to abolish the consensus splice-acceptor site from exon 6; c.1744G>T, p.E582X]. In a second unrelated patient exhibiting microcephaly, we identified two CENPF mutations [c.1744G>T, p.E582X; c.8692 C>T, p.R2898X] by whole exome sequencing. We found that CENP-F colocalised with Ninein at the subdistal appendages of the mother centriole in mouse inner medullary collecting duct cells. Intraflagellar transport protein-88 (IFT-88) colocalised with CENP-F along the ciliary axonemes of renal epithelial cells in age-matched control human fetuses but did not in truncated cilia of mutant CENPF kidneys. Pairwise co-immunoprecipitation assays of mitotic and serum-starved HEKT293 cells confirmed that IFT88 precipitates with endogenous CENP-F.

Conclusions

Our data identify CENPF as a new centriolar disease gene implicated in severe human ciliopathy and microcephaly related phenotypes. CENP-F has a novel putative function in ciliogenesis and cortical neurogenesis.

Mutations in SNX14 cause a distinctive autosomal-recessive cerebellar ataxia and intellectual disability syndrome

Intellectual disability and cerebellar atrophy occur together in a large number of genetic conditions and are frequently associated with microcephaly and/or epilepsy. Here we report the identification of causal mutations in Sorting Nexin 14 (SNX14) found in seven affected individuals from three unrelated consanguineous families who presented with recessively inherited moderate-severe intellectual disability, cerebellar ataxia, early-onset cerebellar atrophy, sensorineural hearing loss, and the distinctive association of progressively coarsening facial features, relative macrocephaly, and the absence of seizures. We used homozygosity mapping and whole-exome sequencing to identify a homozygous nonsense mutation and an in-frame multiexon deletion in two families. A homozygous splice site mutation was identified by Sanger sequencing of SNX14 in a third family, selected purely by phenotypic similarity. This discovery confirms that these characteristic features represent a distinct and recognizable syndrome. SNX14 encodes a cellular protein containing Phox (PX) and regulator of G protein signaling (RGS) domains. Weighted gene coexpression network analysis predicts that SNX14 is highly coexpressed with genes involved in cellular protein metabolism and vesicle-mediated transport. All three mutations either directly affected the PX domain or diminished SNX14 levels, implicating a loss of normal cellular function. This manifested as increased cytoplasmic vacuolation as observed in cultured fibroblasts. Our findings indicate an essential role for SNX14 in neural development and function, particularly in development and maturation of the cerebellum.

Mutation of SALL2 causes recessive ocular coloboma in humans and mice

Ocular coloboma is a congenital defect resulting from failure of normal closure of the optic fissure during embryonic eye development. This birth defect causes childhood blindness worldwide, yet the genetic etiology is poorly understood. Here, we identified a novel homozygous mutation in the SALL2 gene in members of a consanguineous family affected with non-syndromic ocular coloboma variably affecting the iris and retina. This mutation, c.85G>T, introduces a premature termination codon (p.Glu29*) predicted to truncate the SALL2 protein so that it lacks three clusters of zinc-finger motifs that are essential for DNA-binding activity. This discovery identifies SALL2 as the third member of the Drosophila homeotic Spalt-like family of developmental transcription factor genes implicated in human disease. SALL2 is expressed in the developing human retina at the time of, and subsequent to, optic fissure closure. Analysis of Sall2-deficient mouse embryos revealed delayed apposition of the optic fissure margins and the persistence of an anterior retinal coloboma phenotype after birth. Sall2-deficient embryos displayed correct posterior closure toward the optic nerve head, and upon contact of the fissure margins, dissolution of the basal lamina occurred and PAX2, known to be critical for this process, was expressed normally. Anterior closure was disrupted with the fissure margins failing to meet, or in some cases misaligning leading to a retinal lesion. These observations demonstrate, for the first time, a role for SALL2 in eye morphogenesis and that loss of function of the gene causes ocular coloboma in humans and mice.

Identification and validation of loss of function variants in clinical contexts

The choice of an appropriate variant calling pipeline for exome sequencing data is becoming increasingly more important in translational medicine projects and clinical contexts. Within GOSgene, which facilitates genetic analysis as part of a joint effort of the University College London and the Great Ormond Street Hospital, we aimed to optimize a variant calling pipeline suitable for our clinical context. We implemented the GATK/Queue framework and evaluated the performance of its two callers: the classical UnifiedGenotyper and the new variant discovery tool HaplotypeCaller. We performed an experimental validation of the loss-of-function (LoF) variants called by the two methods using Sequenom technology. UnifiedGenotyper showed a total validation rate of 97.6% for LoF single-nucleotide polymorphisms (SNPs) and 92.0% for insertions or deletions (INDELs), whereas HaplotypeCaller was 91.7% for SNPs and 55.9% for INDELs. We confirm that GATK/Queue is a reliable pipeline in translational medicine and clinical context. We conclude that in our working environment, UnifiedGenotyper is the caller of choice, being an accurate method, with a high validation rate of error-prone calls like LoF variants. We finally highlight the importance of experimental validation, especially for INDELs, as part of a standard pipeline in clinical environments.

ARNT2 mutation causes hypopituitarism, post-natal microcephaly, visual and renal anomalies

We describe a previously unreported syndrome characterized by secondary (post-natal) microcephaly with fronto-temporal lobe hypoplasia, multiple pituitary hormone deficiency, seizures, severe visual impairment and abnormalities of the kidneys and urinary tract in a highly consanguineous family with six affected children. Homozygosity mapping and exome sequencing revealed a novel homozygous frameshift mutation in the basic helix-loop-helix transcription factor gene ARNT2 (c.1373_1374dupTC) in affected individuals. This mutation results in absence of detectable levels of ARNT2 transcript and protein from patient fibroblasts compared with controls, consistent with nonsense-mediated decay of the mutant transcript and loss of ARNT2 function. We also show expression of ARNT2 within the central nervous system, including the hypothalamus, as well as the renal tract during human embryonic development. The progressive neurological abnormalities, congenital hypopituitarism and post-retinal visual pathway dysfunction in affected individuals demonstrates for the first time the essential role of ARNT2 in the development of the hypothalamo-pituitary axis, post-natal brain growth, and visual and renal function in humans.

Characterisation and validation of insertions and deletions in 173 patient exomes

Recent advances in genomics technologies have spurred unprecedented efforts in genome and exome re-sequencing aiming to unravel the genetic component of rare and complex disorders. While in rare disorders this allowed the identification of novel causal genes, the missing heritability paradox in complex diseases remains so far elusive. Despite rapid advances of next-generation sequencing, both the technology and the analysis of the data it produces are in its infancy. At present there is abundant knowledge pertaining to the role of rare single nucleotide variants (SNVs) in rare disorders and of common SNVs in common disorders. Although the 1,000 genome project has clearly highlighted the prevalence of rare variants and more complex variants (e.g. insertions, deletions), their role in disease is as yet far from elucidated.We set out to analyse the properties of sequence variants identified in a comprehensive collection of exome re-sequencing studies performed on samples from patients affected by a broad range of complex and rare diseases (N = 173). Given the known potential for Loss of Function (LoF) variants to be false positive, we performed an extensive validation of the common, rare and private LoF variants identified, which indicated that most of the private and rare variants identified were indeed true, while common novel variants had a significantly higher false positive rate. Our results indicated a strong enrichment of very low-frequency insertion/deletion variants, so far under-investigated, which might be difficult to capture with low coverage and imputation approaches and for which most of study designs would be under-powered. These insertions and deletions might play a significant role in disease genetics, contributing specifically to the underlining rare and private variation predicted to be discovered through next generation sequencing.

Identification of germline susceptibility loci in ETV6-RUNX1-rearranged childhood acute lymphoblastic leukemia

Acute lymphoblastic leukemia (ALL) is a malignant disease of the white blood cells. The etiology of ALL is believed to be multifactorial and likely to involve an interplay of environmental and genetic variables. We performed a genome-wide association study of 355 750 single-nucleotide polymorphisms (SNPs) in 474 controls and 419 childhood ALL cases characterized by a t(12;21)(p13;q22) - the most common chromosomal translocation observed in childhood ALL - which leads to an ETV6-RUNX1 gene fusion. The eight most strongly associated SNPs were followed-up in 951 ETV6-RUNX1-positive cases and 3061 controls from Germany/Austria and Italy, respectively. We identified a novel, genome-wide significant risk locus at 3q28 (TP63, rs17505102, P(CMH)=8.94 × 10(-9), OR=0.65). The separate analysis of the combined German/Austrian sample only, revealed additional genome-wide significant associations at 11q11 (OR8U8, rs1945213, P=9.14 × 10(-11), OR=0.69) and 8p21.3 (near INTS10, rs920590, P=6.12 × 10(-9), OR=1.36). These associations and another association at 11p11.2 (PTPRJ, rs3942852, P=4.95 × 10(-7), OR=0.72) remained significant in the German/Austrian replication panel after correction for multiple testing. Our findings demonstrate that germline genetic variation can specifically contribute to the risk of ETV6-RUNX1-positive childhood ALL. The identification of TP63 and PTPRJ as susceptibility genes emphasize the role of the TP53 gene family and the importance of proteins regulating cellular processes in connection with tumorigenesis.

An APOE haplotype associated with decreased ε4 expression increases the risk of late onset Alzheimer's disease

This paper addresses a tenet of the literature on APOE, i.e., the relationship between the effects of the ε4, one of the established genetic risk factor for Alzheimer's disease (AD), and its expression levels as determined by APOE promoter polymorphisms. Five polymorphisms (-491 rs449647, -427 rs769446, -219 rs405509, and ε rs429358-rs7412) were studied in 1308 AD patients and 1082 control individuals from the Central-Northern Italy. Major findings of the present study are the following: 1) the variants -219T and ε4 increase the risk for late onset AD (LOAD) when they are both present in cis on the same chromosome (in phase); 2) the correlation between the haplotype (-219T/ε4) and AD risk persists when the data are stratified by age; 3) this haplotype likely anticipates the age of onset of the disease. These data, while confirming the association between -219T and AD, highlight the importance of the phase of the alleles for the observed effects on AD risk, suggesting that this information has to be taken into account when assessing the AD genetic risk. Moreover, the data help to clarify the apparent discrepancy that emerges from the genetic analysis where an SNP characterizing the haplotype responsible for an increased risk for LOAD is coherently associated with a reduced expression of ApoE levels. Our data are compatible with the hypothesis of a complex role of ApoE in the AD pathogenesis, with positive and negative effects occurring concomitantly according to its expression levels and its protein-protein interactions largely unclarified.

ParkDB: a Parkinson's disease gene expression database

Parkinson’s disease (PD) is a common, adult-onset, neuro-degenerative disorder characterized by the degeneration of cardinal motor signs mainly due to the loss of dopaminergic neurons in the substantia nigra. To date, researchers still have limited understanding of the key molecular events that provoke neurodegeneration in this disease. Here, we present ParkDB, the first queryable database dedicated to gene expression in PD. ParkDB contains a complete set of re-analyzed, curated and annotated microarray datasets. This resource enables scientists to identify and compare expression signatures involved in PD and dopaminergic neuron differentiation under different biological conditions and across species.

Database URL:http://www2.cancer.ucl.ac.uk/Parkinson_Db2/

Failure to replicate an association of rs5984894 SNP in the PCDH11X gene in a collection of 1,222 Alzheimer's disease affected patients

A recent genome-wide study on late-onset Alzheimer's disease identified a SNP (rs5984894) on Xq21.3 in the PCDH11X gene strongly associated with LOAD individuals of European descent from the United States. We genotyped the same polymorphism in 1222 cases and 938 controls from central-northern Italy and could not confirm the association on the Italian population: multivariate logistic regression adjusted for gender and APOE epsilon4 allele resulted in a global p value of 0.56.

Systems biology and longevity: an emerging approach to identify innovative anti-aging targets and strategies

Human aging and longevity are complex and multi-factorial traits that result from a combination of environmental, genetic, epigenetic and stochastic factors, each contributing to the overall phenotype. The multi-factorial process of aging acts at different levels of complexity, from molecule to cell, from organ to organ systems and finally to organism, giving rise to the dynamic "aging mosaic". At present, an increasing amount of experimental data on genetics, genomics, proteomics and other -omics are available thanks to new high-throughput technologies but a comprehensive model for the study of human aging and longevity is still lacking. Systems biology represents a strategy to integrate and quantify the existing knowledge from different sources into predictive models, to be later tested and then implemented with new experimental data for validation and refinement in a recursive process. The ultimate goal is to compact the new acquired knowledge into a single picture, ideally able to characterize the phenotype at systemic/organism level. In this review we will briefly discuss the aging phenotype in a systems biology perspective, showing four specific examples at different levels of complexity, from a systemic process (inflammation) to a cascade-process pathways (coagulation) and from cellular organelle (proteasome) to single gene-network (PON-1), which could also represent targets for anti-aging strategies.

The impact of phenocopy on the genetic analysis of complex traits

A consistent debate is ongoing on genome-wide association studies (GWAs). A key point is the capability to identify low-penetrance variations across the human genome. Among the phenomena reducing the power of these analyses, phenocopy level (PE) hampers very seriously the investigation of complex diseases, as well known in neurological disorders, cancer, and likely of primary importance in human ageing. PE seems to be the norm, rather than the exception, especially when considering the role of epigenetics and environmental factors towards phenotype. Despite some attempts, no recognized solution has been proposed, particularly to estimate the effects of phenocopies on the study planning or its analysis design. We present a simulation, where we attempt to define more precisely how phenocopy impacts on different analytical methods under different scenarios. With our approach the critical role of phenocopy emerges, and the more the PE level increases the more the initial difficulty in detecting gene-gene interactions is amplified. In particular, our results show that strong main effects are not hampered by the presence of an increasing amount of phenocopy in the study sample, despite progressively reducing the significance of the association, if the study is sufficiently powered. On the opposite, when purely epistatic effects are simulated, the capability of identifying the association depends on several parameters, such as the strength of the interaction between the polymorphic variants, the penetrance of the polymorphism and the alleles (minor or major) which produce the combined effect and their frequency in the population. We conclude that the neglect of the possible presence of phenocopies in complex traits heavily affects the analysis of their genetic data.

PON1 is a longevity gene: results of a meta-analysis

Paraoxonase 1 (PON1) is one of the most studied genes regarding cardiovascular risk, oxidative stress and inflammation. Several lines of evidence suggests that PON1 promotes an atheroprotective effect. Patients carrying PON1 codon 192 QQ genotype display a higher risk of cardiovascular events, the major cause of mortality in the elderly: it can be predicted that gene variants increasing the risk of mortality will be under-represented in long-living individuals. We first reported that PON1 R allele (R+) carriers are significantly more represented in Italian centenarians; subsequently this topic has been addressed by many other groups, and here we report a meta-analysis on 11 studies in different populations selected by a review of the literature available in PubMed and testing the effect of the Q192R polymorphism on human ageing. QUORUM guidelines for meta-analysis have been followed, and a total number of 5962 subjects have been included: 2795 young controls (<65 years of age) and 3167 old subjects (>65 years of age). The Mantel-Haenszel weighting for pooling in presence of a fixed effects model has been applied. The meta-analysis of R carriers showed a significant result with an overall OR of 1.16 (1.04-1.30, 95% CI, p=0.006). The meta-analysis of QR genotype also showed a significant result, with an overall OR of 1.14 (1.02-1.27, 95% CI, p=0.016). The results show that PON1 gene variants at codon 192 impact on the probability of attaining longevity, and that subjects carrying RR and QR genotypes (R+ carriers) are favoured in reaching extreme ages. These results likely represent the counterpart of the effects observed on cardiovascular diseases (CVD), as centenarians and nonagenarians escaped or delayed the onset of the major age-related diseases, including CVD.

Human longevity and 11p15.5: a study in 1321 centenarians

The 11p15.5 chromosomal region (2.8 Mb) is of particular interest as it encloses five genes (HRAS1, SIRT3, TH, INS and IGF2), the variability of which was found to be associated with life extension by association studies. Mostly important, the above genes are homologous of genes that modulate lifespan in model organisms. We scanned the area in four European sample groups for a total of 1321 centenarians and 1140 younger subjects, who shared with centenarians ethnicity and geographical origin, with a set of 239 SNPs. No significant results (P<0.05) have been found on the earlier associated loci (ie, TH, IGF2, INS and HRAS1), and this study could not confirm the earlier findings on each of those genes. A meta-analysis was carried out on the SIRT3 SNP data; a total number of 2461 samples were included, but no positive association was found except for one SNP having a significant effect (rs939915). The same meta-analysis approach has been applied to the other 229 markers, and six SNPs have been found significant for the frequent genotype (rs4073591, DEAF1-rs4073590, KRTAP5-6-rs11040489, rs4930001, TSPAN32-rs800140 and rs16928120). This experience, although unable to confirm the earlier findings of the literature, highlights all the common difficulties of such studies in human longevity. Despite the rather negative findings presented here, the results derived from unprecedented studies involving such a large number of centenarians should be disseminated, thus contributing to set up adequate strategies to disentangle complex and likely heterogeneous phenotypes.

Human models of aging and longevity

BACKGROUND

The aging phenotype in humans is very heterogeneous and can be described as a complex mosaic resulting from the interaction of a variety of environmental, stochastic and genetic-epigenetic variables. Therefore, each old person must be considered as a singleton, and consequently the definition of 'aging phenotype' is very difficult.

OBJECTIVE

We discuss the phenotype of centenarians, the best example of successful aging, as well as other models exploited to study human aging and longevity, such as families enriched in long-living subjects, twins and cohorts of unrelated subjects.

METHODS

A critical review of literature available until March 2008.

CONCLUSIONS

No single model can be considered the gold standard for the study of aging and longevity, instead the combination of results obtained from different models must be considered in order to better understand these complex phenomena. We propose that a systems biology concept such as that of 'bow-tie' architecture, useful for managing information flow, could help in this demanding task.

Association between the interleukin-1beta polymorphisms and Alzheimer's disease: a systematic review and meta-analysis

The pro-inflammatory cytokine interleukin(IL)-1beta is a main component in inflammatory pathways and is overexpressed in the brain of Alzheimer's disease (AD) patients. Several studies report associations between IL-1beta polymorphisms and AD, but findings from different studies are controversial. Our aim was to verify the correlation between the single nucleotide polymorphisms (SNPs) of the IL-1beta, at sites -511 and +3953, and AD by meta-analysis. Computerized bibliographic searches of PUBMED and AlzGene database (http://www.alzgene.org) were supplemented with manual searches of reference lists. There is evidence for association between IL-1beta +3953 SNP and AD, with an OR=1.60 (95% C.I.: 1.16-2.22; Z=2.83 p=0.005) for TT genotype. No significant difference in genotype distribution of the IL-1beta -511 SNP in AD was obtained, but high between-study heterogeneity was found. To reduce heterogeneity, subgroup analyses were performed using, as stratifying variables, characteristics of the population under study (age, gender, type of AD diagnosis, Mini Mental State Examination of the controls) and characteristics related to the study design (statistical power of individual studies). The frequency of the IL-1beta -511 TT genotype resulted significantly higher than other genotypes only when the Caucasian studies with the highest statistical power were included in the subgroup analysis (OR=1.32; 95% C.I.: 1.03-1.69; p=0.03), with no evidence of between-study heterogeneity. Our data support an association between the TT genotype of IL-1beta +3953 SNP and AD, and suggest a possible association of the -511 TT genotype. Unreplicability of the results seems to be due mainly to the lack of statistical power of the individual studies.

A genetic-demographic approach reveals male-specific association between survival and tumor necrosis factor (A/G)-308 polymorphism

The (A/G)-308 polymorphism of the tumor necrosis factor alpha gene (TNF) is associated with age-related diseases, but its influence on longevity is controversial. We genotyped for this polymorphism 747 Italian volunteers (401 women and 346 men, age 19-110 years). By applying a genetic-demographic (GD) approach we found that, in men, the survival function of allele A carriers is lower than that of noncarriers at all the ages (p =.044). After defining (by exploiting again demographic information) three age classes, we found that the frequency of men carrying the A allele decreases with age (p =.019), thus confirming the GD analysis results. The same analyses gave negative results in women. Therefore, allele A has a detrimental effect on life expectancy, and this effect is specific to men. A haplotype analysis carried out in men by screening the TNFa, TNFc, and TNFe microsatellite polymorphisms (spanning about 20 kb) confirmed the association of the TNF region with life expectancy.

Role of mitochondrial DNA in longevity, aging and age-related diseases in humans: a reappraisal

The genetic variability of H. sapiens mitochondrial DNA (mtDNA) can be either germ-line inherited or somatically acquired, and its effect on aging and longevity as well as on the pathogenesis of complex age-related diseases is a hot topic. Here we illustrate the complexity of such studies, related to the large genetic variability of mtDNA in different populations and the fact that the rate of the aging process is different in different cells, tissues and organs. As far as concern Alzheimer's disease, the accumulation of somatic mutations in several tissues have been investigated, as well as the inherited mtDNA variability. However, the issue is still controversial and further studies are needed to clarify the role of mtDNA variants in Alzheimer's disease. This review is aimed to summarize the most recent advances in this field. By high throughput mtDNA sequencing and the study of large cohorts of ethnically homogeneous subjects/patients, it is now possible to perform high dimensionality studies in order to clarify the genetic associations among several inherited mtDNA variants and longevity or age-associated diseases in humans.

Genes, ageing and longevity in humans: problems, advantages and perspectives

Many epidemiological data indicate the presence of a strong familial component of longevity that is largely determined by genetics, and a number of possible associations between longevity and allelic variants of genes have been described. A breakthrough strategy to get insight into the genetics of longevity is the study of centenarians, the best example of successful ageing. We review the main results regarding nuclear genes as well as the mitochondrial genome, focusing on the investigations performed on Italian centenarians, compared to those from other countries. These studies produced interesting results on many putative "longevity genes". Nevertheless, many discrepancies are reported, likely due to the population-specific interactions between gene pools and environment. New approaches, including large-scale studies using high-throughput techniques, are urgently needed to overcome the limits of traditional association studies performed on a limited number of polymorphisms in order to make substantial progress to disentangle the genetics of a trait as complex as human longevity.

Complexity of Anti-immunosenescence Strategies in Humans

Immunosenescence is characterized by three main aspects: (i) the shrinkage of the T cell repertoire and the accumulation of oligoclonal expansions (megaclones) of memory/effector cells directed toward ubiquitary infectious agents; (ii) the involution of the thymus and the exhaustion of naïve T cells; and (iii) a chronic inflammatory status called inflamm-aging. We present here possible strategies to counteract these main aspects of immunosenescence in humans with particular attention to the reduction of antigenic load by pathogens, such as CMV, and the normalization of intestinal microflora, the possible utilization of IL-7 to reverse thymic involution, the purging of megaclones, the forced expression of CD28 on T lymphocytes, the reduction of inflamm-aging and the administration of nutrients such as vitamin D. Possible drawbacks of all these strategies are discussed. Finally, the complexity of a rejuvenation approach is stressed, with particular attention to the inhibitory role played by the "old microenvironment" on the performance of progenitor cells, the best candidate to counteract the decline in regenerative potential characteristic of organs and tissues from old organisms.

Genotype of inflammatory cytokines in limbal stem cell graft in Italian patients

We tested the hypothesis that the genetic capability to mount an inflammatory response might contribute to the inter-individual variability of limbal stem cell graft (LSCG) outcome. Two functional polymorphisms in the IL-6 and TNF-alpha promoter regions were genotyped in 35 patients. A new score system (clinical assessment score, CAS) was set up in order to classify patients' clinical profile, and the main parameters relevant for LSCG as well as for the follow-up of the patients. Patients carrying at both loci a genotype associated with a lower production of both cytokines were classified as "low producers" (LP), while all the others were classified as "intermediate or high producers" (HP). LP patients did not show any difference in CAS before and after transplantation while a significant difference was present in HP patients. A similar trend was evident in the 35 months of follow-up. Polymorphisms of IL-6 and TNF-alpha can be used to identify subgroups of patients with higher risk of unsuccessful outcome.

Neuroinflammation and the genetics of Alzheimer's disease: the search for a pro-inflammatory phenotype

The role of interleukin 1 (IL-1) and interleukin 6 (IL-6) in the pathogenesis of Alzheimer's disease (AD) is reviewed within the framework of "inflamm-aging", i.e., the characteristic chronic pro-inflammatory status which develops in old age, and neuroinflammation, i.e., the peculiar inflammatory process which is present in the brain of AD patients. In particular, the data suggesting that several IL-1 and IL-6 gene polymorphisms can contribute to the risk of developing AD are reviewed. The possibility as well as the difficulty in identifying a pro-inflammatory phenotype, and its importance for the prevention, diagnosis and therapy of AD and other age-related pathologies are discussed.

Do men and women follow different trajectories to reach extreme longevity? Italian Multicenter Study on Centenarians (IMUSCE)

Gender accounts for important differences in the incidence and prevalence of a variety of age-related diseases. Considering people of far advanced age, demographic data document a clear-cut prevalence of females compared to males, suggesting that sex-specific mortality rates follow different trajectories during aging. In the present investigation, we report data from a nationwide study on Italian centenarians (a total of 1162 subjects), and from two studies on centenarians living in two distinct zones of Italy, i.e., the island of Sardinia (a total of 222 subjects) and the Mantova province (Northern Italy) (a total of 43 subjects). The female/male ratio was about 2:1 in Sardinia, 4:1 in the whole of Italy, and about 7:1 in the Mantova province. Thus, a complex interaction of environmental, historical and genetic factors, differently characterizing the various parts of Italy, likely plays an important role in determining the gender-specific probability of achieving longevity. Gender differences in the health status of centenarians are also reported, and an innovative score method to classify long-lived people in different health categories, according to clinical and functional parameters, is proposed. Our data indicate that not only is this selected group of people, as a whole, highly heterogeneous, but also that a marked gender difference exists, since male centenarians are less heterogeneous and more healthy than female centenarians. Immunological factors regarding the age-related increase in pro-inflammatory status, and the frequency of HLA ancestral haplotypes also show gender differences that likely contribute to the different strategies that men and women seem to follow to achieve longevity. Concerning the different impact of genetic factors on the probability of reaching the extreme limits of the human life-span, emerging evidence (regarding mtDNA haplogroups, Thyrosine Hydroxilase, and IL-6 genes) suggests that female longevity is less dependent on genetics than male longevity, and that female centenarians likely exploited a healthier life-style and more favorable environmental conditions, owing to gender-specific cultural and anthropological characteristics of the Italian society in the last 100 years.