Evaluation of sensitivity and specificity in RNA-Seq-based detection of grapevine viral pathogens.

Virus detection is a crucial step for the implementation of clean stock programs that preserve healthy crop species. Viral infections in grapevine, a vegetatively propagated perennial crop, cannot be eradicated from the vineyards by the application of agrochemicals and must be curtailed at the stage of nursery production during the propagation of planting material. Viral detection is routinely performed using enzyme-linked immunosorbent assays (ELISA) or Reverse Transcription-quantitative Polymerase Chain Reactions (RT-qPCR). High throughput sequencing (HTS) approaches have the potential to detect all viral pathogens in a plant specimen. However, to date, no published HTS-based study has used threshold selection based on ROC curves for discriminating positive from negative samples. To fill this gap, we assessed the specificity and sensitivity of different sequencing and bioinformatics approaches for nine common viruses, which were tested in the same specimens using ELISA and/or RT-qPCR. The normalized detection thresholds giving the best results were 19.28 Fragments Per Kilobase of transcript per Million mapped reads (FPKM) for alignment-based total RNA-Seq approaches, 386 Reads Per Million mapped reads (RPM) for metagenomics-based total RNA-Seq, 1572 FPKM for alignment-based small RNA-Seq analysis and 0.97 % of contigs for de novo analysis of small RNA-Seq data. Validation of the proposed thresholds using independent specimens collected over time from the same stocks and other specimens collected from nearby stocks that had derived from the same propagating material showed that HTS approaches are accurate, with RNA-Seq approaches showing better performance than small RNA-Seq.

Resistance to Sharka in Apricot: Comparison of Phase-Reconstructed Resistant and Susceptible Haplotypes of 'Lito' Chromosome 1 and Analysis of Candidate Genes.

Sharka, a common disease among most stone fruit crops, is caused by the Plum Pox Virus (PPV). Resistant genotypes have been found in apricot (Prunus armeniaca L.), one of which-the cultivar 'Lito' heterozygous for the resistance-has been used to map a major quantitative trait locus (QTL) on linkage group 1, following a pseudo-test-cross mating design with 231 individuals. In addition, 19 SNP markers were selected from among the hundreds previously developed, which allowed the region to be limited to 236 kb on chromosome 1. A 'Lito' bacterial artificial chromosome (BAC) library was produced, screened with markers of the region, and positive BAC clones were sequenced. Resistant (R) and susceptible (S) haplotypes were assembled independently. To refine the assembly, the whole genome of 'Lito' was sequenced to high coverage (98×) using PacBio technology, enabling the development of a detailed assembly of the region that was able to predict and annotate the genes in the QTL region. The selected cultivar 'Lito' allowed not only to discriminate structural variants between the two haplotypic regions but also to distinguish specific allele expression, contributing towards mining the PPVres locus. In light of these findings, genes previously indicated (i.e., MATHd genes) to have a possible role in PPV resistance were further analyzed, and new candidates were discussed. Although the results are not conclusive, the accurate and independent assembly of R and S haplotypes of 'Lito' is a valuable resource to predict and test alternative transcription and regulation mechanisms underpinning PPV resistance.

Metagenomic profiles of different types of Italian high-moisture Mozzarella cheese.

The microbiota of different types of Italian high-moisture Mozzarella cheese produced using cow or buffalo milk, acidified with natural or selected cultures, and sampled at the dairy or at the mass market, was evaluated using a Next Generation Sequencing approach, in order to identify possible drivers of the bacterial diversity. Cow Mozzarella and buffalo Mozzarella acidified with commercial cultures were dominated by Streptococcus thermophilus, while buffalo samples acidified with natural whey cultures showed similar prevalence of L. delbrueckii subsp. bulgaricus, L. helveticus and S. thermophilus. Moreover, several species of non-starter lactic acid bacteria were frequently detected. The diversity in cow Mozzarella microbiota was much higher than that of water buffalo samples. Cluster analysis clearly separated cow's cheeses from buffalo's ones, the former having a higher prevalence of psychrophilic taxa, and the latter of Lactobacillus and Streptococcus. A higher prevalence of psychrophilic species and potential spoilers was observed in samples collected at the mass retail, suggesting that longer exposures to cooling temperatures and longer production-to-consumption times could significantly affect microbiota diversity. Our results could help in detecting some kind of thermal abuse during the production or storage of mozzarella cheese.

