Recurrent neural network for predicting absence of heterozygosity from low pass WGS with ultra-low depth.

BACKGROUND: The absence of heterozygosity (AOH) is a kind of genomic change characterized by a long contiguous region of homozygous alleles in a chromosome, which may cause human genetic disorders. However, no method of low-pass whole genome sequencing (LP-WGS) has been reported for the detection of AOH in a low-pass setting of less than onefold. We developed a method, termed CNVseq-AOH, for predicting the absence of heterozygosity using LP-WGS with ultra-low sequencing data, which overcomes the sparse nature of typical LP-WGS data by combing population-based haplotype information, adjustable sliding windows, and recurrent neural network (RNN). We tested the feasibility of CNVseq-AOH for the detection of AOH in 409 cases (11 AOH regions for model training and 863 AOH regions for validation) from the 1000 Genomes Project (1KGP). AOH detection using CNVseq-AOH was also performed on 6 clinical cases with previously ascertained AOHs by whole exome sequencing (WES). RESULTS: Using SNP-based microarray results as reference (AOHs detected by CNVseq-AOH with at least a 50% overlap with the AOHs detected by chromosomal microarray analysis), 409 samples (863 AOH regions) in the 1KGP were used for concordant analysis. For 784 AOHs on autosomes and 79 AOHs on the X chromosome, CNVseq-AOH can predict AOHs with a concordant rate of 96.23% and 59.49% respectively based on the analysis of 0.1-fold LP-WGS data, which is far lower than the current standard in the field. Using 0.1-fold LP-WGS data, CNVseq-AOH revealed 5 additional AOHs (larger than 10 Mb in size) in the 409 samples. We further analyzed AOHs larger than 10 Mb, which is recommended for reporting the possibility of UPD. For the 291 AOH regions larger than 10 Mb, CNVseq-AOH can predict AOHs with a concordant rate of 99.66% with only 0.1-fold LP-WGS data. In the 6 clinical cases, CNVseq-AOH revealed all 15 known AOH regions. CONCLUSIONS: Here we reported a method for analyzing LP-WGS data to accurately identify regions of AOH, which possesses great potential to improve genetic testing of AOH.

Validation and depth evaluation of low-pass genome sequencing in prenatal diagnosis using 387 amniotic fluid samples.

BACKGROUND: Low-pass genome sequencing (LP GS) is an alternative to chromosomal microarray analysis (CMA). However, validations of LP GS as a prenatal diagnostic test for amniotic fluid are rare. Moreover, sequencing depth of LP GS in prenatal diagnosis has not been evaluated. OBJECTIVE: The diagnostic performance of LP GS was compared with CMA using 375 amniotic fluid samples. Then, sequencing depth was evaluated by downsampling. RESULTS: CMA and LP GS had the same diagnostic yield (8.3%, 31/375). LP GS showed all copy number variations (CNVs) detected by CMA and six additional variant of uncertain significance CNVs (>100 kb) in samples with negative CMA results; CNV size influenced LP GS detection sensitivity. CNV detection was greatly influenced by sequencing depth when the CNV size was small or the CNV was located in the azoospermia factor c (AZFc) region of the Y chromosome. Large CNVs were less affected by sequencing depth and more stably detected. There were 155 CNVs detected by LP GS with at least a 50% reciprocal overlap with CNVs detected by CMA. With 25 M uniquely aligned high-quality reads (UAHRs), the detection sensitivity for the 155 CNVs was 99.14%. LP GS using samples with 25 M UAHRs showed the same performance as LP GS using total UAHRs. Considering the detection sensitivity, cost and interpretation workload, 25 M UAHRs are optimal for detecting most aneuploidies and microdeletions/microduplications. CONCLUSION: LP GS is a promising, robust alternative to CMA in clinical settings. A total of 25 M UAHRs are sufficient for detecting aneuploidies and most microdeletions/microduplications.

Anthropogenic impacts on Late Holocene land-cover change and floristic biodiversity loss in tropical southeastern Asia.

