Intron Sequences Research Articles

The functions, oncogenic roles, and clinical significance of circular RNAs in renal cell carcinoma.

Renal cell carcinoma (RCC) is the most common form of malignancy affecting the kidneys. Circular RNAs (circRNAs) are non-coding RNAs that are derived from exonic or intronic sequences through a selective shearing process. There is growing evidence that these circRNAs can influence a range of biological pathways by serving as protein decoys, microRNA sponges, regulators of transcriptional activity, or templates for protein translation. The dysregulation of circRNA expression patterns is a hallmark of RCC and other cancer types, and there is strong evidence that these RNA species can play central roles in the onset and progression of RCC tumors. In the present review, we summarized recent findings on the functional roles and clinical impacts of circRNAs in RCC. Further, we discussed their potential utility as diagnostic biomarkers or targets for therapeutic intervention.

Characterization and genomic analyses of a novel alphabaculovirus isolated from the black armyworm, Spodoptera cosmioides (Lepidoptera: Noctuidae)

The black armyworm Spodoptera cosmioides is a pest of increasing importance in Cry1Ac-Bt toxin crops and non-Bt crops of soybean and cotton in Brazil. Here we characterized a baculovirus isolated from extracts of S. cosmioides that died with symptoms of nuclear polyhedrosis. The putative novel virus exhibited polyhedral occlusion bodies (OBs) with virions containing multiple rod-shaped nucleocapsids, characteristic of alphabaculoviruses. The virus isolate was named Spodoptera cosmioides nucleopolyhedrovirus isolate CNPSo-72 (SpcoNPV-CNPSo-72). SpcoNPV-CNPSo-72 was lethal to third-instar S. cosmioides caterpillars but not to S. frugiperda under the tested viral concentrations. Moreover, SpcoNPV-CNPSo-72 contained a circular 147,763 bp long genome and a G + C content of 44.8% with 151 annotated ORFs (10 unique for baculovirus) and five homologous regions (hrs). The 38 currently defined baculovirus core genes were found in the SpcoNPV-CNPSo-72 genome. After phylogenetic analysis, the novel virus was found to be closely related to other members of Alphabaculovirus, especially to the Spodoptera-infecting viruses, which included Spodoptera eridania nucleopolyhedrovirus isolate 251, Spodoptera litura nucleopolyhedrovirus isolate II, Spodoptera exigua multiple nucleopolyhedrovirus isolate US-1, Spodoptera eridania nucleopolyhedrovirus isolate CNPSo-165, and Spodoptera frugiperda multiple nucleopolyhedrovirus isolate 19. Surprisingly, the new baculoviral genome was found to code for a putative arginine-associated tRNA gene with a predicted intronic sequence of 105 nt. The gene was found inside the bjdp CDS. Overall, baculoviruses are pathogens that lethally infect insect larvae and their study allows a better understanding of large DNA virus evolution, which provides important insights for the development and improvement of biological control agents.

RNA analysis of the GALNS transcript reveals novel pathogenic mechanisms associated with Morquio syndrome A

Morquio syndrome A (Mucopolysaccharidosis IVA, MPS IVA) is an autosomal recessive lysosomal storage disorder caused by deficiency of N-acetyl-galactosamine-6-sulfatase (GALNS) which catabolizes the glycosaminoglycans (GAG), keratan sulfate and chondroitin-6-sulfate. Homozygous or compound heterozygous pathogenic variants in the GALNS result in the deficiency of the enzyme and consequent GAG accumulations. DNA sequence and copy number analysis of the GALNS coding region fails to identify biallelic causative pathogenic variants in up to 15% of patients with Morquio syndrome A. RNA transcript analysis was performed to identify pathogenic alterations in two unrelated families with Morquio syndrome A in whom a single heterozygous or no pathogenic alteration was detected by standard analysis of the GALNS gene. RNA sequencing and quantitative expression analysis identified the overabundance of an aberrant GALNS transcript isoform and a reduction of the clinically relevant isoform (NM_000512.4) in the Morquio syndrome A patients from both families. The aberrant isoform (ENST00000568613.1) was produced by alternative splicing and contained intronic sequence that was likely a cryptic exon predicted to result in a reading frame shift and generation of a premature termination codon. These findings indicated that the aberrant splicing is likely the novel molecular defect in our patients. RNA transcript analysis could be useful to identify pathogenic alterations and increase the yield of molecular diagnosis in patients with Morquio syndrome A whose genetic variants are not found by standard sequencing or gene dosage analysis.

