Bacterial droplet-based single-cell RNA-seq reveals antibiotic-associated heterogeneous cellular states.

We introduce BacDrop, a highly scalable technology for bacterial single-cell RNA sequencing that has overcome many challenges hindering the development of scRNA-seq in bacteria. BacDrop can be applied to thousands to millions of cells from both gram-negative and gram-positive species. It features universal ribosomal RNA depletion and combinatorial barcodes that enable multiplexing and massively parallel sequencing. We applied BacDrop to study Klebsiella pneumoniae clinical isolates and to elucidate their heterogeneous responses to antibiotic stress. In an unperturbed population presumed to be homogeneous, we found within-population heterogeneity largely driven by the expression of mobile genetic elements that promote the evolution of antibiotic resistance. Under antibiotic perturbation, BacDrop revealed transcriptionally distinct subpopulations associated with different phenotypic outcomes including antibiotic persistence. BacDrop thus can capture cellular states that cannot be detected by bulk RNA-seq, which will unlock new microbiological insights into bacterial responses to perturbations and larger bacterial communities such as the microbiome.

Comprehensive genetic and phenotype analysis of 95 individuals with mosaic tuberous sclerosis complex.

Tuberous sclerosis complex (TSC) is a neurogenetic disorder due to loss-of-function TSC1 or TSC2 variants, characterized by tumors affecting multiple organs, including skin, brain, heart, lung, and kidney. Mosaicism for TSC1 or TSC2 variants occurs in 10%-15% of individuals diagnosed with TSC. Here, we report comprehensive characterization of TSC mosaicism by using massively parallel sequencing (MPS) of 330 TSC samples from a variety of tissues and fluids from a cohort of 95 individuals with mosaic TSC. TSC1 variants in individuals with mosaic TSC are much less common (9%) than in germline TSC overall (26%) (p < 0.0001). The mosaic variant allele frequency (VAF) is significantly higher in TSC1 than in TSC2, in both blood and saliva (median VAF: TSC1, 4.91%; TSC2, 1.93%; p = 0.036) and facial angiofibromas (median VAF: TSC1, 7.7%; TSC2 3.7%; p = 0.004), while the number of TSC clinical features in individuals with TSC1 and TSC2 mosaicism was similar. The distribution of mosaic variants across TSC1 and TSC2 is similar to that for pathogenic germline variants in general TSC. The systemic mosaic variant was not present in blood in 14 of 76 (18%) individuals with TSC, highlighting the value of analysis of multiple samples from each individual. A detailed comparison revealed that nearly all TSC clinical features are less common in individuals with mosaic versus germline TSC. A large number of previously unreported TSC1 and TSC2 variants, including intronic and large rearrangements (n = 11), were also identified.

Genetics of Kidney Disease: The Unexpected Role of Rare Disorders.

Hundreds of different genetic causes of chronic kidney disease are now recognized, and while individually rare, taken together they are significant contributors to both adult and pediatric diseases. Traditional genetics approaches relied heavily on the identification of large families with multiple affected members and have been fundamental to the identification of genetic kidney diseases. With the increased utilization of massively parallel sequencing and improvements to genotype imputation, we can analyze rare variants in large cohorts of unrelated individuals, leading to personalized care for patients and significant research advancements. This review evaluates the contribution of rare disorders to patient care and the study of genetic kidney diseases and highlights key advancements that utilize new techniques to improve our ability to identify new gene-disease associations.

Integrating massive RNA-seq data to elucidate transcriptome dynamics in Drosophila melanogaster.

The volume of ribonucleic acid (RNA)-seq data has increased exponentially, providing numerous new insights into various biological processes. However, due to significant practical challenges, such as data heterogeneity, it is still difficult to ensure the quality of these data when integrated. Although some quality control methods have been developed, sample consistency is rarely considered and these methods are susceptible to artificial factors. Here, we developed MassiveQC, an unsupervised machine learning-based approach, to automatically download and filter large-scale high-throughput data. In addition to the read quality used in other tools, MassiveQC also uses the alignment and expression quality as model features. Meanwhile, it is user-friendly since the cutoff is generated from self-reporting and is applicable to multimodal data. To explore its value, we applied MassiveQC to Drosophila RNA-seq data and generated a comprehensive transcriptome atlas across 28 tissues from embryogenesis to adulthood. We systematically characterized fly gene expression dynamics and found that genes with high expression dynamics were likely to be evolutionarily young and expressed at late developmental stages, exhibiting high nonsynonymous substitution rates and low phenotypic severity, and they were involved in simple regulatory programs. We also discovered that human and Drosophila had strong positive correlations in gene expression in orthologous organs, revealing the great potential of the Drosophila system for studying human development and disease.

