Establishing national reference materials for genetic testing of cytochrome P450.

OBJECTIVES: Reference materials for in-vitro diagnostic reagents play a critical role in determining the quality of reagents and ensuring the accuracy of clinical test results. This study aimed to establish a national reference material (NRM) for detecting cytochrome P450 (CYP) genes related to drug metabolism by screening databases on the Chinese population to identify CYP gene polymorphism characteristics. METHODS: To prepare the NRM, we used DNA extracted from healthy human immortalized B lymphoblastoid cell lines as the raw material. Samples of these cell lines were obtained from the Chinese Population PGx Gene Polymorphism Biobank. Further, we used Sanger sequencing, next-generation sequencing, and commercial assay kits to validate the polymorphic genotypes. RESULTS: Among the CYP superfamily genes, we confirmed 24 riboswitch loci related to drug metabolism, with evidence levels of 1A, 2A, 3, and 4. We confirmed the polymorphic loci and validated their genotypes using various sequencing techniques. Our results were consistent with the polymorphism information of samples obtained from the biobank, thus demonstrating high precision and stability of the established NRM. CONCLUSION: An NRM (360 056-202 201) for CYP genetic testing covering 24 loci related to drug metabolism was established and approved to assess in-vitro diagnostic reagents containing CYP family gene polymorphisms and perform clinical inter-room quality evaluations.

Development of preimplantation genetic testing for monogenic reference materials using next-generation sequencing.

OBJECTIVE: Preimplantation genetic testing for monogenic disorders (PGT-M) has been used for over 20 years to detect many serious genetic conditions. However, there is still a lack of reference materials (RMs) to validate the test performance during the development and quality control of PGT-M. METHOD: Sixteen thalassemia cell lines from four thalassemia families were selected to establish the RMs. Each family consisted of parents with heterozygous mutations for α- and/or ß-thalassemia and two children, at least one of whom carried a homozygous thalassemia mutation (proband). The RM panel consisted of 12 DNA samples (parents and probands in 4 families) and 4 simulated embryos (cell lines constructed from blood samples from the four nonproband children). Four accredited genetics laboratories that offer verification of thalassemia samples were invited to evaluate the performance of the RM panel. Furthermore, the stability of the RMs was determined by testing after freeze‒thaw cycles and long-term storage. RESULTS: PGT-M reference materials containing 12 genome DNA (gDNA) reference materials and 4 simulated embryo reference materials for thalassemia testing were successfully established. Next-generation sequencing was performed on the samples. The genotypes and haplotypes of all 16 PGT-M reference materials were concordant across the four labs, which used various testing workflows. These well-characterized PGT-M reference materials retained their stability even after 3 years of storage. CONCLUSION: The establishment of PGT-M reference materials for thalassemia will help with the standardization and accuracy of PGT-M in clinical use.

Harmonization of distributed multi-center analysis based on dried blood spot reference materials supporting the screening of neonatal inherited metabolic disorders.

BACKGROUND: The standardization of quantification data is critical for ensuring the reliability and measurement traceability in the screening of neonatal inherited metabolic disorders. However, the availability of national certified reference materials is limited in China. METHODS: In this study, we developed a series of dried blood spot (DBS) reference materials containing 9 amino acids (AA) and 10 acylcarnitines (AC) for neonatal screening. Four levels of the reference materials were measured with tandem mass spectrometry (MS/MS) by seven laboratories using different commercial In Vitro Diagnostic Device (IVD) kits. Then, 100 clinical samples were measured using both derivatization and non-derivatization methods by the same laboratory. RESULTS: We found high homogeneity and stability at all levels of the reference materials, with the coefficient of variation (CV) of the analytes less than 15%. These reference materials can be used to assess the testing capabilities of different laboratories. Our test also revealed that the correction factors (CF) calculated by the reference materials, along with clinical samples, could increase the consistency for different kits. CONCLUSION: The DBS reference materials proposed in this study provide reliability for the harmonization in multi-center analysis for the screening of neonatal inherited metabolic disorders. And applying our correction method for the screening could improve the data consistency of the DBS samples prepared by different methods.

Nanopore Third-Generation Sequencing for Comprehensive Analysis of Hemoglobinopathy Variants.

