Novel copy number variations and phenotypes of infantile epileptic spasms syndrome.

We summarize the copy number variations (CNVs) and phenotype spectrum of infantile epileptic spasms syndrome (IESS) in a Chinese cohort. The CNVs were identified by genomic copy number variation sequencing. The CNVs and clinical data were analyzed. 74 IESS children with CNVs were enrolled. 35 kinds of CNVs were identified. There were 11 deletions and 5 duplications not reported previously in IESS, including 2 CNVs not reported in epilepsy. 87.8% were de novo, 9.5% were inherited from mother and 2.7% from father. Mosaicism occurred in one patient with Xq21.31q25 duplication. 16.2% (12/74) were 1p36 deletion, and 20.3% (15/74) were 15q11-q13 duplication. The age of seizure onset ranged from 17 days to 24 months. Seizure types included epileptic spasms, focal seizures, tonic seizures, and myoclonic seizures. All patients displayed developmental delay. Additional features included craniofacial anomaly, microcephaly, congenital heart defects, and hemangioma. 29.7% of patients were seizure-free for more than 12 months, and 70.3% still had seizures after trying 2 or more anti-seizure medications. In conclusion, CNVs is a prominent etiology of IESS. 1p36 deletion and 15q duplication occurred most frequently. CNV detection should be performed in patients with IESS of unknown causes, especially in children with craniofacial anomalies and microcephaly.

A capture-based method of prenatal cell-free DNA screening for autosomal recessive non-syndromic hearing loss.

OBJECTIVE: This study aimed to develop and validate a prenatal cell-free DNA (cfDNA) screening method that uses capture-based enrichment to genotype fetal autosomal recessive disorders. This method was applied in pregnancies at high risk of autosomal recessive non-syndromic hearing loss (ARNSHL) to assess its accuracy and effectiveness. METHODS: This assay measured the allele counts in both white blood cell DNA and cfDNA from the blood samples of pregnant women using a capture-based next-generation sequencing method. It then applied a binomial model to infer the fetal genotypes with the maximum likelihood. Ninety-four pregnant couples that were carriers of variants of ARNSHL in GJB2 or SLC26A4 were enrolled. The fetal genotypes deduced using this screening method were compared with the results of genetic diagnosis using amniocentesis. RESULTS: Of the 94 couples, 65 carried more than one variant, resulting in 170 single-nucleotide polymorphism (SNP) loci to be inferred in the fetuses. Of the 170 fetal SNP genotypes, 150 (88.2%) had high confidence calls and 139 (92.7%) of these matched the genotypes obtained by amniocentesis result. Out of the remaining 20 (11.8%) cases with low-confidence calls, only 14 (70.0%) were concordant with genetic diagnosis using amniocentesis. The concordance rate was 100% for sites where the maternal genotype was wild-type homozygous. The discordance was site-biased, with each locus showing a consistent direction of discordance. Genetic diagnosis identified a total of 19 wild-type homozygotes, 46 heterozygotes, 19 compound heterozygotes, and 10 pathogenic homozygotes. This screening method correctly genotyped 81.9% (77/94) of fetuses and demonstrated a sensitivity of 89.7% and a specificity of 89.2% for correctly identifying ARNSHL. CONCLUSION: This capture-based method of prenatal screening by cfDNA demonstrated strong potential for fetal genotyping of autosomal recessive disorders.

A new and improved method of library preparation for non-invasive prenatal testing: plasma to library express technology.

OBJECTIVES: This study aims to develop a novel library preparation method, plasma to library express technology (PLET), to construct next-generation sequencing (NGS) libraries directly from plasma without cell-free DNA (cfDNA) isolation. METHODS: Peripheral blood samples (600) were obtained from a retrospective cohort of 300 pregnant women prior to invasive diagnostic testing. The samples were subsequently distributed between library preparation methodologies, with 300 samples prepared by PLET and 300 by conventional methods for non-invasive prenatal testing (NIPT) to screen for common trisomies using low-pass whole genome next generation sequencing. RESULTS: NIPT conducted on PLET libraries demonstrated comparable metrics to libraries prepared using conventional methods, including 100% sensitivity and specificity. CONCLUSIONS: Our study demonstrates the potential utility of PLET in the clinical setting and highlights its significant advantages, including dramatically reduced process complexity and markedly decreased turnaround time.

Clinical spectrum and genetic variations of LMNA-related muscular dystrophies in a large cohort of Chinese patients.

