Search

Background As a popular and valuable technique, grafting is widely used to protect against soil-borne diseases and nematodes in vegetable production. Growing evidences have revealed that long intergenic ncRNAs (lincRNAs) are strictly regulated and play essential roles in plants development and stress responses. Nevertheless, genome-wide identification and function deciphering of pepper lincRNAs, especially for their roles in improving grafting pepper resistance to Phytophthora capsici is largely unknown. Results In this study, RNA-seq data of grafting and control pepper plants with or without P. capsici inoculation were used to identify lincRNAs. In total, 2,388 reliable lincRNAs were identified. They were relatively longer and contained few exons than protein-coding genes. Similar to coding genes, lincRNAs had higher densities in euchromatin regions; and longer chromosome transcribed more lincRNAs. Expression pattern profiling suggested that lincRNAs commonly had lower expression than mRNAs. Totally, 607 differentially expressed lincRNAs (DE-lincRANs) were identified, of which 172 were found between P. capsici resistance grafting pepper sample GR and susceptible sample LDS. The neighboring genes of DE-lincRNAs and miRNAs competitively sponged by DE-lincRNAs were identified. Subsequently, the expression level of DE-lincRNAs was further confirmed by qRT-PCR and regulation patterns between DE-lincRNAs and neighboring mRNAs were also validated. Function annotation revealed that DE-lincRNAs increased the resistance of grafting prepper to P. capsici by modulating the expression of disease-defense related genes through cis-regulating and/or lincRNA-miRNA-mRNA interaction networks. Conclusions This study identified pepper lincRNAs and suggested their potential roles in increasing the resistance level of grafting pepper to P. capsici.

Background Macular corneal dystrophy (MCD) is a rare corneal stromal dystrophy with bilateral progressive vision loss. The pathogenic gene of MCD is carbohydrate sulfotransferase 6 (CHST6). Herein, we report a novel missense mutation and a rare exon deletion mutation in the CHST6 gene in a Chinese family with MCD. Methods Genomic DNA was extracted from the peripheral blood, and next generation sequencing was used to analyse the gene sequence. The pathogenic mutations were identified in all affected family members. The proband successively received binocular penetrating keratoplasty (PKP), and the corneas were examined by histopathology and colloidal iron staining to prove the diagnosis. A long-term follow-up was made to observe the changes after PKP. Results Genetic analysis demonstrated hemizygous mutations in the proband, including a novel c.520A>C (p.K174Q) missense mutation and a rarely reported exon 3 deletion mutation, which were co-segregated with the MCD phenotypes in the pedigree. The positive colloidal iron staining confirmed the diagnosis of MCD in the proband. However, the clinical phenotype and pathological manifestation of both eyes were different from each other because of complicated keratitis in the left eye. During the nine years of follow-up, visual acuity was improved significantly, and the cornea was transparent without rejection and postoperative recurrence in both eyes. Conclusions The novel hemizygous mutations were thought to contribute to the loss of CHST6 function, which induced typical clinical and pathological features of MCD. PKP was an effective treatment for MCD.

Background Amyotrophic lateral sclerosis (ALS) is a devastating neurodegenerative disease with neuronal cell inclusions composed of neurofilaments and other abnormal aggregative proteins as pathological hallmarks. Approximately 90% of patients have sporadic cases (sALS), and at least 4 genes, i.e. C9orf72, SOD1, FUS and TARDBP, have been identified as the main causative genes, while many others have been proposed as potential risk genes. However, these mutations could explain only ~ 10% of sALS cases. The neurofilament polypeptides encoded by NEFH, NEFM, and NEFL are promising protein biomarkers for ALS and other degenerative diseases. However, whether the genetic variants of these genes were associated with ALS remain ambiguous. Methods Here, we used PCR-Sanger to sequence the exons of these three genes in a cohort of 371 sALS patients and 711 healthy controls (Phase I) and validated the risk variant in another 300 sALS patients and 1076 controls (Phase II). Results A total of 92 variants were identified, including 36 rare heterozygous variants in NEFH, 27 in NEFM, and 16 in NEFL, and only rs568759161 (p.Ser787Arg) in NEFH reached nominal statistical power (P = 0.02 at Phase I, P = 0.009 at Phase II) in the case–control comparison. Together, the Phase I and II studies showed the significantly higher frequency of the variant in cases (9/1342, 0.67%) than in controls (2/3574, 0.07%) (OR 12.06; 95% CI 2.60–55.88; P = 0.0003). No variants passed multiple testing in the discovery cohort, but rs568759161 was associated with ALS in a replication cohort. Conclusions Our results confirmed that NEFH Ser787Arg is a novel sALS risk variant in Chinese subjects, but NEFM and NEFL were not associated with sALS. These data may have implications for genetic counselling and for understanding the pathogenesis of sALS.

The CACNA1H gene encodes the α1 subunit of the low voltage-activated Cav3.2 T-type calcium channel, an important regulator of neuronal excitability. Alternative mRNA splicing can generate multiple channel variants with distinct biophysical properties and expression patterns. Two major splice variants, containing or lacking exon 26 (± 26) have been found in different human tissues. In this study, we report splice variant specific effects of a Cav3.2 mutation found in patients with autosomal dominant writer’s cramp, a specific type of focal dystonia. We had previously reported that the R481C missense mutation caused a gain of function effect when expressed in Cav3.2 (+ 26) by accelerating its recovery from inactivation. Here, we show that when the mutation is expressed in the short variant of the channel (− 26), we observe a significant increase in current density when compared to wild-type Cav3.2 (− 26) but the effect on the recovery from inactivation is lost. Our data add to growing evidence that the functional expression of calcium channel mutations depends on which splice variant is being examined.

BackgroundThe characteristics of Waardenburg syndrome (WS) as a scarce heritable disorder are sensorineural hearing loss and deficits of pigmentation in the skin, hair, and eye. Here, clinical features and detection of the mutation in theMITFgene of WS2 patients are reported in a sizable Iranian family.MethodsA man aged 28-years represented with symptoms of mild unilateral hearing loss (right ear), complete heterochromia iridis, premature graying prior to 30 years of age, and synophrys. In this research, there was a sizable family in Iran comprising three generations with seven WS patients and two healthy members. Whole exome sequencing was applied for proband for the identification of the candidate genetic mutations associated with the disease. The detected mutation in proband and investigated family members was validated by PCR-Sanger sequencing.ResultsA novel heterozygous mutation, NM_198159.3:c.1026dup p.(Asn343Glufs*27), in exon 9 of theMITFgene co-segregated with WS2 in the affected family members. The variant was forecasted as a disease-causing variant by the Mutation Taster. According to the UniProt database, this variant has been located in basic helix-loop-helix (bHLH) domain of the protein with critical role in DNA binding.ConclusionsA frameshift was caused by a nucleotide insertion, c.1026dup, in exon 9 of theMITFgene. This mutation is able to induce an early termination, resulting in forming a truncated protein capable of affecting the normal function of the MITF protein. Helpful information is provided through an exactly described mutations involved in WS to clarify the molecular cause of clinical characteristics of WS and have a contribution to better genetic counseling of WS patients.

We develop a novel computational method, NucHMM, to identify functional nucleosome states associated with cell type-specific combinatorial histone marks and nucleosome organization features such as phasing, spacing and positioning. We test it on publicly available MNase-seq and ChIP-seq data in MCF7, H1, and IMR90 cells and identify 11 distinct functional nucleosome states. We demonstrate these nucleosome states are distinctly associated with the splicing potentiality of skipping exons. This advances our understanding of the chromatin function at the nucleosome level and offers insights into the interplay between nucleosome organization and splicing processes. Supplementary Information The online version contains supplementary material available at 10.1186/s13059-021-02465-1.

Background In a previous study, using a molecular approach, we reported the presence of P. vivax in Namibia. Here, we have extended our investigation to the Duffy antigen genetic profile of individuals of the same cohort with and without Plasmodium infections. Methods Participants with P. vivax (n = 3), P. falciparum (n = 23) mono-infections and co-infections of P. vivax/P. falciparum (n = 4), and P. falciparum/P. ovale (n = 3) were selected from seven regions. Participants with similar age but without any Plasmodium infections (n = 47) were also selected from all the regions. Duffy allelic profile was examined using standard PCR followed by sequencing of amplified products. Sequenced samples were also examined for the presence or absence of G125A mutation in codon 42, exon 2. Results All individuals tested carried the − 67 T > C mutation. However, while all P. vivax infected participants carried the c.G125A mutation, 7/28 P. falciparum infected participants and 9/41 of uninfected participants did not have the c.G125A mutation. The exon 2 region surrounding codon 42, had a c.136G > A mutation that was present in all P. vivax infections. The odds ratio for lack of this mutation with P. vivax infections was (OR 0.015, 95% CI 0.001–0.176; p = 0.001). Conclusion We conclude that P. vivax infections previously reported in Namibia, occurred in Duffy negative participants, carrying the G125A mutation in codon 42. The role of the additional mutation c.136 G > A in exon 2 in P. vivax infections, will require further investigations.

