Autosomal recessive non-syndromic hearing loss (ARNSHL) is a highly heterogeneous disease involving more than 70 pathogenic genes. However, most ARNSHL families have small-sized pedigrees with limited genetic information, rendering challenges for the molecular diagnosis of these patients. Therefore, we attempted to establish a strategy for identifying deleterious variants associated with ARNSHL by applying proband whole-exome sequencing (proband-WES). Aside from desiring to improve molecular diagnostic rates, we also aimed to search for novel deafness genes shared by patients with similar phenotype, making up for the deficiency of small ARNSHL families. In this study, 48.5% (16/33) families were detected the pathogenic variants in eight known deafness genes, including 10 novel variants identified in TMPRSS3 (MIM 605551), MYO15A (MIM 602666), TMC1 (MIM 606706), ADGRV1 (MIM 602851), and PTPRQ (MIM 603317). Apart from six novel variants with a truncating effect (nonsense, deletion, insertion, and splice-site), four novel missense variants were not found in 200 unrelated control population by using Sanger sequencing. It is important to note that none of novel genes were shared across different pedigrees, indicating that a larger sample size might be needed. Proband-WES is a cost-effective and precise way of identifying causative variants in nuclear families with ARNSHL. This economical strategy may be appropriated as a clinical application to provide molecular diagnostics, genetic counseling, and individualized health maintenance measures for patients with ARNSHL at hearing clinics.

Background: Nasal polyp (NP) is a chronic upper airway inflammatory disease that is frequently triggered by defective host-defense. However, the mechanisms underlying the impaired barrier function such as cilia-mediated mucociliary clearance remain poorly understood. Objective: To assess ciliary ultrastructural and ciliogenesis marker expression and the phenotypes of ciliated cells in NP. Methods: NP biopsy samples were obtained from 97 NP patients and inferior turbinate from 32 healthy controls. Immunofluorescence staining, quantitative polymerase chain reaction, and single-cell cytospin staining were performed. We classified the patterns of radial spoke head protein (RSPH) 1, 4A (RSPH4A), 9 (RSPH9), and dynein axonemal heavy chain 5 (DNAH5) localization. A semi-quantitative scoring system was developed to assess their expression patterns and associations with ciliogenesis markers [centrosomal protein 110 (CP110) and forkhead box j1 (FOXJ1)]. Results: Median scores of RSPH1, RSPH4A, RSPH9, and DNAH5 were significantly higher in NP than in healthy controls, particularly in eosinophilic NPs. Expression pattern scores of RSPH1, RSPH4A, RSPH9, and DNAH5 correlated positively with each other in both groups. In primary-cell specimens, abnormal expression patterns were significantly more common in NP. The total fluorescence intensity of CP110 and FOXJ1 was significantly higher in NPs and correlated positively with expression pattern scores of RSPH1, RSPH4A, RSPH9, and DNAH5. A trend towards lengthened cilia was observed in NP. Conclusion: In the chronic airway inflammatory milieu, the up-regulated ciliogenesis correlates with the abnormal expression of ciliary ultrastructural markers (i.e., DNAH5) in NP (particularly eosinophilic NP).

Tibetans residing in the high-altitude inhospitable environment have undergone significant natural selection of their genetic architecture. Recently, highly mutational autosomal short tandem repeats were widely used not only in the anthropology and population genetics to investigate the genetic structure and relationships, but also in the medical genetics to explore the pathogenesis of multiple genetic diseases and in the forensic science to identify individual and parentage relatedness. However, genetic variants and forensic efficiency of DNATyperTM 19 amplification system and genetic background of Kham Tibetan remain uncharacterized. Thus, we genotyped 19 forensic genetic markers in 11,402 Kham Tibetans to gain insight into the genetic diversity of Chinese high-altitude adaptive population. Highly discriminating and polymorphic forensic measures were observed, which indicated that this new-developed DNATyper 19 PCR amplification is suitable for routine forensic identification purposes and Chinese national DNA database establishment. Pairwise genetic distances among the comprehensive population comparisons suggested that this high-altitude adaptive Kham Tibetan has genetically closer relationships with lowlanders of Tibeto-Burman-speaking populations (Chengdu Tibetan, Liangshan Tibetan, and Liangshan Yi). Genetic substructure analyses via phylogenetic reconstruction, principal component analysis, and multidimensional scaling analysis in both nationwide and worldwide contexts suggested that the genetic proximity exists along the linguistic, ethnic, and continental geographical boundary. Further studies with whole-genome sequencing of modern or archaic Kham Tibetans would be useful in reconstructing the Tibetan population history.