Southern China and Southeast Asia witnessed some of their most significant economic and social changes relevant to human land use during the Late Holocene, including the intensification and spread of rice agriculture. Despite rice growth being associated with a number of earth systems impacts, how these changes transformed tropical vegetation in this region of immense endemic biodiversity remains poorly understood. Here, we compile a pollen dataset incorporating ∼150,000 identifications and 233 pollen taxa to examine past changes in floral biodiversity, together with a compilation of records of forest decline across the region using 14 pollen records spanning lowland to mountain sites. Our results demonstrate that the rise of intensive rice agriculture from approximately 2,000 y ago led not only to extensive deforestation but also to remarkable changes of vegetation composition and a reduction in arboreal diversity. Focusing specifically on the Tertiary relic tree species, the freshwater wetland conifer Glyptostrobus (Glyptostrobus pensilis), we demonstrate how key species that had survived changing environmental conditions across millions of years shrank in the face of paddy rice farming and human disturbance.

Noninvasive preeclampsia prediction using plasma cell-free RNA signatures.

BACKGROUND: Preeclampsia, especially preterm preeclampsia and early-onset preeclampsia, is a life-threating pregnancy disorder, and the heterogeneity and complexity of preeclampsia make it difficult to predict risk and to develop treatments. Plasma cell-free RNA carries unique information from human tissue and may be useful for noninvasive monitoring of maternal, placental, and fetal dynamics during pregnancy. OBJECTIVE: This study aimed to investigate various RNA biotypes associated with preeclampsia in plasma and to develop classifiers to predict preterm preeclampsia and early-onset preeclampsia before diagnosis. STUDY DESIGN: We performed a novel, cell-free RNA sequencing method termed polyadenylation ligation-mediated sequencing to investigate the cell-free RNA characteristics of 715 healthy pregnancies and 202 pregnancies affected by preeclampsia before symptom onset. We explored differences in the abundance of different RNA biotypes in plasma between healthy and preeclampsia samples and built preterm preeclampsia and early-onset preeclampsia prediction classifiers using machine learning methods. Furthermore, we validated the performance of the classifiers using the external and internal validation cohorts and assessed the area under the curve and positive predictive value. RESULTS: We detected 77 genes, including messenger RNA (44%) and microRNA (26%), that were differentially expressed in healthy mothers and mothers with preterm preeclampsia before symptom onset, which could separate participants with preterm preeclampsia from healthy samples and that played critical functional roles in preeclampsia physiology. We developed 2 classifiers for predicting preterm preeclampsia and early-onset preeclampsia before diagnosis based on 13 cell-free RNA signatures and 2 clinical features (in vitro fertilization and mean arterial pressure), respectively. Notably, both classifiers showed enhanced performance when compared with the existing methods. The preterm preeclampsia prediction model achieved 81% area under the curve and 68% positive predictive value in an independent validation cohort (preterm, n=46; control, n=151); the early-onset preeclampsia prediction model had an area under the curve of 88% and a positive predictive value of 73% in an external validation cohort (early-onset preeclampsia, n=28; control, n=234). Furthermore, we demonstrated that downregulation of microRNAs may play vital roles in preeclampsia through the upregulation of preeclampsia-relevant target genes. CONCLUSION: In this cohort study, a comprehensive transcriptomic landscape of different RNA biotypes in preeclampsia was presented and 2 advanced classifiers with substantial clinical importance for preterm preeclampsia and early-onset preeclampsia prediction before symptom onset were developed. We demonstrated that messenger RNA, microRNA, and long noncoding RNA can simultaneously serve as potential biomarkers of preeclampsia, holding the promise of prevention of preeclampsia in the future. Abnormal cell-free messenger RNA, microRNA, and long noncoding RNA molecular changes may help to elucidate the pathogenic determinants of preeclampsia and open new therapeutic windows to effectively reduce pregnancy complications and fetal morbidity.

Genetic characteristics of human papillomavirus type 16, 18, 52 and 58 in southern China.

