Replicative senescence and high glucose induce the accrual of self-derived cytosolic nucleic acids in human endothelial cells.

Recent literature shows that loss of replicative ability and acquisition of a proinflammatory secretory phenotype in senescent cells is coupled with the build-in of nucleic acids in the cytoplasm. Its implication in human age-related diseases is under scrutiny. In human endothelial cells (ECs), we assessed the accumulation of intracellular nucleic acids during in vitro replicative senescence and after exposure to high glucose concentrations, which mimic an in vivo condition of hyperglycemia. We showed that exposure to high glucose induces senescent-like features in ECs, including telomere shortening and proinflammatory cytokine release, coupled with the accrual in the cytoplasm of telomeres, double-stranded DNA and RNA (dsDNA, dsRNA), as well as RNA:DNA hybrid molecules. Senescent ECs showed an activation of the dsRNA sensors RIG-I and MDA5 and of the DNA sensor TLR9, which was not paralleled by the involvement of the canonical (cGAS) and non-canonical (IFI16) activation of the STING pathway. Under high glucose conditions, only a sustained activation of TLR9 was observed. Notably, senescent cells exhibit increased proinflammatory cytokine (IL-1ß, IL-6, IL-8) production without a detectable secretion of type I interferon (IFN), a phenomenon that can be explained, at least in part, by the accumulation of methyl-adenosine containing RNAs. At variance, exposure to exogenous nucleic acids enhances both IL-6 and IFN-ß1 expression in senescent cells. This study highlights the accrual of cytoplasmic nucleic acids as a marker of senescence-related endothelial dysfunction, that may play a role in dysmetabolic age-related diseases.

High frequency of heterozygous rare variants of the SLC34A1 and SLC9A3R1 genes in patients with atypical femur fracture.

OBJECTIVE: Atypical femur fractures (AFFs) are rare fragility fractures originating at the lateral cortex of the femur, affecting the subtrochanteric or diaphyseal area of thebone with a transverse morphology. Occurrence of AFF is specifically associated with a small number of rare monogenic congenital metabolic bone disorders, such as hypophosphatasia, and with long-term treatment with antiresorptiondrugs. The exact pathogenesis of these fractures remains poorly understood and, except for cases of diagnosed HPP or other AFF-causing bone diseases, it is not possible to assess which patients are at higher riskof developing AFFs as a consequence of anti-resorption therapy. DESIGN: We genetically screened 25 unrelated patients who had developed at least one AFF. INTERVENTION: Genetic screening was performed through a nextgeneration sequencing analysis with a customized panel containing 76 human genes involved in the regulation of the mineralization processWe genetically screened 25 unrelated patients who had developed at least one AFF. RESULTS: We found a relatively high frequency (32.0%) of heterozygous rare variants inthe SLC34A1 and SLC9A3R1 genes, two genes whose heterozygous inactivating mutations have been respectively associated with autosomal dominant hypophosphatemic nephrolithiasis/osteoporosis types 1 and 2 (NPHLOP1and NPHLOP2). Other heterozygous rare variants were found in the BMPR1B, CYP27B1, FBN1, MEPE, PIGO, and PHOSPHO1 genes, each in a single AFF case (4.0%). CONCLUSIONS AND RELEVANCE: Our findings suggest that rarevariants of SLC34A1 and SLC9A3R1 could represent a possible genetic risk factor for the occurrence of AFFs. On the other hand, AFFs could represent an unsuspected clinical manifestation and/or an anti-resorption therapycorrelatedadverse event in patients with NPHLOP disorders.

Genome-Wide Methylation Changes Associated with Replicative Senescence and Differentiation in Endothelial and Bone Marrow Mesenchymal Stromal Cells.