BACKGROUND: Oxford Nanopore Technology (ONT) third-generation sequencing (TGS) is a versatile genetic diagnostic platform. However, it is nonetheless challenging to prepare long-template libraries for long-read TGS, particularly the ONT method for analysis of hemoglobinopathy variants involving complex structures and occurring in GC-rich and/or homologous regions. METHODS: A multiplex long PCR was designed to prepare library templates, including the whole-gene amplicons for HBA2/1, HBG2/1, HBD, and HBB, as well as the allelic amplicons for targeted deletions and special structural variations. Library construction was performed using long-PCR products, and sequencing was conducted on an Oxford Nanopore MinION instrument. Genotypes were identified based on integrative genomics viewer (IGV) plots. RESULTS: This novel long-read TGS method distinguished all single nucleotide variants and structural variants within HBA2/1, HBG2/1, HBD, and HBB based on the whole-gene sequence reads. Targeted deletions and special structural variations were also identified according to the specific allelic reads. The result of 158 α-/ß-thalassemia samples showed 100% concordance with previously known genotypes. CONCLUSIONS: This ONT TGS method is high-throughput, which can be used for molecular screening and genetic diagnosis of hemoglobinopathies. The strategy of multiplex long PCR is an efficient strategy for library preparation, providing a practical reference for TGS assay development.

Development and characterization of reference materials for EGFR, KRAS, NRAS, BRAF, PIK3CA, ALK, and MET genetic testing.

BACKGROUND: Along with the dramatic development of molecular diagnostic testing for the detection of oncogene variations, reference materials (RMs) have become increasingly important in performance evaluation of genetic testing. OBJECTIVE: In this study, we built a set of RMs for genetic testing based on next-generation sequencing (NGS). METHOD: Solid tumor tissues were selected as the samples of RMs for preparation. NGS was used to determine and validate the variants and the mutation frequency in DNA samples. Digital PCR was used to determine the copy numbers of RNA samples. The performance of the RMs was validated by six laboratories. RESULTS: Thirty common genetic alterations were designed based on these RMs. RMs consisted of a positive reference, a limit of detection reference, and a negative reference. The validation results confirmed the performance of the RMs. CONCLUSION: These RMs may be an attractive tool for the development, validation, and quality monitoring of molecular genetic testing.

Integrative analyses of maternal plasma cell-free DNA nucleosome footprint differences reveal chromosomal aneuploidy fetuses gene expression profile.

BACKGROUND: Chromosomal aneuploidy is the most common birth defect. However, the developmental mechanism and gene expression profile of fetuses with chromosomal aneuploidy are relatively unknown, and the maternal immune changes induced by fetal aneuploidy remain unclear. The inability to obtain the placenta multiple times in real-time is a bottleneck in research on aneuploid pregnancies. Plasma cell-free DNA (cfDNA) carries the gene expression profile information of its source cells and may be used to evaluate the development of fetuses with aneuploidy and the immune changes induced in the mother owing to fetal aneuploidy. METHODS: Here, we carried out whole-genome sequencing of the plasma cfDNA of 101 pregnant women carrying a fetus with trisomy (trisomy 21, n = 42; trisomy 18, n = 28; trisomy 13, n = 31) based on non-invasive prenatal testing (NIPT) screening and 140 normal pregnant women to identify differential genes according to the cfDNA nucleosome profile in the region around the transcription start sites (TSSs). RESULTS: The plasma cfDNA promoter profiles were found to differ between aneuploid and euploid pregnancies. A total of 158 genes with significant differences were identified, of which 43 genes were upregulated and 98 genes were downregulated. Functional enrichment and signaling pathway analysis were performed based on Gene Ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) databases found that these signal pathways were mainly related to the coordination of developmental signals during embryonic development, the control of cell growth and development, regulation of neuronal survival, and immune regulation, such as the MAPK, Hippo, TGF-ß, and Rap1 signaling pathways, which play important roles in the development of embryonic tissues and organs. Furthermore, based on the results of differential gene analysis, a total of 14 immune-related genes with significant differences from the ImmPort database were collected and analyzed. These significantly different immune genes were mainly associated with the maintenance of embryonic homeostasis and normal development. CONCLUSIONS: These results suggest that the distribution characteristics of cfDNA nucleosomes in maternal plasma can be used to reflect the status of fetal development and changes of the immune responses in trisomic pregnancies. Overall, our findings may provide research ideas for non-invasive detection of the physiological and pathological states of other diseases.

Simulated confined placental mosaicism proportion (SCPMP) based on cell-free fetal DNA fraction enrichment can reduce false-positive results in non-invasive prenatal testing.