Improving the Yield of Genetic Diagnosis through Additional Genetic Panel Testing in Hereditary Ophthalmic Diseases.

Numerous hereditary ophthalmic diseases display significant genetic diversity. Consequently, the utilization of gene panel sequencing allows a greater number of patients to receive a genetic diagnosis for their clinical manifestations. We investigated how to improve the yield of genetic diagnosis through additional gene panel sequencing in hereditary ophthalmic diseases. A gene panel sequencing consisting of a customized hereditary retinopathy panel or hereditary retinitis pigmentosa (RP) panel was prescribed and referred to a CAP-accredited clinical laboratory. If no significant mutations associated with hereditary retinopathy and RP were detected in either panel, additional gene panel sequencing was requested for research use, utilizing the remaining panel. After additional gene panel sequencing, a total of 16 heterozygous or homozygous variants were identified in 15 different genes associated with hereditary ophthalmic diseases. Of 15 patients carrying any candidate variants, the clinical symptoms could be tentatively accounted for by genetic mutations in seven patients. However, in the remaining eight patients, given the in silico mutation predictive analysis, variant allele frequency in gnomAD, inheritance pattern, and genotype-phenotype correlation, fully elucidating the clinical manifestations with the identified rare variant was challenging. Our study highlights the utility of gene panel sequencing in achieving accurate diagnoses for hereditary ophthalmic diseases and enhancing the diagnostic yield through additional gene panel sequencing. Thus, gene panel sequencing can serve as a primary tool for the genetic diagnosis of hereditary ophthalmic diseases, even in cases where a single genetic cause is suspected. With a deeper comprehension of the genetic mechanisms underlying these diseases, it becomes feasible.

Population genetic analyses of Eastern Chinese Han nationality using ForenSeq™ DNA Signature Prep Kit.

Currently, the most commonly used method for human identification and kinship analysis in forensic genetics is the detection of length polymorphism in short tandem repeats (STRs) using polymerase chain reaction (PCR) and capillary electrophoresis (CE). However, numerous studies have shown that considerable sequence variations exist in the repeat and flanking regions of the STR loci, which cannot be identified by CE detection. Comparatively, massively parallel sequencing (MPS) technology can capture these sequence differences, thereby enhancing the identification capability of certain STRs. In this study, we used the ForenSeq™ DNA Signature Prep Kit to sequence 58 STRs and 94 individual identification SNPs (iiSNPs) in a sample of 220 unrelated individuals from the Eastern Chinese Han population. Our aim is to obtain MPS-based STR and SNP data, providing further evidence for the study of population genetics and forensic applications. The results showed that the MPS method, utilizing sequence information, identified a total of 486 alleles on autosomal STRs (A-STRs), 97 alleles on X-chromosome STRs (X-STRs), and 218 alleles on Y-chromosome STRs (Y-STRs). Compared with length polymorphism, we observed an increase of 260 alleles (157, 31, and 72 alleles on A-STRs, X-STRs, and Y-STRs, respectively) across 36 STRs. The most substantial increments were observed in DYF387S1 and DYS389II, with increases of 287.5% and 250%, respectively. The most increment in the number of alleles was found at DYF387S1 and DYS389II (287.5% and 250%, respectively). The length-based (LB) and sequence-based (SB) combined random match probability (RMP) of 27 A-STRs were 6.05E-31 and 1.53E-34, respectively. Furthermore, other forensic parameters such as total discrimination power (TDP), cumulative probability of exclusion of trios (CPEtrio), and duos (CPEduo) were significantly improved when using the SB data, and informative data were obtained for the 94 iiSNPs. Collectively, these findings highlight the advantages of MPS technology in forensic genetics, and the Eastern Chinese Han genetic data generated in this study could be used as a valuable reference for future research in this field.

Comparison of massively parallel sequencing to capillary electrophoresis for short tandem repeat genotyping of trace DNA.