Bone marrow mesenchymal stromal cells (BMSCs) are multipotent cells able to self-renew and differentiate, depending on the microenvironment, into adipocytes and osteoblasts. These cells have a limited number of replications and enter replicative senescence during in vitro expansion. The role of DNA methylation (DNAm) assumes importance in cell function and commitment; however, its exact contribution to BMSC differentiation and replicative senescence is still unclear. We performed a genome-wide DNAm analysis on BMSCs cultured in vitro at early passages and induced to differentiate into adipocytes and osteoblasts, and on replicative senescent BMSCs and HUVECs, to identify DNAm patterns of senescence and differentiation. We also compared BMSCs and HUVECs in replicative senescence and found that, in both cellular systems, genome-wide hypomethylation was accompanied by a higher-than-expected overlap of differentially methylated positions (DMPs) and concordance in terms of direction of the change. A Kyoto Encyclopedia of Genes and Genomes (KEGG) pathway enrichment analysis on lineage-independent senescence-associated DMPs revealed 16 common pathways, including Insulin resistance, Molecule adhesion, and Wnt/ß-catenin signaling. In both adipogenesis and osteogenesis, we observed a general demethylation of CpG sites compared with undifferentiated BMSCs with a higher number of DMPs in osteogenesis. KEGG analysis resulted in 30 pathways enriched in osteoblasts and only 2 in adipocytes when compared to undifferentiated cells. When comparing differentiated BMSCs with senescent ones, osteogenesis exhibited a greater overlap with senescence in terms of number of DMPs and direction of methylation change compared to adipogenesis. In conclusion, this study may be useful for future research on general mechanisms that occur in replicative senescence and furthermore to identify trajectories of BMSC differentiation and common aspects of differentiated and senescent cells.

Early downregulation of hsa-miR-144-3p in serum from drug-naïve Parkinson's disease patients.

Advanced age represents one of the major risk factors for Parkinson's Disease. Recent biomedical studies posit a role for microRNAs, also known to be remodelled during ageing. However, the relationship between microRNA remodelling and ageing in Parkinson's Disease, has not been fully elucidated. Therefore, the aim of the present study is to unravel the relevance of microRNAs as biomarkers of Parkinson's Disease within the ageing framework. We employed Next Generation Sequencing to profile serum microRNAs from samples informative for Parkinson's Disease (recently diagnosed, drug-naïve) and healthy ageing (centenarians) plus healthy controls, age-matched with Parkinson's Disease patients. Potential microRNA candidates markers, emerging from the combination of differential expression and network analyses, were further validated in an independent cohort including both drug-naïve and advanced Parkinson's Disease patients, and healthy siblings of Parkinson's Disease patients at higher genetic risk for developing the disease. While we did not find evidences of microRNAs co-regulated in Parkinson's Disease and ageing, we report that hsa-miR-144-3p is consistently down-regulated in early Parkinson's Disease patients. Moreover, interestingly, functional analysis revealed that hsa-miR-144-3p is involved in the regulation of coagulation, a process known to be altered in Parkinson's Disease. Our results consistently show the down-regulation of hsa-mir144-3p in early Parkinson's Disease, robustly confirmed across a variety of analytical and experimental analyses. These promising results ask for further research to unveil the functional details of the involvement of hsa-mir144-3p in Parkinson's Disease.

Clinical, Radiometabolic and Immunologic Effects of Olaparib in Locally Advanced Triple Negative Breast Cancer: The OLTRE Window of Opportunity Trial.

INTRODUCTION: Olaparib is effective in metastatic triple negative breast cancer (TNBC) carrying germline mutations in DNA damage repair (DDR) genes BRCA1/2 (gBRCA-mut). The OLTRE window-of-opportunity trial preliminarily investigated potential pathologic, radiometabolic and immune biomarkers of early-response to olaparib in gBRCA-wild-type (wt) TNBC and, as proof-of-concept in gBRCA-mut HER2-negative BC. METHODS: Patients received olaparib for 3 weeks (3w) before standard neoadjuvant chemotherapy and underwent multiple FDG18-PET/CT scan (basal, after olaparib), clinical assessments (basal, every 3w), tumor biopsies and blood samplings (baseline, after olaparib). Clinical and radiometabolic responses were evaluated according to RECIST1.1 and PERCIST criteria. RESULTS: 27 patients with gBRCA-wt TNBC and 8 with gBRCA-mut BC (6 TNBC, 2 HR+/HER2-negative) were enrolled. Three (11.1%) patients showed mutations in non-BRCA1/2 DDR genes and 4 (14.8%) in other genes. 3w olaparib induced 16/35 and 15/27 partial clinical and radiometabolic responses, including in 40.7% and 50.0% gBRCA-wt patients. gBRCA-mut tumors presented numerically higher tumor-infiltrating lymphocytes (TILs) levels and PD-L1 positive tumors. Clinical responders experienced a reduction in T-regs/T-eff ratio (p=0.05), B and NK lymphocytes (p=0.003 both), with an average increase in T-helpers rate (p<0.001) and CD4/CD8 ratio (p=0.02). Ki67% and TILs did not vary significantly (p=0.67 and p=0.77). A numerical increase in PD-L1 positive cases after olaparib was observed, though non-significant (p=0.134). No differences were observed according to gBRCA status and type of response. CONCLUSIONS: Early-stage TNBC might be a target population for olaparib, irrespective of gBRCA mutations. Future trials should combine TILs, PD-L1 and gBRCA status to better identify candidates for escalated/de-escalated treatment strategies including olaparib.

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.

