Background: The radial spoke head component 4A (RSPH4A) is involved in the assembly of radial spokes, which is essential for motile cilia function. Asthenoteratozoospermia in primary ciliary dyskinesia (PCD) related to RSPH4A variants has not been reported. Materials and Methods: RSPH4A variants were identified and validated using whole-exome and Sanger sequencing in three unrelated Chinese families. High-speed video microscopy analysis (HSVA) was performed to measure the beating frequency and pattern of nasal cilia of the patients and healthy control. Papanicolaou staining and computer-aided sperm analysis were performed to analyze the morphology and motility of the sperm in patient 1. Immunofluorescence was adopted to confirm the structure deficiency of sperm and nasal cilia. Results: Patient 1 from family 1 is a 22-year-old unmarried male presented with bronchiectasis. Semen analysis and sperm Papanicolaou staining confirmed asthenoteratozoospermia. Novel compound heterozygous RSPH4A variants c.2T>C, p.(Met1Thr) and c.1774_1775del, p.(Leu592Aspfs*5) were detected in this patient. Patients 2 and 3 are from two unrelated consanguineous families; they are both females and exhibited bronchiectasis and infertility. Two homozygous RSPH4A variants c.2T>C, p.(Met1Thr) and c.351dupT, p.(Pro118Serfs*2) were detected, respectively. HSVA showed that most of the cilia in patients 1 and 3 were with abnormal rotational movement. The absence of RSPH4A and RSPH1 in patient 1's sperm and patient 3's respiratory cilia was indicated by immunofluorescence. Patient 2 died of pulmonary infection and respiratory failure at the age of 35 during follow-up. Conclusion: Dysfunctional sperm flagellum and motile cilia in the respiratory tract and the fallopian tube were found in patients with RSPH4A variants. Our study enriches the genetic spectrum and clinical phenotypes of RSPH4A variants in PCD, and c.2T>C, p.(Met1Thr) detected in our patients may be a hotspot RSPH4A variant in Chinese.

Type I interferons (IFNs) initiate immune responses to viral infections. Their effects are mediated by the type I IFN receptor, IFNAR, comprised of two subunits: IFNAR1 and IFNAR2. One or both chains of the sheep IFNAR were disrupted in fetal fibroblast lines using CRISPR/Cas9 and 12 lambs were produced by somatic cell nuclear transfer (SCNT). Quantitative reverse transcription-polymerase chain reaction for IFN-stimulated gene expression showed that IFNAR deficient sheep fail to respond to IFN-alpha. Furthermore, fibroblast cells from an IFNAR2 -/- fetus supported significantly higher levels of Zika virus (ZIKV) replication than wild-type fetal fibroblast cells. Although many lambs have died from SCNT related problems or infections, one fertile IFNAR2 -/- ram lived to over 4 years of age, remained healthy, and produced more than 80 offspring. Interestingly, ZIKV infection studies failed to demonstrate a high level of susceptibility. Presumably, these sheep compensated for a lack of type I IFN signaling using the type II, IFN-gamma and type III, IFN-lambda pathways. These sheep constitute a unique model for studying the pathogenesis of viral infection. Historical data supports the concept that ruminants utilize a novel type I IFN, IFN-tau, for pregnancy recognition. Consequently, IFNAR deficient ewes are likely to be infertile, making IFNAR knockout sheep a valuable model for studying pregnancy recognition. A breeding herd of 32 IFNAR2 +/- ewes, which are fertile, has been developed for production of IFNAR2 -/- sheep for both infection and reproduction studies.

A reducing sugar reacts with the protein, resulting in advanced glycation end-products (AGEs), which have been implicated in diabetes-related complications. Recently, it has been found that both type 1 and type 2 diabetic patients suffer from not only glucose but also ribose dysmetabolism. Here, we compared the effects of ribose and glucose glycation on epigenetics, such as histone methylation and demethylation. To prepare ribose-glycated (riboglycated) proteins, we incubated 150 μM bovine serum albumin (BSA) with 1 M ribose at different time periods, and we evaluated the samples by ELISAs, Western blot analysis, and cellular experiments. Riboglycated BSA, which was incubated with ribose for approximately 7 days, showed the strongest cytotoxicity, leading to a significant decrease in the viability of SH-SY5Y cells cultured for 24 h (IC50 = 1.5 μM). A global demethylation of histone 3 (H3K4) was observed in SH-SY5Y cells accompanied with significant increases in lysine-specific demethylase-1 (LSD1) and plant homeodomain finger protein 8 (PHF8) after treatment with riboglycated BSA (1.5 μM), but demethylation did not occur after treatment with glucose-glycated (glucoglycated) proteins or the ribose, glucose, BSA, and Tris-HCl controls. Moreover, a significant demethylation of H3K4, H3K4me3, and H3K4me2, but not H3K4me1, occurred in the presence of riboglycated proteins. A significant increase of formaldehyde was also detected in the medium of SH-SY5Y cells cultured with riboglycated BSA, further indicating the occurrence of histone demethylation. The present study provides a new insight into understanding an epigenetic mechanism of diabetes mellitus (DM) related to ribose metabolic disorders.