Deletion of Sf3b4 in Neural Crest Cells Encapsulates the Variable Expressivity of Craniofacial Abnormalities Observed in Nager and Rodriguez Syndrome Patients

Recently, both Nager syndrome and Rodriguez syndrome patients were identified to carry rare mutations in SF3B4. Both disorders affect the face as well as the limb including both hands and feet. Nager syndrome patients survive and exhibit smaller jaw bones, smaller cheek bones, cleft palate, hearing problems as well as downward slanted eyelids. However, Rodriguez syndrome is much more severe, and patients often die before or soon after birth. We hypothesized that craniofacial abnormalities associated with reduced SF3B4 levels is due to tissue‐specific expression and requirement of this gene in neural crest cells. To generate a conditional mutant mouse line for Sf3b4, we used CRISPR/Cas9 to insert LoxP sequences in intron 1 and 3 of this gene. To test our hypothesis, we mated Sf3b4conditional mutant mice to Wnt1‐Cre2 transgenic mice to delete exons 2 and 3 specifically in neural crest cells. Heterozygous mutant embryos from these matings were normal. Sf3b4L/L;Wnt1‐Cre2 embryos have brain and craniofacial abnormalities, including abnormal forebrain, midbrain and hindbrain, hypoplasia of the maxillomandibular region and midline cleft. Majority of these embryos die before birth and the few mutants that are born are indistinguishable from their wildtype littermates. Additionally, we generated Sf3b4L/‐;Wnt1‐Cre2 mutants with homozygous deletion of Sf3b4 in neural crest cells and heterozygous deletion of Sf3b4 in the body. These mutant embryos also had unfused frontonasal prominence and unfused and severely hypoplastic mandibular region. Sf3b4L/‐;Wnt1‐Cre2begin to die at E13.5 and none survive to birth. Thus, neural crest specific deletion of Sf3b4 in mice is able to model the craniofacial abnormalities presented in patients as well as capture the variable expressivity observed. Our model can further be used to uncover the splicing and molecular events regulated by SF3B4 during craniofacial development.

Promoter Proximity Defines Mutation Window for VH and VΚ Genes Rearranged to Different J Genes.

Somatic hypermutation induced by activation-induced deaminase (AID) occurs at high densities between the Ig V gene promoter and intronic enhancer, which encompasses DNA encoding the rearranged V gene exon and J intron. It has been proposed that proximity between the promoter and enhancer defines the boundaries of mutation in V regions. However, depending on the J gene used, the distance between the promoter and enhancer is quite variable and may result in differential targeting around the V gene. To examine the effect of distance in mutation accumulation, we sequenced 320 clones containing different endogenous rearranged V genes in the IgH and Igκ loci from Peyer's patch B cells of mice. Clones were grouped by their use of different J genes. Distances between the V gene and enhancer ranged from ∼2.3 kb of intron DNA for rearrangements using J1, ∼2.0 kb for rearrangements using J2, ∼1.6 kb for rearrangements using J3 (H) or 4 (κ), and 1.1 kb for rearrangements using J4 (H) or 5 (κ). Strikingly, >90% of intron mutations occurred within 1 kb downstream of the J gene for both H and κ clones, regardless of which J gene was used. Thus, there is no evidence that the intron sequence or enhancer plays a role in determining the extent of mutation. The results indicate that V region intron mutations are targeted by their proximity to the promoter, suggesting they result from AID interactions with RNA polymerase II over a 1-kb region.

From the comparative study of a circRNA originating from an mammalian ATXN2L intron to understanding the genesis of intron lariat-derived circRNAs

