Luca Marcolungo

Background/Objectives: Dilated cardiomyopathy (DCM) is a genetically heterogeneous myocardial disorder. Emerging evidence suggests that some genes implicated in DCM may also be associated with neurological disorders, supporting the concept of genetic pleiotropy. This study explored the intersection between cardiac and neurological genetics, with the aim of identifying candidate genes that may contribute to shared pathogenic pathways linking these clinically distinct conditions. Methods: We performed exome sequencing in 149 patients with echocardiographically confirmed DCM and subsequently applied an in silico filter to a predefined list of 211 genes associated with inherited cardiomyopathies. Variants were classified according to the American College of Medical Genetics and Genomics (ACMG) criteria. Genes with validated evidence for DCM according to the Clinical Genome Resource (ClinGen) were further investigated through the Human Gene Mutation Database (HGMD) and a focused literature review to identify reported associations with neurological disorders. Results: Genetic variants in DCM-associated genes were identified in 105 patients. Overall, 137 variants were detected, including pathogenic variants and variants of uncertain significance. The most frequently involved genes were TTN, FLNC, and MYH6. Several DCM-associated genes also showed reported associations with neurological disorders, including autism spectrum disorder, Alzheimer’s disease, Parkinson’s disease, epilepsy, and schizophrenia. Among them, TTN, FLNC, RYR2, and SCN5A displayed the broadest overlap between cardiac and neurological phenotypes. Conclusions: These descriptive findings show that variants were most frequently observed in TTN, FLNC, and MYH6 and that several genes included in the ClinGen DCM curation framework have also been independently reported in neurological disorders. Because most identified variants were VUS and no control group or systematic neurological phenotyping was available, the findings indicate gene-level co-annotation only and do not establish variant enrichment, shared pathogenic mechanisms, or clinical overlap.

Background/Objectives: Parkinson’s disease (PD) is an adult-onset neurodegenerative disorder whose pathogenesis is still not completely understood. Several lines of evidence suggest that alterations in epigenetic architecture may contribute to the development of this condition. Here, we present a pilot DNA methylation study from peripheral blood in a cohort of Sicilian PD patients and matched controls. Peripheral tissue analysis has previously been shown to reflect molecular and functional profiles relevant to neurological diseases, supporting their validity as a proxy for studying brain-related epigenetic mechanisms. Methods: We analyzed 20 PD patients and 20 healthy controls (19 males and 21 females overall), matched for sex, with an age range of 60–87 years (mean 72.3 years). Peripheral blood DNA was extracted and processed using the Illumina Infinium MethylationEPIC v2.0 BeadChip, which interrogates over 935,000 CpG sites across the genome, including promoters, enhancers, CpG islands, and other regulatory elements. The assay relies on sodium bisulfite conversion of DNA to detect methylation status at single-base resolution. Results: Epigenome-wide association study (EWAS) data allowed for multiple levels of analysis, including immune cell-type deconvolution, estimation of biological age (epigenetic clocks), quantification of stochastic epigenetic mutations (SEMs) as a measure of epigenomic stability, and differential methylation profiling. Immune cell-type inference revealed an increased but not significant proportion of monocytes in PD patients, consistent with previous reports. In contrast, epigenetic clock analysis did not reveal significant differences in biological age acceleration between cases and controls, partially at odds with earlier studies—likely due to the limited sample size. SEMs burden did not differ significantly between groups. Epivariations reveal genes involved in pathways known to be altered in dopaminergic neuron dysfunction and α-synuclein toxicity. Differential methylation analysis, however, yielded 167 CpG sites, of which 55 were located within genes, corresponding to 54 unique loci. Gene Ontology enrichment analysis highlighted significant overrepresentation of pathways with neurological relevance, including regulation of synapse structure and activity, axonogenesis, neuron migration, and synapse organization. Notably, alterations in KIAA0319, a gene involved in neuronal migration, synaptic formation, and cortical development, have previously been associated with Parkinson’s disease at the gene expression level, while methylation changes in FAM50B have been reported in neurotoxic and cognitive contexts; our data suggest, for the first time, a potential epigenetic involvement of both genes in Parkinson’s disease. Conclusions: This pilot study on a Sicilian population provides further evidence that DNA methylation profiling can yield valuable molecular insights into PD. Despite the small sample size, our results confirm previously reported findings and highlight biological pathways relevant to neuronal structure and function that may contribute to disease pathogenesis. These data support the potential of epigenetic profiling of peripheral blood as a tool to advance the understanding of PD and generate hypotheses for future large-scale studies.