BACKGROUND: LMNA-related muscular dystrophy is caused by mutations in LMNA gene. We aimed to identify genetic variations and clinical features in a large cohort of Chinese patients with LMNA mutations in an attempt to establish genotype-phenotype correlation. METHODS: The clinical presentations of patients with LMNA-related muscular dystrophy were recorded using retrospective and prospective cohort study. LMNA mutation analysis was performed by Sanger sequencing or next-generation sequencing. Mosaicism was detected by personal genome machine amplicon deep sequencing for mosaicism. RESULTS: Eighty-four patients were identified to harbour LMNA mutations. Forty-one of those were diagnosed with LMNA-related congenital muscular dystrophy (L-CMD), 32 with Emery-Dreifuss muscular dystrophy (EDMD) and 11 with limb-girdle muscular dystrophy type 1B (LGMD1B). We identified 21 novel and 29 known LMNA mutations. Two frequent mutations were identified: c.745C>T and c.1357C>T. A correlation between the location of mutation and the clinical phenotype was observed: mutations affecting the head and coil 2A domains mainly occurred in L-CMD, while the coil 2B and Ig-like domains mainly related to EDMD and LGMD1B. We found somatic mosaicism in one parent of four probands. Muscle biopsies revealed 11 of 20 biopsied L-CMD exhibited inflammatory changes, and muscle cell ultrastructure showed abnormal nuclear morphology. CONCLUSIONS: Our detailed clinical and genetic analysis of 84 patients with LMNA-related muscular dystrophy expands clinical spectrum and broadens genetic variations caused by LMNA mutations. We identified 21 novel and 29 known LMNA mutations and found two frequent mutations. A correlation between the location of mutation and the clinical severity was observed. Preliminary data suggested that low-dose corticosteroid treatment may be effective.

Somatic LINE-1 retrotransposition in cortical neurons and non-brain tissues of Rett patients and healthy individuals.

Mounting evidence supports that LINE-1 (L1) retrotransposition can occur postzygotically in healthy and diseased human tissues, contributing to genomic mosaicism in the brain and other somatic tissues of an individual. However, the genomic distribution of somatic human-specific LINE-1 (L1Hs) insertions and their potential impact on carrier cells remain unclear. Here, using a PCR-based targeted bulk sequencing approach, we profiled 9,181 somatic insertions from 20 postmortem tissues from five Rett patients and their matched healthy controls. We identified and validated somatic L1Hs insertions in both cortical neurons and non-brain tissues. In Rett patients, somatic insertions were significantly depleted in exons-mainly contributed by long genes-than healthy controls, implying that cells carrying MECP2 mutations might be defenseless against a second exonic L1Hs insertion. We observed a significant increase of somatic L1Hs insertions in the brain compared with non-brain tissues from the same individual. Compared to germline insertions, somatic insertions were less sense-depleted to transcripts, indicating that they underwent weaker selective pressure on the orientation of insertion. Our observations demonstrate that somatic L1Hs insertions contribute to genomic diversity and MeCP2 dysfunction alters their genomic patterns in Rett patients.

Distinctive types of postzygotic single-nucleotide mosaicisms in healthy individuals revealed by genome-wide profiling of multiple organs.

Postzygotic single-nucleotide mosaicisms (pSNMs) have been extensively studied in tumors and are known to play critical roles in tumorigenesis. However, the patterns and origin of pSNMs in normal organs of healthy humans remain largely unknown. Using whole-genome sequencing and ultra-deep amplicon re-sequencing, we identified and validated 164 pSNMs from 27 postmortem organ samples obtained from five healthy donors. The mutant allele fractions ranged from 1.0% to 29.7%. Inter- and intra-organ comparison revealed two distinctive types of pSNMs, with about half originating during early embryogenesis (embryonic pSNMs) and the remaining more likely to result from clonal expansion events that had occurred more recently (clonal expansion pSNMs). Compared to clonal expansion pSNMs, embryonic pSNMs had higher proportion of C>T mutations with elevated mutation rate at CpG sites. We observed differences in replication timing between these two types of pSNMs, with embryonic and clonal expansion pSNMs enriched in early- and late-replicating regions, respectively. An increased number of embryonic pSNMs were located in open chromatin states and topologically associating domains that transcribed embryonically. Our findings provide new insights into the origin and spatial distribution of postzygotic mosaicism during normal human development.

[Consensus on technological standards for non-invasive prenatal screening of pathogenic copy number variations by high-throughput sequencing of maternal plasma cell-free DNA].