Background To analyze and evaluate EGFR, KRAS, and PIK3CA gene mutation rates and clinical distribution in patients with different types of lung cancer Method A total of 221 lung cancer patients treated in our hospital between January 2016 and June 2019 were enrolled. Tissue and whole blood samples were collected and analyzed to determine the mutation status of EGFR, KRAS, and PIK3CA genes. The gene exon mutation rates were determined. Relevant clinical data, such as age, gender, tumor sample type, treatment method, pathologic type, and lung cancer stage were recorded and statistically analyzed. Results The EGFR gene mutation rates in exons E18-E21 were 2.3%, 17.6%, 3.6%, and 20.4%, respectively. E18, E19, and E20 mutations were commonly detected in adenosquamous carcinoma, and E21 mutations were commonly detected in adenocarcinoma. Mutations in exons E18-E21 were frequently detected in patients with lung cancer stages IA, IB, IIA, or IIB, respectively. The KRAS gene mutation rate in lung cancer patients in exon E2 was higher in whole blood and tissue samples than other exon mutations, while the KRAS gene mutation rate in exons E2 and E3 was significantly higher in patients with lung cancer stages IIB and IA, respectively. PIK3CA gene mutations in exons E9 and E20 occurred in patients < 60 years of age. Exon E9-positive mutations were more common in men or patients with squamous cell carcinoma, while exon E20-positive mutations were more common in females. Conclusion The EGFR, KRAS, and PIK3CA gene exon mutation rates differ and were shown to be correlated with different clinical indicators, which have significance in clinical treatment.

Background Dermatofibrosarcoma protuberans (DFSP) is a superficial fibroblastic tumor characterized by high rate of local recurrence and low metastatic potential. Fibrosarcomatous transformation can rarely arise in DFSP either de novo or as recurrent, which represents a form of tumor progression and carries an increased risk of metastasis over classic DFSP. Cytogenetically, DFSP is characterized by a recurrent unbalanced chromosome translocation t (17;22)(q22;q13), leading to the formation of COL1A1-PDGFB fusion transcript that is present in more than 90% of cases. Alternative fusions involving the PDGFD with partners of COL6A3 or EMILIN2 have recently been documented in less than 2% of cases. Herein, we report a DFSP with fibrosarcomtous morphology harboring a novel TNC-PDGFD fusion. Case presentation A 54-year-old female presented with a slowly growing mass in the right thigh. Excision demonstrated a 2-cm ovoid, well-circumscribed, gray-white, mass. Microscopic examination revealed a partially encapsulated subcutaneous nodule without dermal connection. The neoplasm was composed of cellular and fairly uniform spindle cells with brisk mitoses, arranged in elongated fascicles and herringbone patterns, with focal collagenized stroma. The neoplastic cells were positive for CD34 and smooth muscle actin. Fluorescence in-situ hybridization analyses showed negative for COL1A1-PDGFB fusion as well as NTRK1/2/3 rearrangements. A subsequent RNA sequencing detected an in-frame fusion between exon 15 of TNC and exon 6 of PDGFD. This fusion was further confirmed by nested reverse transcription polymerase chain reaction amplification followed by Sanger sequencing. A diagnosis of fibrosarcomatous DFSP was rendered and the patient was in good status at a follow-up of 12 months after the operation. Conclusions We report a fibrosarcomatous DFSP with novel TNC-PDGFD fusion, which adds to the pathologic and genetic spectrum of PDGFD-rearranged DFSP.

Background Lipid phosphate phosphatases (LPP) are critical for regulating the production and degradation of phosphatidic acid (PA), an essential signaling molecule under stress conditions. Thus far, the LPP family genes have not been reported in rapeseed (Brassica napus L.). Results In this study, a genome-wide analysis was carried out to identify LPP family genes in rapeseed that respond to different stress conditions. Eleven BnLPPs genes were identified in the rapeseed genome. Based on phylogenetic and synteny analysis, BnLPPs were classified into four groups (Group I-Group IV). Gene structure and conserved motif analysis showed that similar intron/exon and motifs patterns occur in the same group. By evaluating cis-elements in the promoters, we recognized six hormone- and seven stress-responsive elements. Further, six putative miRNAs were identified targeting three BnLPP genes. Gene ontology analysis disclosed that BnLPP genes were closely associated with phosphatase/hydrolase activity, membrane parts, phosphorus metabolic process, and dephosphorylation. The qRT-PCR based expression profiles of BnLPP genes varied in different tissues/organs. Likewise, several gene expression were significantly up-regulated under NaCl, PEG, cold, ABA, GA, IAA, and KT treatments. Conclusions This is the first report to describe the comprehensive genome-wide analysis of the rapeseed LPP gene family. We identified different phytohormones and abiotic stress-associated genes that could help in enlightening the plant tolerance against phytohormones and abiotic stresses. The findings unlocked new gaps for the functional verification of the BnLPP gene family during stresses, leading to rapeseed improvement.

Background The sequence content of the 3′ UTRs of many mRNA transcripts is regulated through alternative polyadenylation (APA). The study of this process using RNAseq data, though, has been historically challenging. Results To combat this problem, we developed LABRAT, an APA isoform quantification method. LABRAT takes advantage of newly developed transcriptome quantification techniques to accurately determine relative APA site usage and how it varies across conditions. Using LABRAT, we found consistent relationships between gene-distal APA and subcellular RNA localization in multiple cell types. We also observed connections between transcription speed and APA site choice as well as tumor-specific transcriptome-wide shifts in APA isoform abundance in hundreds of patient-derived tumor samples that were associated with patient prognosis. We investigated the effects of APA on transcript expression and found a weak overall relationship, although many individual genes showed strong correlations between relative APA isoform abundance and overall gene expression. We interrogated the roles of 191 RNA-binding proteins in the regulation of APA isoforms, finding that dozens promote broad, directional shifts in relative APA isoform abundance both in vitro and in patient-derived samples. Finally, we find that APA site shifts in the two classes of APA, tandem UTRs and alternative last exons, are strongly correlated across many contexts, suggesting that they are coregulated. Conclusions We conclude that LABRAT has the ability to accurately quantify APA isoform ratios from RNAseq data across a variety of sample types. Further, LABRAT is able to derive biologically meaningful insights that connect APA isoform regulation to cellular and molecular phenotypes.

Background Primary adrenal insufficiency (PAI) is life-threatening, and a definitive aetiological diagnosis is essential for management and prognostication. We conducted this study to investigate the genetic aetiologies of PAI in South China and explore their clinical features. Methods Seventy children were enrolled in this cross-sectional study. Clinical information was collected, and combined genetic tests were performed according to the children’s manifestations. Statistical analysis was performed among the different groups. In silico or in vitro experiments were applied to determine the pathogenicity of novel variants. Results Among the 70 children, 84.3% (59/70) were diagnosed with congenital adrenal hyperplasia (CAH), and 21-hydroxylase deficiency (21-OHD) was genetically confirmed in 91.5% of these cases. Salt wasting (SW), simple virilization (SV), and non-classic (NC) CAH accounted for 66.1% (39/59), 30.5% (18/59), and 3.4% (2/59) of the cases, respectively. The 17-hydroxyprogesterone (17-OHP) and testosterone (TES) levels were significantly higher in children with SW than with SV. The 17-OHP and cortisol levels in female SW patients were significantly higher than those in males. The 17-OHP, cortisol, dehydroepiandrosterone (DHEAS) and TES levels in female SW patients were significantly higher than those in female SV patients. Additionally, 72.7% (8/11) of uncharacterized PAI patients had positive genetic findings. Among all the patients, two novel variants in the CYP21A2 gene (c.833dupT and c.651 + 2T > G) were found. A microdeletion (Xp21.2–21.3) and five novel variants, including 2 in the NR0B1 gene (c.323–324CG > GA and c.1231_1234delCTCA), 2 in the AAAS gene (c.399 + 1G > A and c.250delT) and 1 in the NNT gene (c.2274delT), were detected. The novel variant c.399 + 1G > A in the AAAS gene was further confirmed to lead to exon 4 skipping during mRNA transcription and produce a truncated ALADIN protein. Conclusions We found ethnicity-based differences in the CYP21A2 gene variant spectrum among different study populations. Female 21-OHD patients tended to have higher 17-OHP and TES levels, which warrants caution in relation to the effects of virilization. Novel gene variants detected in the CYP21A2, NR0B1, AAAS and NNT genes expanded the genetic spectrum of PAI, however, further improvement of genetic testing tools beyond our protocol are still needed to uncover the complete aetiology of PAI in children.

