Large-scale DNA sequencing identifies rare variants associated with Systemic Lupus Erythematosus susceptibility in known risk genes.

The identification of rare genetic variants associated to Systemic Lupus Erythematosus (SLE) could also help to understand the pathogenic mechanisms at the basis of the disease. In this study we have analyzed a cohort of 200 Italian SLE patients in order to explore the rare protein-coding variants in five genes (TNFAIP3, STAT4, IL10, TRAF3IP2, and HCP5) already investigated for commons variants found associated in our previous studies. Genomic DNA of 200 SLE patients was sequenced by whole exome sequencing. The identified variants were filtered by frequency and evaluated by in silico predictions. Allelic association analysis was performed with standard Fisher's exact test. Introducing a cutoff at MAF < 0.01, a total of 19 rare variants were identified. Seven of these variants were ultra-rare (MAF < 0.001) and six were absent in the GnomAD database. For TNFAIP3 gene, the variant c.A1939C was observed in 4 SLE patients and it is located in a region enriched in phosphorylation sites and affects the predict affinity of specific kinases. In TRAF3IP2 gene, we observed 5 different rare variants, including the novel variant c.G410A, located in the region that mediates interaction with TRAF6, and therefore a possible risk factor for SLE development. In STAT4 gene, we identified 6 different rare variants. Among these, three missense variants decrease the stability of this protein. Moreover, 3 novel rare variants were detected in 3 SLE patients. In particular, c.A767T variant was predicted as damaging by six prediction tools. Concluding, we have observed that even in genes whose common variability is associated with SLE susceptibility, it is possible to identify rare variants that could have a strong effect in the disease development and could therefore allow a better understanding of the functional domain involved.

The complete diploid reference genome of RPE-1 identifies human phased epigenetic landscapes.

Comparative analysis of recent human genome assemblies highlights profound sequence divergence that peaks within polymorphic loci such as centromeres. This raises the question about the adequacy of relying on human reference genomes to accurately analyze sequencing data derived from experimental cell lines. Here, we generated the complete diploid genome assembly for the human retinal epithelial cells (RPE-1), a widely used non-cancer laboratory cell line with a stable karyotype, to use as matched reference for multi-omics sequencing data analysis. Our RPE1v1.0 assembly presents completely phased haplotypes and chromosome-level scaffolds that span centromeres with ultra-high base accuracy (>QV60). We mapped the haplotype-specific genomic variation specific to this cell line including t(Xq;10q), a stable 73.18 Mb duplication of chromosome 10 translocated onto the microdeleted chromosome X telomere t(Xq;10q). Polymorphisms between haplotypes of the same genome reveals genetic and epigenetic variation for all chromosomes, especially at centromeres. The RPE-1 assembly as matched reference genome improves mapping quality of multi-omics reads originating from RPE-1 cells with drastic reduction in alignments mismatches compared to using the most complete human reference to date (CHM13). Leveraging the accuracy achieved using a matched reference, we were able to identify the kinetochore sites at base pair resolution and show unprecedented variation between haplotypes. This work showcases the use of matched reference genomes for multiomics analyses and serves as the foundation for a call to comprehensively assemble experimentally relevant cell lines for widespread application.