Genomic disorders caused by pathogenic copy number variation (pCNV) have proven to underlie a significant proportion of birth defects. With technological advance, improvement of bioinformatics analysis procedure, and accumulation of clinical data, non-invasive prenatal screening of pCNV (NIPS-pCNV) by high-throughput sequencing of maternal plasma cell-free DNA has been put to use in clinical settings. Specialized standards for clinical application of NIPS-pCNV are required. Based on the discussion, 10 pCNV-associated diseases with well-defined conditions and 5 common chromosomal aneuploidy syndromes are recommended as the target of screening in this consensus. Meanwhile, a standardized procedure for NIPS-pCNV is also provided, which may facilitate propagation of this technique in clinical settings.

A novel NGS-based microsatellite instability (MSI) status classifier with 9 loci for colorectal cancer patients.

BACKGROUND: With the recent emergence of immune checkpoint inhibitors, microsatellite instability (MSI) status has become an important biomarker for immune checkpoint blockade therapy. There are growing technical demands for the integration of different genomic alterations profiling including MSI analysis in a single assay for full use of the limited tissues. METHODS: Tumor and paired control samples from 64 patients with primary colorectal cancer were enrolled in this study, including 14 MSI-high (MSI-H) cases and 50 microsatellite stable (MSS) cases determined by MSI-PCR. All the samples were sequenced by a customized NGS panel covering 2.2 MB. A training dataset of 28 samples was used for selection of microsatellite loci and a novel NGS-based MSI status classifier, USCI-msi, was developed. NGS-based MSI status, single nucleotide variant (SNV) and tumor mutation burden (TMB) were detected for all patients. Most of the patients were also independently detected by immunohistochemistry (IHC) staining. RESULTS: A 9-loci model for detecting microsatellite instability was able to correctly predict MSI status with 100% sensitivity and specificity compared with MSI-PCR, and 84.3% overall concordance with IHC staining. Mutations in cancer driver genes (APC, TP53, and KRAS) were dispersed in MSI-H and MSS cases, while BRAF p.V600E and frameshifts in TCF7L2 gene occurred only in MSI-H cases. Mismatch repair (MMR)-related genes are highly mutated in MSI-H samples. CONCLUSION: We established a new NGS-based MSI classifier, USCI-msi, with as few as 9 microsatellite loci for detecting MSI status in CRC cases. This approach possesses 100% sensitivity and specificity, and performed robustly in samples with low tumor purity.

Inferring fetal fractions from read heterozygosity empowers the noninvasive prenatal screening.

PURPOSE: Fetal fraction (FF) is the percent of cell-free DNA (cfDNA) in the mother's peripheral blood that is of fetal origin, which plays a pivotal role in noninvasive prenatal screening (NIPS). We present a method that can reliably estimate FFs by examining autosome single-nucleotide polymorphisms (SNPs). METHODS: Even at a very low sequencing depth, there are plenty of SNPs covered by more than one read. At those SNPs, we define read heterozygosity and demonstrate that the percent of read heterozygosity is a function of FF, which allows FF to be inferred. RESULTS: We first demonstrated the effectiveness of our method in inferring FF. Then we used the inferred FF as an informative alternative prior to computing Bayes factors to test for aneuploidy, and observed better power than the Z-test. In analysis of clinical samples, we were able to identify female-male twins thanks to the accurate FF inference. CONCLUSION: Knowing FF improves efficacy of NIPS. It brings a powerful Bayesian method, allows "no call" for samples with small FFs, renders screening for XXY syndrome simpler, and permits an adaptive design to sequence at a higher depth for samples with small FFs.

Correction: Inferring fetal fractions from read heterozygosity empowers the noninvasive prenatal screening.

An amendment to this paper has been published and can be accessed via a link at the top of the paper.

Mosaicism and incomplete penetrance of PCDH19 mutations.