Persistent infections of high-risk human papillomaviruses (HPVs) are the leading cause of cervical cancers. We collected cervical exfoliated cell samples from females in Changsha city, Hunan Province and obtained 338 viral genomes of four major HPV types, including HPV 16 (n = 82), 18 (n = 35), 52 (n = 121) and 58 (n = 100). The lineage/sublineage distribution of the four HPVs confirmed previous epidemiological reports, with the predominant prevailing sublineage as A4 (50%), A1 (37%) and A3 (13%) for HPV16, A1 (83%) for HPV18, B2 (86%) for HPV52 and A1 (65%), A3 (19%) and A2 (12%) for HPV58. We also identified two potentially novel HPV18 sublineages, i.e. A6 and A7. Virus mutation analysis further revealed the presence of HPV16 and HPV58 sublineages associated with potentially high oncogenicity. These findings expanded our knowledge of the HPV genetic diversity in China, providing valuable evidence to facilitate HPV DNA screening, vaccine effectiveness evaluation and control strategy development.

[Economic burden of Down syndrome patients and psychological and social discrimination to female caregivers in Changsha, China].

OBJECTIVE: To estimate the social and economic burden of Downs syndrome for patients and their families residing in Changsha, China. METHODS: An 160-item self-administered questionnaire was designed and distributed to the primary caregivers of the patients in March 2020. A total of 81 eligible participants had completed the questionnaire, among which 20 were excluded for incomplete data. A patient perspective was taken to estimate the economic burden of the disease. The social impact of the disease on the patient's family was evaluated through questions adapted from the Stanford Psychological Wellbeing (PWB) Scale. RESULTS: The estimated life-course cost of a Downs syndrome patient in Changsha is 4 985 659 RMB, with the patient and caregiver's loss of income taking the greater proportion. In addition, as the majority of the patients' primary caregivers, female caregivers experienced not only considerable financial hardship caused by the care provision, but also a significant amount of psychological pressure and social discrimination. CONCLUSION: Increased level of social welfare for the patients and social support for their female caregivers are essential for reducing economic burden and improving their quality of life in the area. In addition, prenatal screening and diagnosis for Downs syndrome are important for reducing both the social and economic burden of the disease by preventing its occurrence.

[Health economic evaluation of four prenatal screening strategies for Down syndrome in Changsha, China].

OBJECTIVE: To compare the clinical application and health economic values of non-invasive prenatal testing (NIPT) and second trimester serum screening (STSS). METHODS: A retrospective analysis was carried out on 54 026 singleton pregnant women undergoing NIPT and STSS from March 1, 2018 to December 31, 2019 in Changsha Maternal and Child Health Care Hospital. For pregnant women with high-risk results of NIPT, prenatal diagnosis and follow-up of pregnancy outcomes were conducted. The data was grouped to 4 screening models, and their cost-benefit was analyzed. RESULTS: The sensitivity, specificity and positive predictive value of NIPT were all higher than STSS. Screening models 1 to 4 have prevented the birth of 71, 29, 52 and 54 patients with Down syndrome, respectively. The safety index of screening models 1 to 4 were 0.0036, 0.3944, 02215 and 0.1281, respectively. When the price of NIPT was decreased to 600 RMB, the cost-benefit of the screening models 1 to 4 was 0.46, 0.65, 0.44 and 0.40 million RMB, respectively. CONCLUSION: NIPT has a better detection performance than STSS. When the price of NIPT is 600 RMB, screening model 1 has the best screening effect and the highest accuracy, safety index and health economical value.

T cell receptor ß repertoires as novel diagnostic markers for systemic lupus erythematosus and rheumatoid arthritis.