Unintended off-target mutations induced by CRISPR/Cas9 nucleases may result in unwanted consequences, which will impede the efficient applicability of this technology for genetic improvement. We have recently edited the goat genome through CRISPR/Cas9 by targeting MSTN and FGF5, which increased muscle fiber diameter and hair fiber length, respectively. Using family trio-based sequencing that allow better discrimination of variant origins, we herein generated offspring from edited goats, and sequenced the members of four family trios (gene-edited goats and their offspring) to an average of ∼36.8× coverage. This data was to systematically examined for mutation profiles using a stringent pipeline that comprehensively analyzed the sequence data for de novo single nucleotide variants, indels, and structural variants from the genome. Our results revealed that the incidence of de novo mutations in the offspring was equivalent to normal populations. We further conducted RNA sequencing using muscle and skin tissues from the offspring and control animals, the differentially expressed genes (DEGs) were related to muscle fiber development in muscles, skin development, and immune responses in skin tissues. Furthermore, in contrast to recently reports of Cas9 triggered p53 expression alterations in cultured cells, we provide primary evidence to show that Cas9-mediated genetic modification does not induce apparent p53 expression changes in animal tissues. This work provides adequate molecular evidence to support the reliability of conducting Cas9-mediated genome editing in large animal models for biomedicine and agriculture.

Background: Glioblastoma is a common malignant neuroepithelial neoplasm with poor clinical outcomes and limited treatment options. It is extremely important to search and confirm diverse hub genes that are effective in the advance and prediction of glioblastoma. Methods: We analyzed GSE50161, GSE4290, and GSE68848, the three microarray datasets retrieved from the GEO database. GO function and KEGG pathway enrichment analyses for differentially expressed genes (DEGs) were performed using DAVID. The PPI network of the DEGs was analyzed using the Search Tool for the Retrieval of Interacting Genes database and visualized by Cytoscape software. Hub genes were identified through the PPI network and a robust rank aggregation method. The Cancer Genome Atlas (TCGA) and the Oncomine database were used to validate the hub genes. In addition, a survival curve analysis was conducted to verify the correlation between the expression of hub genes and patient prognosis. Human glioblastoma cells and normal cells were collected, and then RT-PCR, Western blot, and immunofluorescence were conducted to validate the expression of the NDC80 gene. A cell proliferation assay was used to detect the proliferation of glioma cells. The effects of NDC80 expression on migration and invasion of GBM cell lines were evaluated by conducting scratch and transwell assays. Results: A total of 716 DEGs were common to all three microarray datasets, which included 188 upregulated DEGs and 528 downregulated DEGs. Furthermore, we found that among the common DEGs, 10 hub genes showed a high degree of connectivity. The expression of the 10 hub genes in TCGA and the Oncomine database was significantly overexpressed in glioblastoma compared with normal genes. Additionally, the survival analysis showed that the patients with low expression of six genes (BIR5C, CDC20, NDC80, CDK1, TOP2A, and MELK) had a significantly favorable prognosis (p < 0.01). We discovered that NDC80, which has been shown to be important in other cancers, also has an important role in malignant gliomas. The RT-PCR, Western blot, and immunofluorescence results showed that the expression level of NDC80 was significantly higher in human glioblastoma cells than in normal cells. Moreover, we identified that NDC80 increased the proliferation and invasion abilities of human glioblastoma cells. Conclusion: The six genes identified here may be utilized to form a panel of disease progression and predictive biomarkers of glioblastoma for clinical purposes. NDC80, one of the six genes, was discovered to have a potentially important role in GBM, a finding that needs to be further studied.