BACKGROUND: Mutations in the PCDH19 gene have mainly been reported in female patients with epilepsy. To date, PCDH19 mutations have been reported in hundreds of females and only in 10 mosaic male epileptic patients with mosaicism. OBJECTIVE: We aimed to investigate the occurrence of mosaic PCDH19 mutations in 42 families comprising at least one patient with PCDH19-related epilepsy. METHODS: Two male patients with mosaic PCDH19 variants were identified using targeted next-generation sequencing. Forty female patients with PCDH19 variants were identified by Sanger sequencing and Multiple Ligation Probe Amplification (MLPA). Microdroplet digital PCR was used to quantify the mutant allelic fractions (MAFs) in 20 families with PCDH19 variants. RESULTS: Five mosaic individuals, four males and one female, were identified in total. Mosaic variant was confirmed in multiple somatic tissues from one male patient and in blood from the other male patient. Among 22 female patients harbouring a newly occurred PCDH19 variant identified by Sanger sequencing and MLPA, Sanger sequencing revealed two mosaic fathers (9%, 2/22), one with two affected daughters and the other with an affected child. Two asymptomatic mosaic fathers were confirmed as gonosomal mosaicism, with MAFs ranging from 4.16% to 37.38% and from 1.27% to 19.13%, respectively. In 11 families with apparent de novo variants, 1 female patient was identified as a mosaic with a blood MAF of 26.72%. CONCLUSION: Our study provides new insights into phenotype-genotype correlations in PCDH19 related epilepsy and the finding of high-frequency mosaicism has important implications for genetic counselling.

Targeted resequencing of 358 candidate genes for autism spectrum disorder in a Chinese cohort reveals diagnostic potential and genotype-phenotype correlations.

Autism spectrum disorder (ASD) is a childhood neuropsychiatric disorder with a complex genetic architecture. The diagnostic potential of a targeted panel of ASD genes has only been evaluated in small cohorts to date and is especially understudied in the Chinese population. Here, we designed a capture panel with 358 genes (111 syndromic and 247 nonsyndromic) for ASD and sequenced a Chinese cohort of 539 cases evaluated with the Autism Diagnostic Interview-Revised (ADI-R) and the Autism Diagnostic Observation Schedule (ADOS) as well as 512 controls. ASD cases were found to carry significantly more ultra-rare functional variants than controls. A subset of 78 syndromic and 54 nonsyndromic genes was the most significantly associated and should be given high priority in the future screening of ASD patients. Pathogenic and likely pathogenic variants were detected in 9.5% of cases. Variants in SHANK3 and SHANK2 were the most frequent, especially in females, and occurred in 1.2% of cases. Duplications of 15q11-13 were detected in 0.8% of cases. Variants in CNTNAP2 and MEF2C were correlated with epilepsy/tics in cases. Our findings reveal the diagnostic potential of ASD genetic panel testing and new insights regarding the variant spectrum. Genotype-phenotype correlations may facilitate the diagnosis and management of ASD.

Genomic mosaicism in the pathogenesis and inheritance of a Rett syndrome cohort.

PURPOSE: To determine the role of mosaicism in the pathogenesis and inheritance of Rett and Rett-like disorders. METHODS: We recruited 471 Rett and Rett-like patients. Panel-sequencing targeting MECP2, CDKL5, and FOXG1 was performed. Mosaicism was quantified in 147 patients by a Bayesian genotyper. Candidates were validated by amplicon sequencing and digital PCR. Germline mosaicism of 21 fathers with daughters carrying pathogenic MECP2 variants was further quantified. RESULTS: Pathogenic variants of MECP2/CDKL5/FOXG1 were found in 324/471 (68.7%) patients. Somatic MECP2 mosaicism was confirmed in 5/471 (1.1%) patients, including 3/18 males (16.7%) and 2/453 females (0.4%). Three of the five patients with somatic MECP2 mosaicism had mosaicism at MECP2-Arg106. Germline MECP2 mosaicism was detected in 5/21 (23.8%) fathers. CONCLUSION: This is the first systematic screening of somatic and paternal germline MECP2 mosaicism at a cohort level. Our findings indicate that somatic MECP2 mosaicism contributes directly to the pathogenicity of Rett syndrome, especially in male patients. MECP2-Arg106 might be a mosaic hotspot. The high proportion of paternal germline MECP2 mosaicism indicates an underestimated mechanism underlying the paternal origin bias of MECP2 variants. Finally, this study provides an empirical foundation for future studies of genetic disorders caused by de novo variations of strong paternal origin.

Correction: Genomic mosaicism in the pathogenesis and inheritance of a Rett syndrome cohort.

The second author Jiarui Li is now listed as a co-first author according to her contribution to this paper. The list of authors who contributed equally now reads: Qingping Zhang, Xiaoxu Yang, Jiaping Wang, and Jiarui Li. This has now been corrected in both the PDF and HTML versions of the Article.

ATP1A3 mosaicism in families with alternating hemiplegia of childhood.