Circular intronic RNAs (ciRNAs) are still unexplored regarding mechanisms for their emergence. We considered the ATXN2L intron lariat-derived circular RNA (ciRNA-ATXN2L) as an opportunity to conduct a cross-species examination of ciRNA genesis. To this end, we investigated 207 datasets from 4 tissues and from 13 mammalian species. While in eight species, ciRNA-ATXN2L was never detected, in pigs and rabbits, ciRNA-ATXN2L was expressed in all tissues and sometimes at very high levels. Bovine tissues were an intermediate case and in macaques and cats, only ciRNA-ATXN2L traces were detected. The pattern of ciRNA-ATXN2L restricted to only five species is not related to a particular evolution of intronic sequences. To empower our analysis, we considered 221 additional introns including 80 introns where a lariat-derived ciRNA was previously described. The primary driver of micro-ciRNA genesis (< 155 nt as ciRNA-ATXN2L) appears to be the absence of a canonical “A” (i.e. a "tnA" located in the usual branching region) to build the lariat around this adenosine. The balance between available "non canonical-A" (no ciRNA genesis) and "non-A" (ciRNA genesis) for use as a branch point to build the lariat could modify the expression level of ciRNA-ATXN2L. In addition, the rare localization of the 2′-5′ bond in an open RNA secondary structure could also negatively affect the lifetime of ciRNAs (macaque ciRNA-ATXN2L). Our analyses suggest that ciRNA-ATXN2L is likely a functionless splice remnant. This study provides a better understanding of the ciRNAs origin, especially drivers for micro ciRNA genesis.

Comparative Study of Starch Phosphorylase Genes and Encoded Proteins in Various Monocots and Dicots with Emphasis on Maize.

Starch phosphorylase (PHO) is a multimeric enzyme with two distinct isoforms: plastidial starch phosphorylase (PHO1) and cytosolic starch phosphorylase (PHO2). PHO1 specifically resides in the plastid, while PHO2 is found in the cytosol. Both play a critical role in the synthesis and degradation of starch. This study aimed to report the detailed structure, function, and evolution of genes encoding PHO1 and PHO2 and their protein ligand-binding sites in eight monocots and four dicots. “True” orthologs of PHO1 and PHO2 of Oryza sativa were identified, and the structure of the enzyme at the protein level was studied. The genes controlling PHO2 were found to be more conserved than those controlling PHO1; the variations were mainly due to the variable sequence and length of introns. Cis-regulatory elements in the promoter region of both genes were identified, and the expression pattern was analyzed. The real-time quantitative polymerase chain reaction indicated that PHO2 was expressed in all tissues with a uniform pattern of transcripts, and the expression pattern of PHO1 indicates that it probably contributes to the starch biosynthesis during seed development in Zea mays. Under abscisic acid (ABA) treatment, PHO1 was found to be downregulated in Arabidopsis and Hordeum vulgare. However, we found that ABA could up-regulate the expression of both PHO1 and PHO2 within 12 h in Zea mays. In all monocots and dicots, the 3D structures were highly similar, and the ligand-binding sites were common yet fluctuating in the position of aa residues.

Case Review: Whole-Exome Sequencing Analyses Identify Carriers of a Known Likely Pathogenic Intronic BRCA1 Variant in Ovarian Cancer Cases Clinically Negative for Pathogenic BRCA1 and BRCA2 Variants.

Background: Detecting pathogenic intronic variants resulting in aberrant splicing remains a challenge in routine genetic testing. We describe germline whole-exome sequencing (WES) analyses and apply in silico predictive tools of familial ovarian cancer (OC) cases reported clinically negative for pathogenic BRCA1 and BRCA2 variants. Methods: WES data from 27 familial OC cases reported clinically negative for pathogenic BRCA1 and BRCA2 variants and 53 sporadic early-onset OC cases were analyzed for pathogenic variants in BRCA1 or BRCA2. WES data from carriers of pathogenic BRCA1 or BRCA2 variants were analyzed for pathogenic variants in 10 other OC predisposing genes. Loss of heterozygosity analysis was performed on tumor DNA from variant carriers. Results: BRCA1 c.5407-25T>A intronic variant, identified in two affected sisters and one sporadic OC case, is predicted to create a new splice effecting transcription of BRCA1. WES data from BRCA1 c.5407-25T>A carriers showed no evidence of pathogenic variants in other OC predisposing genes. Sequencing the tumor DNA from the variant carrier showed complete loss of the wild-type allele. Conclusions: The findings support BRCA1 c.5407-25T>A as a likely pathogenic variant and highlight the importance of investigating intronic sequences as causal variants in OC families where the involvement of BRCA1 is highly suggestive.