A geroscience approach for Parkinson's disease: Conceptual framework and design of PROPAG-AGEING project.

Advanced age is the major risk factor for idiopathic Parkinson's disease (PD), but to date the biological relationship between PD and ageing remains elusive. Here we describe the rationale and the design of the H2020 funded project "PROPAG-AGEING", whose aim is to characterize the contribution of the ageing process to PD development. We summarize current evidences that support the existence of a continuum between ageing and PD and justify the use of a Geroscience approach to study PD. We focus in particular on the role of inflammaging, the chronic, low-grade inflammation characteristic of elderly physiology, which can propagate and transmit both locally and systemically. We then describe PROPAG-AGEING design, which is based on the multi-omic characterization of peripheral samples from clinically characterized drug-naïve and advanced PD, PD discordant twins, healthy controls and "super-controls", i.e. centenarians, who never showed clinical signs of motor disability, and their offspring. Omic results are then validated in a large number of samples, including in vitro models of dopaminergic neurons and healthy siblings of PD patients, who are at higher risk of developing PD, with the final aim of identifying the molecular perturbations that can deviate the trajectories of healthy ageing towards PD development.

Shedding light on dark genes: enhanced targeted resequencing by optimizing the combination of enrichment technology and DNA fragment length.

The exome contains many obscure regions difficult to explore with current short-read sequencing methods. Repetitious genomic regions prevent the unique alignment of reads, which is essential for the identification of clinically-relevant genetic variants. Long-read technologies attempt to resolve multiple-mapping regions, but they still produce many sequencing errors. Thus, a new approach is required to enlighten the obscure regions of the genome and rescue variants that would be otherwise neglected. This work aims to improve the alignment of multiple-mapping reads through the extension of the standard DNA fragment size. As Illumina can sequence fragments up to 550 bp, we tested different DNA fragment lengths using four major commercial WES platforms and found that longer DNA fragments achieved a higher genotypability. This metric, which indicates base calling calculated by combining depth of coverage with the confidence of read alignment, increased from hundreds to thousands of genes, including several associated with clinical phenotypes. While depth of coverage has been considered crucial for the assessment of WES performance, we demonstrated that genotypability has a greater impact in revealing obscure regions, with ~1% increase in variant calling in respect to shorter DNA fragments. Results confirmed that this approach enlightened many regions previously not explored.

Genomic history of the Italian population recapitulates key evolutionary dynamics of both Continental and Southern Europeans.