Alternating hemiplegia of childhood (AHC) is a rare and severe neurodevelopmental disorder characterized by recurrent hemiplegic episodes. Most AHC cases are sporadic and caused by de novo ATP1A3 pathogenic variants. In this study, the aim was to identify the origin of ATP1A3 pathogenic variants in a Chinese cohort. In 105 probands including 101 sporadic and 4 familial cases, 98 patients with ATP1A3 pathogenic variants were identified, and 96.8% were confirmed as de novo. Micro-droplet digital polymerase chain reaction was applied for detecting ATP1A3 mosaicism in 80 available families. In blood samples, four asymptomatic parents, including two paternal and two maternal, and one proband with a milder phenotype were identified as mosaicism. Six (7.5%) parental mosaicisms were identified in multiple tissues, including four previously identified in blood and two additional cases identified from paternal sperms. Mosaicism was identified in multiple tissues with varied mutant allele fractions (MAFs, 0.03%-33.03%). The results suggested that MAF of mosaicism may be related to phenotype severity. This is the first systematic report of ATP1A3 mosaicism in AHC and showed mosaicism as an unrecognized source of previously considered "de novo" AHC. Identifying ATP1A3 mosaicism provides more evidence for estimating recurrence risk and has implications in genetic counseling of AHC.

Postzygotic single-nucleotide mosaicisms contribute to the etiology of autism spectrum disorder and autistic traits and the origin of mutations.

The roles and characteristics of postzygotic single-nucleotide mosaicisms (pSNMs) in autism spectrum disorders (ASDs) remain unclear. In this study of the whole exomes of 2,361 families in the Simons Simplex Collection, we identified 1,248 putative pSNMs in children and 285 de novo SNPs in children with detectable parental mosaicism. Ultra-deep amplicon resequencing suggested a validation rate of 51%. Analyses of validated pSNMs revealed that missense/loss-of-function (LoF) pSNMs with a high mutant allele fraction (MAF≥ 0.2) contributed to ASD diagnoses (P = 0.022, odds ratio [OR] = 5.25), whereas missense/LoF pSNMs with a low MAF (MAF<0.2) contributed to autistic traits in male non-ASD siblings (P = 0.033). LoF pSNMs in parents were less likely to be transmitted to offspring than neutral pSNMs (P = 0.037), and missense/LoF pSNMs in parents with a low MAF were transmitted more to probands than to siblings (P = 0.016, OR = 1.45). We estimated that pSNMs in probands or de novo mutations inherited from parental pSNMs increased the risk of ASD by approximately 6%. Adding pSNMs into the transmission and de novo association test model revealed 13 new ASD risk genes. These results expand the existing repertoire of genes involved in ASD and shed new light on the contribution of genomic mosaicisms to ASD diagnoses and autistic traits.

Amplicon Resequencing Identified Parental Mosaicism for Approximately 10% of "de novo" SCN1A Mutations in Children with Dravet Syndrome.

The majority of children with Dravet syndrome (DS) are caused by de novo SCN1A mutations. To investigate the origin of the mutations, we developed and applied a new method that combined deep amplicon resequencing with a Bayesian model to detect and quantify allelic fractions with improved sensitivity. Of 174 SCN1A mutations in DS probands which were considered "de novo" by Sanger sequencing, we identified 15 cases (8.6%) of parental mosaicism. We identified another five cases of parental mosaicism that were also detectable by Sanger sequencing. Fraction of mutant alleles in the 20 cases of parental mosaicism ranged from 1.1% to 32.6%. Thirteen (65% of 20) mutations originated paternally and seven (35% of 20) maternally. Twelve (60% of 20) mosaic parents did not have any epileptic symptoms. Their mutant allelic fractions were significantly lower than those in mosaic parents with epileptic symptoms (P = 0.016). We identified mosaicism with varied allelic fractions in blood, saliva, urine, hair follicle, oral epithelium, and semen, demonstrating that postzygotic mutations could affect multiple somatic cells as well as germ cells. Our results suggest that more sensitive tools for detecting low-level mosaicism in parents of families with seemingly "de novo" mutations will allow for better informed genetic counseling.

Cep57 protein is required for cytokinesis by facilitating central spindle microtubule organization.

