A pipeline for sample tagging of whole genome bisulfite sequencing data using genotypes of whole genome sequencing.

BACKGROUND: In large-scale high-throughput sequencing projects and biobank construction, sample tagging is essential to prevent sample mix-ups. Despite the availability of fingerprint panels for DNA data, little research has been conducted on sample tagging of whole genome bisulfite sequencing (WGBS) data. This study aims to construct a pipeline and identify applicable fingerprint panels to address this problem. RESULTS: Using autosome-wide A/T polymorphic single nucleotide variants (SNVs) obtained from whole genome sequencing (WGS) and WGBS of individuals from the Third China National Stroke Registry, we designed a fingerprint panel and constructed an optimized pipeline for tagging WGBS data. This pipeline used Bis-SNP to call genotypes from the WGBS data, and optimized genotype comparison by eliminating wildtype homozygous and missing genotypes, and retaining variants with identical genomic coordinates and reference/alternative alleles. WGS-based and WGBS-based genotypes called from identical or different samples were extensively compared using hap.py. In the first batch of 94 samples, the genotype consistency rates were between 71.01%-84.23% and 51.43%-60.50% for the matched and mismatched WGS and WGBS data using the autosome-wide A/T polymorphic SNV panel. This capability to tag WGBS data was validated among the second batch of 240 samples, with genotype consistency rates ranging from 70.61%-84.65% to 49.58%-61.42% for the matched and mismatched data, respectively. We also determined that the number of genetic variants required to correctly tag WGBS data was on the order of thousands through testing six fingerprint panels with different orders for the number of variants. Additionally, we affirmed this result with two self-designed panels of 1351 and 1278 SNVs, respectively. Furthermore, this study confirmed that using the number of genetic variants with identical coordinates and ref/alt alleles, or identical genotypes could not correctly tag WGBS data. CONCLUSION: This study proposed an optimized pipeline, applicable fingerprint panels, and a lower boundary for the number of fingerprint genetic variants needed for correct sample tagging of WGBS data, which are valuable for tagging WGBS data and integrating multi-omics data for biobanks.

Construction of competing endogenous RNA interaction network as prognostic markers in metastatic melanoma.

Malignant melanoma (MM) is a malignant tumor originating from melanocytes, with high aggressiveness, high metastasis and extremely poor prognosis. MM accounts for 4% of skin cancers and 80% of mortality, and the median survival of patients with metastatic melanoma is only about 6 months, with a five-year survival rate of less than 10%. In recent years, the incidence of melanoma has gradually increased and has become one of the serious diseases that endanger human health. Competitive endogenous RNA (ceRNA) is the main model of the mechanism by which long chain non-coding RNAs (lncRNAs) play a regulatory role in the disease. LncRNAs can act as a "sponge", competitively attracting small RNAs (micoRNAs; miRNAs), thus interfering with miRNA function, and affect the expression of target gene messenger RNAs (mRNAs), ultimately promoting tumorigenesis and progression. Bioinformatics analysis can identify potentially prognostic and therapeutically relevant differentially expressed genes in MM, finding lncRNAs, miRNAs and mRNAs that are interconnected through the ceRNA network, providing further insight into gene regulation and prognosis of metastatic melanoma. Weighted co-expression networks were used to identify lncRNA and mRNA modules associated with the metastatic phenotype, as well as the co-expression genes contained in the modules. A total of 17 lncRNAs, six miRNAs, and 11 mRNAs were used to construct a ceRNA interaction network that plays a regulatory role in metastatic melanoma patients. The prognostic risk model was used as a sorter to classify the survival prognosis of melanoma patients. Four groups of ceRNA interaction triplets were finally obtained, which miR-3662 might has potential implication for the treatment of metaststic melanoma patients, and futher experiments confirmed the regulating relationship and phenotype of this assumption. This study provides new targets to regulate metastatic process, predict metastatic potential and indicates that the miR-3662 can be used in the treatment of melanoma.