BACKGROUND: The cline of human genetic diversity observable across Europe is recapitulated at a micro-geographic scale by variation within the Italian population. Besides resulting from extensive gene flow, this might be ascribable also to local adaptations to diverse ecological contexts evolved by people who anciently spread along the Italian Peninsula. Dissecting the evolutionary history of the ancestors of present-day Italians may thus improve the understanding of demographic and biological processes that contributed to shape the gene pool of European populations. However, previous SNP array-based studies failed to investigate the full spectrum of Italian variation, generally neglecting low-frequency genetic variants and examining a limited set of small effect size alleles, which may represent important determinants of population structure and complex adaptive traits. To overcome these issues, we analyzed 38 high-coverage whole-genome sequences representative of population clusters at the opposite ends of the cline of Italian variation, along with a large panel of modern and ancient Euro-Mediterranean genomes. RESULTS: We provided evidence for the early divergence of Italian groups dating back to the Late Glacial and for Neolithic and distinct Bronze Age migrations having further differentiated their gene pools. We inferred adaptive evolution at insulin-related loci in people from Italian regions with a temperate climate, while possible adaptations to pathogens and ultraviolet radiation were observed in Mediterranean Italians. Some of these adaptive events may also have secondarily modulated population disease or longevity predisposition. CONCLUSIONS: We disentangled the contribution of multiple migratory and adaptive events in shaping the heterogeneous Italian genomic background, which exemplify population dynamics and gene-environment interactions that played significant roles also in the formation of the Continental and Southern European genomic landscapes.

Intronic ATTTC repeat expansions in STARD7 in familial adult myoclonic epilepsy linked to chromosome 2.

Familial Adult Myoclonic Epilepsy (FAME) is characterised by cortical myoclonic tremor usually from the second decade of life and overt myoclonic or generalised tonic-clonic seizures. Four independent loci have been implicated in FAME on chromosomes (chr) 2, 3, 5 and 8. Using whole genome sequencing and repeat primed PCR, we provide evidence that chr2-linked FAME (FAME2) is caused by an expansion of an ATTTC pentamer within the first intron of STARD7. The ATTTC expansions segregate in 158/158 individuals typically affected by FAME from 22 pedigrees including 16 previously reported families recruited worldwide. RNA sequencing from patient derived fibroblasts shows no accumulation of the AUUUU or AUUUC repeat sequences and STARD7 gene expression is not affected. These data, in combination with other genes bearing similar mutations that have been implicated in FAME, suggest ATTTC expansions may cause this disorder, irrespective of the genomic locus involved.

Small supernumerary marker chromosomes: A legacy of trisomy rescue?

We studied by a whole genomic approach and trios genotyping, 12 de novo, nonrecurrent small supernumerary marker chromosomes (sSMC), detected as mosaics during pre- or postnatal diagnosis and associated with increased maternal age. Four sSMCs contained pericentromeric portions only, whereas eight had additional non-contiguous portions of the same chromosome, assembled together in a disordered fashion by repair-based mechanisms in a chromothriptic event. Maternal hetero/isodisomy was detected with a paternal origin of the sSMC in some cases, whereas in others two maternal alleles in the sSMC region and biparental haplotypes of the homologs were detected. In other cases, the homologs were biparental while the sSMC had the same haplotype of the maternally inherited chromosome. These findings strongly suggest that most sSMCs are the result of a multiple-step mechanism, initiated by maternal meiotic nondisjunction followed by postzygotic anaphase lagging of the supernumerary chromosome and its subsequent chromothripsis.

Insertional translocation involving an additional nonchromothriptic chromosome in constitutional chromothripsis: Rule or exception?

BACKGROUND: Chromothripsis, which is the local massive shattering of one or more chromosomes and their reassembly in a disordered array with frequent loss of some fragments, has been mainly reported in association with abnormal phenotypes. We report three unrelated healthy persons, two of which parenting a child with some degree of intellectual disability, carrying a chromothripsis involving respectively one, two, and three chromosomes, which was detected only after whole-genome sequencing. Unexpectedly, in all three cases a fragment from one of the chromothripsed chromosomes resulted to be inserted within a nonchromothripsed one. METHODS: Conventional cytogenetic techniques, paired-end whole-genome sequencing, polymerase chain reaction, and Sanger sequencing were used to characterize complex rearrangements, copy-number variations, and breakpoint sequences in all three families. RESULTS: In two families, one parent was carrier of a balanced chromothripsis causing in the index case a deletion and a noncontiguous duplication at 3q in case 1, and a t(6;14) translocation associated with interstitial 14q deletion in case 2. In the third family, an unbalanced chromothripsis involving chromosomes 6, 7, and 15 was inherited to the proband by the mosaic parent. In all three parents, the chromothripsis was concurrent with an insertional translocation of a portion of one of the chromothriptic chromosomes within a further chromosome that was not involved in the chromothripsis event. CONCLUSION: Our findings show that (a) both simple and complex unbalanced rearrangements may result by the recombination of a cryptic parental balanced chromothripsis and that (b) insertional translocations are the spy of more complex rearrangements and not simply a three-breakpoint event.