Expansion of Cancer Risk Profile for BRCA1 and BRCA2 Pathogenic Variants

The clinical importance of genetic testing of BRCA1 and BRCA2 in breast, ovarian, prostate, and pancreatic cancers is widely recognized. However, there is insufficient evidence to include other cancer types that are potentially associated with BRCA1 and BRCA2 in clinical management guidelines. To evaluate the association of BRCA1 and BRCA2 pathogenic variants with additional cancer types and their clinical characteristics in 100 914 individuals across 14 cancer types. This case-control analysis to identify cancer types and clinical characteristics associated with pathogenic variants in BRCA1 and BRCA2 included DNA samples and clinical information from 63 828 patients with 14 common cancer types and 37 086 controls that were sourced from a multi-institutional hospital-based registry, BioBank Japan, between April 2003 and March 2018. The data were analyzed between August 2019 and October 2021. Germline pathogenic variants in coding regions and 2 bp flanking intronic sequences in BRCA1 and BRCA2 were identified by a multiplex polymerase chain reaction-based target sequence method. Associations of (likely) pathogenic variants with each cancer type were assessed by comparing pathogenic variant carrier frequency between patients in each cancer type and controls. A total of 65 108 patients (mean [SD] age at diagnosis, 64.1 [11.6] years; 27 531 [42.3%] female) and 38 153 controls (mean [SD] age at registration, 61.8 [14.6] years; 17 911 [46.9%] female) were included in this study. A total of 315 unique pathogenic variants were identified. Pathogenic variants were associated with P < 1 × 10-4 with an odds ratio (OR) of greater than 4.0 in biliary tract cancer (OR, 17.4; 95% CI, 5.8-51.9) in BRCA1, esophageal cancer (OR, 5.6; 95% CI, 2.9-11.0) in BRCA2, and gastric cancer (OR, 5.2; 95% CI, 2.6-10.5) in BRCA1, and (OR, 4.7; 95% CI, 3.1-7.1) in BRCA2 in addition to the 4 established cancer types. We also observed an association with 2 and 4 other cancer types in BRCA1 and BRCA2, respectively. Biliary tract, female breast, ovarian, and prostate cancers showed enrichment of carrier patients according to the increased number of reported cancer types in relatives. The results of this large-scale registry-based case-control study suggest that pathogenic variants in BRCA1 and BRCA2 were associated with the risk of 7 cancer types. These results indicate broader clinical relevance of BRCA1 and BRCA2 genetic testing.

Using iRNA-seq analysis to predict gene expression regulatory level and activity in Zea mays tissues.

Plants regulate gene expression at the transcriptional and post-transcriptional levels to produce a variety of functionally diverse cells and tissues that ensure normal growth, development, and environmental response. Although distinct gene expression patterns have been characterized between different plant tissues, the specific role of transcriptional regulation of tissue-specific expression is not well-characterized in plants. RNA-seq, while widely used to assay for changes in transcript abundance, does not discriminate between differential expression caused by mRNA degradation and active transcription. Recently, the presence of intron sequences in RNA-seq analysis of libraries constructed with total RNA has been found to coincide with genes undergoing active transcription. We have adapted the intron RNA-sequencing analysis to determine genome-wide transcriptional activity in 2 different maize (Zea mays) tissues: husk and V2-inner stem tissue. A total of 5,341 genes were predicted to be transcriptionally differentially expressed between the 2 tissues, including many genes expected to have biological activity relevant to the functional and developmental identity of each tissue. Correlations with transcriptional enhancer and transcription factor activity support the validity of intron RNA-sequencing predictions of transcriptional regulation. A subset of transcription factors was further analyzed using gene regulatory network analysis to determine the possible impact of their activation. The predicted regulatory patterns between these genes were used to model a potential gene regulatory network of transcription factors and regulatory targets.