Background: 5-Methylcytidine (m5C) is the most common RNA modification and plays an important role in multiple tumors including cervical cancer (CC). We aimed to develop a novel gene signature by identifying m5C modification subtypes of CC to better predict the prognosis of patients. Methods: We obtained the expression of 13 m5C regulatory factors from The Cancer Genome Atlas (TCGA all set, 257 patients) to determine m5C modification subtypes by the "nonnegative matrix factorization" (NMF). Then the "limma" package was used to identify differentially expressed genes (DEGs) between different subtypes. According to these DEGs, we performed Cox regression and Kaplan-Meier (KM) survival analysis to establish a novel gene signature in TCGA training set (128 patients). We also verified the risk prediction effect of gene signature in TCGA test set (129 patients), TCGA all set (257 patients) and GSE44001 (300 patients). Furthermore, a nomogram including this gene signature and clinicopathological parameters was established to predict the individual survival rate. Finally, the expression and function of these signature genes were explored by qRT-PCR, immunohistochemistry (IHC) and proliferation, colony formation, migration and invasion assays. Results: Based on consistent clustering of 13 m5C-modified genes, CC was divided into two subtypes (C1 and C2) and the C1 subtype had a worse prognosis. The 4-gene signature comprising FNDC3A, VEGFA, OPN3 and CPE was constructed. In TCGA training set and three validation sets, we found the prognosis of patients in the low-risk group was much better than that in the high-risk group. A nomogram incorporating the gene signature and T stage was constructed, and the calibration plot suggested that it could accurately predict the survival rate. The expression levels of FNDC3A, VEGFA, OPN3 and CPE were all high in cervical cancer tissues. Downregulation of FNDC3A, VEGFA or CPE expression suppressed the proliferation, migration and invasion of SiHa cells. Conclusions: Two m5C modification subtypes of CC were identified and then a 4-gene signature was established, which provide new feasible methods for clinical risk assessment and targeted therapies for CC.

Some studies have proven that autophagy and lncRNA play important roles in AML. Several autophagy related lncRNA signatures have been shown to affect the survival of patients in some other cancers. However, the role of autophagy related lncRNA in AML has not been explored yet. Hence, this study aims to find an autophagy related lncRNA signature that can affect survival for AML patients.

The emerging evidence has demonstrated the critical roles of long non-coding RNAs (lncRNAs) as regulators in the tumor immune microenvironment (TIME). However, the tumor immune infiltration-associated lncRNAs and their clinical significance in colon cancer have not yet been thoroughly investigated. This study performed an integrative analysis of lncRNA expression profiles and immune cell infiltration profiles and identified 258 immune infiltration-associated lncRNAs. Of them, four lncRNAs (AC008494.3, LINC00926, AC022034.1, and SNHG26) were significantly and independently associated with the patient's overall survival. Finally, we developed a tumor immune infiltration-associated lncRNA signature (TIILncSig) comprising of these four lncRNAs, which can divide colon cancer patients of The Cancer Genome Atlas (TCGA) into high-risk and low-risk groups with a significantly different outcome [Hazard ratio (HR) = 2.718, 95% CI = 1.955-3.779, p < 0.001]. Prognostic performance of the TIILncSig was further validated in another independent colon cancer cohort (HR = 1.832, 95% CI = 1.045-3.21, p = 0.034). Results of multivariate Cox regression and stratification analysis demonstrated that the TIILncSig is an independent predictive factor from other clinical features (HR = 2.687, 95% CI = 1.912-3.776, p < 0.001 for TCGA cohort and HR = 1.837, 95% CI = 1.047-3.223, p = 0.034 for GSE17538 cohort). Literature analysis provided experimental evidence supporting roles of the TIILncSig in cancer carcinogenesis and progression and immune regulation. Summary, our study will help to understand the mechanisms of lncRNAs in immune regulation in the tumor microenvironment and provide novel biomarkers or targets for prognosis prediction and therapy decision-making for patients with colon cancer.

Oral squamous cell carcinoma (OSCC) is a solid tumor, which originates from squamous epithelium, with about 400,000 new-cases/year worldwidely. Presently, chemoradiotherapy is the most important adjuvant treatment for OSCC, mostly in advanced tumors. However, clinical resistance to chemotherapy still leads to poor prognosis of OSCC patients. Via high-throughput analysis of gene expression database of OSCC, we investigated the molecular mechanisms underlying cisplatin resistance in OSCC, analyzing the differentially expressed genes (DEGs) and their regulatory relationship, to clarify the molecular basis of OSCC chemotherapy resistance and provide a theoretical foundation for the treatment of patients with OSCC and individualized therapeutic targets accurately.