Background: Despite the significant survival benefits of anti-PD-1/PD-L1 immunotherapy, non-small cell lung cancer (NSCLC) remains one of the most common tumors and major causes of cancer-related deaths worldwide. Thus, there is an urgent need to identify new therapeutic targets for this refractory disease. Methods: In this study, microarray datasets GSE27262, GSE75037, GSE102287, and GSE21933 were integrated by Venn diagram. We performed functional clustering and pathway enrichment analyses using R. Through the STRING database and Cytoscape, we conducted protein-protein interaction (PPI) network analysis and identified the key genes, which were verified by the GEPIA2 and UALCAN portal. Validation of actin-binding protein anillin (ANLN) was performed by quantitative real-time polymerase chain reaction and Western blotting. Additionally, Kaplan-Meier methods were used to compute the survival analyses. Results: In total, 126 differentially expressed genes were identified, which were enriched in mitotic nuclear division, mitotic cell cycle G2/M transition, vasculogenesis, spindle, and peroxisome proliferator-activated receptor signaling pathway. 12 central node genes were identified in the PPI network complex. The survival analysis revealed that high transcriptional levels were associated with inferior survival in NSCLC patients. The clinical implication of ANLN was further explored; its protein expression showed a gradually increasing trend from grade I to III. Conclusion: These Key genes may be involved in the carcinogenesis and progression of NSCLC, which may serve as useful targets for NSCLC diagnosis and treatment.

Folate deficiency is associated with a broad range of human disorders, including anemia, fetal neural tube defects, age-associated dementia and several types of cancer. It is well established that a subgroup of rare fragile sites (RFSs) containing expanded CGG trinucleotide repeat (TNR) sequences display instability when cells are deprived of folate. However, given that folate sensitive RFSs exist in a very small percentage of the population, they are unlikely to be the cause of the widespread health problems associated with folate deficiency. We hypothesized that folate deficiency could specifically affect DNA replication at regions containing CG-rich repeat sequences. For this, we identified a region on human chromosome 2 (Chr2) comprising more than 300 CG-rich TNRs (termed "FOLD1") by examining the human genome database. Via the analysis of chromosome shape and segregation in mitosis, we demonstrate that, when human cells are cultured under folate stress conditions, Chr2 is prone to display a "kink" or "bending" at FOLD1 in metaphase and nondisjunction in anaphase. Furthermore, long-term folate deprivation causes Chr2 aneuploidy. Our results provide new evidence on the abnormalities folate deficiency could cause in human cells. This could facilitate future studies on the deleterious health conditions associated with folate deficiency.

Familial chylomicronemia syndrome (FCS) is a rare monogenic autosomal recessive disease caused by loss-of-function mutations in genes involved in chylomicron breakdown through hydrolysis of triglycerides into free fatty acids. Patients are often diagnosed in early childhood with extremely high triglyceride levels and symptoms including abdominal pain, eruptive cutaneous xanthomata, hepatosplenomegaly, and significant cognitive, psychological, and social impairment. The most serious medical condition suffered by FCS patients is recurrent acute pancreatitis. Lipoprotein lipase (LPL) gene mutation accounts for majority of the known pathogenic mutations. Early diagnosis and strict low-fat diet are critical for successful management of the triglyceride concentration to lower the risk of pancreatitis. The true prevalence of FCS in China is unknown and here we report a Chinese female preterm neonate presented with an extremely high triglyceride level of 22.11 mmol/L on day 13 after birth. Clinical and laboratory workup including whole-exome sequencing revealed two novel compound heterozygous LPL mutations (c.406G > C and c.829G > C) that are co-segregated with her non-consanguineous parents, consistent with autosomal recessive inheritance. A diagnosis of FCS based on clinical, biochemical, and genetic ground was made to guide her management.