CircRNA perspective: new strategies for RNA therapy.

Circular RNA (circRNA) is generated via intronic complementary sequence (ICS)-driven, RNA binding protein (RBP)-driven, or lariat-mediated back-splicing (BS) biogenesis models. The physiological and pathological functions of circRNAs occur through miRNA sponging, interacting with circRNA binding proteins (cRBPs), protein translation, or transcriptional regulation. CircRNAs play roles in regulation of cancer development, cardiovascular diseases, neurogenesis and neurological disorders, and metabolic diseases. CircRNA perspective: new strategies for RNA therapy: (Trends in Molecular Medicine , 343–344; 2022)Wu et al.Trends in Molecular MedicineMay 1, 2022In BriefIn the published version of Figure 2, the term 'circRNA' is mispresented as 'ciRNA' due to an editing error. This has been corrected in the revised figure. We apologize for any confusion this caused. Full-Text PDF

Identification of a Novel Deep Intronic Variant by Whole Genome Sequencing Combined With RNA Sequencing in a Chinese Patient With Menkes Disease

Background: Menkes disease (MD) is a rare X-linked connective tissue disorder of copper metabolism caused by pathogenic variant(s) in ATP7A gene. The aim of the present study is to determine the clinical characteristics and molecular basis of one patient with MD. Methods: One 10-month-old Chinese boy who met the clinical manifestations of MD was enrolled in this study. Whole genome sequencing (WGS) was performed in the patient in order to identify the variant(s), followed by Sanger sequencing. RNA sequencing (RNA-seq) from whole blood was subsequently applied to assess the effect of variant on transcription levels, and reverse transcriptase-polymerase chain reaction (RT-PCR) was performed for further validation. In addition, X chromosome inactivation (XCI) status of the patient’s mother at the DNA level was measured by capillary electrophoresis. Results: The patient suffered from intermittent convulsions for more than 6 months, with psychomoto retardation and neurodegenerations. The patient also had curly hair, hypopigmented skin, cutis laxa, decreased muscle strength and hypotonia. MRI showed the intracranial arteries were tortuous with some “spiral” changes. The patient’s serum ceruloplasmin level was low. WGS revealed one novel hemizygous variant, c.2627-501C > T (NM_000,052.7), located in the deep intronic sequence of ATP7A gene. Sanger sequencing confirmed that the variant was inherited from his mother. RNA-seq confirmed the variant itself, and identified a pseudo-exon inserted between exons 12 and 13 in mRNA of ATP7A. The sequencing results of RT-PCR from the patient confirmed this finding, while neither of his parents detected aberrant splicing. The Capillary electrophoresis results showed that the patient’s mother had a skewed XCI. Conclusion: Our finding of the variant enlarges the variant spectrum in the ATP7A gene. This is a novel deep intronic variant which leads to the activation of a pseudo-exons in the ATP7A gene, and it demonstrates the usefulness of WGS combined with RNA-seq, in terms of revealing disease-causing variants in non-coding regions. Furthermore, the fact that the deep intronic variants cause disease by the activation of pseudo-exon inclusion indicates that in MD this might be an important mechanism.

Methylome Analysis in Nonfunctioning and GH-Secreting Pituitary Adenomas.

Pituitary adenomas (PAs), usually benign lesions, can sometimes present with “aggressive” features (rapid growth, local invasiveness, scarce response to conventional treatments). Despite the fact that a few genetic alterations have been associated to this clinical behavior, the role of epigenetic modifications, mainly methylation and miRNAs activity, is now opening new frontiers in this field. We evaluated the methylation profile of 21 PA (11 GH-omas, 10 nonfunctioning tumors—NFPAs) samples from TNS surgery and 5 normal pituitaries, collected at our neurosurgery between 2015 and 2017. DNA was extracted and sequenced, selecting 184,841 target regions. Moreover, methylation profiles were correlated with demographic, radiological, and clinicopathological features. NFPAs showed higher methylation levels vs. GH-omas, with 178 differentially methylated regions (DMRs) mainly consisting of noncoding and intronic sequences, and mostly localized in the open sea regions. We also found three hypermethylated genes (C7orf50, GNG7, and BAHCC1) involved in tumorigenesis processes and potentially influencing pituitary tumor pathophysiology. Among the clinicopathological features, only the maximum diameter resulted significantly higher in NFPAs. Our data provide further evidence of the complex epigenetic background of pituitary tumors. In line with the current literature, we confirmed a significant prevalence of hypermethylation in NFPAs vs. GH-omas, whose pathophysiological consequence is yet to be defined.