The evolutionary and admixture history of Han Chinese have been widely discussed via traditional autosomal and uniparental genetic markers [e.g., short tandem repeats, low-density single nucleotide polymorphisms). However, their fine-scale genetic landscapes (admixture scenarios and natural selection signatures) based on the high-density allele/haplotype sharing patterns have not been deeply characterized. Here, we collected and generated genome-wide data of 50 Han Chinese individuals from four populations in Guizhou Province, one of the most ethnolinguistically diverse regions, and merged it with over 3,000 publicly available modern and ancient Eurasians to describe the genetic origin and population admixture history of Guizhou Hans and their neighbors. PCA and ADMIXTURE results showed that the studied four populations were homogeneous and grouped closely to central East Asians. Genetic homogeneity within Guizhou populations was further confirmed via the observed strong genetic affinity with inland Hmong-Mien people through the observed genetic clade in Fst and outgroup f 3 /f 4-statistics. qpGraph-based phylogenies and f 4-based demographic models illuminated that Guizhou Hans were well fitted via the admixture of ancient Yellow River Millet farmers related to Lajia people and southern Yangtze River farmers related to Hanben people. Further ChromoPainter-based chromosome painting profiles and GLOBETROTTER-based admixture signatures confirmed the two best source matches for southwestern Hans, respectively, from northern Shaanxi Hans and southern indigenes with variable mixture proportions in the historical period. Further three-way admixture models revealed larger genetic contributions from coastal southern East Asians into Guizhou Hans compared with the proposed inland ancient source from mainland Southeast Asia. We also identified candidate loci (e.g., MTUS2, NOTCH4, EDAR, ADH1B, and ABCG2) with strong natural selection signatures in Guizhou Hans via iHS, nSL, and ihh, which were associated with the susceptibility of the multiple complex diseases, morphology formation, alcohol and lipid metabolism. Generally, we provided a case and ideal strategy to reconstruct the detailed demographic evolutionary history of Guizhou Hans, which provided new insights into the fine-scale genomic formation of one ethnolinguistically specific targeted population from the comprehensive perspectives of the shared unlinked alleles, linked haplotypes, and paternal and maternal lineages.

Cornelia de Lange syndrome (CdLS) is a genetic disorder characterized by multisystemic malformations. Mutation in the NIPBL gene accounts for nearly 60% of the cases. This study reports the clinical and genetic findings of three cases of CdLS from unrelated Chinese families. Clinically, all the three cases were classified as classic CdLS based on the cardinal (distinctive facial features and limb malformations) and suggestive (developmental delay, growth retardation, microcephaly, hirsutism, etc.) manifestations. SNP array detected a novel de novo heterozygous microdeletion of 0.2 Mb [arr[GRCh37]5p13.2(36848530_37052821) × 1] that spans the first 43 exons of NIPBL in the fetus with nuchal translucency thickening in case 1. Whole-exome sequencing in family trios plus Sanger sequencing validation identified a de novo heterozygous NIPBL c.5566G>A (p.R1856G) mutation in the fetus with intrauterine growth retardation in case 2 and a novel de novo heterozygous NIPBL c.448dupA (p.S150Kfs*23) mutation in the proband (an 8-month-old girl) in case 3. The cases presented in this study may serve as references for increasing our understanding of the mutation spectrum of NIPBL in association with CdLS.

Background: Skeletal dysplasias (SDs) are a heterogeneous group of genetic disorders that primarily affect bone and cartilage. This study aims to identify the genetic causes for fetal SDs, and evaluates the diagnostic yield of prenatal whole-exome sequencing (WES) for this disorder. Methods: WES was performed on 38 fetuses with sonographically identified SDs and normal results of karyotype and single nucleotide polymorphism (SNP) analysis. Candidate variants were selected by bioinformatics analysis, and verified by Sanger sequencing. Results: WES revealed pathogenic or likely pathogenic variants associated with SDs in 65.79% (25/38) of fetuses, variants of uncertain significance (VUS) in SDs-related genes in 10.53% (4/38) cases, and incidental findings in 31.58% (12/38) fetuses. The SDs-associated variants identified in the present study affected 10 genes, and 35.71% (10/28) of the variants were novel. Conclusion: WES has a high diagnostic rate for prenatal SDs, which improves pregnancy management, prenatal counseling and recurrence risk assessment for future pregnancies. The newly identified variants expanded mutation spectrum of this disorder.

