Towards a better understanding of the low recall of insertion variants with short-read based variant callers.

BACKGROUND: Since 2009, numerous tools have been developed to detect structural variants using short read technologies. Insertions >50 bp are one of the hardest type to discover and are drastically underrepresented in gold standard variant callsets. The advent of long read technologies has completely changed the situation. In 2019, two independent cross technologies studies have published the most complete variant callsets with sequence resolved insertions in human individuals. Among the reported insertions, only 17 to 28% could be discovered with short-read based tools. RESULTS: In this work, we performed an in-depth analysis of these unprecedented insertion callsets in order to investigate the causes of such failures. We have first established a precise classification of insertion variants according to four layers of characterization: the nature and size of the inserted sequence, the genomic context of the insertion site and the breakpoint junction complexity. Because these levels are intertwined, we then used simulations to characterize the impact of each complexity factor on the recall of several structural variant callers. We showed that most reported insertions exhibited characteristics that may interfere with their discovery: 63% were tandem repeat expansions, 38% contained homology larger than 10 bp within their breakpoint junctions and 70% were located in simple repeats. Consequently, the recall of short-read based variant callers was significantly lower for such insertions (6% for tandem repeats vs 56% for mobile element insertions). Simulations showed that the most impacting factor was the insertion type rather than the genomic context, with various difficulties being handled differently among the tested structural variant callers, and they highlighted the lack of sequence resolution for most insertion calls. CONCLUSIONS: Our results explain the low recall by pointing out several difficulty factors among the observed insertion features and provide avenues for improving SV caller algorithms and their combinations.

MinYS: mine your symbiont by targeted genome assembly in symbiotic communities.

Most metazoans are associated with symbionts. Characterizing the effect of a particular symbiont often requires getting access to its genome, which is usually done by sequencing the whole community. We present MinYS, a targeted assembly approach to assemble a particular genome of interest from such metagenomic data. First, taking advantage of a reference genome, a subset of the reads is assembled into a set of backbone contigs. Then, this draft assembly is completed using the whole metagenomic readset in a de novo manner. The resulting assembly is output as a genome graph, enabling different strains with potential structural variants coexisting in the sample to be distinguished. MinYS was applied to 50 pea aphid resequencing samples, with variable diversity in symbiont communities, in order to recover the genome sequence of its obligatory bacterial symbiont, Buchnera aphidicola. It was able to return high-quality assemblies (one contig assembly in 90% of the samples), even when using increasingly distant reference genomes, and to retrieve large structural variations in the samples. Because of its targeted essence, it outperformed standard metagenomic assemblers in terms of both time and assembly quality.

Role of the acquisition of a type 3 secretion system in the emergence of novel pathogenic strains of Xanthomonas.

Cases of emergence of novel plant-pathogenic strains are regularly reported that reduce the yields of crops and trees. However, the molecular mechanisms underlying such emergence are still poorly understood. The acquisition by environmental non-pathogenic strains of novel virulence genes by horizontal gene transfer has been suggested as a driver for the emergence of novel pathogenic strains. In this study, we tested such an hypothesis by transferring a plasmid encoding the type 3 secretion system (T3SS) and four associated type 3 secreted proteins (T3SPs) to the non-pathogenic strains of Xanthomonas CFBP 7698 and CFBP 7700, which lack genes encoding T3SS and any previously known T3SPs. The resulting strains were phenotyped on Nicotiana benthamiana using chlorophyll fluorescence imaging and image analysis. Wild-type, non-pathogenic strains induced a hypersensitive response (HR)-like necrosis, whereas strains complemented with T3SS and T3SPs suppressed this response. Such suppression depends on a functional T3SS. Amongst the T3SPs encoded on the plasmid, Hpa2, Hpa1 and, to a lesser extent, XopF1 collectively participate in suppression. Monitoring of the population sizes in planta showed that the sole acquisition of a functional T3SS by non-pathogenic strains impairs growth inside leaf tissues. These results provide functional evidence that the acquisition via horizontal gene transfer of a T3SS and four T3SPs by environmental non-pathogenic strains is not sufficient to make strains pathogenic. In the absence of a canonical effector, the sole acquisition of a T3SS seems to be counter-selective, and further acquisition of type 3 effectors is probably needed to allow the emergence of novel pathogenic strains.