High-throughput genotyping enables the large-scale analysis of genetic diversity in population genomics and genome-wide association studies that combine the genotypic and phenotypic characterization of large collections of accessions. Sequencing-based approaches for genotyping are progressively replacing traditional genotyping methods because of the lower ascertainment bias. However, genome-wide genotyping based on sequencing becomes expensive in species with large genomes and a high proportion of repetitive DNA. Here we describe the use of CRISPR-Cas9 technology to deplete repetitive elements in the 3.76-Gb genome of lentil (Lens culinaris), 84% consisting of repeats, thus concentrating the sequencing data on coding and regulatory regions (single-copy regions). We designed a custom set of 566,766 gRNAs targeting 2.9 Gbp of repeats and excluding repetitive regions overlapping annotated genes and putative regulatory elements based on ATAC-seq data. The novel depletion method removed ∼40% of reads mapping to repeats, increasing those mapping to single-copy regions by ∼2.6-fold. When analyzing 25 million fragments, this repeat-to-single-copy shift in the sequencing data increased the number of genotyped bases of ∼10-fold compared to nondepleted libraries. In the same condition, we were also able to identify ∼12-fold more genetic variants in the single-copy regions and increased the genotyping accuracy by rescuing thousands of heterozygous variants that otherwise would be missed because of low coverage. The method performed similarly regardless of the multiplexing level, type of library or genotypes, including different cultivars and a closely related species (L. orientalis). Our results showed that CRISPR-Cas9-driven repeat depletion focuses sequencing data on single-copy regions, thus improving high-density and genome-wide genotyping in large and repetitive genomes.

High-throughput chromosome conformation capture (Hi-C) is widely used for scaffolding in de novo assembly because it produces highly contiguous genomes, but its indirect statistical approach can introduce connection errors. We employed optical mapping (Bionano Genomics) as an orthogonal scaffolding technology to assess the structural solidity of Hi-C reconstructed scaffolds. Optical maps were used to assess the correctness of five de novo genome assemblies based on long-read sequencing for contig generation and Hi-C for scaffolding. Hundreds of inconsistencies were found between the reconstructions generated using the Hi-C and optical mapping approaches. Manual inspection, exploiting raw long-read sequencing data and optical maps, confirmed that several of these conflicts were derived from Hi-C joining errors. Such misjoins were widespread, involved the connection of both small and large contigs, and even overlapped annotated genes. We conclude that the integration of optical mapping data after, not before, Hi-C-based scaffolding, improves the quality of the assembly and limits reconstruction errors by highlighting misjoins that can then be subjected to further investigation.

Rapid and sensitive assays for the identification of plant pathogens are necessary for the effective management of crop diseases. The main limitation of current diagnostic testing is the inability to combine broad and sensitive pathogen detection with the identification of key strains, pathovars, and subspecies. Such discrimination is necessary for quarantine pathogens, whose management is strictly dependent on genotype identification. To address these needs, we have established and evaluated a novel all-in-one diagnostic assay based on nanopore sequencing for the detection and simultaneous characterization of quarantine pathogens, using Xylella fastidiosa as a case study. The assay proved to be at least as sensitive as standard diagnostic tests and the quantitative results agreed closely with qPCR-based analysis. The same sequencing results also allowed discrimination between subspecies when present either individually or in combination. Pathogen detection and typing were achieved within 13 min of sequencing owing to the use of an internal control that allowed to stop sequencing when sufficient data had accumulated. These advantages, combined with the use of portable equipment, will facilitate the development of next-generation diagnostic assays for the efficient monitoring of other plant pathogens.