In forensic genetics, massively parallel sequencing (MPS) offers several advantages over the current golden standard, capillary electrophoresis (CE): additional sequence information, shorter amplicon lengths, and the simultaneous analysis of many markers. These benefits result in a reduced number of reactions necessary while improving the amount of data obtained, thereby conserving valuable sample extracts. This proves particularly advantageous for the analysis of trace DNA. This study assessed the suitability of MPS for short tandem repeat (STR) typing of low template samples compared with results obtained through CE. The MPS genotypes showed higher concordance to reference genotypes, with donor alleles being more frequently assigned to be the major contributor, meeting the requirements for database entry. However, the MPS workflow is more time-consuming and associated with higher costs.

Analysis of homozygous allele mismatches in paternity tests with massively parallel sequencing.

In paternity testing, short tandem repeats (STRs) allele mismatches are often detected. Nowadays, polymerase chain reaction- and capillary electrophoresis (CE)-based STR genotyping is the most commonly used method to distinguish alleles based on their length. However, it could not detect alleles of the same size with sequence differences. Massively parallel sequencing (MPS) can determine not only allele sizes but also sequences, which could explain the causes of allele mismatches. Additionally, more types of genetic markers can be detected in a single assay, which increases the discriminatory power and facilitates the analysis of paternity tests. In this study, we analyzed 11 cases with homozygous allele mismatches from routine DNA trio paternity tests using the CE platform. Samples were sequenced using the ForenSeq DNA Signature Prep Kit and the MiSeq FGx Sequencing System. The results show that of the eight father-child mismatch cases and three mother-child mismatch cases, five cases with D5S818 and D8S1179 and one case at D13S317 were classified as non-amplification. The other three cases and two cases could be defined as mutations. This study suggests that MPS-based STR genotyping can provide additional information that allows more accurate interpretation of allelic mismatches in paternity testing.

Sequenced-based Rwanda population provides insights into demographic history.

Rwanda is known as the heart of Africa, reflecting the history of the world. Colonization and genocide have led to Rwanda's existing genetic structure. Herein, we used massively parallel sequencing to analyze 296 loci in 185 Rwandans and constructed a database for Rwandan forensic data for the first time. We found the following results: First, forensic parameters demonstrated that all loci were highly informative and could be used for forensic identification and paternity tests in Rwandans. Second, we found that the differences in genetic background between Rwandans and other African populations were similar but slight, as indicated by the massively parallel sequencing panel. Rwandans belonged to the African population and were inseparable from populations from neighboring countries. Also, Rwandans were closer to the European and American populations because of colonization, war, and other reasons. There was no scientific basis for racial classification established by colonization. Further research still needs to be carried out on more loci and larger Rwandan samples.

Integrative lncRNA, circRNA, and mRNA analysis reveals expression profiles of six forensic body fluids/tissue.

RNAs have attracted much attention in forensic body fluid/tissue identification (BFID) due to their tissue-specific expression characteristics. Among RNAs, long RNAs (e.g., mRNA) have a higher probability of containing more polymorphic sites that can be used to assign the specific donor of the body fluid/tissue. However, few studies have characterized their overall profiles in forensic science. In this study, we sequenced the transcriptomes of 30 samples from venous blood, menstrual blood, semen, saliva, vaginal secretion, and skin tissue, obtaining a comprehensive picture of mRNA, lncRNA, and circRNA profiles. A total of 90,305 mRNAs, 102,906 lncRNAs (including 19,549 novel lncRNAs), and 40,204 circRNAs were detected. RNA type distribution, length distribution, and expression distribution were presented according to their annotation and expression level, and many novel body fluid/tissue-specific RNA markers were identified. Furthermore, the cognate relations among the three RNAs were analyzed according to gene annotations. Finally, SNPs and InDels from RNA transcripts were genotyped, and 21,611 multi-SNP and 4,471 multi-InDel transcriptomic microhaplotypes (tMHs) were identified. These results provide a comprehensive understanding of transcriptome profiles, which could provide new avenues for tracing the origin of the body fluid/tissue and identifying an individual.

Developmental validation of the STRSeqTyper122 kit for massively parallel sequencing of forensic STRs.