The Tibetan Plateau (TP) is considered to be one of the last terrestrial environments conquered by the anatomically modern human. Understanding of the genetic background of highland Tibetans plays a pivotal role in archeology, anthropology, genetics, and forensic investigations. Here, we genotyped 22 forensic genetic markers in 1,089 Tibetans residing in Nagqu Prefecture and collected 1,233,013 single nucleotide polymorphisms (SNPs) in the highland East Asians (Sherpa and Tibetan) from the Simons Genome Diversity Project and ancient Tibetans from Nepal and Neolithic farmers from northeastern Qinghai-Tibetan Plateau from public databases. We subsequently merged our two datasets with other worldwide reference populations or eastern ancient Eurasians to gain new insights into the genetic diversity, population movements, and admixtures of high-altitude East Asians via comprehensive population genetic statistical tools [principal component analysis (PCA), multidimensional scaling plot (MDS), STRUCTURE/ADMIXTURE, f3 , f4 , qpWave/qpAdm, and qpGraph]. Besides, we also explored their forensic characteristics and extended the Chinese National Database based on STR data. We identified 231 alleles with the corresponding allele frequencies spanning from 0.0005 to 0.5624 in the forensic low-density dataset, in which the combined powers of discrimination and the probability of exclusion were 1-1.22E-24 and 0.999999998, respectively. Additionally, comprehensive population comparisons in our low-density data among 57 worldwide populations via the Nei's genetic distance, PCA, MDS, NJ tree, and STRUCTURE analysis indicated that the highland Tibeto-Burman speakers kept the close genetic relationship with ethnically close populations. Findings from the 1240K high-density dataset not only confirmed the close genetic connection between modern Highlanders, Nepal ancients (Samdzong, Mebrak, and Chokhopani), and the upper Yellow River Qijia people, suggesting the northeastern edge of the TP served as a geographical corridor for ancient population migrations and interactions between highland and lowland regions, but also evidenced that late Neolithic farmers permanently colonized into the TP by adopting cold-tolerant barley agriculture that was mediated via the acculturation of idea via the millet farmer and not via the movement of barley agriculturalist as no obvious western Eurasian admixture signals were identified in our analyzed modern and ancient populations. Besides, results from the qpAdm-based admixture proportion estimation and qpGraph-based phylogenetic relationship reconstruction consistently demonstrated that all ancient and modern highland East Asians harbored and shared the deeply diverged Onge/Hoabinhian-related eastern Eurasian lineage, suggesting a common Paleolithic genetic legacy existed in high-altitude East Asians as the first layer of their gene pool.

There is still a considerable proportion of patients with inherited peripheral neuropathy (IPN) whose pathogenic genes are unknown. This study was intended to investigate whether the GGC repeat expansion in the NOTCH2NLC is presented in some patients with IPN. A total of 142 unrelated mainland Chinese patients with highly suspected diagnosis of IPN without any known causative gene were recruited. Repeat-primed polymerase chain reaction (RP-PCR) was performed to screen GGC repeat expansion in NOTCH2NLC, followed by fluorescence amplicon length analysis-PCR (AL-PCR) to determine the GGC repeat size. Detailed clinical data as well as nerve, muscle, and skin biopsy were reviewed and analyzed in the NOTCH2NLC-related IPN patients. In total, five of the 142 patients (3.52%) were found to have pathogenic GGC expansion in NOTCH2NLC, with repeat size ranging from 126 to 206 repeats. All the NOTCH2NLC-related IPN patients presented with adult-onset motor-sensory and autonomic neuropathy that predominantly affected the motor component of peripheral nerves. While tremor and irritating dry cough were noted in four-fifths of the patients, no other signs of the central nervous system were presented. Electrophysiological studies revealed both demyelinating and axonal changes of polyneuropathy that were more severe in lower limbs and asymmetrically in upper limbs. Sural nerve pathology was characterized by multiple fibers with thin myelination, indicating a predominant demyelinating process. Muscle pathology was consistent with neuropathic changes. P62-positive intranuclear inclusions were observed in nerve, skin, and muscle tissues. Our study has demonstrated that GGC expansion in NOTCH2NLC is associated with IPN presenting as predominant motor-sensory and autonomic neuropathy, which expands the phenotype of the NOTCH2NLC-related repeat expansion spectrum. Screening of GGC repeat expansions in the NOTCH2NLC should be considered in patients presenting with peripheral neuropathy with tremor and irritating dry cough.