Do you cov me? Effect of coverage reduction on metagenome shotgun sequencing studies.

Shotgun metagenomics sequencing is a powerful tool for the characterization of complex biological matrices, enabling analysis of prokaryotic and eukaryotic organisms and viruses in a single experiment, with the possibility of reconstructing de novo the whole metagenome or a set of genes of interest. One of the main factors limiting the use of shotgun metagenomics on wide scale projects is the high cost associated with the approach. However, we demonstrate that-for some applications-it is possible to use shallow shotgun metagenomics to characterize complex biological matrices while reducing costs. We measured the variation of several summary statistics simulating a decrease in sequencing depth by randomly subsampling a number of reads. The main statistics that were compared are alpha diversity estimates, species abundance, detection threshold, and ability of reconstructing the metagenome in terms of length and completeness. Our results show that a classification of prokaryotic, eukaryotic and viral communities can be accurately performed even using very low number of reads, both in mock communities and in real complex matrices. With samples of 100,000 reads, the alpha diversity estimates were in most cases comparable to those obtained with the full sample, and the estimation of the abundance of all the present species was in excellent agreement with those obtained with the full sample. On the contrary, any task involving the reconstruction of the metagenome performed poorly, even with the largest simulated subsample (1M reads). The length of the reconstructed assembly was smaller than the length obtained with the full dataset, and the proportion of conserved genes that were identified in the meta-genome was drastically reduced compared to the full sample. Shallow shotgun metagenomics can be a useful tool to describe the structure of complex matrices, but it is not adequate to reconstruct-even partially-the metagenome.

Exome Capture with Heterologous Enrichment in Pig (Sus scrofa).

The discovery of new protein-coding DNA variants related to carcass traits is very important for the Italian pig industry, which requires heavy pigs with higher thickness of subcutaneous fat for Protected Designation of Origin (PDO) productions. Exome capture techniques offer the opportunity to focus on the regions of DNA potentially related to the gene and protein expression. In this research a human commercial target enrichment kit was used to evaluate its performances for pig exome capture and for the identification of DNA variants suitable for comparative analysis. Two pools of 30 pigs each, crosses of Italian Duroc X Large White (DU) and Commercial hybrid X Large White (HY), were used and NGS libraries were prepared with the SureSelectXT Target Enrichment System for Illumina Paired-End Sequencing Library (Agilent). A total of 140.2 M and 162.5 M of raw reads were generated for DU and HY, respectively. Average coverage of all the exonic regions for Sus scrofa (ENSEMBL Sus_scrofa.Sscrofa10.2.73.gtf) was 89.33X for DU and 97.56X for HY; and 35% of aligned bases uniquely mapped to off-target regions. Comparison of sequencing data with the Sscrofa10.2 reference genome, after applying hard filtering criteria, revealed a total of 232,530 single nucleotide variants (SNVs) of which 20.6% mapped in exonic regions and 49.5% within intronic regions. The comparison of allele frequencies of 213 randomly selected SNVs from exome sequencing and the same SNVs analyzed with a Sequenom MassARRAY® system confirms that this "human-on-pig" approach offers new potentiality for the identification of DNA variants in protein-coding genes.

Heterozygous reelin mutations cause autosomal-dominant lateral temporal epilepsy.

Autosomal-dominant lateral temporal epilepsy (ADLTE) is a genetic epilepsy syndrome clinically characterized by focal seizures with prominent auditory symptoms. ADLTE is genetically heterogeneous, and mutations in LGI1 account for fewer than 50% of affected families. Here, we report the identification of causal mutations in reelin (RELN) in seven ADLTE-affected families without LGI1 mutations. We initially investigated 13 ADLTE-affected families by performing SNP-array linkage analysis and whole-exome sequencing and identified three heterozygous missense mutations co-segregating with the syndrome. Subsequent analysis of 15 small ADLTE-affected families revealed four additional missense mutations. 3D modeling predicted that all mutations have structural effects on protein-domain folding. Overall, RELN mutations occurred in 7/40 (17.5%) ADLTE-affected families. RELN encodes a secreted protein, Reelin, which has important functions in both the developing and adult brain and is also found in the blood serum. We show that ADLTE-related mutations significantly decrease serum levels of Reelin, suggesting an inhibitory effect of mutations on protein secretion. We also show that Reelin and LGI1 co-localize in a subset of rat brain neurons, supporting an involvement of both proteins in a common molecular pathway underlying ADLTE. Homozygous RELN mutations are known to cause lissencephaly with cerebellar hypoplasia. Our findings extend the spectrum of neurological disorders associated with RELN mutations and establish a link between RELN and LGI1, which play key regulatory roles in both the developing and adult brain.

