Targeted next-generation sequencing based on PCR amplification involves pooling of hundreds to thousands of primers, for preamplification and subsequent parallel single/multiplex PCR. It is often necessary to allocate the set of primers into subpools, a common issue being potential cross-hybridization. For smaller numbers of primers, pool division can be done manually using trial and error to minimize potential hybridization, but this becomes inefficient and time consuming with increasing numbers of primer pairs. We developed PrimerPooler that automates swapping of primer pairs between any user-defined number of subpools to obtain combinations with low-potential interactions. PrimerPooler performs inter-/intra-primer hybridization analysis to identify the adverse interactions, as well as simultaneous mapping of all primers onto a genome sequence in a single run without requiring a prior index of the genome. This allows detection of overlapping primer pairs and allocation of these primer pairs into separate subpools where tiling approaches are used. Using PrimerPooler, 1153 primer pairs were assigned to three preamplification pools (388, 389 and 376 primer pairs each), then 144 subpools of six- to nine-plex PCR for Fluidigm Access Array PCR, followed by Illumina MiSeq sequencing. With optimized experimental protocols, an average of 3269 reads was achieved for the targeted regions, with 95% of targets covered by at least 50 reads, the minimal depth of reads for confident variant calling. PrimerPooler provides a fast and highly efficient stratification of primer pairs for targeted enrichment, thus ensuring representative amplification of the targeted sequences. PrimerPooler is also able to analyse degenerate primers, and is thus also useful for microbiological identification and related target sequencing.
Genome-wide association studies (GWASs) have identified a number of genetic risk loci associated with systemic sclerosis (SSc) and Crohn’s disease (CD), some of which confer susceptibility to both diseases. In order to identify new risk loci shared between these two immune-mediated disorders, we performed a cross-disease meta-analysis including GWAS data from 5,734 SSc patients, 4,588 CD patients and 14,568 controls of European origin. We identified 4 new loci shared between SSc and CD, IL12RB2, IRF1/SLC22A5, STAT3 and an intergenic locus at 6p21.31. Pleiotropic variants within these loci showed opposite allelic effects in the two analysed diseases and all of them showed a significant effect on gene expression. In addition, an enrichment in the IL-12 family and type I interferon signaling pathways was observed among the set of SSc-CD common genetic risk loci. In conclusion, through the first cross-disease meta-analysis of SSc and CD, we identified genetic variants with pleiotropic effects on two clinically distinct immune-mediated disorders. The fact that all these pleiotropic SNPs have opposite allelic effects in SSc and CD reveals the complexity of the molecular mechanisms by which polymorphisms affect diseases.
The identification of the genetic risk factors that could discriminate non- thrombotic from thrombotic antiphospholipid antibodies (aPLA) carriers will improve prognosis of these patients. Several human studies have shown the presence of aPLAs associated with atherosclerotic plaque, which is a known risk factor for thrombosis. Hence, in order to determine the implication of atherosclerosis in the risk of developing thrombosis in aPLA positive patients, we performed a genetic association study with 3 candidate genes, APOH, LDLR and PCSK9. For genetic association study we analyzed 190 aPLA carriers -100 with non-thrombotic events and 90 with thrombotic events- and 557 healthy controls. Analyses were performed by χ2 test and were corrected by false discovery rate. To evaluate the functional implication of the newly established susceptibility loci, we performed expression analyses in 86 aPLA carrier individuals (43 with thrombotic manifestations and 43 without it) and in 45 healthy controls. Our results revealed significant associations after correction in SNPs located in LDLR gene with aPLA carriers and thrombotic aPLA carriers, when compared with healthy controls. The most significant association in LDLR gene was found between SNP rs129083082 and aPLA carriers in recessive model (adjusted P-value = 2.55 x 10−3; OR = 2.18; 95%CI = 1.49–3.21). Furthermore, our work detected significant allelic association after correction between thrombotic aPLA carriers and healthy controls in SNP rs562556 located in PCSK9 gene (adjusted P-value = 1.03 x 10−2; OR = 1.60; 95%CI = 1.24–2.06). Expression level study showed significantly decreased expression level of LDLR gene in aPLA carriers (P-value <0.0001; 95%CI 0.16–2.10; SE 0.38–1.27) in comparison to the control group. Our work has identified LDLR gene as a new susceptibility gene associated with the development of thrombosis in aPLA carriers, describing for the first time the deregulation of LDLR expression in individuals with aPLAs. Besides, thrombotic aPLA carriers also showed significant association with PCSK9 gene, a regulator of LDLR plasma levels. These results highlight the importance of atherosclerotic processes in the development of thrombosis in patients with aPLA.