OBJECTIVE: T cell receptor (TCR) diversity determines the autoimmune responses in systemic lupus erythematosus (SLE) and rheumatoid arthritis (RA) and is closely associated with autoimmune diseases prognosis and prevention. However, the characteristics of variations in TCR diversity and their clinical significance is still unknown. Large series of patients must be studied in order to elucidate the effects of these variations. METHODS: Peripheral blood from 877 SLE patients, 206 RA patients and 439 healthy controls (HC) were amplified for the TCR repertoire and sequenced using a high-throughput sequencer. We have developed a statistical model to identify disease-associated TCR clones and diagnose autoimmune diseases. RESULTS: Significant differences were identified in variable (V), joining (J) and V-J pairing between the SLE or RA and HC groups. These differences can be utilised to discriminate the three groups with perfect accuracy (V: area under receiver operating curve > 0.99). One hundred ninety-eight SLE-associated and 53 RA-associated TCRs were identified and used for diseases classification by cross validation with high specificity and sensitivity. Disease-associated clones showed common features and high similarity between both autoimmune diseases. SLE displayed higher TCR heterogeneity than RA with several organ specific properties. Furthermore, the association between clonal expansion and the concentration of disease-associated clones with disease severity were identified, and pathogen-related TCRs were enriched in both diseases. CONCLUSIONS: These characteristics of the TCR repertoire, particularly the disease-associated clones, can potentially serve as biomarkers and provide novel insights for disease status and therapeutical targets in autoimmune diseases.

Modern pollen distribution in the northeastern Indian Ocean and its significance.

In order to provide a reference for reconstructing the paleoclimate of the northeastern Indian Ocean, 36 airborne pollen samples were analyzed using methods for airborne pollen, and 26 surface water samples were analyzed using a lab method for surface water. We found that little pollen is airborne over the Indian Ocean in spring, but airborne pollen types and concentrations can help to deduce paleomonsoon strength and direction. The conclusions included the following: (1) Pollen in the sediment was transported mainly via ocean currents instead of the early summer or spring wind. (2) Airborne pollen types and concentrations are proportional to the wind speed and inversely proportional to the pollen distance transported and depend on whether the wind is from the land or from the sea. If the wind is from the land, the pollen concentration is proportional to the angle between the wind direction and the coastline. (3) The pollen concentration in the sample collected from a water depth of 30-45 m is higher than in the samples collected from a depth of 5 m. The pollen concentration and salinity are higher in the equatorial area than in the Northern Hemisphere.

Nicotine-Induced Airway Smooth Muscle Cell Proliferation Involves TRPC6-Dependent Calcium Influx Via α7 nAChR.

BACKGROUND/AIMS: The proliferation of human bronchial smooth muscle cells (HBSMCs) is a key pathophysiological component of airway remodeling in chronic obstructive pulmonary disease (COPD) for which pharmacotherapy is limited, and only slight improvements in survival have been achieved in recent decades. Cigarette smoke is a well-recognized risk factor for COPD; however, the pathogenesis of cigarette smoke-induced COPD remains incompletely understood. This study aimed to investigate the mechanisms by which nicotine affects HBSMC proliferation. METHODS: Cell viability was assessed with a CCK-8 assay. Proliferation was measured by cell counting and EdU immunostaining. Fluorescence calcium imaging was performed to measure intracellular Ca2+ concentration ([Ca2+]i). RESULTS: The results showed that nicotine promotes HBSMC proliferation, which is accompanied by elevated store-operated calcium entry (SOCE), receptor-operated calcium entry (ROCE) and basal [Ca2+]i in HBSMCs. Moreover, we also confirmed that canonical transient receptor potential protein 6 (TRPC6) and α7 nicotinic acetylcholine receptor (α7 nAChR) are involved in nicotine-induced upregulation of cell proliferation. Furthermore, we verified that activation of the PI3K/Akt signaling pathway plays a pivotal role in nicotine-enhanced proliferation and calcium influx in HBSMCs. Inhibition of α7 nAChR significantly decreased Akt phosphorylation levels, and LY294002 inhibited the protein expression levels of TRPC6. CONCLUSION: Herein, these data provide compelling evidence that calcium entry via the α7 nAChR-PI3K/Akt-TRPC6 signaling pathway plays an important role in the physiological regulation of airway smooth muscle cell proliferation, representing an important target for augmenting airway remodeling.