Genetic diversity and phylogenetic analyses of Vitex agnus castus L. populations using ISSR-PCR and chloroplast DNA trnL intron and trnL-F sequences

In this study, genetic diversity and phylogenetic analysis of some Vitex agnus-castus L. populations were conducted based on ISSR-PCR and chloroplast DNA trnL intron, and trnL-F sequence analyses. Vitex agnus-castus populations were detected in Aydın province, Turkey. Fresh leaf samples from the populations were collected and brought to the laboratory for genomic DNA isolation. 15 ISSR primers were used to determine the genetic diversity ofe Vitex agnus-castus populations. A total of 138 bands were obtained in ISSR analysis, 85 of which were polymorphic and 53 were monomorphic. Polymorphism rate was determined as 61.59 %. trnC, trnD, trnE, and trnF primers were used for PCR amplification of the chloroplast trnL intron and trnL-F region. A total of 138 bands were obtained by ISSR analysis. For trnL intron analyses, nucleotide lengths of 13 populations were between 508 and 516. The average nucleotide composition consisted of 38.5 % T, 18.3 % C, 27.5 % A and 15.7 % G. In trnL-F assays, the nucleotide lengths of the 13 populations ranged from 330 to 353. The average nucleotide composition consisted of 29.4 % T, 18.1 % C, 32.9 % A and 19.6 % G. The results of the phylogenetic trees constructed using some trnL intron and trnL-F sequences of Vitex doniana, Vitex trifolia, Vitex triflora, Vitex turczaninowii, Vitex queenslandica, Vitex axillariflora, Vitex rotundifolia and Vitex negundo species obtained from NCBI were compared. As a result of the study, polymorphisms were obtained at a rate of 61.59% from the ISSR analysis. In addition, the phylogenetic relationship between chloroplast trnL intron and trnL-F sequences of Vitex agnus-castus populations along with the other species was revealed.

Gene loss, genome rearrangement, and accelerated substitution rates in plastid genome of Hypericum ascyron (Hypericaceae)

BackgroundComparative genomic analysis exhibits dynamic evolution of plastid genome (plastome) in the clusioid clade of Malpighiales, which comprise five families, including multiple inversions and gene losses. Little is known about the plastome evolution in Hypericaceae, a large family in the clade. Only the plastome of one species, Cratoxylum cochinchinense, has been published.ResultsWe generated a complete plastome sequence for Hypericum ascyron, providing the first complete plastome from the tribe Hypericeae (Hypericaceae). The H. ascyron plastome exhibits dynamic changes in gene and intron content, structure, and sequence divergence compared to the C. cochinchinense plastome from the tribe Cratoxyleae (Hypericaceae). Transcriptome data determined the evolutionary fate of the missing plastid genes infA, rps7, rps16, rpl23, and rpl32 in H. ascyron. Putative functional transfers of infA, rps7, and rpl32 were detected to the nucleus, whereas rps16 and rpl23 were substituted by nuclear-encoded homologs. The plastid rpl32 was integrated into the nuclear-encoded SODcp gene. Our findings suggested that the transferred rpl32 had undergone subfunctionalization by duplication rather than alternative splicing. The H. ascyron plastome rearrangements involved seven inversions, at least three inverted repeat (IR) boundary shifts, which generated gene relocations and duplications. Accelerated substitution rates of plastid genes were observed in the H. ascyron plastome compared with that of C. cochinchinense plastid genes. The higher substitution rates in the accD and clpP were correlated with structural change, including a large insertion of amino acids and losses of two introns, respectively. In addition, we found evidence of positive selection of the clpP, matK, and rps3 genes in the three branches related to H. ascyron. In particular, the matK gene was repeatedly under selection within the family Hypericaceae. Selective pressure in the H. ascyron matK gene was associated with the loss of trnK-UUU and relocation into the IR region.ConclusionsThe Hypericum ascyron plastome sequence provides valuable information for improving the understanding of plastome evolution among the clusioid of the Malpighiales. Evidence for intracellular gene transfer from the plastid to the nucleus was detected in the nuclear transcriptome, providing insight into the evolutionary fate of plastid genes in Hypericaceae.