Massively parallel sequencing allows for integrated genotyping of different types of forensic markers, which reduces DNA consumption, simplifies experimental processes, and provides additional sequence-based genetic information. The STRseqTyper122 kit genotypes 63 autosomal STRs, 16 X-STRs, 42 Y-STRs, and the Amelogenin locus. Amplicon sizes of 117 loci were below 300 bp. In this study, MiSeq FGx sequencing metrics for STRseqTyper122 were presented. The genotyping accuracy of this kit was examined by comparing to certified genotypes of NIST standard reference materials and results from five capillary electrophoresis-based kits. The sensitivity of STRseqTyper122 reached 125 pg, and > 80% of the loci were correctly called with 62.5 pg and 31.25 pg input genomic DNA. Repeatability, species specificity, and tolerance for DNA degradation and PCR inhibitors of this kit were also evaluated. STRseqTyper122 demonstrated reliable performance with routine case-work samples and provided a powerful tool for forensic applications.

Genetic Evaluation for Monogenic Disorders of Low Bone Mass and Increased Bone Fragility: What Clinicians Need to Know.

PURPOSE OF REVIEW: The purpose of this review is to outline the principles of clinical genetic testing and to provide practical guidance to clinicians in navigating genetic testing for patients with suspected monogenic forms of osteoporosis. RECENT FINDINGS: Heritability assessments and genome-wide association studies have clearly shown the significant contributions of genetic variations to the pathogenesis of osteoporosis. Currently, over 50 monogenic disorders that present primarily with low bone mass and increased risk of fractures have been described. The widespread availability of clinical genetic testing offers a valuable opportunity to correctly diagnose individuals with monogenic forms of osteoporosis, thus instituting appropriate surveillance and treatment. Clinical genetic testing may identify the appropriate diagnosis in a subset of patients with low bone mass, multiple or unusual fractures, and severe or early-onset osteoporosis, and thus clinicians should be aware of how to incorporate such testing into their clinical practices.

The Development and Characterisation of ssDNA Aptamers via a Modified Cell-SELEX Methodology for the Detection of Human Red Blood Cells.

Blood is one of the most commonly found biological fluids at crime scenes, with the detection and identification of blood holding a high degree of evidential value. It can provide not only information about the nature of the crime but can also lead to identification via DNA profiling. Presumptive tests for blood are usually sensitive but not specific, so small amounts of the substrate can be detected, but false-positive results are often encountered, which can be misleading. Novel methods for the detection of red blood cells based on aptamer-target interactions may be able to overcome these issues. Aptamers are single-stranded DNA or RNA sequences capable of undergoing selective antigen association due to three-dimensional structure formation. The use of aptamers as a target-specific moiety poses several advantages and has the potential to replace antibodies within immunoassays. Aptamers are cheaper to produce, display no batch-to-batch variation and can allow for a wide range of chemical modifications. They can help limit cross-reactivity, which is a hindrance to current forensic testing methods. Within this study, a modified Systematic Evolution of Ligands by Exponential Enrichment (SELEX) process was used to generate aptamers against whole red blood cells. Obtained aptamer pools were analysed via massively parallel sequencing to identify viable sequences that demonstrate a high affinity for the target. Using bioinformatics platforms, aptamer candidates were identified via their enrichment profiles. Binding characterisation was also conducted on two selected aptamer candidates via fluorescent microscopy and qPCR to visualise and quantify aptamer binding. The potential for these aptamers is broad as they can be utilised within a range of bioassays for not only forensic applications but also other analytical science and medical applications. Potential future work includes the incorporation of developed aptamers into a biosensing platform that can be used at crime scenes for the real-time detection of human blood.

Peak Scores Significantly Depend on the Relationships between Contextual Signals in ChIP-Seq Peaks.

Chromatin immunoprecipitation followed by massively parallel DNA sequencing (ChIP-seq) is a central genome-wide method for in vivo analyses of DNA-protein interactions in various cellular conditions. Numerous studies have demonstrated the complex contextual organization of ChIP-seq peak sequences and the presence of binding sites for transcription factors in them. We assessed the dependence of the ChIP-seq peak score on the presence of different contextual signals in the peak sequences by analyzing these sequences from several ChIP-seq experiments using our fully enumerative GPU-based de novo motif discovery method, Argo_CUDA. Analysis revealed sets of significant IUPAC motifs corresponding to the binding sites of the target and partner transcription factors. For these ChIP-seq experiments, multiple regression models were constructed, demonstrating a significant dependence of the peak scores on the presence in the peak sequences of not only highly significant target motifs but also less significant motifs corresponding to the binding sites of the partner transcription factors. A significant correlation was shown between the presence of the target motifs FOXA2 and the partner motifs HNF4G, which found experimental confirmation in the scientific literature, demonstrating the important contribution of the partner transcription factors to the binding of the target transcription factor to DNA and, consequently, their important contribution to the peak score.