[Endocytosis pathway and intracellular distribution of heparosan].

In this study, the endocytosis pathway of heparosan and its intracellular distribution were investigated in MCF-7 tumor cells and COS7 normal cells. The endocytosis inhibition and cellular probe location experiments showed that MCF-7 tumor cells took heparosan more efficiently and selectively than COS7 cells. The cellular uptake of heparosan was energy-dependent in both MCF-7 tumor cells and COS7 normal cells. Moreover, the major endocytosis pathway of heparosan into MCF-7 tumor cells was caveolin-mediated endocytosis and macropinocytosis. The internalized heparosan was mainly located in lysosomes of the cells.

Noninvasive prenatal testing (NIPT) in twin pregnancies with treatment of assisted reproductive techniques (ART) in a single center.

OBJECTIVE: The objective of the study is to report the performance of noninvasive prenatal testing (NIPT) in twin pregnancies after the treatment of assisted reproductive technology (ART). METHOD: In two years period, 565 pregnant women with ART twin pregnancies were prospectively tested by NIPT for screening for trisomy 21 (T21), 18 (T18), and 13 (T13) by sequencing cell-free DNA in maternal plasma. Positive NIPT results were confirmed by karyotyping, while negative results were interviewed after delivery. Pregnant decision based on NIPT and confirmation results was discussed during post-test counseling. RESULTS: In total of 565 cases, NIPT had a failure rate of 0.9% (5/565). Four cases of T21 were identified by NIPT and confirmed by karyotyping, resulting in 100% (95%CI 39.8%-100%) positive predictive value. Among 556 cases with NIPT negative results, 506 cases (91.0%) were confirmed by follow-up of postnatal phenotypes, while 33 cases (5.9%) had adverse pregnant outcomes with unconfirmed reasons because of the lack of cytogenetic samples. The remaining 17 cases (3.1%) refused follow-up. No false negative result was reported. CONCLUSION: With apparently high positive predictive value and low false positive rate, NIPT has the potential to be used as a good alternative approach of conventional prenatal screening at the first trimester in ART twin pregnancy. © 2016 John Wiley & Sons, Ltd.

Non-invasive prediction of preeclampsia using the maternal plasma cell-free DNA profile and clinical risk factors.

Background: Preeclampsia (PE) is a pregnancy complication defined by new onset hypertension and proteinuria or other maternal organ damage after 20 weeks of gestation. Although non-invasive prenatal testing (NIPT) has been widely used to detect fetal chromosomal abnormalities during pregnancy, its performance in combination with maternal risk factors to screen for PE has not been extensively validated. Our aim was to develop and validate classifiers that predict early- or late-onset PE using the maternal plasma cell-free DNA (cfDNA) profile and clinical risk factors. Methods: We retrospectively collected and analyzed NIPT data of 2,727 pregnant women aged 24-45 years from four hospitals in China, which had previously been used to screen for fetal aneuploidy at 12 + 0 ~ 22 + 6 weeks of gestation. According to the diagnostic criteria for PE and the time of diagnosis (34 weeks of gestation), a total of 143 early-, 580 late-onset PE samples and 2,004 healthy controls were included. The wilcoxon rank sum test was used to identify the cfDNA profile for PE prediction. The Fisher's exact test and Mann-Whitney U-test were used to compare categorical and continuous variables of clinical risk factors between PE samples and healthy controls, respectively. Machine learning methods were performed to develop and validate PE classifiers based on the cfDNA profile and clinical risk factors. Results: By using NIPT data to analyze cfDNA coverages in promoter regions, we found the cfDNA profile, which was differential cfDNA coverages in gene promoter regions between PE and healthy controls, could be used to predict early- and late-onset PE. Maternal age, body mass index, parity, past medical histories and method of conception were significantly differential between PE and healthy pregnant women. With a false positive rate of 10%, the classifiers based on the combination of the cfDNA profile and clinical risk factors predicted early- and late-onset PE in four datasets with an average accuracy of 89 and 80% and an average sensitivity of 63 and 48%, respectively. Conclusion: Incorporating cfDNA profiles in classifiers might reduce performance variations in PE models based only on clinical risk factors, potentially expanding the application of NIPT in PE screening in the future.