CRISPR-Mediated Enzyme Fragment Complementation Assay for Quantification of the Stability of Splice Isoforms.

Small-molecule splicing modulators exemplified by an FDA-approved drug, risdiplam, are a new pharmacological modality for regulating the expression and stability of splice isoforms. We report a CRISPR-mediated enzyme fragment complementation (EFC) assay to quantify the splice isoform stability. The EFC assay harnessed a 42 amino acid split of a β-galactosidase (designate α-tag), which could be fused at the termini of the target genes using CRISPR/cas9. The α-tagged splice isoform would be quantified by measuring the enzymatic activity upon complementation with the rest of β-galactosidase. This EFC assay retained all the sequences of introns and exons of the target gene in the native genomic environment that recapitulates the cell biology of the diseases of interest. For a proof-of-concept, we developed a CRISPR-mediated EFC assay targeting the exon 7 of the survival of motor neuron 2 (SMN2) gene. The EFC assay is compatible with 384-well plates and robustly quantified the splicing modulation activity of small molecules. In this study, we also discovered that a coumarin derivative, compound 4, potently modulated SMN2 exon 7 splicing at as low as 1.1 nM.

The evolution, evolvability and engineering of gene regulatory DNA.

Mutations in non-coding regulatory DNA sequences can alter gene expression, organismal phenotype and fitness1-3. Constructing complete fitness landscapes, in which DNA sequences are mapped to fitness, is a long-standing goal in biology, but has remained elusive because it is challenging to generalize reliably to vast sequence spaces4-6. Here we build sequence-to-expression models that capture fitness landscapes and use them to decipher principles of regulatory evolution. Using millions of randomly sampled promoter DNA sequences and their measured expression levels in the yeast Saccharomyces cerevisiae, we learn deep neural network models that generalize with excellent prediction performance, and enable sequence design for expression engineering. Using our models, we study expression divergence under genetic drift and strong-selection weak-mutation regimes to find that regulatory evolution is rapid and subject to diminishing returns epistasis; that conflicting expression objectives in different environments constrain expression adaptation; and that stabilizing selection on gene expression leads to the moderation of regulatory complexity. We present an approach for using suchmodels to detect signatures of selection on expression from natural variation in regulatory sequences and use it to discover an instance of convergent regulatory evolution. We assess mutational robustness, finding that regulatory mutation effect sizes follow a power law, characterize regulatory evolvability, visualize promoter fitness landscapes, discover evolvability archetypes and illustrate the mutational robustness of natural regulatory sequence populations. Our work provides a general framework for designing regulatorysequences andaddressing fundamental questions in regulatory evolution.

Genomic features, aroma profiles, and probiotic potential of the Debaryomyces hansenii species complex strains isolated from Korean soybean fermented food

Fermented soybean products are gaining attention in the food industry owing to their nutritive value and health benefits. In this study, we performed genomic analysis and physiological characterization of two Debaryomyces spp. yeast isolates obtained from a Korean traditional fermented soy sauce “ganjang”. Both Debaryomyces hansenii ganjang isolates KD2 and C11 showed halotolerance to concentrations of up to 15% NaCl and improved growth in the presence of salt. Ploidy and whole-genome sequencing analyses indicated that the KD2 genome is haploid, whereas the C11 genome is heterozygous diploid with two distinctive subgenomes. Interestingly, phylogenetic analysis using intron sequences indicated that the C11 strain was generated via hybridization between D. hansenii and D. tyrocola ancestor strains. The D. hansenii KD2 and D. hansenii-hybrid C11 produced various volatile flavor compounds associated with butter, caramel, cheese, and fruits, and showed high bioconversion activity from ferulic acid to 4-vinylguaiacol, a characteristic flavor compound of soybean products. Both KD2 and C11 exhibited viability in the presence of bile salts and at low pH and showed immunomodulatory activity to induce high levels of the anti-inflammatory cytokine IL-10. The safety of the yeast isolates was confirmed by analyzing virulence and acute oral toxicity. Together, the D. hansenii ganjang isolates possess physiological properties beneficial for improving the flavor and nutritional value of fermented products.