Fast all versus all genotype comparison using DNA/RNA sequencing data: method and workflow.

BACKGROUND: Massively parallel sequencing includes many liquid handling steps which introduce the possibility of sample swaps, mixing, and duplication. The unique profile of inherited variants in human genomes allows for comparison of sample identity using sequence data. A comparison of all samples vs. each other (all vs. all) provides both identification of mismatched samples and the possibility of resolving swapped samples. However, all vs. all comparison complexity grows as the square of the number of samples, so efficiency becomes essential. RESULTS: We have developed a tool for fast all vs. all genotype comparison using low level bitwise operations built into the Perl programming language. Importantly, we have also developed a complete workflow allowing users to start with either raw FASTQ sequence files, aligned BAM files, or genotype VCF files and automatically generate comparison metrics and summary plots. The tool is freely available at https://github.com/teerjk/TimeAttackGenComp/ . CONCLUSIONS: A fast and easy to use method for genotype comparison as described here is an important tool to ensure high quality and robust results in sequencing studies.

Spectrum of variants in a large Chinese Gitelman syndrome cohort.

Gitelman syndrome (GS) is caused by SLC12A3 biallelic variants. A previous study showed that large rearrangements (LRGs) of SLC12A3 accounted for the low sensitivity of genetic testing. However, a systematic screening for LRGs in Chinese GS patients is lacking. Massively parallel sequencing (MPS) and multiplex ligation-dependent probe amplification (MLPA) were performed to sequence the genomic DNA of patients with clinically diagnosed GS. Of 165 index cases, MPS identified 151 cases with two or more affected alleles and 14 cases with one variant allele. LRGs were detected by MLPA in 20 out of 27 cases, including 15 cases with suspected LRGs by MPS. Among these 20 cases with LRGs, the results of MPS and MLPA were identical in only 8 cases. Additional LRGs in 6 cases were detected by MLPA alone. In 6 cases, E4_E6del was identified by MPS, while E4_E5del and Intron6del were identified by MLPA. Among the 102 distinct variants, 30 are novel. LRGs were found in 20 cases (12.1%). LRGs were found in 12.1% of our Chinese GS patients cohort. We show that MPS and MLPA are two complementary techniques with the ability to improve the diagnostic yield of GS.

Development and validation of a custom panel including 114 InDels using massively parallel sequencing for forensic application.

Insertion/deletion polymorphisms (InDels) have particular characteristics, such as a relatively low mutation rate, small amplicon size, and no stutter artifacts when genotyped via the capillary electrophoresis platform. It would be an important complementary tool for individual identification and certain kinship analyses. At present, massively parallel sequencing (MPS) has shown excellent application value in forensic studies. Therefore, in this study, we developed a custom MPS InDel panel that contains 114 InDels [77 autosomal InDels (A-InDels), 32 X-chromosomal InDels (X-InDels), and 5 Y-chromosomal InDels) based on previous studies. To assess this panel's performance, several validation experiments were performed, including sensitivity, inhibitor, degraded DNA testing, species specificity, concordance, repeatability, case-type samples, and population studies. The results showed that the lowest DNA input was 0.25 ng. All genotypes were obtained in the presence of 80 ng/µL humic acid, 2000 µmol/L calcium, 3000 µmol/L EDTA and indigo. In degraded DNA testing, 90% of loci could be detected for 16-day-old formalin-fixed hearts. In addition, this panel has good species specificity. The values of combined power of discrimination and the combined power of exclusion for 77 A-InDels were 1-3.9951 × 10-32 and 1-4.2956 × 10-7 , respectively. The combined mean exclusion chance for 32 X-InDels was 0.99999 in trios and 0.99904 in duos. The validation results indicate that this newly developed MPS multiplex system is a robust tool for forensic applications.

Application of SNPs with low minor allele frequencies in missing person identification (MPI) through kinship analysis of DNA mixtures.