Massively parallel sequencing for chromosomal abnormality testing in trophectoderm cells of human blastocysts.

Preimplantation genetic diagnosis and screening are widely accepted for chromosomal abnormality identification to avoid transferring embryos with genetic defects. Massively parallel sequencing (MPS) is a rapidly developing approach for genome analysis with increasing application in clinical practice. The purpose of this study was to use MPS for identification of aneuploidies and unbalanced chromosomal rearrangements after blastocyst biopsy. Trophectoderm (TE) samples of 38 blastocysts from 16 in vitro fertilization cycles were subjected to analysis. Low-coverage whole genome sequencing was performed using the Illumina HiSeq2000 platform with a novel algorithm purposely created for chromosomal analysis. The efficiency of this MPS approach was estimated by comparing results obtained by an Affymetrix single-nucleotide polymorphism (SNP) array. Whole genome amplification (WGA) products of TE cells were detected by MPS, with an average of 0.07× depth and 5.5% coverage of the human genome. Twenty-six embryos (68.4%) were detected as euploid, while six embryos (15.8%) contained uniform aneuploidies. Four of these (10.5%) were with solely unbalanced chromosomal rearrangements, whereas the remaining two embryos (5.3%) showed both aneuploidies and unbalanced rearrangements. Almost all these results were confirmed by the SNP array, with the exception of one sample, where different sizes of unbalanced rearrangements were detected, possibly due to chromosomal GC bias in array analysis. Our study demonstrated MPS could be applied to accurately detect embryonic chromosomal abnormality with a flexible and cost-effective strategy and higher potential accuracy.

Molecular characterization of cold-responsive basic helix-loop-helix transcription factors MabHLHs that interact with MaICE1 in banana fruit.

Basic helix-loop-helix (bHLH) transcription factors (TFs) are ubiquitously involved in the response of higher plants to various abiotic stresses. However, little is known about bHLH TFs involved in the cold stress response in economically important fruits. Here, five novel full-length bHLH genes, designated as MabHLH1-MabHLH5, were isolated and characterized from banana fruit. Gene expression profiles revealed that MabHLH1/2/4 were induced by cold stress and methyl jasmonate (MeJA) treatment. Transient assays in tobacco BY2 protoplasts showed that MabHLH1/2/4 promoters were activated by cold stress and MeJA treatments. Moreover, protein-protein interaction analysis demonstrated that MabHLH1/2/4 not only physically interacted with each other to form hetero-dimers in the nucleus, but also interacted with an important upstream component of cold signaling MaICE1, with different interaction domains at their N-terminus. These results indicate that banana fruit cold-responsive MabHLHs may form a big protein complex in the nucleus with MaICE1. Taken together, our findings advance our understanding of the possible involvement of bHLH TFs in the regulatory network of ICE-CBF cold signaling pathway.

The U-shaped association between hospitalization time and fall incidence in inpatients using publicly available data: A cross-sectional study in Japan.

AIM: To examine the association between the hospitalization time and fall incidence. DESIGN: A secondary analysis using the Dryad Digital Repository public database. METHODS: Data were extracted from the Fukushima Medical University Hospital cohort study between August 2008 and September 2009. The final analytic sample included 8,598 participants, 156 of who fell. The risk of fall incidents according to hospitalization time was estimated using logistic proportional hazards models, and restricted cubic splines with four knots model were developed. RESULTS: The median hospitalization time was 9.00 (4.00, 17.00) days. The incidence of falls was 1.81% (N = 156). A U-shaped association between the hospitalization time and fall incidence, with an inflextion point of 8 days. We found a decreasing fall incidence as the hospitalization time increased from 0 to 8 days (OR 0.72 [0.62 ~ 0.83], p < .001); beyond 8 days, the fall incidence increased as the hospitalization time increased (OR 1.06 [1.04 ~ 1.09]).