Background Intensive evidence has highlighted the effect of aberrant alternative splicing (AS) events on cancer progression when triggered by dysregulation of the SR protein family. Nonetheless, the underlying mechanism in breast cancer (BRCA) remains elusive. Here we sought to explore the molecular function of SRSF1 and identify the key AS events regulated by SRSF1 in BRCA. Methods We conducted a comprehensive analysis of the expression and clinical correlation of SRSF1 in BRCA based on the TCGA dataset, Metabric database and clinical tissue samples. Functional analysis of SRSF1 in BRCA was conducted in vitro and in vivo. SRSF1-mediated AS events and their binding motifs were identified by RNA-seq, RNA immunoprecipitation-PCR (RIP-PCR) and in vivo crosslinking followed by immunoprecipitation (CLIP), which was further validated by the minigene reporter assay. PTPMT1 exon 3 (E3) AS was identified to partially mediate the oncogenic role of SRSF1 by the P-AKT/C-MYC axis. Finally, the expression and clinical significance of these AS events were validated in clinical samples and using the TCGA database. Results SRSF1 expression was consistently upregulated in BRCA samples, positively associated with tumor grade and the Ki-67 index, and correlated with poor prognosis in a hormone receptor-positive (HR+) cohort, which facilitated proliferation, cell migration and inhibited apoptosis in vitro and in vivo. We identified SRSF1-mediated AS events and discovered the SRSF1 binding motif in the regulation of splice switching of PTPMT1. Furthermore, PTPMT1 splice switching was regulated by SRSF1 by binding directly to its motif in E3 which partially mediated the oncogenic role of SRSF1 by the AKT/C-MYC axis. Additionally, PTPMT1 splice switching was validated in tissue samples of BRCA patients and using the TCGA database. The high-risk group, identified by AS of PTPMT1 and expression of SRSF1, possessed poorer prognosis in the stage I/II TCGA BRCA cohort. Conclusions SRSF1 exerts oncogenic roles in BRCA partially by regulating the AS of PTPMT1, which could be a therapeutic target candidate in BRCA and a prognostic factor in HR+ BRCA patient.

Background Defects in interleukin 10 (IL10) and its receptors are particularly involved in very early onset inflammatory bowel disease (VEOIBD). However, large fragment deletions of IL10 receptor A (IL10RA) are rare. Methods VEOIBD patients with confirmed mutations in the IL10RA gene were enrolled from January 1, 2019 to June 30, 2020. The clinical features and endoscopic-radiological findings of the patients with large fragment deletions of the IL10RA gene were determined and followed up. Results Thirty-five patients with IL10RA gene mutations, namely, 28 compound heterozygous mutations and 7 homozygote mutations, were enrolled in this study. Six patients carried the reported point mutation c.301C > T (p. R101RW) or c.537 G > A (p. T179T) in one locus and a large fragment deletion in exon 1 in another locus, which were novel mutations in this gene. A 333-bp deletion of exon 1 (117857034–11857366 del) was the main mutation in this locus in 85.7% of the patients with large fragment deletions. The time of disease onset ranged from birth to 4 years, and diarrhea was the main initial symptom. In total, 6/7 patients had perianal complications, including perianal abscess, fistula and skin tags. Six patients accepted thalidomide treatment, 5/7 accepted mesalamine, 3/7 accepted hematopoietic stem cell transplantation (HSCT), and 3/7 were waiting for HSCT. Conclusions We identified a novel large deletion of exon 1 involving the IL10RA gene for the first time and showed the characteristics of VEOIBD patients. This study expands the spectrum of Chinese VEOIBD patients with IL0RA gene mutations.

Background Apoptosis-inducing factor (AIF), as a mitochondrial flavoprotein, plays a fundamental role in mitochondrial bioenergetics that is critical for cell survival and also mediates caspase-independent cell death once it is released from mitochondria and translocated to the nucleus under ischemic stroke or neurodegenerative diseases. Although alternative splicing regulation of AIF has been implicated, it remains unknown which AIF splicing isoform will be induced under pathological conditions and how it impacts mitochondrial functions and neurodegeneration in adult brain. Methods AIF splicing induction in brain was determined by multiple approaches including 5′ RACE, Sanger sequencing, splicing-specific PCR assay and bottom-up proteomic analysis. The role of AIF splicing in mitochondria and neurodegeneration was determined by its biochemical properties, cell death analysis, morphological and functional alterations and animal behavior. Three animal models, including loss-of-function harlequin model, gain-of-function AIF3 knockin model and conditional inducible AIF splicing model established using either Cre-loxp recombination or CRISPR/Cas9 techniques, were applied to explore underlying mechanisms of AIF splicing-induced neurodegeneration. Results We identified a nature splicing AIF isoform lacking exons 2 and 3 named as AIF3. AIF3 was undetectable under physiological conditions but its expression was increased in mouse and human postmortem brain after stroke. AIF3 splicing in mouse brain caused enlarged ventricles and severe neurodegeneration in the forebrain regions. These AIF3 splicing mice died 2–4 months after birth. AIF3 splicing-triggered neurodegeneration involves both mitochondrial dysfunction and AIF3 nuclear translocation. We showed that AIF3 inhibited NADH oxidase activity, ATP production, oxygen consumption, and mitochondrial biogenesis. In addition, expression of AIF3 significantly increased chromatin condensation and nuclear shrinkage leading to neuronal cell death. However, loss-of-AIF alone in harlequin or gain-of-AIF3 alone in AIF3 knockin mice did not cause robust neurodegeneration as that observed in AIF3 splicing mice. Conclusions We identified AIF3 as a disease-inducible isoform and established AIF3 splicing mouse model. The molecular mechanism underlying AIF3 splicing-induced neurodegeneration involves mitochondrial dysfunction and AIF3 nuclear translocation resulting from the synergistic effect of loss-of-AIF and gain-of-AIF3. Our study provides a valuable tool to understand the role of AIF3 splicing in brain and a potential therapeutic target to prevent/delay the progress of neurodegenerative diseases.

BackgroundTrichomes are excellent model systems for the analysis of cell differentiation and play essential roles in plant protection. From cucumber inbred line ‘WD1’, we identified an EMS-induced trichome abnormally developing mutant,nps, which exhibited smaller, denser and no pyramid-shaped head trichomes.ResultsUsing F2and BC1populations constructed from a cross betweennpsand ‘9930’, the genetic analysis showed that thenpstrait is controlled by a single recessive nuclear gene. We identifiedCsNpsby map-based cloning with 576 individuals of the F2population generated from the cross ofnpsand inbred line ‘9930’. TheCsNpswas located at a 13.4-kb genomic region on chromosome 3, which region contains three predicted genes. Sequence analysis showed that only one single nucleotide mutation (C → T) between 9930 andnpswas found in the second exon ofCsa3G748220, a plant-specific class I HD-Zip gene. The result of allelism test also indicated thatnpsis a novel allelic mutant ofMict(Micro-trichome). Thus,npswas renamedmict-L130F. By comparing the transcriptome ofmict-L130Fvs WD1 and 06–2 (mict) vs 06–1 (wildtype, near-isogenic line of 06–2), several potential target genes that may be related to trichome development were identified.ConclusionsOur results demonstrate thatMict-L130Fis involved in the morphogenesis of trichomes. Map-based cloning of theMict-L130Fgene could promote the study of trichome development in cucumber.