Suggestive linkage of familial mesial temporal lobe epilepsy to chromosome 3q26.

PURPOSE: To describe the clinical findings in a family with a benign form of mesial temporal lobe epilepsy and to identify the causative genetic factors. METHODS: All participants were personally interviewed and underwent neurologic examination. The affected subjects underwent EEG and most of them neuroradiological examinations (MRI). All family members were genotyped with the HumanCytoSNP-12 v1.0 beadchip and linkage analysis was performed with Merlin and Simwalk2 programs. Exome sequencing was performed on HiSeq2000, after exome capture with SureSelect 50 Mb kit v2.0. RESULTS: The family had 6 members with temporal lobe epilepsy. Age at seizure onset ranged from 8 to 13 years. Five patients had epigastric auras often associated to oro-alimentary automatic activity, 3 patients presented loss of contact, and 2 experienced secondary generalizations. Febrile seizures occurred in 2 family members, 1 of whom also had temporal lobe epilepsy. EEG showed focal slow waves and epileptic abnormalities on temporal regions in 1 patient and was normal in the other affected individuals. MRI was normal in all temporal lobe epilepsy patients. We performed single nucleotide polymorphism-array linkage analysis of the family and found suggestive evidence of linkage (LOD score=2.106) to a region on chromosome 3q26. Haplotype reconstruction supported the linkage data and showed that the majority of unaffected family members carried the haplotype at risk. Whole exome sequencing failed to identify pathogenic mutations in genes of the candidate region. CONCLUSIONS: Our data suggest the existence of a novel locus for benign familial mesial temporal lobe epilepsy on chromosome 3q26. Our failure to identify pathogenic mutations in genes of this region may be due to limitations of the exome sequencing technology.

The SSR-based molecular profile of 1005 grapevine (Vitis vinifera L.) accessions uncovers new synonymy and parentages, and reveals a large admixture amongst varieties of different geographic origin.

A collection of 1005 grapevine accessions was genotyped at 34 microsatellite loci (SSR) with the aim of analysing genetic diversity and exploring parentages. The comparison of molecular profiles revealed 200 groups of synonymy. The removal of perfect synonyms reduced the database to 745 unique genotypes, on which population genetic parameters were calculated. The analysis of kinship uncovered 74 complete pedigrees, with both parents identified. Many of these parentages were not previously known and are of considerable historical interest, e.g. Chenin blanc (Sauvignon × Traminer rot), Covè (Harslevelu selfed), Incrocio Manzoni 2-14 and 2-15 (Cabernet franc × Prosecco), Lagrein (Schiava gentile × Teroldego), Malvasia nera of Bolzano (Perera × Schiava gentile), Manzoni moscato (Raboso veronese × Moscato d'Amburgo), Moscato violetto (Moscato bianco × Duraguzza), Muscat of Alexandria (Muscat blanc à petit grain × Axina de tres bias) and others. Statistical robustness of unexpected pedigrees was reinforced with the analysis of an additional 7-30 SSRs. Grouping the accessions by profile resulted in a weak correlation with their geographical origin and/or current area of cultivation, revealing a large admixture of local varieties with those most widely cultivated, as a result of ancient commerce and population flow. The SSRs with tri- to penta-nucleotide repeats adopted for the present study showed a great capacity for discriminating amongst accessions, with probabilities of identity by chance as low as 1.45 × 10(-27) and 9.35 × 10(-12) for unrelated and full sib individuals, respectively. A database of allele frequencies and SSR profiles of 32 reference cultivars are provided.