Populational pan-ethnic screening panel enabled by deep whole genome sequencing.

Birth defect is a global threat to the public health systems. Mitigating neonatal anomalies is hampered by elusive molecular mechanisms of pathogenic mutations and poor subsequent translation into preventative measures. Applying appropriate strategies in China to promote reproductive health is particularly challenging, as the Chinese population compromises complex genomic diversity due to the inclusion of many ethnic groups with distinct genetic backgrounds. To investigate and evaluate the feasibility of implementing a pan-ethnic screening strategy, and guide future reproductive counselling, high-quality variants associated with autosome recessive (AR) diseases derived from the largest publicly available cohort of the Chinese population were re-analysed using a bottom-up approach. The analyses of gene carrier rates (GCRs) across distinct ethnic groups revealed that substantial heterogeneity existed potentially due to diverse evolutionary selection. The sampling population, sequencing coverage and underlying population structure contributed to the differential variants observed between ChinaMAP and the East Asian group in gnomAD. Beyond characteristics of GCR, potential druggable targets were additionally explored according to genomic features and functional roles of investigated genes, demonstrating that phase separation could be a therapeutic target for autosomal recessive diseases. A further examination of estimated GCR across ethnic groups indicated that most genes shared by at least two populations could be utilised to direct the design of a pan-ethnic screening application once sequencing and interpreting costs become negligible. To this end, a list of autosomal recessive disease genes is proposed based on the prioritised rank of GCR to formulate a tiered screening strategy.

Accuracy and depth evaluation of clinical low pass genome sequencing in the detection of mosaic aneuploidies and CNVs.

BACKGROUND: Low-pass genome sequencing (LP GS) has shown distinct advantages over traditional methods for the detection of mosaicism. However, no study has systematically evaluated the accuracy of LP GS in the detection of mosaic aneuploidies and copy number variants (CNVs) in prenatal diagnosis. Moreover, the influence of sequencing depth on mosaicism detection of LP GS has not been fully evaluated. METHODS: To evaluate the accuracy of LP GS in the detection of mosaic aneuploidies and mosaic CNVs, 27 samples with known aneuploidies and CNVs and 1 negative female sample were used to generate 6 simulated samples and 21 virtual samples, each sample contained 9 different mosaic levels. Mosaic levels were simulated by pooling reads or DNA from each positive sample and the negative sample according to a series of percentages (ranging from 3 to 40%). Then, the influence of sequencing depth on LP GS in the detection of mosaic aneuploidies and CNVs was evaluated by downsampling. RESULTS: To evaluate the accuracy of LP GS in the detection of mosaic aneuploidies and CNVs, a comparative analysis of mosaic levels was performed using 6 simulated samples and 21 virtual samples with 35 M million (M) uniquely aligned high-quality reads (UAHRs). For mosaic levels > 30%, the average difference (detected mosaic levels vs. theoretical mosaic levels) of 6 mosaic CNVs in simulated samples was 4.0%, and the average difference (detected mosaic levels vs. mosaic levels of Y chromosome) of 6 mosaic aneuploidies and 15 mosaic CNVs in virtual samples was 2.7%. Furthermore, LP GS had a higher detection rate and accuracy for the detection of mosaic aneuploidies and CNVs of larger sizes, especially mosaic aneuploidies. For depth evaluation, the results of LP GS in downsampling samples were compared with those of LP GS using 35 M UAHRs. The detection sensitivity of LP GS for 6 mosaic aneuploidies and 15 mosaic CNVs in virtual samples increased with UAHR. For mosaic levels > 30%, the total detection sensitivity reached a plateau at 30 M UAHRs. With 30 M UAHRs, the total detection sensitivity was 99.2% for virtual samples. CONCLUSIONS: We demonstrated the accuracy of LP GS in mosaicism detection using simulated data and virtual samples, respectively. Thirty M UAHRs (single-end 35 bp) were optimal for LP GS in the detection of mosaic aneuploidies and most mosaic CNVs larger than 1.48 Mb (Megabases) with mosaic levels > 30%. These results could provide a reference for laboratories that perform clinical LP GS in the detection of mosaic aneuploidies and CNVs.