Background Mutations in the DMD gene encoding dystrophin—a critical structural element in muscle cells—cause Duchenne muscular dystrophy (DMD), which is the most common fatal genetic disease. Clustered regularly interspaced short palindromic repeat (CRISPR)-mediated gene editing is a promising strategy for permanently curing DMD. Methods In this study, we developed a novel strategy for reframing DMD mutations via CRISPR-mediated large-scale excision of exons 46–54. We compared this approach with other DMD rescue strategies by using DMD patient-derived primary muscle-derived stem cells (DMD-MDSCs). Furthermore, a patient-derived xenograft (PDX) DMD mouse model was established by transplanting DMD-MDSCs into immunodeficient mice. CRISPR gene editing components were intramuscularly delivered into the mouse model by adeno-associated virus vectors. Results Results demonstrated that the large-scale excision of mutant DMD exons showed high efficiency in restoring dystrophin protein expression. We also confirmed that CRISPR from Prevotella and Francisella 1(Cas12a)-mediated genome editing could correct DMD mutation with the same efficiency as CRISPR-associated protein 9 (Cas9). In addition, more than 10% human DMD muscle fibers expressed dystrophin in the PDX DMD mouse model after treated by the large-scale excision strategies. The restored dystrophin in vivo was functional as demonstrated by the expression of the dystrophin glycoprotein complex member β-dystroglycan. Conclusions We demonstrated that the clinically relevant CRISPR/Cas9 could restore dystrophin in human muscle cells in vivo in the PDX DMD mouse model. This study demonstrated an approach for the application of gene therapy to other genetic diseases.

Background Cannabis, an important industrial crop, has a high sensitivity to photoperiods. The flowering time of cannabis is one of its important agronomic traits, and has a significant effect on its yield and quality. The CONSTANS-like (COL) gene plays a key role in the regulation of flowering in this plant. However, the specific roles of the COL gene family in cannabis are still unknown. Results In this study, 13 CsCOL genes were identified in the cannabis genome. Phylogenetic analysis implied that the CsCOL proteins were divided into three subgroups, and each subgroup included conserved intron/exon structures and motifs. Chromosome distribution analysis showed that 13 CsCOL genes were unevenly distributed on 7 chromosomes, with chromosome 10 having the most CsCOL members. Collinearity analysis showed that two syntenic gene pairs of CsCOL4 and CsCOL11 were found in both rice and Gossypium raimondii. Of the 13 CsCOL genes, CsCOL6 and CsCOL12 were a pair of tandem duplicated genes, whereas CsCOL8 and CsCOL11 may have resulted from segmental duplication. Furthermore, tissue-specific expression showed that 10 CsCOL genes were preferentially expressed in the leaves, 1 CsCOL in the stem, and 2 CsCOL in the female flower. Most CsCOL exhibited a diurnal oscillation pattern under different light treatment. Additionally, sequence analysis showed that CsCOL3 and CsCOL7 exhibited amino acid differences among the early-flowering and late flowering cultivars. Conclusion This study provided insight into the potential functions of CsCOL genes, and highlighted their roles in the regulation of flowering time in cannabis. Our results laid a foundation for the further elucidation of the functions of COL genes in cannabis.

SORL1 is strongly associated with both sporadic and familial forms of Alzheimer’s disease (AD), but a lack of information about alternatively spliced transcripts currently limits our understanding of the role of SORL1 in AD. Here, we describe a SORL1 transcript (SORL1-38b) characterized by inclusion of a novel exon (E38b) that encodes a truncated protein. We identified E38b-containing transcripts in several brain regions, with the highest expression in the cerebellum and showed that SORL1-38b is largely located in neuronal dendrites, which is in contrast to the somatic distribution of transcripts encoding the full-length SORLA protein (SORL1-fl). SORL1-38b transcript levels were significantly reduced in AD cerebellum in three independent cohorts of postmortem brains, whereas no changes were observed for SORL1-fl. A trend of lower 38b transcript level in cerebellum was found for individuals carrying the risk variant at rs2282649 (known as SNP24), although not reaching statistical significance. These findings suggest synaptic functions for SORL1-38b in the brain, uncovering novel aspects of SORL1 that can be further explored in AD research.

Background Extensive molecular differences exist between proliferative and differentiated cells. Here, we conduct a meta-analysis of publicly available transcriptomic datasets from preimplantation and differentiation stages examining the architectural properties and content of genes whose abundance changes significantly across developmental time points. Results Analysis of preimplantation embryos from human and mouse showed that short genes whose introns are enriched in Alu (human) and B (mouse) elements, respectively, have higher abundance in the blastocyst compared to the zygote. These highly expressed genes encode ribosomal proteins or metabolic enzymes. On the other hand, long genes whose introns are depleted in repetitive elements have lower abundance in the blastocyst and include genes from signaling pathways. Additionally, the sequences of the genes that are differentially expressed between the blastocyst and the zygote contain distinct collections of pyknon motifs that differ between up- and down-regulated genes. Further examination of the genes that participate in the stem cell-specific protein interaction network shows that their introns are short and enriched in Alu (human) and B (mouse) elements. As organogenesis progresses, in both human and mouse, we find that the primarily short and repeat-rich expressed genes make way for primarily longer, repeat-poor genes. With that in mind, we used a machine learning-based approach to identify gene signatures able to classify human adult tissues: we find that the most discriminatory genes comprising these signatures have long introns that are repeat-poor and include transcription factors and signaling-cascade genes. The introns of widely expressed genes across human tissues, on the other hand, are short and repeat-rich, and coincide with those with the highest expression at the blastocyst stage. Conclusions Protein-coding genes that are characteristic of each trajectory, i.e., proliferation/pluripotency or differentiation, exhibit antithetical biases in their intronic and exonic lengths and in their repetitive-element content. While the respective human and mouse gene signatures are functionally and evolutionarily conserved, their introns and exons are enriched or depleted in organism-specific repetitive elements. We posit that these organism-specific repetitive sequences found in exons and introns are used to effect the corresponding genes’ regulation.

Background Our previous study have shown that the PSMD11 protein was an important survival factor for cancer cells except for its key role in regulation of assembly and activity of the 26S proteasome. To further investigate the role of PSMD11 in carcinogenesis, we constructed a conditional exon 5 floxed allele of PSMD11 (PSMD11flx) in mice. Results It was found that homozygous PSMD11 flx/flx mice showed normal and exhibited a normal life span and fertility, and showed roughly equivalent expression of PSMD11 in various tissues, suggesting that the floxed allele maintained the wild-type function. Cre recombinase could induce efficient knockout of the floxed PSMD11 allele both in vitro and in vivo. Mice with constitutive single allele deletion of PSMD11 derived from intercrossing between PSMD11flx/flx and CMV-Cre mice were all viable and fertile, and showed apparent growth retardation, suggesting that PSMD11 played a significant role in the development of mice pre- or postnatally. No whole-body PSMD11 deficient embryos (PSMD11−/−) were identified in E7.5–8.5 embryos in uteros, indicating that double allele knockout of PSMD11 leads to early embryonic lethality. To avoid embryonic lethality produced by whole-body PSMD11 deletion, we further developed conditional PSMD11 global knockout mice with genotype Flp;FSF-R26CAG − CreERT2/+; PSMD11flx/flx, and demonstrated that PSMD11 could be depleted in a temporal and tissue-specific manner. Meanwhile, it was found that depletion of PSMD11 could induce massive apoptosis in MEFs. Conclusions In summary, our data demonstrated that we have successfully generated a conditional knockout allele of PSMD11 in mice, and found that PSMD11 played a key role in early and postnatal development in mice, the PSMD11 flx/flx mice will be an invaluable tool to explore the functions of PSMD11 in development and diseases.