Lung cancer is the most commonly diagnosed type of cancer and the leading cause of cancer-associated death worldwide. MicroRNAs (miRNAs) are non-coding single-stranded RNA molecules of ∼20-25 nucleotides in length. Single nucleotide polymorphisms are a class of genetic variation in the human genome, which when present in miRNA genes are associated with the risk of developing cancer. This study aimed to identify whether the miRNA (miR)-608 polymorphism rs4919510 influenced the incidence of lung cancer, and to explore the underlying mechanisms of miR-608 in the pathogenesis of the disease. A total of 37 patients with non-small cell lung cancer (NSCLC) were selected to determine the expression levels of miR-608; 96 NSCLC patients and 136 cancer-free healthy controls were recruited to determine the incidence of miR-608 rs4919510 in lung cancer patients. Additionally, the impact of miR-608 on the expression of predicted target genes, cell migration, viability, proliferation, and apoptosis was also assessed. We found that the presence of miR-608 rs4919510 did not affect the susceptibility of patients to NSCLC or the maturation of miR-608. miR-608 expression levels were found to be downregulated in NSCLC tissues. Furthermore, overexpression of miR-608 promoted doxorubicin-induced apoptosis in NSCLC cell lines A549 and HCC4006 by inhibiting the expression of transcription factor activating enhancer-binding protein 4 (TFAP4), and high expression levels of TFAP4 were observed in NSCLC tissues. Therefore, our results may provide valuable insights for the chemotherapeutical treatment of NSCLC.

Various methods have been proposed for genomic prediction (GP) in livestock. These methods have mainly focused on statistical considerations and did not include genome annotation information. In this study, to improve the predictive performance of carcass traits in Chinese Simmental beef cattle, we incorporated the genome annotation information into GP. Single nucleotide polymorphisms (SNPs) were annotated to five genomic classes: intergenic, gene, exon, protein coding sequences, and 3'/5' untranslated region. Haploblocks were constructed for all markers and these five genomic classes by defining a biologically functional unit, and haplotype effects were modeled in both numerical dosage and categorical coding strategies. The first-order epistatic effects among SNPs and haplotypes were modeled using a categorical epistasis model. For all makers, the extension from the SNP-based model to a haplotype-based model improved the accuracy by 5.4-9.8% for carcass weight (CW), live weight (LW), and striploin (SI). For the five genomic classes using the haplotype-based prediction model, the incorporation of gene class information into the model improved the accuracies by an average of 1.4, 2.1, and 1.3% for CW, LW, and SI, respectively, compared with their corresponding results for all markers. Including the first-order epistatic effects into the prediction models improved the accuracies in some traits and genomic classes. Therefore, for traits with moderate-to-high heritability, incorporating genome annotation information of gene class into haplotype-based prediction models could be considered as a promising tool for GP in Chinese Simmental beef cattle, and modeling epistasis in prediction can further increase the accuracy to some degree.

Mesocentrotus nudus is an important commercially aquatic species because of its high edible and medicinal values. However, wild stocks have dramatically decreased in recent decades. Understanding the population structure and genetic diversity can provide vital information for genetic conservation and improvement. In the present study, the genotyping-by-sequencing (GBS) approach was adopted to identify the genome-wide single-nucleotide polymorphisms (SNPs) from a collection of 80 individuals consisting of five geographical populations (16 individuals from each population), covering the natural habitats of M. nudus in China seas. An average of 0.96-Gb clean reads per sample were sequenced, and a total of 51,738 biallelic SNPs were identified. Based on these SNPs, diversity index analysis showed that all populations have a similar pattern with positive Fis (0.136) and low Ne (724.3). Low genetic differentiation and high genetic connectivity among five geographical populations were detected by pairwise Fst, principal component analysis (PCA), admixture, and phylogenetic analysis. Besides, two YWL individuals originating from an isolated ancestor may imply that there is a genetically differentiated population in the adjacent sea. Overall, the results showed that GBS is an effective method to detect genome-wide SNPs for M. nudus and suggested that the protective measures and the investigation with larger spatial scale and sample size for M. nudus should be carried out in the future.