OBJECTIVE: To decrease the false-positive rate of NIPT using cell-free fetal DNA (cffDNA) fraction enrichment and the simulated confined placental mosaicism proportion (SCPMP) threshold application via cffDNA quantification. METHOD: Using a cffDNA enrichment method, 303 plasma samples with positive NIPT results (Z-score > 3.0; 200 true-positive and 103 false-positive cases) were re-sequenced. A method to calculate the SCPMP based on the quantified cffDNA fraction was developed; the SCPMP threshold between true- and false-positive NIPT results was determined and used for re-analyses. RESULTS: With enrichment, the fetal fraction of the 303 samples was 26.9 ± 8.4%, compared to 11.0 ± 3.2% without enrichment. The optimized threshold method with double determination using the Z-value-defined SCPMP can reduce the false-positive rates for trisomies 21, 18, and 13 by 87%, 80%, and 88.9%, respectively. CONCLUSION: Our optimized method can decrease the false-positive rate of NIPT results.

Comprehensive evaluation of BRCA1/2 variant interpretation ability among laboratories in China.

High-quality interpretation of BRCA1/2 variants plays a critical role in the clinical practice of precision medicine. However, a comprehensive system to evaluate the quality and accuracy of variant interpretation has yet to be established. This study investigates the performance of an interpretation system in evaluating the capacities of BRCA1/2 interpretation among distinct laboratories in China. The evaluation system is based on a reference database that contains 750 different variants in BRCA1/2 Evaluation was performed among 41 laboratories in China. We classified their performance into five levels. Only level A was considered qualified. This level allows for a 0.3% error rate for clinical decision-related misinterpretation; 26 of 41 laboratories (63%) met the qualified standard, while 7 laboratories were at levels D and E, which indicated egregious mistakes and systemic problems in variant interpretation. Due to strict quality demands, the interpretation of several variants was amended, which largely influenced the quality rate. The number of qualified laboratories would decrease from 26 to 17 if those incorrect recommended interpretations were not corrected. This evaluation system provides a potential approach for standardisation of variant interpretation and lowers the discordance of variant interpretation between different laboratories. A well-designed interpretation ability evaluation is essential to evaluate the interpretation level of laboratories before they provide service in real-world clinical settings.

Robust Benchmark Structural Variant Calls of An Asian Using State-of-the-art Long-read Sequencing Technologies.

The importance of structural variants (SVs) for human phenotypes and diseases is now recognized. Although a variety of SV detection platforms and strategies that vary in sensitivity and specificity have been developed, few benchmarking procedures are available to confidently assess their performances in biological and clinical research. To facilitate the validation and application of these SV detection approaches, we established an Asian reference material by characterizing the genome of an Epstein-Barr virus (EBV)-immortalized B lymphocyte line along with identified benchmark regions and high-confidence SV calls. We established a high-confidence SV callset with 8938 SVs by integrating four alignment-based SV callers, including 109× Pacific Biosciences (PacBio) continuous long reads (CLRs), 22× PacBio circular consensus sequencing (CCS) reads, 104× Oxford Nanopore Technologies (ONT) long reads, and 114× Bionano optical mapping platform, and one de novo assembly-based SV caller using CCS reads. A total of 544 randomly selected SVs were validated by PCR amplification and Sanger sequencing, demonstrating the robustness of our SV calls. Combining trio-binning-based haplotype assemblies, we established an SV benchmark for identifying false negatives and false positives by constructing the continuous high-confidence regions (CHCRs), which covered 1.46 gigabase pairs (Gb) and 6882 SVs supported by at least one diploid haplotype assembly. Establishing high-confidence SV calls for a benchmark sample that has been characterized by multiple technologies provides a valuable resource for investigating SVs in human biology, disease, and clinical research.

A Deep-Learning Pipeline for TSS Coverage Imputation From Shallow Cell-Free DNA Sequencing.

Cell-free DNA (cfDNA) serves as a footprint of the nucleosome occupancy status of transcription start sites (TSSs), and has been subject to wide development for use in noninvasive health monitoring and disease detection. However, the requirement for high sequencing depth limits its clinical use. Here, we introduce a deep-learning pipeline designed for TSS coverage profiles generated from shallow cfDNA sequencing called the Autoencoder of cfDNA TSS (AECT) coverage profile. AECT outperformed existing single-cell sequencing imputation algorithms in terms of improvements to TSS coverage accuracy and the capture of latent biological features that distinguish sex or tumor status. We built classifiers for the detection of breast and rectal cancer using AECT-imputed shallow sequencing data, and their performance was close to that achieved by high-depth sequencing, suggesting that AECT could provide a broadly applicable noninvasive screening approach with high accuracy and at a moderate cost.

