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.

TLR4 activation initiates a signaling cascade leading to the production of type I IFNs and of the downstream IFN-stimulated genes (ISGs). Recently, a number of IFN-induced long non-coding RNAs (lncRNAs) that feed-back regulate the IFN response have been identified. Dysregulation of this process, collectively known as the “Interferon (IFN) Response,” represents a common molecular basis in the development of autoimmune and autoinflammatory disorders. Concurrently, alteration of lncRNA profile has been described in several type I IFN-driven autoimmune diseases. In particular, both TLR activation and the upregulation of ISGs in peripheral blood mononuclear cells have been identified as possible contributors to the pathogenesis of systemic sclerosis (SSc), a connective tissue disease characterized by vascular abnormalities, immune activation, and fibrosis. However, hitherto, a potential link between specific lncRNA and the presence of a type I IFN signature remains unclear in SSc. In this study, we identified, by RNA sequencing, a group of lncRNAs related to the IFN and anti-viral response consistently modulated in a type I IFN-dependent manner in human monocytes in response to TLR4 activation by LPS. Remarkably, these lncRNAs were concurrently upregulated in a total of 46 SSc patients in different stages of their disease as compared to 18 healthy controls enrolled in this study. Among these lncRNAs, Negative Regulator of the IFN Response (NRIR) was found significantly upregulated in vivo in SSc monocytes, strongly correlating with the IFN score of SSc patients. Weighted Gene Co-expression Network Analysis showed that NRIR-specific modules, identified in the two datasets, were enriched in “type I IFN” and “viral response” biological processes. Protein coding genes common to the two distinct NRIR modules were selected as putative NRIR target genes. Fifteen in silico-predicted NRIR target genes were experimentally validated in NRIR-silenced monocytes. Remarkably, induction of CXCL10 and CXCL11, two IFN-related chemokines associated with SSc pathogenesis, was reduced in NRIR-knockdown monocytes, while their plasmatic level was increased in SSc patients. Collectively, our data show that NRIR affects the expression of ISGs and that dysregulation of NRIR in SSc monocytes may account, at least in part, for the type I IFN signature present in SSc patients.

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.

The common bean (Phaseolus vulgaris L.) is a crucial legume crop and an ideal evolutionary model to study adaptive diversity in wild and domesticated populations. Here, we present a common bean pan-genome based on five high-quality genomes and whole-genome reads representing 339 genotypes. It reveals ~234 Mb of additional sequences containing 6,905 protein-coding genes missing from the reference, constituting 49% of all presence/absence variants (PAVs). More non-synonymous mutations are found in PAVs than core genes, probably reflecting the lower effective population size of PAVs and fitness advantages due to the purging effect of gene loss. Our results suggest pan-genome shrinkage occurred during wild range expansion. Selection signatures provide evidence that partial or complete gene loss was a key adaptive genetic change in common bean populations with major implications for plant adaptation. The pan-genome is a valuable resource for food legume research and breeding for climate change mitigation and sustainable agriculture.

Myotonic dystrophy type 2 (DM2) is caused by CCTG repeat expansions in the CNBP gene, comprising 75 to >11,000 units and featuring extensive mosaicism, making it challenging to sequence fully expanded alleles. To overcome these limitations, we used PCR-free Cas9-mediated nanopore sequencing to characterize CNBP repeat expansions at the single-nucleotide level in nine DM2 patients. The length of normal and expanded alleles can be assessed precisely using this strategy, agreeing with traditional methods, and revealing the degree of mosaicism. We also sequenced an entire ~50 kbp expansion, which has not been achieved previously for DM2 or any other repeat-expansion disorders. Our approach precisely counted the repeats and identified the repeat pattern for both short interrupted and uninterrupted alleles. Interestingly, in the expanded alleles, only two DM2 samples featured the expected pure CCTG repeat pattern, while the other seven presented also TCTG blocks at the 3′ end, which have not been reported before in DM2 patients, but confirmed hereby with orthogonal methods. The demonstrated approach simultaneously determines repeat length, structure/motif, and the extent of somatic mosaicism, promising to improve the molecular diagnosis of DM2 and achieve more accurate genotype–phenotype correlations for the better stratification of DM2 patients in clinical trials.