Background Heat stress (HS) is a major stress event in the life of an animal, with detrimental upshots in production and health. Long-non-coding RNAs (lncRNAs) play an important role in many biological processes by transcriptional regulation. However, no research has been reported on the characterization and functionality of lncRNAs in heat-stressed rats. Results We studied expression levels of lncRNAs in rats during HS, using strand-specific RNA sequencing. Six rats, three in each of Control (22 ± 1 °C) and H120 (42 °C for 120 min) experimental groups, were used to screen for lncRNAs in their liver and adrenal glands. Totally, 4498 and 7627 putative lncRNAs were identified in liver and adrenal glands of the Control and H120 groups, respectively. The majority of lncRNAs were relatively shorter and contained fewer exons than protein-coding transcripts. In total, 482 (174 up-regulated and 308 down-regulated) and 271 (126 up-regulated and 145 down-regulated) differentially-expressed lncRNAs (DElncRNAs, P trans−/cis- and sequence similarity regulatory modes. Functional annotation analyses indicated that these DEGs were mostly significantly enriched in insulin signalling, myeloid leukaemia, and glucagon signalling pathways. Similarly, 437, 73 and 41 target DEGs in the adrenal glands were mostly significantly enriched in the cell cycle (trans-prediction) and lysosome pathways (cis-prediction). The DElncRNAs interacting with DEGs that encode heat shock proteins (HSPs) may play an important role in HS response, which include Hsf4, Dnaja1, Dnajb4, Hsph1 and Hspb1 in the liver, and Dnajb13 and Hspb8 in the adrenal glands. The strand-specific RNA sequencing findings were also further verified through RT-qPCR. Conclusions This study is the first to provide a detailed characterization and functional analysis of expression levels of lncRNAs in liver and adrenal glands of heat-stressed rats, which provides basis for further studies on the biological functions of lncRNAs under heat stress in rats and other mammalian species.

Background Rapid diagnostic tests (RDTs) targeting histidine rich protein 2(HRP2) are widely used for diagnosis of Plasmodium falciparum infections. Besides PfHRP2, the PfHRP3 antigen contributes to the detection of P. falciparum infections in PfHRP2 RDTs. However, the performance HRP2-based RDT is affected by pfhrp2/3 gene deletions resulting in false-negative test results. The objective of this study was to determine the presence and prevalence of pfhrp2/3 gene deletions including the respective flanking regions among symptomatic patients in Assosa zone, Northwest Ethiopia. Methods A health-facility based cross-sectional study was conducted in febrile patients seeking a malaria diagnosis in 2018. Blood samples were collected by finger-prick for microscopic examination of blood smears, malaria RDT, and molecular analysis using dried blood spots (DBS) prepared on Whatman filter paper. A total of 218 P. falciparum positive samples confirmed by quantitative PCR were included for molecular assay of pfhrp2/3 target gene. Results Of 218 P. falciparum positive samples, exon 2 deletions were observed in 17.9% of pfhrp2 gene and in 9.2% of pfhrp3 gene. A high proportion of deletions in short segments of pfhrp2 exon1-2 (50%) was also detected while the deletions of the pfhrp3 exon1-2 gene were 4.1%. The deletions were extended to the downstream and upstream of the flanking regions in pfhrp2/3 gene (above 30%). Of eighty-six PfHRP2 RDT negative samples, thirty-six lacked pfhrp2 exon 2. Five PfHRP2 RDT negative samples had double deletions in pfhrp2 exon 2 and pfhrp3 exon2. Of these double deletions, only two of the samples with a parasite density above 2000 parasite/µl were positive by the microscopy. Three samples with intact pfhrp3 exon2 in the pfhrp2 exon2 deleted parasite isolates were found to be positive by PfHRP2 RDT and microscopy with a parasite density above 10,000/µl. Conclusion This study confirms the presence of deletions of pfhrp2/3 gene including the flanking regions. Pfhrp2/3 gene deletions results in false-negative results undoubtedly affect the current malaria control and elimination effort in the country. However, further countrywide investigations are required to determine the magnitude of pfhrp2/3 gene deletions and its consequences on routine malaria diagnosis.

Background Variants that regulate transcription, such as expression quantitative trait loci (eQTL), have shown enrichment in genome-wide association studies (GWAS) for mammalian complex traits. However, no study has reported eQTL in sheep, although it is an important agricultural species for which many GWAS of complex meat traits have been conducted. Using RNA sequence data produced from liver and muscle from 149 sheep and imputed whole-genome single nucleotide polymorphisms (SNPs), our aim was to dissect the genetic architecture of the transcriptome by associating sheep genotypes with three major molecular phenotypes including gene expression (geQTL), exon expression (eeQTL) and RNA splicing (sQTL). We also examined these three types of eQTL for their enrichment in GWAS of multi-meat traits and fatty acid profiles. Results Whereas a relatively small number of molecular phenotypes were significantly heritable (h2 > 0, P P P FAM184B gene, hundreds of sQTL in liver and muscle mapped to the CAST gene, and hundreds of sQTL in liver mapped to the C6 gene. These three genes are associated with body composition or fatty acid profiles. Conclusions We detected a large number of significant eQTL and found that the overlap of variants between eQTL types and tissues was prevalent. Many eQTL were also QTL for meat traits. Our study fills a gap in the knowledge on the regulatory variants and their role in complex traits for the sheep model.

Background Microexons, exons that are ≤ 30 nucleotides, are a highly conserved and dynamically regulated set of cassette exons. They have key roles in nervous system development and function, as evidenced by recent results demonstrating the impact of microexons on behaviour and cognition. However, microexons are often overlooked due to the difficulty of detecting them using standard RNA-seq aligners. Results Here, we present MicroExonator, a novel pipeline for reproducible de novo discovery and quantification of microexons. We process 289 RNA-seq datasets from eighteen mouse tissues corresponding to nine embryonic and postnatal stages, providing the most comprehensive survey of microexons available for mice. We detect 2984 microexons, 332 of which are differentially spliced throughout mouse embryonic brain development, including 29 that are not present in mouse transcript annotation databases. Unsupervised clustering of microexons based on their inclusion patterns segregates brain tissues by developmental time, and further analysis suggests a key function for microexons in axon growth and synapse formation. Finally, we analyse single-cell RNA-seq data from the mouse visual cortex, and for the first time, we report differential inclusion between neuronal subpopulations, suggesting that some microexons could be cell type-specific. Conclusions MicroExonator facilitates the investigation of microexons in transcriptome studies, particularly when analysing large volumes of data. As a proof of principle, we use MicroExonator to analyse a large collection of both mouse bulk and single-cell RNA-seq datasets. The analyses enabled the discovery of previously uncharacterized microexons, and our study provides a comprehensive microexon inclusion catalogue during mouse development.

BackgroundFamilial progressive hyper- and hypopigmentation (FPHH, MIM 145250) is a rare hereditary skin disorder that is predominantly characterized by progressive, diffuse, partly blotchy hyperpigmented lesions intermingled with scattered hypopigmented spots, lentigines and sometimes Cafe-au-lait spots (CALs). Heterozygous mutations of the KIT ligand (KITLG, MIM 184745) gene are responsible for FPHH. To date, only eightKITLGmutations have been reported to be associated with FPHH, and no clear genotype–phenotype correlations have been established. This study aimed to identify the causative mutations in theKITLGgene in two Chinese FPHH patients.MethodsDirect sequencing of the coding regions ofKITLGwas performed. Pathogenicity prediction was performed using bioinformatics tools, including SIFT, Polyphen2, and SWISS-MODEL, and the results were further evaluated according to the 2015 American College of Medical Genetics and Genomics (ACMG) guidelines.ResultsThe novel mutation c.104A > T (p.Asn35Ile) and the recurrent mutation c.101C > T (p.Thr34Ile) inKITLGwere identified. As shown using SIFT and Polyphen-2 software, both mutations identified in this study were predicted to be detrimental variations. Three-dimensional protein structure modeling indicated that the mutant KITLG proteins might affect the affinity of KITLG for its receptor, c-KIT. According to the 2015 ACMG guidelines, the novel mutation c.104A > T was ‘likely pathogenic’.ConclusionsTo date, most of the identifiedKITLGmutations have been clustered within the conserved VTNNV motif (amino acids 33–37) in exon 2. The known mutations are only involved in 33 V, 34 T, 36 N, and 37 V but not 35 N. We have now identified a novel mutation inKITLG, c.104A > T, that was first reported in FPHH within the conserved 35 N motif. These results strengthen our understanding of FPHH and expand the mutational spectrum of theKITLGgene.

De novo and acquired resistance, which are mainly mediated by genetic alterations, are barriers to effective routine chemotherapy. However, the mechanisms underlying gastric cancer (GC) resistance to chemotherapy are still unclear. We showed that the long noncoding RNA CRNDE was related to the chemosensitivity of GC in clinical samples and a PDX model. CRNDE was decreased and inhibited autophagy flux in chemoresistant GC cells. CRNDE directly bound to splicing protein SRSF6 to reduce its protein stability and thus regulate alternative splicing (AS) events. We determined that SRSF6 regulated the PICALM exon 14 skip splice variant and triggered a significant S-to-L isoform switch, which contributed to the expression of the long isoform of PICALM (encoding PICALML). Collectively, our findings reveal the key role of CRNDE in autophagy regulation, highlighting the significance of CRNDE as a potential prognostic marker and therapeutic target against chemoresistance in GC. Supplementary Information The online version contains supplementary material available at 10.1186/s12943-020-01299-y.

Background Patterns of methylation influence lifespan, but methylation and lifespan may also depend on diet, or differ between genotypes. Prior to this study, interactions between diet and genotype have not been explored together to determine their influence on methylation. The invertebrate Daphnia magna is an excellent choice for testing the epigenetic response to the environment: parthenogenetic offspring are identical to their siblings (making for powerful genetic comparisons), they are relatively short lived and have well-characterised inter-strain life-history trait differences. We performed a survival analysis in response to caloric restriction and then undertook a 47-replicate experiment testing the DNA methylation response to ageing and caloric restriction of two strains of D. magna. Results Methylated cytosines (CpGs) were most prevalent in exons two to five of gene bodies. One strain exhibited a significantly increased lifespan in response to caloric restriction, but there was no effect of food-level CpG methylation status. Inter-strain differences dominated the methylation experiment with over 15,000 differently methylated CpGs. One gene, Me31b, was hypermethylated extensively in one strain and is a key regulator of embryonic expression. Sixty-one CpGs were differentially methylated between young and old individuals, including multiple CpGs within the histone H3 gene, which were hypermethylated in old individuals. Across all age-related CpGs, we identified a set that are highly correlated with chronological age. Conclusions Methylated cytosines are concentrated in early exons of gene sequences indicative of a directed, non-random, process despite the low overall DNA methylation percentage in this species. We identify no effect of caloric restriction on DNA methylation, contrary to our previous results, and established impacts of caloric restriction on phenotype and gene expression. We propose our approach here is more robust in invertebrates given genome-wide CpG distributions. For both strain and ageing, a single gene emerges as differentially methylated that for each factor could have widespread phenotypic effects. Our data showed the potential for an epigenetic clock at a subset of age positions, which is exciting but requires confirmation.

BackgroundSerum alpha-fetoprotein (AFP) is the approved serum marker for hepatocellular carcinoma (HCC) screening. However, not all HCC patients show high (≥ 20 ng/mL) serum AFP, and the molecular mechanisms of HCCs with normal (MethodsWe performed a genome-wide survey of AS events in 249 HCCs with normal (n = 131) and high (n = 118) serum AFP levels using RNA-sequencing data obtained from The Cancer Genome Atlas.ResultsIn group comparisons of RNA-seq profiles from HCCs with normal and high serum AFP levels, 161 differential AS events (125 genes; ΔPSI > 0.05, FDR FN1andFAM20A) directly interact with AFP; these relate to the regulation of IGF transport and post-translational protein phosphorylation. Interestingly, 42 genes and 27 genes were associated with gender and vascular invasion (VI), respectively, but only eighteen genes were significant in survival analysis. We especially highlight thatFN1exhibited increased differential expression of AS events (ΔPSI > 0.05), in which exons 25 and 33 were more frequently skipped in HCCs with normal (low) serum AFP compared to those with high serum AFP. Moreover, these events were gender and VI dependent.ConclusionWe found that AS may influence the regulation of transcriptional differences inherent in the occurrence of HCC maintaining normal rather than elevated serum AFP levels.

Background Alternative splicing isoforms have been reported as a new and robust class of diagnostic biomarkers. Over 95% of human genes are estimated to be alternatively spliced as a powerful means of producing functionally diverse proteins from a single gene. The emergence of next-generation sequencing technologies, especially RNA-seq, provides novel insights into large-scale detection and analysis of alternative splicing at the transcriptional level. Advances in Proteomic Technologies such as liquid chromatography coupled tandem mass spectrometry (LC–MS/MS), have shown tremendous power for the parallel characterization of large amount of proteins in biological samples. Although poor correspondence has been generally found from previous qualitative comparative analysis between proteomics and microarray data, significantly higher degrees of correlation have been observed at the level of exon. Combining protein and RNA data by searching LC–MS/MS data against a customized protein database from RNA-Seq may produce a subset of alternatively spliced protein isoform candidates that have higher confidence. Results We developed a bioinformatics workflow to discover alternative splicing biomarkers from LC–MS/MS using RNA-Seq. First, we retrieved high confident, novel alternative splicing biomarkers from the breast cancer RNA-Seq database. Then, we translated these sequences into in silico Isoform Junction Peptides, and created a customized alternative splicing database for MS searching. Lastly, we ran the Open Mass spectrometry Search Algorithm against the customized alternative splicing database with breast cancer plasma proteome. Twenty six alternative splicing biomarker peptides with one single intron event and one exon skipping event were identified. Further interpretation of biological pathways with our Integrated Pathway Analysis Database showed that these 26 peptides are associated with Cancer, Signaling, Metabolism, Regulation, Immune System and Hemostasis pathways, which are consistent with the 256 alternative splicing biomarkers from the RNA-Seq. Conclusions This paper presents a bioinformatics workflow for using RNA-seq data to discover novel alternative splicing biomarkers from the breast cancer proteome. As a complement to synthetic alternative splicing database technique for alternative splicing identification, this method combines the advantages of two platforms: mass spectrometry and next generation sequencing and can help identify potentially highly sample-specific alternative splicing isoform biomarkers at early-stage of cancer.

Background Eukaryotic protein-coding genes consist of exons and introns. Exon–intron borders are conserved between species and thus their changes might be observed only on quite long evolutionary distances. One of the rarest types of change, in which intron relocates over a short distance, is called "intron sliding", but the reality of this event has been debated for a long time. The main idea of a search for intron sliding is to use the most accurate genome annotation and genome sequence, as well as high-quality transcriptome data. We applied them in a search for sliding introns in mammals in order to widen knowledge about the presence or absence of such phenomena in this group. Results We didn’t find any significant evidence of intron sliding in the primate group (human, chimpanzee, rhesus macaque, crab-eating macaque, green monkey, marmoset). Only one possible intron sliding event supported by a set of high quality transcriptomes was observed between EIF1AX human and sheep gene orthologs. Also, we checked a list of previously observed intron sliding events in mammals and showed that most likely they are artifacts of genome annotations and are not shown in subsequent annotation versions as well as are not supported by transcriptomic data. Conclusions We assume that intron sliding is indeed a very rare evolutionary event if it exists at all. Every case of intron sliding needs a lot of supportive data for detection and confirmation.

BackgroundGerminomas (IG) account for up to 50% of all intracranial germ cell tumours. These tumours are reputed to be more prevalent in Oriental populations in comparison to Western cohorts. Biological characteristics of IG in other ethnic groups are unknown. Singapore is a multi-ethnic country with diverse cultures. Owing to inter-racial heterogeneity, the authors hypothesize there are molecular differences between paediatric IG patients in our local population. The aims of this study are exploratory: firstly, to identify molecular characteristics in this tumour type and circulating CSF unique to different racial cohorts; and next, to corroborate our findings with published literature.MethodsThis is a single-institution, retrospective study of prospectively collected data. Inclusion criteria encompass all paediatric patients with histologically confirmed IG. Excess CSF and brain tumour tissues are collected for molecular analysis. Tumour tissues are subjected to a next generation sequencing (NGS) targeted panel forKIT and PDGRA. All CSF samples are profiled via a high-throughput miRNA multiplexed workflow. Results are then corroborated with existing literature and public databases.ResultsIn our cohort of 14 patients, there areKITexon variants in the tumour tissues and CSF miRNAs corroborative with published studies. Separately, there are alsoKITexon variants and miRNAs not previously highlighted in IG. A subgroup analysis demonstrates differential CSF miRNAs between Chinese and Malay IG patients.ConclusionThis is the first in-depth molecular study of a mixed ethnic population of paediatric IGs from a Southeast Asian cohort. Validation studies are required to assess the relevance of novel findings in our study.

Background Central nervous system (CNS) hemangioblastomas are the most frequent cause of mortality in patients with Von Hippel–Lindau (VHL) disease, an autosomal dominant genetic disease resulting from germline mutations in the VHL tumor suppressor gene, with most mutations occurring in the exons. To date, there have been no reports of CNS hemangioblastoma cases related to pathogenic variants in intron 2 of VHL, which encodes a tumor suppressor protein (i.e., pVHL) that regulates hypoxia-inducible factor proteins. Case presentation We report the presence of a base substitution of c.464-1G > C and c.464-2A > G in the intron 2 of VHL causing CNS hemangioblastomas in six patients with VHL from two Chinese families. The clinical information about the two pathogentic variants has been submitted to ClinVar database. The ClinVar accession for NM_000551.3(VHL):c.464-1G > C was SCV001371687. This finding may provide a new approach for diagnosing and researching VHL-associated hemangioblastomas. Conclusions This is the first report of a pathogenic variant at intron 2 in VHL-associated hemangioblastomas. Gene sequencing showed that not only exonic but also intronic mutations can lead to the development of CNS hemangioblastomas.

Background The naked caryopsis character in barley is a domestication-associated trait defined by loss-of-function of the NUD gene. The functional NUD gene encodes an Apetala 2/Ethylene-Response Factor (AP2/ERF) controlling the formation of a cementing layer between pericarp and both lemma and palea. The downstream genes regulated by the NUD transcription factor and molecular mechanism of a cementing layer formation are still not sufficiently described. A naturally occurring 17-kb deletion in the nud locus is associated with the emergence of naked barley. Naked barley has been traditionally used for food and nowadays is considered as a dietary component for functional nutrition. Results In the present study, we demonstrate that targeted knockout of the NUD gene using RNA-guided Cas9 endonuclease leads to the phenotype conversion from hulled to naked barley. Using in vivo pre-testing systems, highly effective guide RNAs targeting the first exon of the NUD gene were selected. Expression cassettes harboring the cas9 and guide RNA genes were used to transform barley cv. Golden Promise via Agrobacterium-mediated DNA transfer. The recessive naked grain phenotype was observed in 57% of primary transformants, which indicates a frequent occurrence of homozygous or biallelic mutations. T-DNA-free homozygous lines with independently generated mutations in the NUD gene were obtained in the T1 generation. At homozygous state, all obtained mutations including one- and two-amino acid losses with the translational reading frame being retained invariably caused the naked grain phenotype. Conclusions The hulled and naked barley isogenic lines generated are a perfect experimental model for further studies on pleiotropic consequences of nud mutations on overall plant performance under particular consideration of yield-determining traits. Due to the high β-glucan content of its grains, naked barley is considered as being of particular dietary value. The possibility to convert hulled into naked barley cultivars by targeted mutagenesis allows breeders to extend the potential utilization of barley by the provision of functional food.

Background Our previous study found that chicken KLF7 was an important regulator in formation of adipose tissue. In the present study, we analyzed the association for DNA methylation in chicken KLF7 with its transcripts of abdominal adipose tissue and blood metabolic indicators. Results The KLF7 transcripts of the adipose tissue of Chinese yellow broilers were associated with age (F = 6.67, P = 0.0035). In addition, the KLF7 transcripts were negatively correlated with blood glucose levels (r = − 0.61841, P = 0.0140). The DNA methylation levels of 26 CpG loci in the chicken KLF7 promoter and Exon 2 were studied by Sequenom MassArray. A total of 22 valid datasets were obtained. None of them was significantly different in relation to age (P > 0.05). However, the DNA methylation levels in the promoter were lower than those in Exon 2 (T = 40.74, P KLF7 transcripts and blood high-density lipoprotein levels, respectively, and many CpG loci were correlated with each other (P KLF7 transcripts (P KLF7 transcripts was built, and the contribution of z1 to the variation on KLF7 transcripts was 34.29%. Conclusions In conclusion, the KLF7 transcripts of chicken abdominal adipose tissue might be inhibited by DNA methylation in the promoter, and it might be related to the DNA methylation level of PCpG6.

Background Stargardt disease (STGD1) is a common recessive hereditary macular dystrophy in early adulthood or childhood, with an estimated prevalence of 1:8000 to 1:10,000. ABCA4 is the causative gene for STGD1. The current study aims at identifying the novel disease-related ABCA4 variants in Han Chinese families with STGD1 using next-generation sequencing (NGS). Methods In the present study, 12 unrelated Han Chinese families (19 males and 17 females) with STGD1 were tested by panel-based NGS. In order to capture the coding exons and the untranslated regions (UTRs) plus 30 bp of intronic flanking sequences of 792 genes, which were closely associated with usual ophthalmic genetic disease, we designed a customized panel, namely, Target_Eye_792_V2 chip. STGD1 patients were clinically diagnosed by experienced ophthalmologists. All the detected variants were filtered and analyzed through the public databases and in silico programs to assess potential pathogenicity. Results Twenty-one ABCA4 mutant variants were detected in 12 unrelated Han Chinese families with STGD1, containing 14 missense, three splicing, two frameshift, one small deletion, and one nonsense variants. Base on the American College of Medical Genetics (ACMG) guidelines, 8 likely pathogenic and 13 pathogenic variants were determined. The functional consequences of these mutant variants were predicted through in silico programs. Of the 21 mutant variants in ABCA4, two novel coding variants c.3017G > A and c.5167 T > C and one novel null variant c.3051-1G > A were detected in three unrelated probands. Conclusions By panel-based NGS, 21 ABCA4 variants were confirmed in 12 unrelated Han Chinese families. Among them, 3 novel mutant variants were found, which further expanded the ABCA4 mutation spectrum in STGD1 patients.

Background Mandibulofacial dysostosis with microcephaly (MFDM) is a rare autosomal dominant genetic disease characterized by intellectual and growth retardations, as well as major microcephaly, induced by missense and splice site variants or microdeletions in the EFTUD2 gene. Case presentation Here, we investigate the case of a young girl with symptoms of MFDM and a normal karyotype. Whole-exome sequencing of the family was performed to identify genetic alterations responsible for this phenotype. We identified a de novo synonymous variant in the EFTUD2 gene. We demonstrated that this synonymous variant disrupts the donor splice-site in intron 9 resulting in the skipping of exon 9 and a frameshift that leads to a premature stop codon. Conclusions We present the first case of MFDM caused by a synonymous variant disrupting the donor splice site, leading to exon skipping.

Background Deregulated circular RNAs (circRNAs) are associated with the development of cancer and therapy resistance. However, functional research of circRNAs mostly focus on potential miRNA or protein binding and more potential regulation of circRNA on host gene DNA in cancers are yet to be inspected. Method We performed total RNA sequencing on clinical breast cancer samples and identified the expression patterns of circRNAs and corresponding host genes in patient blood, tumor and adjacent normal tissues. qPCR, northern blot and in situ hybridization were used to validate the dysregulation of circRNA circSMARCA5. A series of procedures including R-loop dot-blotting, DNA-RNA immunoprecipitation and mass spectrum, etc. were conducted to explore the regulation of circSMARCA5 on the transcription of exon 15 of SMARCA5. Moreover, immunofluorescence and in vivo experiments were executed to investigate the overexpression of circSMARCA5 with drug sensitivities. Results We found that circRNAs has average higher expression over its host linear genes in peripheral blood. Compared to adjacent normal tissues, circSMARCA5 is decreased in breast cancer tissues, contrary to host gene SMARCA5. The enforced expression of circSMARCA5 induced drug sensitivity of breast cancer cell lines in vitro and in vivo. Furthermore, we demonstrated that circSMARCA5 can bind to its parent gene locus, forming an R-loop, which results in transcriptional pausing at exon 15 of SMARCA5. CircSMARCA5 expression resulted in the downregulation of SMARCA5 and the production of a truncated nonfunctional protein, and the overexpression of circSMARCA5 was sufficient to improve sensitivity to cytotoxic drugs. Conclusion Our results revealed a new regulatory mechanism for circRNA on its host gene and provided evidence that circSMARCA5 may serve as a therapeutic target for drug-resistant breast cancer patients.

BackgroundDe novo RNA-Seq assembly is a powerful method for analysing transcriptomes when the reference genome is not available or poorly annotated. However, due to the short length of Illumina reads it is usually impossible to reconstruct complete sequences of complex genes and alternative isoforms. Recently emerged possibility to generate long RNA reads, such as PacBio and Oxford Nanopores, may dramatically improve the assembly quality, and thus the consecutive analysis. While reference-based tools for analysing long RNA reads were recently developed, there is no established pipeline for de novo assembly of such data.ResultsIn this work we present a novel method that allows to perform high-quality de novo transcriptome assemblies by combining accuracy and reliability of short reads with exon structure information carried out from long error-prone reads. The algorithm is designed by incorporating existing hybridSPAdes approach into rnaSPAdes pipeline and adapting it for transcriptomic data.ConclusionTo evaluate the benefit of using long RNA reads we selected several datasets containing both Illumina and Iso-seq or Oxford Nanopore Technologies (ONT) reads. Using an existing quality assessment software, we show that hybrid assemblies performed with rnaSPAdes contain more full-length genes and alternative isoforms comparing to the case when only short-read data is used.Availability and implementationrnaSPAdes is implemented in C++ and Python and is freely available for Linux and MacOS under GPLv2 license at cab.spbu.ru/software/rnaspades/ and github.com/ablab/spades.

Background Holocarboxylase synthetase (HLCS) deficiency is a rare inborn disorder of biotin metabolism, which results in defects in several biotin-dependent carboxylases and presents with metabolic ketoacidosis and skin lesions. Case presentation In this paper, we report a Chinese Han pedigree with HLCS deficiency diagnosed by using next-generation sequencing and validated with Sanger sequencing of the HLCS and BTD genes. The Chinese proband carries the common missense mutation c.1522C > T (p.Arg508Trp) in exon 9 of the HLCS gene, which generates an increased Km value for biotin. A novel frameshift mutation c.1006_1007delGA (p.Glu336Thrfs*15) in exon 6 of the HLCS gene is predicted to be deleterious through PROVEAN and MutationTaster. A novel heterozygous mutation, c.638_642delAACAC (p.His213Profs*4), in the BTD gene is also identified. Conclusions The Chinese proband carries the reported Arg508Trp variant, the novel 2-bp frameshift mutation c.1006_1007delGA (p.Glu336Thrfs*15), which expands the mutational spectrum of the HLCS gene, and the novel heterozygous mutation c.638_642delAACAC (p.His213Profs*4), which expands the mutational spectrum of the BTD gene. Furthermore, reversible hearing damage is rarely reported in patients with HLCS deficiency, which deserves further discussion.

Background The major histocompatibility complex (MHC) is responsible for detecting and addressing foreign pathogens inside the body. While the general structure of MHC genes is relatively well conserved among mammalian species, it is notably different among ruminants due to a chromosomal inversion that splits MHC type II genes into two subregions (IIa, IIb). Recombination rates are reportedly high between these subregions, and a lack of linkage has been documented in domestic ruminants. However, no study has yet examined the degree of linkage between these subregions in a wild ruminant. The white-tailed deer (Odocoileus virginianus), a popular ruminant of the Cervidae family, is habitually plagued by pathogens in its natural environment (e.g. Haemonchus contortus, Elaeophora). Due to the association between MHC haplotypes and disease susceptibility, a deeper understanding of MHC polymorphism and linkage between MHC genes can further aid in this species’ successful management. We sequenced MHC-DRB exon 2 (IIa) and MHC-DOB exon 2 (IIb) on the MiSeq platform from an enclosed white-tailed deer population located in Alabama. Results We identified 12 new MHC-DRB alleles, and resampled 7 alleles, which along with other published alleles brings the total number of documented alleles in white-tailed deer to 30 for MHC-DRB exon 2. The first examination of MHC-DOB in white-tailed deer found significantly less polymorphism (11 alleles), as was expected of a non-classical MHC gene. While MHC-DRB was found to be under positive, diversifying selection, MHC-DOB was found to be under purifying selection for white-tailed deer. We found no significant linkage disequilibrium between MHC-DRB and MHC-DOB, suggesting that these loci are unlikely to be closely linked. Conclusions Overall, this study identified 12 new MHC-DRB exon 2 alleles and characterized a new, non-classical, MHC II gene (MHC-DOB) for white-tailed deer. We also found a lack of significant linkage between these two loci, which supports previous findings of a chromosomal inversion within the MHC type II gene region in ruminants, and suggests that white-tailed deer may have a recombination hotspot between these MHC regions similar to that found for Bos taurus.

Background The disease gene of fragile X syndrome, FMR1 gene, encodes fragile X mental retardation protein (FMRP). The alternative splicing (AS) of FMR1 can affect the structure and function of FMRP. However, the biological functions of alternatively spliced isoforms remain elusive. In a previous study, we identified a new 140bp exon from the intron 9 of human FMR1 gene. In this study, we further examined the biological functions of this new exon and its underlying signaling pathways. Results qRT-PCR results showed that this novel exon is commonly expressed in the peripheral blood of normal individuals. Comparative genomics showed that sequences paralogous to the 140 bp sequence only exist in the genomes of primates. To explore the biological functions of the new transcript, we constructed recombinant eukaryotic expression vectors and lentiviral overexpression vectors. Results showed that the spliced transcript encoded a truncated protein which was expressed mainly in the cell nucleus. Additionally, several genes, including the BEX1 gene involved in mGluR-LTP or mGluR-LTD signaling pathways were significantly influenced when the truncated FMRP was overexpressed. Conclusions our work identified a new exon from amid intron 9 of human FMR1 gene with wide expression in normal healthy individuals, which emphasizes the notion that the AS of FMR1 gene is complex and may in a large part account for the multiple functions of FMRP.

Hepatitis B virus (HBV), along with Hepatitis C virus chronic infection, represents a major risk factor for hepatocellular carcinoma (HCC) development. However, molecular mechanisms involved in the development of HCC are not yet completely understood. Recent studies have indicated that mutations in CTNNB1 gene encoding for β-catenin protein lead to aberrant activation of the Wnt/ β-catenin pathway. The mutations in turn activate several downstream genes, including c-Myc, promoting the neoplastic process. The present study evaluated the mutational profile of the CTNNB1 gene and expression levels of CTNNB1 and c-Myc genes in HBV-related HCC, as well as in cirrhotic and control tissues. Mutational analysis of the β-catenin gene and HBV genotyping were conducted by direct sequencing. Expression of β-catenin and c-Myc genes was assessed using real-time PCR. Among the HCC cases, 18.1% showed missense point mutation in exon 3 of CTNNB1, more frequently in codons 32, 33, 38 and 45. The frequency of mutation in the hotspots of exon 3 was significantly higher in non-viral HCCs (29.4%) rather than HBV-related cases (12.7%, P = 0.021). The expression of β-catenin and c-Myc genes was found upregulated in cirrhotic tissues in association with HBV infection. Mutations at both phosphorylation and neighboring sites were associated with increased activity of the Wnt pathway. The results demonstrated that mutated β-catenin caused activation of the Wnt pathway, but the rate of CTNNB1 gene mutations was not related to HBV infection. HBV factors may deregulate the Wnt pathway by causing epigenetic alterations in the HBV-related HCC.

Molecular biomarkers provide both diagnostic and prognostic results for patients with diffuse glioma, the most common primary brain tumor in adults. Here, we used a long-read nanopore-based sequencing technique to simultaneously assess IDH mutation status and MGMT methylation level in 4 human cell lines and 8 fresh human brain tumor biopsies. Currently, these biomarkers are assayed separately, and results can take days to weeks. We demonstrated the use of nanopore Cas9-targeted sequencing (nCATS) to identify IDH1 and IDH2 mutations within 36 h and compared this approach against currently used clinical methods. nCATS was also able to simultaneously provide high-resolution evaluation of MGMT methylation levels not only at the promoter region, as with currently used methods, but also at CpGs across the proximal promoter region, the entirety of exon 1, and a portion of intron 1. We compared the methylation levels of all CpGs to MGMT expression in all cell lines and tumors and observed a positive correlation between intron 1 methylation and MGMT expression. Finally, we identified single nucleotide variants in 3 target loci. This pilot study demonstrates the feasibility of using nCATS as a clinical tool for cancer precision medicine.

Background Primary myoepithelial carcinoma of the lung is a rare subtype in lung cancer. Comprehensive molecular profiling of myoepithelial carcinoma of the lung is absent, neither was clinical evidence of targeted therapy available for this disease. Therefore, the optimal treatment regimen of this tumor needs to be established. Case presentation Here we present a case of a 68-year-old patient with stage IVB primary myoepithelial carcinoma of the lung who harbored EGFR exon 19 deletion and KRAS mutation and underwent icotinib targeted therapy, achieving partial response (PR) with progression free survival (PFS) of 3 months. Conclusion To our knowledge, this study describes the first documented case of primary myoepithelial carcinoma lung cancer patient harboring EGFR exon 19 deletion and KRAS mutation, and showed clinical efficacy of epidermal growth factor receptor tyrosine kinase inhibitors (EGFR-TKI) treatment in this patient.