Introduction: Systematic pan-cancer analysis of the roles and regulatory mechanisms for PRR7-AS1 is currently not available. Methods: In the present study, a comprehensive bioinformatic approach was used to mine the underlying oncogenic effects of PRR7-AS1, including expression status, prognostic value and immune characteristics. Results: We discovered that PRR7-AS1 expression was remarkably upregulated in most cancer types and exhibited a negative correlation with the prognosis. Furthermore, PRR7-AS1 expression was inversely connected with the majority of tumor-infiltrating immune cells, immune scores and immune checkpoint gene expression in pancancer. There was also a significant correlation between PRR7-AS1 expression status and tumor mutational burden, microsatellite instability, and neoantigens in certain tumors. PRR7-AS1 had the best predictive power for immune checkpoint blockade efficacy compared to other well-recognized biomarkers. PRR7-AS1 overexpression could affect cytotoxic T cells-mediated antitumor responses. Functional enrichment analysis revealed that PRR7-AS1 might be involved in the metabolic pathways. Super enhancer activity might have participated in the regulation of PRR7-AS1 expression. And we constructed the competitive endogenous RNA networks for PRR7-AS1. Discussion: In general, PRR7-AS1 had the potential to be a diagnostic, prognostic and immune biomarker for pan cancer. PRR7-AS1 was correlated with an immunosuppressive microenvironment and was a new potential target for immunotherapy. Epigenetic factors were the driving forces for PRR7-AS1 overexpression in tumors.

Non-small-cell lung cancer (NSCLC) is divided into three major histological types, namely, lung adenocarcinoma (LUAD), lung squamous cell carcinoma (LUSC), and large-cell lung carcinoma (LCLC). We previously identified that 4.1N/EPB41L1 acts as a tumor suppressor and is reduced in NSCLC patients. In the current study, we explored the underlying epigenetic mechanisms of 4.1N/EPB41L1 reduction in NSCLC. The 4.1N/EPB41L1 gene promoter region was highly methylated in LUAD and LUSC patients. LUAD patients with higher methylation level in the 4.1N/EPB41L1 gene promoter (TSS1500, cg13399773 or TSS200, cg20993403) had a shorter overall survival time (Log-rank p = 0.02 HR = 1.509 or Log-rank p = 0.016 HR = 1.509), whereas LUSC patients with higher methylation level in the 4.1N/EPB41L1 gene promoter (TSS1500 cg13399773, TSS1500 cg07030373 or TSS200 cg20993403) had a longer overall survival time (Log-rank p = 0.045 HR = 0.5709, Log-rank p = 0.018 HR = 0.68 or Log-rank p = 0.014 HR = 0.639, respectively). High methylation of the 4.1N/EPB41L1 gene promoter appeared to be a relatively early event in LUAD and LUSC. DNA methyltransferase inhibitor 5-Aza-2'-deoxycytidine restored the 4.1N/EPB41L1 expression at both the mRNA and protein levels. MiR-454-3p was abnormally highly expressed in NSCLC and directly targeted 4.1N/EPB41L1 mRNA. MiR-454-3p expression was significantly correlated with 4.1N/EPB41L1 expression in NSCLC patients (r = -0.63, p < 0.0001). Therefore, we concluded that promoter hypermethylation of the 4.1N/EPB41L1 gene and abnormally high expressed miR-454-3p work at different regulation levels but in concert to restrict 4.1N/EPB41L1 expression in NSCLC. Taken together, this work contributes to elucidate the underlying epigenetic disruptions of 4.1N/EPB41L1 deficiency in NSCLC.