Centrosome Linker-induced Tetraploid Segregation Errors Link Rhabdoid Phenotypes and Lethal Colorectal Cancers.

Centrosome anomalies contribute to tumorigenesis, but it remains unclear how they are generated in lethal cancer phenotypes. Here, it is demonstrated that human microsatellite instable (MSI) and BRAFV600E-mutant colorectal cancers with a lethal rhabdoid phenotype are characterized by inactivation of centrosomal functions. A splice site mutation that causes an unbalanced dosage of rootletin (CROCC), a centrosome linker component required for centrosome cohesion and separation at the chromosome 1p36.13 locus, resulted in abnormally shaped centrosomes in rhabdoid cells from human colon tissues. Notably, deleterious deletions at 1p36.13 were recurrent in a subgroup of BRAFV600E-mutant and microsatellite stable (MSS) rhabdoid colorectal cancers, but not in classical colorectal cancer or pediatric rhabdoid tumors. Interfering with CROCC expression in near-diploid BRAFV600E-mutant/MSI colon cancer cells disrupts bipolar mitotic spindle architecture, promotes tetraploid segregation errors, resulting in a highly aggressive rhabdoid-like phenotype in vitro Restoring near-wild-type levels of CROCC in a metastatic model harboring 1p36.13 deletion results in correction of centrosome segregation errors and cell death, revealing a mechanism of tolerance to mitotic errors and tetraploidization promoted by deleterious 1p36.13 loss. Accordingly, cancer cells lacking 1p36.13 display far greater sensitivity to centrosome spindle pole stabilizing agents in vitro These data shed light on a previously unknown link between centrosome cohesion defects and lethal cancer phenotypes providing new insight into pathways underlying genome instability.Implications: Mis-segregation of chromosomes is a prominent feature of chromosome instability and intratumoral heterogeneity recurrent in metastatic tumors for which the molecular basis is unknown. This study provides insight into the mechanism by which defects in rootletin, a centrosome linker component causes tetraploid segregation errors and phenotypic transition to a clinically devastating form of malignant rhabdoid tumor. Mol Cancer Res; 16(9); 1385-95. ©2018 AACR.

Chromothripsis and ring chromosome 22: a paradigm of genomic complexity in the Phelan-McDermid syndrome (22q13 deletion syndrome).

INTRODUCTION: Phelan-McDermid syndrome (PMS) is caused by SHANK3 haploinsufficiency. Its wide phenotypic variation is attributed partly to the type and size of 22q13 genomic lesion (deletion, unbalanced translocation, ring chromosome), partly to additional undefined factors. We investigated a child with severe global neurodevelopmental delay (NDD) compatible with her distal 22q13 deletion, complicated by bilateral perisylvian polymicrogyria (BPP) and urticarial rashes, unreported in PMS. METHODS: Following the cytogenetic and array-comparative genomic hybridization (CGH) detection of a r(22) with SHANK3 deletion and two upstream duplications, whole-genome sequencing (WGS) in blood and whole-exome sequencing (WES) in blood and saliva were performed to highlight potential chromothripsis/chromoanagenesis events and any possible BPP-associated variants, even in low-level mosaicism. RESULTS: WGS confirmed the deletion and highlighted inversion and displaced order of eight fragments, three of them duplicated. The microhomology-mediated insertion of partial Alu-elements at one breakpoint junction disrupted the topological associating domain joining NFAM1 to the transcriptional coregulator TCF20. WES failed to detect BPP-associated variants. CONCLUSIONS: Although we were unable to highlight the molecular basis of BPP, our data suggest that SHANK3 haploinsufficiency and TCF20 misregulation, both associated with intellectual disability, contributed to the patient's NDD, while NFAM1 interruption likely caused her skin rashes, as previously reported. We provide the first example of chromoanasynthesis in a constitutional ring chromosome and reinforce the growing evidence that chromosomal rearrangements may be more complex than estimated by conventional diagnostic approaches and affect the phenotype by global alteration of the topological chromatin organisation rather than simply by deletion or duplication of dosage-sensitive genes.