Silicosis is a fatal occupational lung disease which currently has no effective clinical cure. Recent studies examining the underlying mechanism of silicosis have primarily examined experimental models, which may not perfectly reflect the nature of human silicosis progression. A comprehensive profiling of the molecular changes in human silicosis lungs is urgently needed. Here, we conducted RNA sequencing (RNA-seq) on the lung tissues of 10 silicosis patients and 7 non-diseased donors. A total of 2,605 differentially expressed genes (DEGs) and critical pathway changes were identified in human silicosis lungs. Further, the DEGs in silicosis were compared with those in idiopathic pulmonary fibrosis (IPF) and chronic obstructive pulmonary diseases (COPD), to extend current knowledge about the disease mechanisms and develop potential drugs. This analysis revealed both common and specific regulations in silicosis, along with several critical genes (e.g., MUC5AC and FGF10), which are potential drug targets for silicosis treatment. Drugs including Plerixafor and Retinoic acid were predicted as potential candidates in treating silicosis. Overall, this study provides the first transcriptomic fingerprint of human silicosis lungs. The comparative transcriptome analyses comprehensively characterize pathological regulations resulting from silicosis, and provide valuable cues for silicosis treatment.

During the carcinogenesis of cervical cancer, the DNA of human papillomavirus (HPV) is frequently integrated into the human genome, which might be a biomarker for the early diagnosis of cervical cancer. Although the detection sensitivity of virus infection status increased significantly through the Illumina sequencing platform, there were still disadvantages remain for further improvement, including the detection accuracy and the complex integrated genome structure identification, etc. Nanopore sequencing has been proven to be a fast yet accurate technique of detecting pathogens in clinical samples with significant longer sequencing length. However, the identification of virus integration sites, especially HPV integration sites was seldom carried out by using nanopore platform. In this study, we evaluated the feasibility of identifying HPV integration sites by nanopore sequencer. Specifically, we re-sequenced the integration sites of a previously published sample by both nanopore and Illumina sequencing. After analyzing the results, three points of conclusions were drawn: first, 13 out of 19 previously published integration sites were found from all three datasets (i.e., nanopore, Illumina, and the published data), indicating a high overlap rate and comparability among the three platforms; second, our pipeline of nanopore and Illumina data identified 66 unique integration sites compared with previous published paper with 13 of them being verified by Sanger sequencing, indicating the higher integration sites detection sensitivity of our results compared with published data; third, we established a pipeline which could be used in HPV integration site detection by nanopore sequencing data without doing error correction analysis. In summary, a new nanopore data analysis method was tested and proved to be reliable in integration sites detection compared with methods of existing Illumina data analysis pipeline with less sequencing data required. It provides a solid evidence and tool to support the potential application of nanopore in virus status identification.

Recent advances in genomics and proteomics generated a large amount of trans regulatory data such as those mediated by RNA binding proteins (RBPs) and microRNAs. Since many trans regulators target 3' UTR of mRNA transcripts, it is likely that there would be interactions, i.e., competitive or cooperative effect, among these trans factors. We compiled the available RBP and microRNA binding sites, mapped them to the mRNA transcripts, and correlated the binding data with mRNA expression data generated by The Cancer Genome Atlas (TCGA). We separated pairs of RBPs and microRNAs into three scenarios: those that have overlapping target sites on the same mRNA transcript (overlapping), those that have target sites on the same mRNA transcript but non-overlapping (neighboring), and those that do not target the same mRNA transcript (independent). Through a regression analysis on expression profiles, we indeed observed interaction effect between RBPs and microRNAs in the majority of the cancer expression data sets. We further discussed implication of such widespread interactions in the context of cancer and diseases.

Postnatal development and maturation of pineal gland is a highly dynamic period of tissue remodeling and phenotype maintenance, which is genetically controlled by programmed gene expression regulations. However, limited molecular characterization, particularly regarding long noncoding RNAs (lncRNA), is available for postnatal pineal at a whole transcriptome level. The present study first characterized the comprehensive pineal transcriptome profiles using strand-specific RNA-seq to illustrate the dynamic mRNA/lncRNA expression at three developmental stages (infancy, puberty, and adulthood). The results showed that 21,448 mRNAs and 8,166 novel lncRNAs were expressed in pig postnatal pineal gland. Among these genes, 3,573 mRNAs and 851 lncRNAs, including the 5-hydroxytryptamine receptors, exhibited significant dynamic regulation along maturation process, while the expression of homeobox genes didn't show significant differences. Gene Ontology analysis revealed that the differentially expressed genes (DEGs) were significantly enriched in ion transport and synaptic transmission, highlighting the critical role of calcium signaling in postnatal pineal development. Additionally, co-expression analysis revealed the DEGs could be grouped into 12 clusters with distinct expression patterns. Many differential lncRNAs were functionally enriched in co-expressed clusters of genes related to ion transport, transcription regulation, DNA binding, and visual perception. Our study first provided an overview of postnatal pineal transcriptome dynamics in pig and demonstrated that dynamic lncRNA regulation of developmental transitions impact pineal physiology.