Genome-wide identification, characterization and relative expression analysis of putative iron homeostasis genes: NAS, NAAT, and DMAS in hexaploid wheat and its progenitors

The present study used different in silico approaches to identify and characterize the putative iron homeostasis genes in hexaploid wheat (Chinese Spring) and its diploid wild progenitors (Triticum urartu, Aegilops speltoides, and Aegilops tauschii). A total of six NAAT genes, three NAAT1 and other three NAAT2 types on long arms of group 1 chromosomes, three DMAS genes on group 4 chromosomes and twenty NAS genes on group 2, 3, 4, 5 and 6 chromosomes were identified and characterized in hexaploid wheat whereas, two NAAT genes from T.urartu (TuNAAT1 and TuNAAT2) and one each from Ae. speltoides (Ae.spelNAAT2) and Ae.tauschii (Ae.tausNAAT1), a single DMAS gene each from T. urartu, Ae. speltoides and Ae.tauschii and five NAS genes each from the diploid wild progenitors was reported on similar group of chromosomes. Structure analysis among the hexaploid wheat and wild progenitors revealed differences in lengths of NAAT and DMAS genes. This variation in the length of homologous genes was found to be associated with variation in the length of intron sequences which may influence the expression pattern of these genes. Further, gene expression analyses (quantitative RT-PCR) of four of the iron homeostasis genes, namely, TaNAAT, TaDMAS, TaIDEF2, and TaNAS9-2AS was performed in hexaploid wheat. Expression of three genes, TaNAAT, TaDMAS, and TaNAS9-2AS in wheat roots was found to be directly regulated by the iron deficiency and age of the plants, whereas TaIDEF2 was appeared to be post-transcriptionally regulated. The identification and characterization of these genes, provides a useful genetic resource for better understanding of iron homeostasis genes in hexaploid wheat. This information could be used in subsequent wheat biofortification programmes using molecular breeding approaches.

Comparative Chloroplast Genome Analysis of Wax Gourd (Benincasa hispida) with Three Benincaseae Species, Revealing Evolutionary Dynamic Patterns and Phylogenetic Implications.

Benincasa hispida (wax gourd) is an important Cucurbitaceae crop, with enormous economic and medicinal importance. Here, we report the de novo assembly and annotation of the complete chloroplast genome of wax gourd with 156,758 bp in total. The quadripartite structure of the chloroplast genome comprises a large single-copy (LSC) region with 86,538 bp and a small single-copy (SSC) region with 18,060 bp, separated by a pair of inverted repeats (IRa and IRb) with 26,080 bp each. Comparison analyses among B. hispida and three other species from Benincaseae presented a significant conversion regarding nucleotide content, genome structure, codon usage, synonymous and non-synonymous substitutions, putative RNA editing sites, microsatellites, and oligonucleotide repeats. The LSC and SSC regions were found to be much more varied than the IR regions through a divergent analysis of the species within Benincaseae. Notable IR contractions and expansions were observed, suggesting a difference in genome size, gene duplication and deletion, and the presence of pseudogenes. Intronic gene sequences, such as trnR-UCU–atpA and atpH–atpI, were observed as highly divergent regions. Two types of phylogenetic analysis based on the complete cp genome and 72 genes suggested sister relationships between B. hispida with the Citrullus, Lagenaria, and Cucumis. Variations and consistency with previous studies regarding phylogenetic relationships are discussed. The cp genome of B. hispida provides valuable genetic information for the detection of molecular markers, research on taxonomic discrepancies, and the inference of the phylogenetic relationships of Cucurbitaceae.