On site DNA barcoding by nanopore sequencing.

Biodiversity research is becoming increasingly dependent on genomics, which allows the unprecedented digitization and understanding of the planet's biological heritage. The use of genetic markers i.e. DNA barcoding, has proved to be a powerful tool in species identification. However, full exploitation of this approach is hampered by the high sequencing costs and the absence of equipped facilities in biodiversity-rich countries. In the present work, we developed a portable sequencing laboratory based on the portable DNA sequencer from Oxford Nanopore Technologies, the MinION. Complementary laboratory equipment and reagents were selected to be used in remote and tough environmental conditions. The performance of the MinION sequencer and the portable laboratory was tested for DNA barcoding in a mimicking tropical environment, as well as in a remote rainforest of Tanzania lacking electricity. Despite the relatively high sequencing error-rate of the MinION, the development of a suitable pipeline for data analysis allowed the accurate identification of different species of vertebrates including amphibians, reptiles and mammals. In situ sequencing of a wild frog allowed us to rapidly identify the species captured, thus confirming that effective DNA barcoding in the field is possible. These results open new perspectives for real-time-on-site DNA sequencing thus potentially increasing opportunities for the understanding of biodiversity in areas lacking conventional laboratory facilities.

Centenarians as extreme phenotypes: An ecological perspective to get insight into the relationship between the genetics of longevity and age-associated diseases.

In this review, we address the genetic continuum between aging and age-related diseases, with particular attention to the ecological perspective. We describe the connections between genes that promote longevity and genes associated with age-related diseases considering tradeoff mechanisms in which the same genetic variants could have different effects according to the tissue considered and could be involved in several biological pathways. Then we describe mechanisms of antagonistic pleiotropy, focusing on the complex interplay between genetic variants and environmental changes (internal or external). We sustain the use of centenarians as "super-controls" for the study of the major age-related diseases, starting from the concept that the maximization of the phenotypic differences in the considered cohort, achieved by selecting the most divergent phenotypes, could be useful for increasing the significant differences observed in the genetic association study. We describe the potential impact of the population genetic variability in the study of human longevity and the possible contribution of the past selective pressures in shaping the current genomic background of individuals. In conclusion, we illustrate recent findings emerged from whole-genome sequencing of long-lived individuals and future perspectives for interpreting the huge amount of genetic data that will be generated in the next future.

Identification of new BMP6 pro-peptide mutations in patients with iron overload.

Hereditary Hemochromatosis (HH) is a genetically heterogeneous disorder caused by mutations in at least five different genes (HFE, HJV, TFR2, SLC40A1, HAMP) involved in the production or activity of the liver hormone hepcidin, a key regulator of systemic iron homeostasis. Nevertheless, patients with an HH-like phenotype that remains completely/partially unexplained despite extensive sequencing of known genes are not infrequently seen at referral centers, suggesting a role of still unknown genetic factors. A compelling candidate is Bone Morphogenetic Protein 6 (BMP6), which acts as a major activator of the BMP-SMAD signaling pathway, ultimately leading to the upregulation of hepcidin gene transcription. A recent seminal study by French authors has described three heterozygous missense mutations in BMP6 associated with mild to moderate late-onset iron overload (IO). Using an updated next-generation sequencing (NGS)-based genetic test in IO patients negative for the classical HFE p.Cys282Tyr mutation, we found three BMP6 heterozygous missense mutations in four patients from three different families. One mutation (p.Leu96Pro) has already been described and proven to be functional. The other two (p.Glu112Gln, p.Arg257His) were novel, and both were located in the pro-peptide domain known to be crucial for appropriate BMP6 processing and secretion. In silico modeling also showed results consistent with their pathogenetic role. The patients' clinical phenotypes were similar to that of other patients with BMP6-related IO recently described. Our results independently add further evidence to the role of BMP6 mutations as likely contributing factors to late-onset moderate IO unrelated to mutations in the established five HH genes.

Identification of novel mutations in hemochromatosis genes by targeted next generation sequencing in Italian patients with unexplained iron overload.