The reconstruction of individual haplotypes can facilitate the interpretation of disease risks; however, high costs and technical challenges still hinder their assessment in clinical settings. Second-generation sequencing is the gold standard for variant discovery but, due to the production of short reads covering small genomic regions, allows only indirect haplotyping based on statistical methods. In contrast, third-generation methods such as the nanopore sequencing platform developed by Oxford Nanopore Technologies (ONT) generate long reads that can be used for direct haplotyping, with fewer drawbacks. However, robust standards for variant phasing in ONT-based target resequencing efforts are not yet available. In this study, we presented a streamlined proof-of-concept workflow for variant calling and phasing based on ONT data in a clinically relevant 12-kb region of the APOE locus, a hotspot for variants and haplotypes associated with aging-related diseases and longevity. Starting with sequencing data from simple amplicons of the target locus, we demonstrated that ONT data allow for reliable single-nucleotide variant (SNV) calling and phasing from as little as 60 reads, although the recognition of indels is less efficient. Even so, we identified the best combination of ONT read sets (600) and software (BWA/Minimap2 and HapCUT2) that enables full haplotype reconstruction when both SNVs and indels have been identified previously using a highly-accurate sequencing platform. In conclusion, we established a rapid and inexpensive workflow for variant phasing based on ONT long reads. This allowed for the analysis of multiple samples in parallel and can easily be implemented in routine clinical practice, including diagnostic testing.

Structural Variants (SVs) are a widely unexplored source of genetic variation, both due to methodological limitations and because they are generally associated to deleterious effects. However, with the advent of long-range genomic platforms, it has become easier to directly detect SVs. In the same direction, clonally propagated crops provide a unique opportunity to study SVs, offering a suitable genomic environment for their accumulation in heterozygosis. In particular, it has been reported that SVs generate drastic levels of heterozygosity in grapevines. ‘Nebbiolo’ ( Vitis vinifera L.) is a grapevine cultivar typical of north-western Italy, appreciated for its use in producing high-quality red wines. Here, we aimed to analyze the frequency of SVs in ‘Nebbiolo’, at three different organizational levels. For this purpose, we generated genomic data based on long-reads, linked-reads and optical mapping. We assembled a reference genome for this cultivar and compared two different clones, including V. vinifera reference genome (PN40024) in our comparisons. Our results indicate that SVs differentially occurring between ‘Nebbiolo’ clones might be rare, while SVs differentiating haplotypes of the same individual are as abundant as those that occur differentially between cultivars.

Whole-exome sequencing (WES) enrichment platforms are usually evaluated by measuring the depth of coverage at target regions. However, variants called in WES are reported in the variant call format (VCF) file, which is filtered by minimum site coverage and mapping quality. Therefore, genotypability (base calling calculated by combining depth of coverage with the confidence of read alignment) should be considered as a more informative parameter to assess the performance of WES. We found that the mapping quality of reads aligned to difficult target regions was improved by increasing the DNA fragment length well above the average exon size. We tested three different DNA fragment lengths using four major commercial WES platforms and found that longer DNA fragments achieved a higher percentage of callable bases in the target regions and thus improved the genotypability of many genes, including several associated with clinical phenotypes. DNA fragment size also affected the uniformity of coverage, which in turn influences genotypability, indicating that different platforms are optimized for different DNA fragment lengths. Finally, we found that although the depth of coverage continued to increase in line with the sequencing depth (overall number of reads), base calling reached saturation at a depth of coverage that depended on the enrichment platform and DNA fragment length. This confirmed that genotypability provides better estimates for the optimal sequencing depth of each fragment size/enrichment platform combination.

Genetic markers (DNA barcodes) are often used to support and confirm species identification. Barcode sequences can be generated in the field using portable systems based on the Oxford Nanopore Technologies (ONT) MinION sequencer. However, to achieve a broader application, current proof-of-principle workflows for on-site barcoding analysis must be standardized to ensure a reliable and robust performance under suboptimal field conditions without increasing costs. Here, we demonstrate the implementation of a new on-site workflow for DNA extraction, PCR-based barcoding, and the generation of consensus sequences. The portable laboratory features inexpensive instruments that can be carried as hand luggage and uses standard molecular biology protocols and reagents that tolerate adverse environmental conditions. Barcodes are sequenced using MinION technology and analyzed with ONTrack, an original de novo assembly pipeline that requires as few as 1000 reads per sample. ONTrack-derived consensus barcodes have a high accuracy, ranging from 99.8 to 100%, despite the presence of homopolymer runs. The ONTrack pipeline has a user-friendly interface and returns consensus sequences in minutes. The remarkable accuracy and low computational demand of the ONTrack pipeline, together with the inexpensive equipment and simple protocols, make the proposed workflow particularly suitable for tracking species under field conditions.