Cytokinesis is the final stage of cell division in which the cytoplasm of a cell is divided into two daughter cells after the segregation of genetic material, and the central spindle and midbody are considered to be the essential structures required for the initiation and completion of cytokinesis. Here, we determined that the centrosome protein Cep57, which is localized to the central spindle and midbody, acts as a spindle organizer and is required for cytokinesis. Depletion of Cep57 disrupted microtubule assembly of the central spindle and further led to abnormal midbody localization of MKLP1, Plk1, and Aurora B, which resulted in cytokinesis failure and the formation of binuclear cells. Furthermore, we found that Cep57 directly recruited Tektin 1 to the midbody matrix to regulate microtubule organization. Thus, our data reveal that Cep57 is essential for cytokinesis via regulation of central spindle assembly and formation of the midbody.

A robust microsatellite instability detection model for unpaired colorectal cancer tissue samples.

BACKGROUND: Microsatellite instability (MSI) is a key biomarker for cancer immunotherapy and prognosis. Integration of MSI testing into a next-generation-sequencing (NGS) panel could save tissue sample, reduce turn-around time and cost, and provide MSI status and comprehensive genomic profiling in single test. We aimed to develop an MSI calling model to detect MSI status along with the NGS panel-based profiling test using tumor-only samples. METHODS: From January 2019 to December 2020, a total of 174 colorectal cancer (CRC) patients were enrolled, including 31 MSI-high (MSI-H) and 143 microsatellite stability (MSS) cases. Among them, 56 paired tumor and normal samples (10 MSI-H and 46 MSS) were used for modeling, and another 118 tumor-only samples were used for validation. MSI polymerase chain reaction (MSI-PCR) was performed as the gold standard. A baseline was built for the selected microsatellite loci using the NGS data of 56 normal blood samples. An MSI detection model was constructed by analyzing the NGS data of tissue samples. The performance of the model was compared with the results of MSI-PCR. RESULTS: We first intersected the target genomic regions of the NGS panels used in this study to select common microsatellite loci. A total of 42 loci including 23 mononucleotide repeat sites and 19 longer repeat sites were candidates for modeling. As mononucleotide repeat sites are more sensitive and specific for detecting MSI status than sites with longer length motif and the mononucleotide repeat sites performed even better than the total sites, a model containing 23 mononucleotide repeat sites was constructed and named Colorectal Cancer Microsatellite Instability test (CRC-MSI). The model achieved 100% sensitivity and 100% specificity when compared with MSI-PCR in both training and validation sets. Furthermore, the CRC-MSI model was robust with the tumor content as low as 6%. In addition, 8 out of 10 MSI-H samples showed alternations in the four mismatch repair genes ( MLH1 , MSH2 , MSH6 , and PMS2 ). CONCLUSION: MSI status can be accurately determined along the targeted NGS panels using only tumor samples. The performance of mononucleotide repeat sites surpasses loci with longer repeat motif in MSI calling.

Establishment and validation of a multigene model to predict the risk of relapse in hormone receptor-positive early-stage Chinese breast cancer patients.

BACKGROUND: Breast cancer patients who are positive for hormone receptor typically exhibit a favorable prognosis. It is controversial whether chemotherapy is necessary for them after surgery. Our study aimed to establish a multigene model to predict the relapse of hormone receptor-positive early-stage Chinese breast cancer after surgery and direct individualized application of chemotherapy in breast cancer patients after surgery. METHODS: In this study, differentially expressed genes (DEGs) were identified between relapse and nonrelapse breast cancer groups based on RNA sequencing. Gene set enrichment analysis (GSEA) was performed to identify potential relapse-relevant pathways. CIBERSORT and Microenvironment Cell Populations-counter algorithms were used to analyze immune infiltration. The least absolute shrinkage and selection operator (LASSO) regression, log-rank tests, and multiple Cox regression were performed to identify prognostic signatures. A predictive model was developed and validated based on Kaplan-Meier analysis, receiver operating characteristic curve (ROC). RESULTS: A total of 234 out of 487 patients were enrolled in this study, and 1588 DEGs were identified between the relapse and nonrelapse groups. GSEA results showed that immune-related pathways were enriched in the nonrelapse group, whereas cell cycle- and metabolism-relevant pathways were enriched in the relapse group. A predictive model was developed using three genes ( CKMT1B , SMR3B , and OR11M1P ) generated from the LASSO regression. The model stratified breast cancer patients into high- and low-risk subgroups with significantly different prognostic statuses, and our model was independent of other clinical factors. Time-dependent ROC showed high predictive performance of the model. CONCLUSIONS: A multigene model was established from RNA-sequencing data to direct risk classification and predict relapse of hormone receptor-positive breast cancer in Chinese patients. Utilization of the model could provide individualized evaluation of chemotherapy after surgery for breast cancer patients.