Hereditary hemochromatosis, one of the commonest genetic disorder in Caucasians, is mainly associated to homozygosity for the C282Y mutation in the HFE gene, which is highly prevalent (allele frequency up to near 10% in Northern Europe) and easily detectable through a widely available "first level" molecular test. However, in certain geographical regions like the Mediterranean area, up to 30% of patients with a HH phenotype has a negative or non-diagnostic (i.e. simple heterozygosity) test, because of a known heterogeneity involving at least four other genes (HAMP, HJV, TFR2, and SLC40A1). Mutations in such genes are generally rare/private, making the diagnosis of atypical HH essentially a matter of exclusion in clinical practice (from here the term of "non-HFE" HH), unless cumbersome traditional sequencing is applied. We developed a Next Generation Sequencing (NGS)-based test targeting the five HH genes, and applied it to patients with clinically relevant iron overload (IO) and a non-diagnostic first level genetic test. We identified several mutations, some of which were novel (i.e. HFE W163X, HAMP R59X, and TFR2 D555N) and allowed molecular reclassification of "non-HFE" HH clinical diagnosis, particularly in some highly selected IO patients without concurring acquired risk factors. This NGS-based "second level" genetic test may represent a useful tool for molecular diagnosis of HH in patients in whom HH phenotype remains unexplained after the search of common HFE mutations.

The Use of Non-Variant Sites to Improve the Clinical Assessment of Whole-Genome Sequence Data.

Genetic testing, which is now a routine part of clinical practice and disease management protocols, is often based on the assessment of small panels of variants or genes. On the other hand, continuous improvements in the speed and per-base costs of sequencing have now made whole exome sequencing (WES) and whole genome sequencing (WGS) viable strategies for targeted or complete genetic analysis, respectively. Standard WGS/WES data analytical workflows generally rely on calling of sequence variants respect to the reference genome sequence. However, the reference genome sequence contains a large number of sites represented by rare alleles, by known pathogenic alleles and by alleles strongly associated to disease by GWAS. It's thus critical, for clinical applications of WGS and WES, to interpret whether non-variant sites are homozygous for the reference allele or if the corresponding genotype cannot be reliably called. Here we show that an alternative analytical approach based on the analysis of both variant and non-variant sites from WGS data allows to genotype more than 92% of sites corresponding to known SNPs compared to 6% genotyped by standard variant analysis. These include homozygous reference sites of clinical interest, thus leading to a broad and comprehensive characterization of variation necessary to an accurate evaluation of disease risk. Altogether, our findings indicate that characterization of both variant and non-variant clinically informative sites in the genome is necessary to allow an accurate clinical assessment of a personal genome. Finally, we propose a highly efficient extended VCF (eVCF) file format which allows to store genotype calls for sites of clinical interest while remaining compatible with current variant interpretation software.

The Solanum commersonii Genome Sequence Provides Insights into Adaptation to Stress Conditions and Genome Evolution of Wild Potato Relatives.

Here, we report the draft genome sequence of Solanum commersonii, which consists of ∼830 megabases with an N50 of 44,303 bp anchored to 12 chromosomes, using the potato (Solanum tuberosum) genome sequence as a reference. Compared with potato, S. commersonii shows a striking reduction in heterozygosity (1.5% versus 53 to 59%), and differences in genome sizes were mainly due to variations in intergenic sequence length. Gene annotation by ab initio prediction supported by RNA-seq data produced a catalog of 1703 predicted microRNAs, 18,882 long noncoding RNAs of which 20% are shown to target cold-responsive genes, and 39,290 protein-coding genes with a significant repertoire of nonredundant nucleotide binding site-encoding genes and 126 cold-related genes that are lacking in S. tuberosum. Phylogenetic analyses indicate that domesticated potato and S. commersonii lineages diverged ∼2.3 million years ago. Three duplication periods corresponding to genome enrichment for particular gene families related to response to salt stress, water transport, growth, and defense response were discovered. The draft genome sequence of S. commersonii substantially increases our understanding of the domesticated germplasm, facilitating translation of acquired knowledge into advances in crop stability in light of global climate and environmental changes.