The need to identify a missing person (MP) through kinship analysis of DNA samples found at a crime scene has become increasingly prevalent. DNA samples from MPs can be severely degraded, contain little DNA and mixed with other contributors, which often makes it difficult to apply conventional methods in practice. This study developed a massively parallel sequencing-based panel that contains 1661 single-nucleotide polymorphisms (SNPs) with low minor allele frequencies (MAFs) (averaged at 0.0613) in the Chinese Han population, and the strategy for relationship inference from DNA mixtures comprising different numbers of contributors (NOCs) and of varying allele dropout probabilities. Based on the simulated dataset and genotyping results of 42 artificial DNA mixtures (NOC = 2-4), it was observed that the present SNP panel was sufficient for balanced mixtures when referenced to the closest relatives (parents/offspring and full siblings). When the mixture profiles suffered from dropout, incorrect assignments were markedly associated with relatedness, NOC and the dropout level. We, therefore, indicate that SNPs with low MAFs could be reliably interpreted for MP identification through the kinship analysis of complex DNA mixtures. Further studies should be extended to more possible scenarios to test the feasibility of this present approach.

Noninvasive fetal genotyping of single nucleotide variants and linkage analysis for prenatal diagnosis of monogenic disorders.

BACKGROUND: High-cost, time-consuming and complex processes of several current approaches limit the use of noninvasive prenatal diagnosis (NIPD) for monogenic disorders in clinical application. Thus, a more cost-effective and easily implementable approach is required. METHODS: We established a low-cost and convenient test to noninvasively deduce fetal genotypes of the mutation and single nucleotide polymorphisms (SNPs) loci by means of targeted amplification combined with deep sequencing of maternal genomic and plasma DNA. The sequential probability ratio test was performed to detect the allelic imbalance in maternal plasma. This method can be employed to directly examine familial pathogenic mutations in the fetal genome, as well as infer the inheritance of parental haplotypes through a group of selected SNPs linked to the pathogenic mutation. RESULTS: The fetal mutations in 17 families with different types of monogenic disorders including hemophilia A, von Willebrand disease type 3, Duchenne muscular dystrophy, hyper-IgM type 1, glutaric acidemia type I, Nagashima-type palmoplantar keratosis, and familial exudative vitreoretinopathy were identified in the study. The mutations included various forms: point mutations, gene inversion, deletions/insertions and duplication. The results of 12 families were verified by sequencing of amniotic fluid samples, the accuracy of the approach in fetal genotyping at the mutation and SNPs loci was 98.85% (172/174 loci), and the no-call rate was 28.98% (71/245 loci). The overall accuracy was 12/12 (100%). Moreover, the approach was successfully applied in plasma samples with a fetal fraction as low as 2.3%. CONCLUSIONS: We have shown in this study that the approach is a cost-effective, less time consuming and accurate method for NIPD of monogenic disorders.

Landscape of chromatin remodeling gene alterations in endometrial carcinoma.

OBJECTIVE: Chromatin remodeling genes (CRGs) encode components of epigenetic regulatory mechanisms and alterations in these genes have been identified in several tumor types, including gynecologic cancers. In this study, we sought to investigate the prevalence and clinicopathological associations of CRG alterations in endometrial carcinoma (EC). METHODS: We performed a retrospective analysis of 660 ECs sequenced using a clinical massively parallel sequencing assay targeting up to 468 genes, including 25 CRGs, and defined the presence of somatic CRG alterations. Clinicopathologic features were obtained for all cases. Immunohistochemical interrogation of ARID1A and PTEN proteins was performed in a subset of samples. RESULTS: Of the 660 ECs sequenced, 438 (66.4%) harbored CRG alterations covered by our panel. The most commonly altered CRG was ARID1A (46%), followed by CTCF (21%), KMT2D (18%), KMT2B (17%), BCOR (16%), ARID1B (12%) and SMARCA4 (11%). We found that ARID1A genetic alterations were preferentially bi-allelic and often corresponded to altered ARID1A protein expression in ECs. We further observed that ARID1A alterations were often subclonal when compared to PTEN alterations, which were primarily clonal in ECs harboring both mutations. Finally, CRG alterations were associated with an increased likelihood of myometrial and lymphovascular invasion in endometrioid ECs. CONCLUSION: CRG alterations are common in EC and are associated with clinicopathologic features and likely play a crucial role in EC.