BDdb: a comprehensive platform for exploration and utilization of birth defect multi-omics data.

BACKGROUND: Birth defects pose a major challenge to infant health. Thus far, however, the causes of most birth defects remain cryptic. Over the past few decades, considerable effort has been expended on disclosing the underlying mechanisms related to birth defects, yielding myriad treatises and data. To meet the increasing requirements for data resources, we developed a freely accessible birth defect multi-omics database (BDdb, http://t21omics.cngb.org ) consisting of multi-omics data and potential disease biomarkers. RESULTS: In total, omics datasets from 136 Gene Expression Omnibus (GEO) Series records, including 5245 samples, as well as 869 biomarkers of 22 birth defects in six different species, were integrated into the BDdb. The database provides a user-friendly interface for searching, browsing, and downloading data of interest. The BDdb also enables users to explore the correlations among different sequencing methods, such as chromatin immunoprecipitation sequencing (ChIP-Seq) and RNA sequencing (RNA-Seq) from different studies, to obtain the information on gene expression patterns from diverse aspects. CONCLUSION: To the best of our knowledge, the BDdb is the first comprehensive database associated with birth defects, which should benefit the diagnosis and prevention of birth defects.

Influence of low tumor content on tumor mutational burden estimation by whole-exome sequencing and targeted panel sequencing.

BACKGROUND: Tumor mutational burden (TMB) is a promising biomarker for stratifying patient subpopulation who would benefit from immune checkpoint blockade (ICB) therapies. Although great efforts have been made for standardizing TMB measurement, mutation calling and TMB quantification can be challenging in samples with low tumor content including liquid biopsies. The effect of varying tumor content on TMB estimation by different assay methods has never been systematically investigated. METHOD: We established a series of reference standard DNA samples derived from 11 pairs of tumor-normal matched human cell lines across different cancer types. Each tumor cell line was mixed with its matched normal at 0% (control), 1%, 2%, 5%, and 10% mass-to-mass ratio to mimic the clinical samples with low tumor content. TMB of these reference standards was evaluated by both ∼1000× whole-exome sequencing (wesTMB) and targeted panel sequencing (psTMB) at four different vendors. Both regression and classification analyses of TMB were performed for theoretical investigation and clinical practice purposes. RESULTS: Linear regression model was established that demonstrated in silico psTMB determined by regions of interest (ROI) as a great representative of wesTMB based on TCGA dataset. It was also true in our reference standard samples as the predicted psTMB interval based on the observed wesTMB captured the intended 90% of the in silico psTMB values. Although ∼1000× deep WES was applied, reference standard samples with less than 5% of tumor proportions are below the assay limit of detection (LoD) of wesTMB quantification. However, predicted wesTMB based on observed psTMB accurately classify (>0.97 AUC) for TMB high and low patient stratification even in samples with 2% of tumor content, which is more clinically relevant, as TMB determination should be a qualitative assay for TMB high and low patient classification. One targeted panel sequencing vendor using an optimized blood psTMB pipeline can further classify TMB status accurately (>0.82 AUC) in samples with only 1% of tumor content. CONCLUSIONS: We developed a linear model to establish the quantitative correlation between wesTMB and psTMB. A set of DNA reference standards was produced in aid to standardize TMB measurements in samples with low tumor content across different targeted sequencing panels. This study is a significant contribution aiming to harmonize TMB estimation and extend its future application in clinical samples with low tumor content including liquid biopsy.

[Study on Proficiency Testing of Determination Potassium in Serum].

OBJECTIVE: To evaluate the capacity of laboratories participated in the proficiency testing (PT) of determination potassium in serum and improve the quality of testing, and put forward technical suggestions for unsatisfied laboratories. METHODS: According to the requirements of CNAS related documents, the homogeneity and stability of the real PT sample were evaluated by one-way ANOVA and t test, respectively. The values of real PT samples were assigned by reference method which was used in PT results assay. It is required that the deviation of value of real PT samples (code:2, 3, 5) between the measured value and the assigned value shall be within ±15.0%. The precision of values for all samples should not be greater than 3.0%. RESULTS: All the laboratories submitted valid data according to the requirements. Only one laboratory did not meet the requirements, and the satisfaction rate was 90.9%. CONCLUSIONS: The ability of most of laboratories are accurate and reliable.

The Setup and Application of Reference Material in Sequencing-Based Noninvasive Prenatal Testing.

INTRODUCTION: The sequencing-based noninvasive prenatal testing (NIPT) has been successfully integrated into clinical practice and facilitated the early detection of fetal chromosomal anomalies. However, a comprehensive reference material to evaluate and quality control NIPT services from different NIPT providers remains unavailable. METHODS: In this study, we established a set of NIPT reference material consisting of 192 simulated samples. Most of the potential factors influencing the accuracy of NIPT, such as fetal fraction, mosaicism, and interfering substances, were included in the reference material. We compared the performance of chromosomal abnormalities detection on 3 widely used sequencers (NextSeq 500, BGISEQ-500, and Ion Proton) based on the reference material. RESULTS: All 3 sequencers provided highly accurate and reliable results to samples with ≥3.5% fetal fractions and high percentage of mosaicism. CONCLUSIONS: The established reference material can serve as a universal standard quality control for the current and new-coming NIPT providers based on various sequencers.

Construction of a reference material panel for detecting KRAS/NRAS/EGFR/BRAF/MET mutations in plasma ctDNA.

BACKGROUND: The absence of high-quality next-generation sequencing (NGS) reference material (RM) has impeded the clinical use of liquid biopsies with plasma cell-free DNA (cfDNA) in China. OBJECTIVE: This study aimed to develop a national RM panel for external quality assessment and performance evaluation during kit registration of non-small-cell lung cancer (NSCLC)-related Kirsten rat sarcoma viral oncogene (KRAS)/neuroblastoma ras oncogene (NRAS)/epidermal growth factor receptor (EGFR)/B-type Raf kinase (BRAF)/mesenchymal-epithelial transition factor (MET) genetic assays using plasma circulating tumor DNA (ctDNA). METHODS: Mutation cell lines detected by NGS and validated by Sanger sequencing were selected to establish the RM. Cell line genomic DNA was sheared and used to spike basal plasma cfDNA at 10% concentration. Then, the calibration accuracy was determined by four sequencing platforms. Average values were adopted and diluted to 0.1%, 0.3%, 1% and 3% concentrations with basal plasma as the RM panel. Then, five manufacturers were invited to evaluate the performance of the RM panel. RESULTS: 20 cell lines with 23 clinically important mutations were selected, including six mutations in KRAS, two mutations in NRAS, three in BRAF, four in phosphatidylinositol-4,5-bisphosphate 3-kinase catalytic subunit alpha (PIK3CA), six in EGFR, one EGFR Gain (4-5 copy) and one MET Gain (2-5 copy). The RM panel consisted of 87 samples, including these 21 mutations at four concentrations (0.1%, 0.3%, 1% and 3%), one MET gain, one EGFR gain and one wild type. The detection rate was 100% for the 3%, 1% and 0.3% samples at all five companies. For the 0.1% concentration, 15 samples had inconsistent results, but at least three companies had correct results for each mutation. CONCLUSION: RM for a KRAS/NRAS/EGFR/BRAF/MET mutation panel for plasma ctDNA was developed, which will be essential for quality control of the performance of independent laboratories.

Mutation Screening of mtDNA Combined Targeted Exon Sequencing in a Cohort With Suspected Hereditary Optic Neuropathy.

Purpose: Leber hereditary optic neuropathy (LHON) and autosomal dominant optic atrophy (ADOA) are the two commonest forms of hereditary optic neuropathy. The aim of this study was to comprehensively investigate the incidence and spectrum of mutations in patients with suspected hereditary optic neuropathy by combining mitochondrial DNA (mtDNA) genome-wide and targeted exon sequencing. Methods: A cohort of 1101 subjects were recruited to participate in the study, comprising 177 families (177 probands and their family members, a total of 537 subjects, including 254 patients) and 164 sporadic cases with suspected hereditary optic neuropathy, and 400 unrelated control subjects for genetic analysis: all subjects (including control subjects) underwent a comprehensive ophthalmologic examination and were subjected to sequencing analysis of mtDNA genome-wide and targeted exon. Overall, targeted exon sequencing was used to screen 792 genes associated with common hereditary eye diseases, and the mtDNA genome-wide were screened by next-generation sequencing. Results: We found variants detected in 168 (40.2%, 168/418) of the 418 patients screened. Among these, 132 cases (78.6%, 132/168) were detected with known LHON disease-causing mtDNA variants; 40 cases (23.8%, 40/168) were detected with nuclear DNA (ntDNA) variants, which included 36 cases (21.4%, 36/168) with detected OPA1 mutations, 4 patients (2.4%, 4/168) with detected OPA3 mutations, and 2 patients (1.2%, 2/168) with detected TMEM126A homozygous mutation. Coexistence variation (mtDNA/mtDNA [n = 16], ntDNA/ntDNA [n = 4], mtDNA/ntDNA [n = 7]) was found in 27 patients (16.4%, 27/165), including mtDNA/ntDNA coexistence variation that was detected in seven patients. Among these ntDNA mutations, 38 distinct disease-causing variants, including autosomal recessive heterozygous mutations, were detected, which included 22 novel variants and two de novo variants. Total haplogroup distribution showed that 34.5% (29/84) and 28.6% (24/84) of the affected subjects with m.11778G>A belonged to haplogroup D and M, with a high frequency of subhaplogroups D4, D5, and M7. Conclusions: The LHON-mtDNA mutations are the commonest genetic defects in this Chinese cohort, followed by the OPA1 mutations. To our knowledge, this is the first comprehensive study of LHON, ADOA, and autosomal recessive optic atrophy combined with mtDNA genome-wide and targeted exon sequencing, as well as haplogroup analysis, in a large cohort of Chinese patients with suspected hereditary optic neuropathy. Our findings provide a powerful basis for genetic counseling in patients with suspected hereditary optic neuropathy. Translational Relevance: We applied mtDNA genome-wide sequencing combined with panel-based targeted exon sequencing to explore the pathogenic variation spectrum and genetic characteristics of patients with suspected hereditary optic neuropathy, providing a comprehensive research strategy for clinical assistant diagnosis, treatment, and genetic counseling.

Detection of chromosomal abnormalities in spontaneous miscarriage by low­coverage next­generation sequencing.

Chromosomal abnormalities (CAs) can cause spontaneous miscarriage and increase the incidence of subsequent pregnancy loss and other complications. Presently, CAs are detected mainly by array comparative genomic hybridization (CGH) and single nucleotide polymorphism microarrays. The present study developed a low­coverage next­generation sequencing method to detect CAs in spontaneous miscarriage and assess its clinical performance. In total, 1,401 patients who had experienced an abortion were enrolled in the present study and divided into two groups. In group I, 437 samples that had been previously validated by array CGH were used to establish a method to detect CAs using a semiconductor sequencing platform. In group II, 964 samples, which were not verified, were assessed using established methods with respect to clinical significance. Copy number variant (CNV)­positive and euploidy samples were verified by array CGH and short tandem repeat profiling, respectively, based on quantitative fluorescent PCR. The low­coverage sequencing method detected CNVs >1 Mb in length and a total of 3.5 million unique reads. Similar results to array CGH were obtained in group I, except for six CNVs <1 Mb long. In group II, there were 341 aneuploidies, 195 CNVs, 25 mosaicisms and 403 euploidies. Overall, among the 1,401 abortion samples, there were 536 aneuploidies, 263 CNVs, 34 mosaicisms, and 568 euploidies. Trisomies were present in all autosomal chromosomes. The most common aneuploidies were T16, monosomy X, T22, T15, T21 and T13. Furthermore, one tetrasomy 21, one CNV associated with Wolf­Hirschhorn syndrome, one associated with DiGeorge syndrome and one associated with both Prader­Willi and Angelman syndromes were identified. These four cases were confirmed by short tandem repeat profiling and array CGH. Quantitative fluorescent PCR revealed nine polyploidy samples. The present method demonstrated equivalent efficacy to that of array CGH in detecting CNVs >1 Mb, with advantages of requiring less input DNA and lower cost.

An integrated Asian human SNV and indel benchmark established using multiple sequencing methods.

Sequencing technologies have been rapidly developed recently, leading to the breakthrough of sequencing-based clinical diagnosis, but accurate and complete genome variation benchmark would be required for further assessment of precision medicine applications. Despite the human cell line of NA12878 has been successfully developed to be a variation benchmark, population-specific variation benchmark is still lacking. Here, we established an Asian human variation benchmark by constructing and sequencing a stabilized cell line of a Chinese Han volunteer. By using seven different sequencing strategies, we obtained ~3.88 Tb clean data from different laboratories, hoping to reach the point of high sequencing depth and accurate variation detection. Through the combination of variations identified from different sequencing strategies and different analysis pipelines, we identified 3.35 million SNVs and 348.65 thousand indels, which were well supported by our sequencing data and passed our strict quality control, thus should be high confidence variation benchmark. Besides, we also detected 5,913 high-quality SNVs which had 969 sites were novel and  located in the high homologous regions supported by long-range information in both the co-barcoding single tube Long Fragment Read (stLFR) data and PacBio HiFi CCS data. Furthermore, by using the long reads data (stLFR and HiFi CCS), we were able to phase more than 99% heterozygous SNVs, which helps to improve the benchmark to be haplotype level. Our study provided comprehensive sequencing data as well as the integrated variation benchmark of an Asian derived cell line, which would be valuable for future sequencing-based clinical development.

Development of a genomic DNA reference material panel for thalassemia genetic testing.

INTRODUCTION: Thalassemia is one of the most common autosomal recessive inherited diseases worldwide, and it is also highly prevalent and variable in southern China. Various types of genetic testing technologies have been developed for diagnosis and screening of thalassemia. Characterized genomic DNA reference materials (RMs) are necessary for assay development, validation, proficiency testing, and quality assurance. However, there are no publicly available RMs for thalassemia genetic testing as yet. METHODS: To address the need for the publicly available DNA RMs for thalassemia genetic testing, the National Institutes for Food and Drug Control and the China National GeneBank established 32 new cell lines with three wild-type genotypes and 29 distinct genotypes of thalassemia which account for approximately 90% thalassemia carriers in China. The genomic DNA of 32 cell lines was characterized by four clinical genetic testing laboratories using different genetic testing methods and technology platforms. RESULTS: The genotyping results are concordant among four laboratories. In addition, the results of stability test demonstrated that the genotypes of these DNA samples are not influenced by preanalytical conditions such as long-term exposure to high-temperature (37°C) environment and repeated freeze-thawing. CONCLUSION: We developed the first national panel of 32 genomic DNA RMs which are renewable and publicly available for the quality assurance of various genetic testing methods and will facilitate research and development in thalassemia genetic testing.

Panel-based targeted exome sequencing reveals novel candidate susceptibility loci for age-related cataracts in Chinese Cohort.

BACKGROUND: Age-related cataracts (ARC) is the most common blinding eye disease worldwide, and its incidence tend to become younger. However, the relationship between genetic factors and mechanisms is not fully understood. The aim of the study was to further clarify the relationship between ARC and genetic mechanisms in East Asian populations and to elucidate the pathogenesis. METHODS: The study collected 191 sporadic cataracts and 208 healthy people from the eastern provinces of China, with an average age of about 60 years. All participants were subjected to a comprehensive ophthalmic clinical examination and peripheral blood samples were collected and their genomic DNA was extracted. Mutations were screened among 792 candidate genes to enhance understanding of the disease through targeted capture and high-throughput sequencing. RESULTS: We identified novel candidate susceptibility gene, which may serve as a potential susceptibility factor leading to an increase in the incidence of age-related cataracts. Three novel loci are associated with age-related cataracts significant significance: rs129882 in DBH (p = 5.27E-07, odds ratio = 3.9), rs1800280 in DMD (p = 2.85E-06, odds ratio = 1.4) and rs2871776 in ATP13A2 (p = 4.18E-05, odds ratio = 0.04). Gene-gene interaction analysis revealed that the most significant interactions between genes include the interaction between DBH and TUB (rs17847537 in TUB, rs129882 in DBH, p-value = 2.12E-14), and the interaction between DBH and DMD (rs1800280 in DMD, rs129882 in DBH, p-value = 2.12E-14). Pathway analysis shows that the most significant processes are concentrated in response to light stimulation (adjusted p-Value = 5.56E-03), response to radiation (adjusted P-Value = 5.56E-03), abiotic stimulus (adjusted p-Value = 5.56E-03). eQTL analysis shows that DBH rs129882 could regulate the expression of DBH mRNA in various tissues including retina. CONCLUSION: Our study indicates rs129882 and rs1800280 loci are associated with age-related cataracts, which enlarge the gene map of age-related cataracts.