Test development, optimization and validation of a WGS pipeline for genetic disorders.

BACKGROUND: With advances in massive parallel sequencing (MPS) technology, whole-genome sequencing (WGS) has gradually evolved into the first-tier diagnostic test for genetic disorders. However, deployment practice and pipeline testing for clinical WGS are lacking. METHODS: In this study, we introduced a whole WGS pipeline for genetic disorders, which included the entire process from obtaining a sample to clinical reporting. All samples that underwent WGS were constructed using polymerase chain reaction (PCR)-free library preparation protocols and sequenced on the MGISEQ-2000 platform. Bioinformatics pipelines were developed for the simultaneous detection of various types of variants, including single nucleotide variants (SNVs), insertions and deletions (indels), copy number variants (CNVs) and balanced rearrangements, mitochondrial (MT) variants, and other complex variants such as repeat expansion, pseudogenes and absence of heterozygosity (AOH). A semiautomatic pipeline was developed for the interpretation of potential SNVs and CNVs. Forty-five samples (including 14 positive commercially available samples, 23 laboratory-held positive cell lines and 8 clinical cases) with known variants were used to validate the whole pipeline. RESULTS: In this study, a whole WGS pipeline for genetic disorders was developed and optimized. Forty-five samples with known variants (6 with SNVs and Indels, 3 with MT variants, 5 with aneuploidies, 1 with triploidy, 23 with CNVs, 5 with balanced rearrangements, 2 with repeat expansions, 1 with AOHs, and 1 with exon 7-8 deletion of SMN1 gene) validated the effectiveness of our pipeline. CONCLUSIONS: This study has been piloted in test development, optimization, and validation of the WGS pipeline for genetic disorders. A set of best practices were recommended using our pipeline, along with a dataset of positive samples for benchmarking.

Genomic Sequencing as a First-Tier Screening Test and Outcomes of Newborn Screening.

Importance: Newborn screening via biochemical tests is in use worldwide. The availability of genetic sequencing has allowed rapid screening for a substantial number of monogenic disorders. However, the outcomes of this strategy have not been evaluated in a general newborn population. Objective: To evaluate the outcomes of applying gene panel sequencing as a first-tier newborn screening test. Design, Setting, and Participants: This cohort study included newborns who were prospectively recruited from 8 screening centers in China between February 21 and December 31, 2021. Neonates with positive results were followed up before July 5, 2022. Exposures: All participants were concurrently screened using dried blood spots. The screen consisted of biochemical screening tests and a targeted gene panel sequencing test for 128 conditions. The biochemical and genomic tests could both detect 43 of the conditions, whereas the other 85 conditions were screened solely by the gene panel. Main Outcomes and Measures: The primary outcomes were the number of patients detected by gene panel sequencing but undetected by the biochemical test. Results: This study prospectively recruited 29 601 newborns (15 357 [51.2%] male). The mean (SD) gestational age was 39.0 (1.5) weeks, and the mean (SD) birth weight was 3273 (457) g. The gene panel sequencing screened 813 infants (2.7%; 95% CI, 2.6%-2.9%) as positive. By the date of follow-up, 402 infants (1.4%; 95% CI, 1.2%-1.5%) had been diagnosed, indicating the positive predictive value was 50.4% (95% CI, 50.0%-53.9%). The gene panel sequencing identified 59 patients undetected by biochemical tests, including 20 patients affected by biochemically and genetically screened disorders and 39 patients affected by solely genetically screened disorders, which translates into 1 out of every 500 newborns (95% CI, 1/385-1/625) benefiting from the implementation of gene panels as a first-tier screening test. Conclusions and Relevance: In this cohort study, the use of gene panel sequencing in a general newborn population as a first-tier screening test improved the detection capability of traditional screening, providing an evidence-based suggestion that it could be considered as a crucial method for first-tier screening.