Multiple horizontal transfer events of a DNA transposon into turtles, fishes, and a frog

Horizontal transfer of transposable elements (HTT) has been reported across many species and the impact of such events on genome structure and function has been well described. However, few studies have focused on reptilian genomes, especially HTT events in Testudines (turtles). Here, as a consequence of investigating the repetitive content of Malaclemys terrapin terrapin (Diamondback turtle) we found a high similarity DNA transposon, annotated in RepBase as hAT-6_XT, shared between other turtle species, ray-finned fishes, and a frog. hAT-6_XT was notably absent in reptilian taxa closely related to turtles, such as crocodiles and birds. Successful invasion of DNA transposons into new genomes requires the conservation of specific residues in the encoded transposase, and through structural analysis, these residues were identified indicating some retention of functional transposition activity. We document six recent independent HTT events of a DNA transposon in turtles, which are known to have a low genomic evolutionary rate and ancient repeats.

Differential allelic representation (DAR) identifies candidate eQTLs and improves transcriptome analysis

In comparisons between mutant and wild-type genotypes, transcriptome analysis can reveal the direct impacts of a mutation, together with the homeostatic responses of the biological system. Recent studies have highlighted that, when the effects of homozygosity for recessive mutations are studied in non-isogenic backgrounds, genes located proximal to the mutation on the same chromosome often appear over-represented among those genes identified as differentially expressed (DE). One hypothesis suggests that DE genes chromosomally linked to a mutation may not reflect functional responses to the mutation but, instead, result from an unequal distribution of expression quantitative trait loci (eQTLs) between sample groups of mutant or wild-type genotypes. This is problematic because eQTL expression differences are difficult to distinguish from genes that are DE due to functional responses to a mutation. Here we show that chromosomally co-located differentially expressed genes (CC-DEGs) are also observed in analyses of dominant mutations in heterozygotes. We define a method and a metric to quantify, in RNA-sequencing data, localised differential allelic representation (DAR) between those sample groups subjected to differential expression analysis. We show how the DAR metric can predict regions prone to eQTL-driven differential expression, and how it can improve functional enrichment analyses through gene exclusion or weighting-based approaches. Advantageously, this improved ability to identify probable eQTLs also reveals examples of CC-DEGs that are likely to be functionally related to a mutant phenotype. This supports a long-standing prediction that selection for advantageous linkage disequilibrium influences chromosome evolution. By comparing the genomes of zebrafish (Danio rerio) and medaka (Oryzias latipes), a teleost with a conserved ancestral karyotype, we find possible examples of chromosomal aggregation of CC-DEGs during evolution of the zebrafish lineage. Our method for DAR analysis requires only RNA-sequencing data, facilitating its application across new and existing datasets.

LncRNA DANA1 promotes drought tolerance and histone deacetylation of drought responsive genes in Arabidopsis

Although many long noncoding RNAs have been discovered in plants, little is known about their biological function and mode of action. Here we show that the drought-induced long intergenic noncoding RNA DANA1 interacts with the L1p/L10e family member protein DANA1-INTERACTING PROTEIN 1 (DIP1) in the cell nucleus of Arabidopsis, and both DANA1 and DIP1 promote plant drought resistance. DANA1 and DIP1 increase histone deacetylase HDA9 binding to the CYP707A1 and CYP707A2 loci. DIP1 further interacts with PWWP3, a member of the PEAT complex that associates with HDA9 and has histone deacetylase activity. Mutation of DANA1 enhances CYP707A1 and CYP707A2 acetylation and expression resulting in impaired drought tolerance, in agreement with dip1 and pwwp3 mutant phenotypes. Our results demonstrate that DANA1 is a positive regulator of drought response and that DANA1 works jointly with the novel chromatin-related factor DIP1 on epigenetic reprogramming of the plant transcriptome during the response to drought.

Fake IDs? Widespread misannotation of DNA transposons as a general transcription factor

Accurate annotation of genes and transposable elements (TEs) is vital for understanding genomes, but current annotation pipelines often misannotate TEs as genes. This study reveals how the general transcription factor II-I repeat domain-containing protein 2 (GTF2IRD2) erroneously annotated DNA transposons in non-mammalian species, as it contains a 3' fused hAT transposase domain. We also demonstrate the generality of this problem by identifying misannotated TEs as genes in other vertebrate genomes. Such misannotations can lead to errors in phylogenetic analyses and wasted time for investigators. The study proposes adding a final TE-check to gene annotation pipelines to mitigate this problem.

Chromosomal level genome assembly of medicinal plant Sophora flavescens

Sophora flavescens is a medicinal plant in the genus Sophora of the Fabaceae family. The root of S. flavescens is known in China as Kushen and has a long history of wide use in multiple formulations of Traditional Chinese Medicine (TCM). In this study, we used third-generation Nanopore long-read sequencing technology combined with Hi-C scaffolding technology to de novo assemble the S. flavescens genome. We obtained a chromosomal level high-quality S. flavescens draft genome. The draft genome size is approximately 2.08 Gb, with more than 80% annotated as Transposable Elements (TEs), which have recently and rapidly proliferated. This genome size is ~5x larger than its closest sequenced relative Lupinus albus L. . We annotated 60,485 genes and examined their expression profiles in leaf, stem and root tissues, and also characterised the genes and pathways involved in the biosynthesis of major bioactive compounds, including alkaloids, flavonoids and isoflavonoids. The assembled genome highlights the very different evolutionary trajectories that have occurred in recently diverged Fabaceae, leading to smaller duplicated genomes.

The long non-coding RNA DANA2 positively regulates drought tolerance by recruiting ERF84 to promote JMJ29-mediated histone demethylation

Tens of thousands of long non-coding RNAs have been uncovered in plants, but few of them have been comprehensively studied for their biological function and molecular mechanism of their mode of action. Here, we show that the Arabidopsis long non-coding RNA DANA2 interacts with an AP2/ERF transcription factor ERF84 in the cell nucleus and then affects the transcription of JMJ29 that encodes a Jumonji C domain-containing histone H3K9 demethylase. Both RNA sequencing (RNA-seq) and genetic analyses demonstrate that DANA2 positively regulates drought stress responses through JMJ29. JMJ29 positively regulates the expression of ERF15 and GOLS2 by modulation of H3K9me2 demethylation. Accordingly, mutation of JMJ29 causes decreased ERF15 and GOLS2 expression, resulting in impaired drought tolerance, in agreement with drought-sensitive phenotypes of dana2 and erf84 mutants. Taken together, these results demonstrate that DANA2 is a positive regulator of drought response and works jointly with the transcriptional activator ERF84 to modulate JMJ29 expression in plant response to drought.

Long noncoding RNA ARTA controls ABA response through MYB7 nuclear trafficking in Arabidopsis

In eukaryotes, transcription factors are a crucial element in the regulation of gene expression, and nuclear translocation is the key to the function of transcription factors. Here, we show that the long intergenic noncoding RNA ARTA interacts with an importin β-like protein, SAD2, through a long noncoding RNA-binding region embedded in the carboxyl terminal, and then it blocks the import of the transcription factor MYB7 into the nucleus. Abscisic acid (ABA)-induced ARTA expression can positively regulate ABI5 expression by fine-tuning MYB7 nuclear trafficking. Therefore, the mutation of arta represses ABI5 expression, resulting in desensitization to ABA, thereby reducing Arabidopsis drought tolerance. Our results demonstrate that lncRNA can hijack a nuclear trafficking receptor to modulate the nuclear import of a transcription factor during plant responses to environmental stimuli.

Detecting Horizontal Transfer of Transposons

Transposable elements (TEs) are prevalent genomic components which can replicate as a function of mobilization in eukaryotes. Not only do they alter genome structure, they also play regulatory functions or organize chromatin structure. In addition to vertical parent-to-offspring inheritance, TEs can also horizontally "jump" between species, known as horizontal transposon transfer (HTT). This can rapidly alter the course of genome evolution. In this chapter, we provide a practical framework to detect HTT events. Our HTT detection framework is based on the use of sequence alignment to determine the divergence/conservation profiles of TE families to determine the history of expansion events. In summary, it includes (a) workflow of HTT detection from Ab initio identified TEs; (b) workflow for detecting HTT for specific, curated TEs; and (c) workflow for validating detected HTT candidates. Our framework covers two common scenarios of HTT detection in the modern omics era, and we believe it will serve as a valuable toolbox for the TE and genomics research community.

Transient Polycomb activity represses developmental genes in growing oocytes

BACKGROUND

Non-genetic disease inheritance and offspring phenotype are substantially influenced by germline epigenetic programming, including genomic imprinting. Loss of Polycomb Repressive Complex 2 (PRC2) function in oocytes causes non-genetically inherited effects on offspring, including embryonic growth restriction followed by post-natal offspring overgrowth. While PRC2-dependent non-canonical imprinting is likely to contribute, less is known about germline epigenetic programming of non-imprinted genes during oocyte growth. In addition, de novo germline mutations in genes encoding PRC2 lead to overgrowth syndromes in human patients, but the extent to which PRC2 activity is conserved in human oocytes is poorly understood.

RESULTS

In this study, we identify a discrete period of early oocyte growth during which PRC2 is expressed in mouse growing oocytes. Deletion of Eed during this window led to the de-repression of 343 genes. A high proportion of these were developmental regulators, and the vast majority were not imprinted genes. Many of the de-repressed genes were also marked by the PRC2-dependent epigenetic modification histone 3 lysine 27 trimethylation (H3K27me3) in primary-secondary mouse oocytes, at a time concurrent with PRC2 expression. In addition, we found H3K27me3 was also enriched on many of these genes by the germinal vesicle (GV) stage in human oocytes, strongly indicating that this PRC2 function is conserved in the human germline. However, while the 343 genes were de-repressed in mouse oocytes lacking EED, they were not de-repressed in pre-implantation embryos and lost H3K27me3 during pre-implantation development. This implies that H3K27me3 is a transient feature that represses a wide range of genes in oocytes.

CONCLUSIONS

Together, these data indicate that EED has spatially and temporally distinct functions in the female germline to repress a wide range of developmentally important genes and that this activity is conserved in the mouse and human germlines.

Compound Kushen injection reduces severity of radiation-induced gastrointestinal mucositis in rats

Mucositis, or damage/injury to mucous membranes of the alimentary, respiratory, or genitourinary tract, is the major side effect associated with anticancer radiotherapies. Because there is no effective treatment for mucositis at present, this is a particular issue as it limits the dose of therapy in cancer patients and significantly affects their quality of life. Gastrointestinal mucositis (GIM) occurs in patients receiving radiotherapies to treat cancers of the stomach, abdomen, and pelvis. It involves inflammation and ulceration of the gastrointestinal (GI) tract causing diarrhea, nausea and vomiting, abdominal pain, and bloating. However, there is currently no effective treatment for this debilitating condition. In this study, we investigated the potential of a type of traditional Chinese medicine (TCM), compound Kushen injection (CKI), as a treatment for GIM. It has previously been shown that major groups of chemical compounds found in CKI have anti-inflammatory effects and are capable of inhibiting the expression of pro-inflammatory cytokines. Intraperitoneal administration of CKI to Sprague Dawley (SD) rats that concurrently received abdominal irradiation over five fractions resulted in reduced severity of GIM symptoms compared to rats administered a vehicle control. Histological examination of the intestinal tissues revealed significantly less damaged villus epithelium in CKI-administered rats that had reduced numbers of apoptotic cells in the crypts. Furthermore, it was also found that CKI treatment led to decreased levels of inflammatory factors including lower levels of interleukin (IL)-1β and IL-6 as well as myeloperoxidase (MPO)-producing cells in the intestinal mucosa. Together, our data indicate a novel effect of CKI to reduce the symptoms of radiation-induced GIM by inhibiting inflammation in the mucosa and apoptosis of epithelial cells.

Horizontal Transposon Transfer and Its Implications for the Ancestral Ecology of Hydrophiine Snakes

Transposable elements (TEs), also known as jumping genes, are sequences able to move or copy themselves within a genome. As TEs move throughout genomes they often act as a source of genetic novelty, hence understanding TE evolution within lineages may help in understanding environmental adaptation. Studies into the TE content of lineages of mammals such as bats have uncovered horizontal transposon transfer (HTT) into these lineages, with squamates often also containing the same TEs. Despite the repeated finding of HTT into squamates, little comparative research has examined the evolution of TEs within squamates. Here we examine a diverse family of Australo-Melanesian snakes (Hydrophiinae) to examine if the previously identified, order-wide pattern of variable TE content and activity holds true on a smaller scale. Hydrophiinae diverged from Asian elapids ~30 Mya and have since rapidly diversified into six amphibious, ~60 marine and ~100 terrestrial species that fill a broad range of ecological niches. We find TE diversity and expansion differs between hydrophiines and their Asian relatives and identify multiple HTTs into Hydrophiinae, including three likely transferred into the ancestral hydrophiine from fish. These HTT events provide the first tangible evidence that Hydrophiinae reached Australia from Asia via a marine route.

Future Directions for the Discovery of Natural Product-Derived Immunomodulating Drugs

Drug discovery from natural sources is going through a renaissance, having spent many decades in the shadow of synthetic molecule drug discovery, despite the fact that natural product-derived compounds occupy a much greater chemical space than those created through synthetic chemistry methods. With this new era comes new possibilities, not least the novel targets that have emerged in recent times and the development of state-of-the-art technologies that can be applied to drug discovery from natural sources. Although progress has been made with some immunomodulating drugs, there remains a pressing need for new agents that can be used to treat the wide variety of conditions that arise from disruption, or over-activation, of the immune system; natural products may therefore be key in filling this gap. Recognising that, at present, there is no authoritative article that details the current state-of-the-art of the immunomodulatory activity of natural products, this in-depth review has arisen from a joint effort between the International Union of Basic and Clinical Pharmacology (IUPHAR) Natural Products and Immunopharmacology, with contributions from a Powered by Editorial Manager® and ProduXion Manager® from Aries Systems Corporation number of world-leading researchers in the field of natural product drug discovery, to provide a "position statement" on what natural products has to offer in the search for new immunomodulatory argents. To this end, we provide a historical look at previous discoveries of naturally occurring immunomodulators, present a picture of the current status of the field and provide insight into the future opportunities and challenges for the discovery of new drugs to treat immune-related diseases.

Genome Stability Is in the Eye of the Beholder: CR1 Retrotransposon Activity Varies Significantly across Avian Diversity

Since the sequencing of the zebra finch genome it has become clear that avian genomes, while largely stable in terms of chromosome number and gene synteny, are more dynamic at an intrachromosomal level. A multitude of intrachromosomal rearrangements and significant variation in transposable element (TE) content have been noted across the avian tree. TEs are a source of genome plasticity, because their high similarity enables chromosomal rearrangements through nonallelic homologous recombination, and they have potential for exaptation as regulatory and coding sequences. Previous studies have investigated the activity of the dominant TE in birds, chicken repeat 1 (CR1) retrotransposons, either focusing on their expansion within single orders, or comparing passerines with nonpasserines. Here, we comprehensively investigate and compare the activity of CR1 expansion across orders of birds, finding levels of CR1 activity vary significantly both between and within orders. We describe high levels of TE expansion in genera which have speciated in the last 10 Myr including kiwis, geese, and Amazon parrots; low levels of TE expansion in songbirds across their diversification, and near inactivity of TEs in the cassowary and emu for millions of years. CR1s have remained active over long periods of time across most orders of neognaths, with activity at any one time dominated by one or two families of CR1s. Our findings of higher TE activity in species-rich clades and dominant families of TEs within lineages mirror past findings in mammals and indicate that genome evolution in amniotes relies on universal TE-driven processes.

Horizontal transfer and subsequent explosive expansion of a DNA transposon in sea kraits (Laticauda)

Transposable elements (TEs) are self-replicating genetic sequences and are often described as important ‘drivers of evolution’. This driving force is because TEs promote genomic novelty by enabling rearrangement, and through exaptation as coding and regulatory elements. However, most TE insertions potentially lead to neutral or harmful outcomes, therefore host genomes have evolved machinery to suppress TE expansion. Through horizontal transposon transfer (HTT) TEs can colonize new genomes, and since new hosts may not be able to regulate subsequent replication, these TEs may proliferate rapidly. Here, we describe HTT of the Harbinger-Snek DNA transposon into sea kraits (Laticauda), and its subsequent explosive expansion within Laticauda genomes. This HTT occurred following the divergence of Laticauda from terrestrial Australian elapids approximately 15–25 Mya. This has resulted in numerous insertions into introns and regulatory regions, with some insertions into exons which appear to have altered UTRs or added sequence to coding exons. Harbinger-Snek has rapidly expanded to make up 8–12% of Laticauda spp. genomes; this is the fastest known expansion of TEs in amniotes following HTT. Genomic changes caused by this rapid expansion may have contributed to adaptation to the amphibious-marine habitat.

The long noncoding RNA FRILAIR regulates strawberry fruit ripening by functioning as a noncanonical target mimic

Long noncoding RNAs (lncRNAs) are emerging as important regulators in plant development, but few of them have been functionally characterized in fruit ripening. Here, we have identified 25,613 lncRNAs from strawberry ripening fruits based on RNA-seq data from poly(A)-depleted libraries and rRNA-depleted libraries, most of which exhibited distinct temporal expression patterns. A novel lncRNA, FRILAIR harbours the miR397 binding site that is highly conserved in diverse strawberry species. FRILAIR overexpression promoted fruit maturation in the Falandi strawberry, which was consistent with the finding from knocking down miR397, which can guide the mRNA cleavage of both FRILAIR and LAC11a (encoding a putative laccase-11-like protein). Moreover, LAC11a mRNA levels were increased in both FRILAIR overexpressing and miR397 knockdown fruits, and accelerated fruit maturation was also found in LAC11a overexpressing fruits. Overall, our study demonstrates that FRILAIR can act as a noncanonical target mimic of miR397 to modulate the expression of LAC11a in the strawberry fruit ripening process.

New Environment, New Invaders-Repeated Horizontal Transfer of LINEs to Sea Snakes

Although numerous studies have found horizontal transposon transfer (HTT) to be widespread across metazoans, few have focused on HTT in marine ecosystems. To investigate potential recent HTTs into marine species, we searched for novel repetitive elements in sea snakes, a group of elapids which transitioned to a marine habitat at most 18 Ma. Our analysis uncovered repeated HTTs into sea snakes following their marine transition. The seven subfamilies of horizontally transferred LINE retrotransposons we identified in the olive sea snake (Aipysurus laevis) are transcribed, and hence are likely still active and expanding across the genome. A search of 600 metazoan genomes found all seven were absent from other amniotes, including terrestrial elapids, with the most similar LINEs present in fish and marine invertebrates. The one exception was a similar LINE found in sea kraits, a lineage of amphibious elapids which independently transitioned to a marine environment 25 Ma. Our finding of repeated horizontal transfer events into marine snakes greatly expands past findings that the marine environment promotes the transfer of transposons. Transposons are drivers of evolution as sources of genomic sequence and hence genomic novelty. We identified 13 candidate genes for HTT-induced adaptive change based on internal or neighboring HTT LINE insertions. One of these, ADCY4, is of particular interest as a part of the KEGG adaptation pathway “Circadian Entrainment.” This provides evidence of the ecological interactions between species influencing evolution of metazoans not only through specific selection pressures, but also by contributing novel genomic material.

The tuatara genome reveals ancient features of amniote evolution

The tuatara (Sphenodon punctatus)-the only living member of the reptilian order Rhynchocephalia (Sphenodontia), once widespread across Gondwana1,2-is an iconic species that is endemic to New Zealand2,3. A key link to the now-extinct stem reptiles (from which dinosaurs, modern reptiles, birds and mammals evolved), the tuatara provides key insights into the ancestral amniotes2,4. Here we analyse the genome of the tuatara, which-at approximately 5 Gb-is among the largest of the vertebrate genomes yet assembled. Our analyses of this genome, along with comparisons with other vertebrate genomes, reinforce the uniqueness of the tuatara. Phylogenetic analyses indicate that the tuatara lineage diverged from that of snakes and lizards around 250 million years ago. This lineage also shows moderate rates of molecular evolution, with instances of punctuated evolution. Our genome sequence analysis identifies expansions of proteins, non-protein-coding RNA families and repeat elements, the latter of which show an amalgam of reptilian and mammalian features. The sequencing of the tuatara genome provides a valuable resource for deep comparative analyses of tetrapods, as well as for tuatara biology and conservation. Our study also provides important insights into both the technical challenges and the cultural obligations that are associated with genome sequencing.

The effect of compound kushen injection on cancer cells: Integrated identification of candidate molecular mechanisms

Traditional Chinese Medicine (TCM) preparations are often extracts of single or multiple herbs containing hundreds of compounds, and hence it has been difficult to study their mechanisms of action. Compound Kushen Injection (CKI) is a complex mixture of compounds extracted from two medicinal plants and has been used in Chinese hospitals to treat cancer for over twenty years. To demonstrate that a systematic analysis of molecular changes resulting from complex mixtures of bioactives from TCM can identify a core set of differentially expressed (DE) genes and a reproducible set of candidate pathways. We used in vitro cancer models to measure the effect of CKI on cell cycle phases and apoptosis, and correlated those phenotypes with CKI induced changes in gene expression. We treated two cancer cell lines with or without CKI and assessed the resulting phenotypes by employing cell viability and proliferation assays. Based on these results, we carried out high-throughput transcriptome data analysis to identify genes and candidate pathways perturbed by CKI. We integrated these differential gene expression results with previously reported results and carried out validation of selected differentially expressed genes. CKI induced cell-cycle arrest and apoptosis in the cancer cell lines tested. In these cells CKI also altered the expression of 363 core candidate genes associated with cell cycle, apoptosis, DNA replication/repair, and various cancer pathways. Of these, 7 are clinically relevant to cancer diagnosis or therapy, 14 are cell cycle regulators, and most of these 21 candidates are downregulated by CKI. Comparison of our core candidate genes to a database of plant medicinal compounds and their effects on gene expression identified one-to-one, one-to-many and many-to-many regulatory relationships between compounds in CKI and DE genes. By identifying genes and promising candidate pathways associated with CKI treatment based on our transcriptome-based analysis, we have shown that this approach is useful for the systematic analysis of molecular changes resulting from complex mixtures of bioactives.

A New Strategy for Identifying Mechanisms of Drug-drug Interaction Using Transcriptome Analysis: Compound Kushen Injection as a Proof of Principle

Drug-drug interactions (DDIs), especially with herbal medicines, are complex, making it difficult to identify potential molecular mechanisms and targets. We introduce a workflow to carry out DDI research using transcriptome analysis and interactions of a complex herbal mixture, Compound Kushen Injection (CKI), with cancer chemotherapy drugs, as a proof of principle. Using CKI combined with doxorubicin or 5-Fu on cancer cells as a model, we found that CKI enhanced the cytotoxic effects of doxorubicin on A431 cells while protecting MDA-MB-231 cells treated with 5-Fu. We generated and analysed transcriptome data from cells treated with single treatments or combined treatments and our analysis showed that opposite directions of regulation for pathways related to DNA synthesis and metabolism which appeared to be the main reason for different effects of CKI when used in combination with chemotherapy drugs. We also found that pathways related to organic biosynthetic and metabolic processes might be potential targets for CKI when interacting with doxorubicin and 5-Fu. Through co-expression analysis correlated with phenotype results, we selected the MYD88 gene as a candidate major regulator for validation as a proof of concept for our approach. Inhibition of MYD88 reduced antagonistic cytotoxic effects between CKI and 5-Fu, indicating that MYD88 is an important gene in the DDI mechanism between CKI and chemotherapy drugs. These findings demonstrate that our pipeline is effective for the application of transcriptome analysis to the study of DDIs in order to identify candidate mechanisms and potential targets.

Fractional Deletion of Compound Kushen Injection Indicates Cytokine Signaling Pathways are Critical for its Perturbation of the Cell Cycle

We used computational and experimental biology approaches to identify candidate mechanisms of action of aTraditional Chinese Medicine, Compound Kushen Injection (CKI), in a breast cancer cell line (MDA-MB-231). Because CKI is a complex mixture of plant secondary metabolites, we used a high-performance liquid chromatography (HPLC) fractionation and reconstitution approach to define chemical fractions required for CKI to induce apoptosis. The initial fractionation separated major from minor compounds, and it showed that major compounds accounted for little of the activity of CKI. Furthermore, removal of no single major compound altered the effect of CKI on cell viability and apoptosis. However, simultaneous removal of two major compounds identified oxymatrine and oxysophocarpine as critical with respect to CKI activity. Transcriptome analysis was used to correlate compound removal with gene expression and phenotype data. Many compounds in CKI are required to trigger apoptosis but significant modulation of its activity is conferred by a small number of compounds. In conclusion, CKI may be typical of many plant based extracts that contain many compounds in that no single compound is responsible for all of the bioactivity of the mixture and that many compounds interact in a complex fashion to influence a network containing many targets.

Yiqi Chutan Tang Reduces Gefitinib-Induced Drug Resistance in Non-Small-Cell Lung Cancer by Targeting Apoptosis and Autophagy

High incidence and mortality rates for non-small-cell lung cancer (NSCLC) lead to low survival rates. Epidermal growth factor receptor-tyrosine kinase inhibitors (EGFR-TKI) are commonly first prescribed for NSCLC patients with EGFR mutations. However, most patients with sensitizing EGFR mutations become resistant to EGFR-TKI after 9-13 months treatment. Yiqi Chutan Tang (YQCT) has been prescribed as a treatment to this issue for over 20 years. In this report, high-performance liquid chromatography (HPLC) analysis, flow cytometry, western blot analysis, and functional annotation analysis were applied to uncover the molecular mechanisms of YQCT. Our results show the application of YQCT reduces gefitinib-induced drug resistance, induces slight cell cycle arrest, enhances gefitinib-induced apoptosis, and activates the autophagy. These results indicate that at the molecular level YQCT can reduce drug resistance and improve anti-cancer effects when associated with gefitinib, which could be a result of enhancement of apoptosis and autophagy in the EGFR-TKI resistant cells of NSCLC. This research provides a new treatment strategy for patients with EGFR-TKI resistance in NSCLC. © 2019 The Authors. Cytometry Part A published by Wiley Periodicals, Inc. on behalf of International Society for Advancement of Cytometry.

An effective drug sensitizing agent increases gefitinib treatment by down regulating PI3K/Akt/mTOR pathway and up regulating autophagy in non-small cell lung cancer

Gefitinib is one of commonly used first-line treatment options for patients with positive EGFR mutation in non-small cell lung cancer (NSCLC). However, most patients with gefitinib treatment relapse over time due to the loss of drug sensitivity. Compound Kushen injection (CKI) has been used to treat lung cancer, including EGFR-mutated NSCLC. In this report, we examined the anti-cancer and drug sensitivity increased activities of CKI in gefitinib less sensitive NSCLC cell lines H1650 and H1975. Bioinformatics analysis was applied to uncover gene regulation and molecular mechanisms of CKI. Our results indicated that when associating with gefitinib in a dose-dependent fashion, CKI demonstrated the ability to inhibit the proliferation and to increase the sensitivity to gefitinib treatment in gefitinib less sensitive cell lines. This could be the results of down regulation of the PI3K/Akt/mTOR pathway and up regulation of autophagy, which were identified as the potential primary targets of CKI to increase gefitinib treatment effect.

Understanding the Mechanistic Contribution of Herbal Extracts in Compound Kushen Injection With Transcriptome Analysis

Herbal compatibility is the knowledge of which herbs to combine in traditional Chinese medicine (TCM) formulations. The lack of understanding of herbal compatibility is one of the key problems for the application and popularization of TCM in western society. Because of the chemical complexity of herbal medicines, it is simpler to begin to conduct compatibility research based on herbs rather than component plant secondary metabolites. We have used transcriptome analysis to explore the effects and interactions of two plant extracts (Kushen and Baituling) combined in Compound Kushen Injection (CKI). Based on shared chemical compounds and in vitro cytotoxicity comparisons, we found that both the major compounds in CKI, and the cytotoxicity effects of CKI were mainly derived from the extract of Kushen (Sophorae flavescentis). We generated and analyzed transcriptome data from MDA-MB-231 cells treated with single-herb extracts or CKI and results showed that Kushen contributed to the perturbation of the majority of cytotoxicity/cancer related pathways in CKI such as cell cycle and DNA replication. We also found that Baituling (Heterosmilax yunnanensis Gagnep) could not only enhance the cytotoxic effects of Kushen in CKI, but also activate immune-related pathways. Our analyses predicted that IL-1β gene expression was upregulated by Baituling in CKI and we confirmed that IL-1β protein expression was increased using an ELISA assay. Altogether, these findings help to explain the rationale for combining Kushen and Baituling in CKI, and show that transcriptome analysis using single herb extracts is an effective method for understanding herbal compatibility in TCM.

Effect of Compound Kushen Injection, a Natural Compound Mixture, and Its Identified Chemical Components on Migration and Invasion of Colon, Brain, and Breast Cancer Cell Lines

Traditional Chinese Medicines are promising sources of new agents for controlling cancer metastasis. Compound Kushen Injection (CKI), prepared from medicinal plants Sophora flavescens and Heterosmilax chinensis, disrupts cell cycle and induces apoptosis in breast cancer; however, effects on migration and invasion remained unknown. CKI, fractionated mixtures, and isolated components were tested in migration assays with colon (HT-29, SW-480, DLD-1), brain (U87-MG, U251-MG), and breast (MDA-MB-231) cancer cell lines. Human embryonic kidney (HEK-293) and human foreskin fibroblast (HFF) served as non-cancerous controls. Wound closure, transwell invasion, and live cell imaging showed CKI reduced motility in all eight lines. Fractionation and reconstitution of CKI demonstrated combinations of compounds were required for activity. Live cell imaging confirmed CKI strongly reduced migration of HT-29 and MDA-MB-231 cells, moderately slowed brain cancer cells, and had a small effect on HEK-293. CKI uniformly blocked invasiveness through extracellular matrix. Apoptosis was increased by CKI in breast cancer but not in non-cancerous lines. Cell viability was unaffected by CKI in all cell lines. Transcriptomic analyses of MDA-MB-231indicated down-regulation of actin cytoskeletal and focal adhesion genes with CKI treatment, consistent with observed impairment of cell migration. The pharmacological complexity of CKI is important for effective blockade of cancer migration and invasion.

Cell cycle, energy metabolism and DNA repair pathways in cancer cells are suppressed by Compound Kushen Injection

BACKGROUND

In this report we examine candidate pathways perturbed by Compound Kushen Injection (CKI), a Traditional Chinese Medicine (TCM) that we have previously shown to alter the gene expression patterns of multiple pathways and induce apoptosis in cancer cells.

METHODS

We have measured protein levels in Hep G2 and MDA-MB-231 cells for genes in the cell cycle pathway, DNA repair pathway and DNA double strand breaks (DSBs) previously shown to have altered expression by CKI. We have also examined energy metabolism by measuring [ADP]/[ATP] ratio (cell energy charge), lactate production and glucose consumption. Our results demonstrate that CKI can suppress protein levels for cell cycle regulatory proteins and DNA repair while increasing the level of DSBs. We also show that energy metabolism is reduced based on reduced glucose consumption and reduced cellular energy charge.

RESULTS

Our results validate these pathways as important targets for CKI. We also examined the effect of the major alkaloid component of CKI, oxymatrine and determined that it had no effect on DSBs, a small effect on the cell cycle and increased the cell energy charge.

CONCLUSIONS

Our results indicate that CKI likely acts through the effect of multiple compounds on multiple targets where the observed phenotype is the integration of these effects and synergistic interactions.

Horizontal transfer of BovB and L1 retrotransposons in eukaryotes

BACKGROUND

Transposable elements (TEs) are mobile DNA sequences, colloquially known as jumping genes because of their ability to replicate to new genomic locations. TEs can jump between organisms or species when given a vector of transfer, such as a tick or virus, in a process known as horizontal transfer. Here, we propose that LINE-1 (L1) and Bovine-B (BovB), the two most abundant TE families in mammals, were initially introduced as foreign DNA via ancient horizontal transfer events.

RESULTS

Using analyses of 759 plant, fungal and animal genomes, we identify multiple possible L1 horizontal transfer events in eukaryotic species, primarily involving Tx-like L1s in marine eukaryotes. We also extend the BovB paradigm by increasing the number of estimated transfer events compared to previous studies, finding new parasite vectors of transfer such as bed bug, leech and locust, and BovB occurrences in new lineages such as bat and frog. Given that these transposable elements have colonised more than half of the genome sequence in today's mammals, our results support a role for horizontal transfer in causing long-term genomic change in new host organisms.

CONCLUSIONS

We describe extensive horizontal transfer of BovB retrotransposons and provide the first evidence that L1 elements can also undergo horizontal transfer. With the advancement of genome sequencing technologies and bioinformatics tools, we anticipate our study to be a valuable resource for inferring horizontal transfer from large-scale genomic data.

Transposable elements and gene expression during the evolution of amniotes

Background

Transposable elements (TEs) are primarily responsible for the DNA losses and gains in genome sequences that occur over time within and between species. TEs themselves evolve, with clade specific LTR/ERV, LINEs and SINEs responsible for the bulk of species-specific genomic features. Because TEs can contain regulatory motifs, they can be exapted as regulators of gene expression. While TE insertions can provide evolutionary novelty for the regulation of gene expression, their overall impact on the evolution of gene expression is unclear. Previous investigators have shown that tissue specific gene expression in amniotes is more similar across species than within species, supporting the existence of conserved developmental gene regulation. In order to understand how species-specific TE insertions might affect the evolution/conservation of gene expression, we have looked at the association of gene expression in six tissues with TE insertions in six representative amniote genomes.

Results

A novel bootstrapping approach has been used to minimise the conflation of effects of repeat types on gene expression. We compared the expression of orthologs containing recent TE insertions to orthologs that contained older TE insertions, and the expression of non-orthologs containing recent TE insertions to non-orthologs with older TE insertions. Both orthologs and non-orthologs showed significant differences in gene expression associated with TE insertions. TEs were found associated with species-specific changes in gene expression, and the magnitude and direction of expression changes were noteworthy. Overall, orthologs containing species-specific TEs were associated with lower gene expression, while in non-orthologs, non-species specific TEs were associated with higher gene expression. Exceptions were SINE elements in human and chicken, which had an opposite association with gene expression compared to other species.

Conclusions

Our observed species-specific associations of TEs with gene expression support a role for TEs in speciation/response to selection by species. TEs do not exhibit consistent associations with gene expression and observed associations can vary depending on the age of TE insertions. Based on these observations, it would be prudent to refrain from extrapolating these and previously reported associations to distantly related species.

Divergent genome evolution caused by regional variation in DNA gain and loss between human and mouse

The forces driving the accumulation and removal of non-coding DNA and ultimately the evolution of genome size in complex organisms are intimately linked to genome structure and organisation. Our analysis provides a novel method for capturing the regional variation of lineage-specific DNA gain and loss events in their respective genomic contexts. To further understand this connection we used comparative genomics to identify genome-wide individual DNA gain and loss events in the human and mouse genomes. Focusing on the distribution of DNA gains and losses, relationships to important structural features and potential impact on biological processes, we found that in autosomes, DNA gains and losses both followed separate lineage-specific accumulation patterns. However, in both species chromosome X was particularly enriched for DNA gain, consistent with its high L1 retrotransposon content required for X inactivation. We found that DNA loss was associated with gene-rich open chromatin regions and DNA gain events with gene-poor closed chromatin regions. Additionally, we found that DNA loss events tended to be smaller than DNA gain events suggesting that they were able to accumulate in gene-rich open chromatin regions due to their reduced capacity to interrupt gene regulatory architecture. GO term enrichment showed that mouse loss hotspots were strongly enriched for terms related to developmental processes. However, these genes were also located in regions with a high density of conserved elements, suggesting that despite high levels of DNA loss, gene regulatory architecture remained conserved. This is consistent with a model in which DNA gain and loss results in turnover or "churning" in regulatory element dense regions of open chromatin, where interruption of regulatory elements is selected against.

Superior ab initio identification, annotation and characterisation of TEs and segmental duplications from genome assemblies

Transposable Elements (TEs) are mobile DNA sequences that make up significant fractions of amniote genomes. However, they are difficult to detect and annotate ab initio because of their variable features, lengths and clade-specific variants. We have addressed this problem by refining and developing a Comprehensive ab initio Repeat Pipeline (CARP) to identify and cluster TEs and other repetitive sequences in genome assemblies. The pipeline begins with a pairwise alignment using krishna, a custom aligner. Single linkage clustering is then carried out to produce families of repetitive elements. Consensus sequences are then filtered for protein coding genes and then annotated using Repbase and a custom library of retrovirus and reverse transcriptase sequences. This process yields three types of family: fully annotated, partially annotated and unannotated. Fully annotated families reflect recently diverged/young known TEs present in Repbase. The remaining two types of families contain a mixture of novel TEs and segmental duplications. These can be resolved by aligning these consensus sequences back to the genome to assess copy number vs. length distribution. Our pipeline has three significant advantages compared to other methods for ab initio repeat identification: 1) we generate not only consensus sequences, but keep the genomic intervals for the original aligned sequences, allowing straightforward analysis of evolutionary dynamics, 2) consensus sequences represent low-divergence, recently/currently active TE families, 3) segmental duplications are annotated as a useful by-product. We have compared our ab initio repeat annotations for 7 genome assemblies to other methods and demonstrate that CARP compares favourably with RepeatModeler, the most widely used repeat annotation package.

A comprehensive genomic history of extinct and living elephants

Elephantids are the world's most iconic megafaunal family, yet there is no comprehensive genomic assessment of their relationships. We report a total of 14 genomes, including 2 from the American mastodon, which is an extinct elephantid relative, and 12 spanning all three extant and three extinct elephantid species including an ∼120,000-y-old straight-tusked elephant, a Columbian mammoth, and woolly mammoths. Earlier genetic studies modeled elephantid evolution via simple bifurcating trees, but here we show that interspecies hybridization has been a recurrent feature of elephantid evolution. We found that the genetic makeup of the straight-tusked elephant, previously placed as a sister group to African forest elephants based on lower coverage data, in fact comprises three major components. Most of the straight-tusked elephant's ancestry derives from a lineage related to the ancestor of African elephants while its remaining ancestry consists of a large contribution from a lineage related to forest elephants and another related to mammoths. Columbian and woolly mammoths also showed evidence of interbreeding, likely following a latitudinal cline across North America. While hybridization events have shaped elephantid history in profound ways, isolation also appears to have played an important role. Our data reveal nearly complete isolation between the ancestors of the African forest and savanna elephants for ∼500,000 y, providing compelling justification for the conservation of forest and savanna elephants as separate species.

Similar Evolutionary Trajectories for Retrotransposon Accumulation in Mammals

The factors guiding retrotransposon insertion site preference are not well understood. Different types of retrotransposons share common replication machinery and yet occupy distinct genomic domains. Autonomous long interspersed elements accumulate in gene-poor domains and their nonautonomous short interspersed elements accumulate in gene-rich domains. To determine genomic factors that contribute to this discrepancy we analyzed the distribution of retrotransposons within the framework of chromosomal domains and regulatory elements. Using comparative genomics, we identified large-scale conserved patterns of retrotransposon accumulation across several mammalian genomes. Importantly, retrotransposons that were active after our sample-species diverged accumulated in orthologous regions. This suggested a similar evolutionary interaction between retrotransposon activity and conserved genome architecture across our species. In addition, we found that retrotransposons accumulated at regulatory element boundaries in open chromatin, where accumulation of particular retrotransposon types depended on insertion size and local regulatory element density. From our results, we propose a model where density and distribution of genes and regulatory elements canalize retrotransposon accumulation. Through conservation of synteny, gene regulation and nuclear organization, mammalian genomes with dissimilar retrotransposons follow similar evolutionary trajectories.

Seed weight differences between wild and domesticated soybeans are associated with specific changes in gene expression

KEY MESSAGE

Our study systematically explored potential genes and molecular pathways as candidates for differences in seed weight resulting from soybean domestication. In addition, potential contributions of lncRNAs to seed weight were also investigated. Soybeans have a long history of domestication in China, and there are several significant phenotypic differences between cultivated and wild soybeans, for example, seeds of cultivars are generally larger and heavier than those from wild accessions. We analyzed seed transcriptomes from thirteen soybean samples, including six landraces and seven wild accessions using strand-specific RNA sequencing. Differentially expressed genes related to seed weight were identified, and some of their homologs were associated with seed development in Arabidopsis. We also identified 1251 long intergenic noncoding RNAs (lincRNAs), 243 intronic RNAs and 81 antisense lncRNAs de novo from these soybean transcriptomes. We then profiled the expression patterns of lncRNAs in cultivated and wild soybean seeds, and found that transcript levels of a number of lncRNAs were sample-specific. Moreover, gene transcript and lincRNA co-expression network analysis showed that some soybean lincRNAs might have functional roles as they were hubs of co-expression modules. In conclusion, this study systematically explored potential genes and molecular pathways as candidates for differences in seed weight resulting from soybean domestication, and will provide a useful future resource for molecular breeding of soybeans.

Expression patterns of HENMT1 and PIWIL1 in human testis: implications for transposon expression

This study aimed to define the expression patterns of HENMT1 and PIWI proteins in human testis and investigate their association with transposon expression, infertility sub-type or development of testicular germ cell tumours (TGCTs). Testis biopsies showing normal spermatogenesis were used to identify normal localisation patterns of HENMT1 and PIWIL1 by immunolocalisation and RT-PCR after laser microdissection. 222 testis biopsies representing normal spermatogenesis, hypospermatogenesis, spermatogenic arrests, Sertoli cell-only (SCO) tumours and TGCTs were analysed by RT-qPCR for expression of HENMT1/PIWIL1/PIWIL2/PIWIL3/PIWIL4 and LINE-1 Additionally, HENMT1-overexpressing TCam2 seminoma cell lines were analysed for the same parameters by RT-qPCR. We found that HENMT1 and PIWIL1 are coexpressed in pachytene spermatocytes and spermatids. Expression of HENMT1, PIWIL1 and PIWIL2 was mainly dependent on germ cell content but low levels of expression were also detected in some SCO samples. Levels of HENMT1, PIWIL1 and PIWIL2 expression were low in TGCT. Samples with HENMT1, PIWIL2 and PIWIL4 expression showed significantly (P < 0.05) lower transposon expression compared to samples without expression in the same histological group. HENMT1-overexpressing TCam2 cells showed lower LINE-1 expression than empty vector-transfected control lines. Our findings support that the transposon-regulating function of the piRNA pathway found in the mouse is conserved in adult human testis. HENMT1 and PIWI proteins are expressed in a germ-cell-specific manner and required for transposon control.

Transposable elements (TEs) contribute to stress-related long intergenic noncoding RNAs in plants

Noncoding RNAs have been extensively described in plant and animal transcriptomes by using high-throughput sequencing technology. Of these noncoding RNAs, a growing number of long intergenic noncoding RNAs (lincRNAs) have been described in multicellular organisms, however the origins and functions of many lincRNAs remain to be explored. In many eukaryotic genomes, transposable elements (TEs) are widely distributed and often account for large fractions of plant and animal genomes yet the contribution of TEs to lincRNAs is largely unknown. By using strand-specific RNA-sequencing, we profiled the expression patterns of lincRNAs in Arabidopsis, rice and maize, and identified 47 611 and 398 TE-associated lincRNAs (TE-lincRNAs), respectively. TE-lincRNAs were more often derived from retrotransposons than DNA transposons and as retrotransposon copy number in both rice and maize genomes so did TE-lincRNAs. We validated the expression of these TE-lincRNAs by strand-specific RT-PCR and also demonstrated tissue-specific transcription and stress-induced TE-lincRNAs either after salt, abscisic acid (ABA) or cold treatments. For Arabidopsis TE-lincRNA11195, mutants had reduced sensitivity to ABA as demonstrated by longer roots and higher shoot biomass when compared to wild-type. Finally, by altering the chromatin state in the Arabidopsis chromatin remodelling mutant ddm1, unique lincRNAs including TE-lincRNAs were generated from the preceding untranscribed regions and interestingly inherited in a wild-type background in subsequent generations. Our findings not only demonstrate that TE-associated lincRNAs play important roles in plant abiotic stress responses but lincRNAs and TE-lincRNAs might act as an adaptive reservoir in eukaryotes.

A novel hypothesis-unbiased method for Gene Ontology enrichment based on transcriptome data

Gene Ontology (GO) classification of statistically significantly differentially expressed genes is commonly used to interpret transcriptomics data as a part of functional genomic analysis. In this approach, all significantly expressed genes contribute equally to the final GO classification regardless of their actual expression levels. Gene expression levels can significantly affect protein production and hence should be reflected in GO term enrichment. Genes with low expression levels can also participate in GO term enrichment through cumulative effects. In this report, we have introduced a new GO enrichment method that is suitable for multiple samples and time series experiments that uses a statistical outlier test to detect GO categories with special patterns of variation that can potentially identify candidate biological mechanisms. To demonstrate the value of our approach, we have performed two case studies. Whole transcriptome expression profiles of Salmonella enteritidis and Alzheimer's disease (AD) were analysed in order to determine GO term enrichment across the entire transcriptome instead of a subset of differentially expressed genes used in traditional GO analysis. Our result highlights the key role of inflammation related functional groups in AD pathology as granulocyte colony-stimulating factor receptor binding, neuromedin U binding, and interleukin were remarkably upregulated in AD brain when all using all of the gene expression data in the transcriptome. Mitochondrial components and the molybdopterin synthase complex were identified as potential key cellular components involved in AD pathology.

LINEs between Species: Evolutionary Dynamics of LINE-1 Retrotransposons across the Eukaryotic Tree of Life

LINE-1 (L1) retrotransposons are dynamic elements. They have the potential to cause great genomic change because of their ability to ‘jump’ around the genome and amplify themselves, resulting in the duplication and rearrangement of regulatory DNA. Active L1, in particular, are often thought of as tightly constrained, homologous and ubiquitous elements with well-characterized domain organization. For the past 30 years, model organisms have been used to define L1s as 6–8 kb sequences containing a 5′-UTR, two open reading frames working harmoniously in cis, and a 3′-UTR with a polyA tail. In this study, we demonstrate the remarkable and overlooked diversity of L1s via a comprehensive phylogenetic analysis of elements from over 500 species from widely divergent branches of the tree of life. The rapid and recent growth of L1 elements in mammalian species is juxtaposed against the diverse lineages found in other metazoans and plants. In fact, some of these previously unexplored mammalian species (e.g. snub-nosed monkey, minke whale) exhibit L1 retrotranspositional ‘hyperactivity’ far surpassing that of human or mouse. In contrast, non-mammalian L1s have become so varied that the current classification system seems to inadequately capture their structural characteristics. Our findings illustrate how both long-term inherited evolutionary patterns and random bursts of activity in individual species can significantly alter genomes, highlighting the importance of L1 dynamics in eukaryotes.

In depth analysis of the Sox4 gene locus that consists of sense and natural antisense transcripts

SRY (Sex Determining Region Y)-Box 4 or Sox4 is an important regulator of the pan-neuronal gene expression during post-mitotic cell differentiation within the mammalian brain. Sox4 gene locus has been previously characterized with multiple sense and overlapping natural antisense transcripts [1], [2]. Here we provide accompanying data on various analyses performed and described in Ling et al. [2]. The data include a detail description of various features found at Sox4 gene locus, additional experimental data derived from RNA-Fluorescence in situ Hybridization (RNA-FISH), Western blotting, strand-specific reverse-transcription quantitative polymerase chain reaction (RT-qPCR), gain-of-function and in situ hybridization (ISH) experiments. All the additional data provided here support the existence of an endogenous small interfering- or PIWI interacting-like small RNA known as Sox4_sir3, which origin was found within the overlapping region consisting of a sense and a natural antisense transcript known as Sox4ot1.

HENMT1 and piRNA Stability Are Required for Adult Male Germ Cell Transposon Repression and to Define the Spermatogenic Program in the Mouse

piRNAs are critical for transposable element (TE) repression and germ cell survival during the early phases of spermatogenesis, however, their role in adult germ cells and the relative importance of piRNA methylation is poorly defined in mammals. Using a mouse model of HEN methyltransferase 1 (HENMT1) loss-of-function, RNA-Seq and a range of RNA assays we show that HENMT1 is required for the 2’ O-methylation of mammalian piRNAs. HENMT1 loss leads to piRNA instability, reduced piRNA bulk and length, and ultimately male sterility characterized by a germ cell arrest at the elongating germ cell phase of spermatogenesis. HENMT1 loss-of-function, and the concomitant loss of piRNAs, resulted in TE de-repression in adult meiotic and haploid germ cells, and the precocious, and selective, expression of many haploid-transcripts in meiotic cells. Precocious expression was associated with a more active chromatin state in meiotic cells, elevated levels of DNA damage and a catastrophic deregulation of the haploid germ cell gene expression. Collectively these results define a critical role for HENMT1 and piRNAs in the maintenance of TE repression in adult germ cells and setting the spermatogenic program.

Piwi-interacting RNAs (piRNAs) are small non-coding RNAs found in great abundance within both embryonic and adult male germ cells. Within embryonic germ cells, piRNAs have a well-recognized role in transposable element (TE) silencing, however, their role in adult cells remains poorly defined. Here we demonstrate that HENMT1 dysfunction and the resultant piRNA instability dramatically impacts multiple aspects of adult germ cell biology. Specifically, pachytene piRNAs are required to maintain TE silencing in adult germ cells and to set the spermatogenic gene expression program. piRNA loss leads to a more active chromatin state in the regulatory regions of numerous normally haploid germ cell genes and their precocious expression during meiosis, followed by a catastrophic deregulation of gene expression in haploid cells and male sterility.

Mammalian genome evolution as a result of epigenetic regulation of transposable elements

Transposable elements (TEs) make up a large proportion of mammalian genomes and are a strong evolutionary force capable of rewiring regulatory networks and causing genome rearrangements. Additionally, there are many eukaryotic epigenetic defense mechanisms able to transcriptionally silence TEs. Furthermore, small RNA molecules that target TE DNA sequences often mediate these epigenetic defense mechanisms. As a result, epigenetic marks associated with TE silencing can be reestablished after epigenetic reprogramming - an event during the mammalian life cycle that results in widespread loss of parental epigenetic marks. Furthermore, targeted epigenetic marks associated with TE silencing may have an impact on nearby gene expression. Therefore, TEs may have driven species evolution via their ability to heritably alter the epigenetic regulation of gene expression in mammals.

Revealing Missing Human Protein Isoforms Based on Ab Initio Prediction, RNA-seq and Proteomics

Biological and biomedical research relies on comprehensive understanding of protein-coding transcripts. However, the total number of human proteins is still unknown due to the prevalence of alternative splicing. In this paper, we detected 31,566 novel transcripts with coding potential by filtering our ab initio predictions with 50 RNA-seq datasets from diverse tissues/cell lines. PCR followed by MiSeq sequencing showed that at least 84.1% of these predicted novel splice sites could be validated. In contrast to known transcripts, the expression of these novel transcripts were highly tissue-specific. Based on these novel transcripts, at least 36 novel proteins were detected from shotgun proteomics data of 41 breast samples. We also showed L1 retrotransposons have a more significant impact on the origin of new transcripts/genes than previously thought. Furthermore, we found that alternative splicing is extraordinarily widespread for genes involved in specific biological functions like protein binding, nucleoside binding, neuron projection, membrane organization and cell adhesion. In the end, the total number of human transcripts with protein-coding potential was estimated to be at least 204,950.

Interspersed repeats in the horse (Equus caballus); spatial correlations highlight conserved chromosomal domains

The interspersed repeat content of mammalian genomes has been best characterized in human, mouse and cow. In this study, we carried out de novo identification of repeated elements in the equine genome and identified previously unknown elements present at low copy number. The equine genome contains typical eutherian mammal repeats, but also has a significant number of hybrid repeats in addition to clade-specific Long Interspersed Nuclear Elements (LINE). Equus caballus clade specific LINE 1 (L1) repeats can be classified into approximately five subfamilies, three of which have undergone significant expansion. There are 1115 full-length copies of these equine L1, but of the 103 presumptive active copies, 93 fall within a single subfamily, indicating a rapid recent expansion of this subfamily. We also analysed both interspersed and simple sequence repeats (SSR) genome-wide, finding that some repeat classes are spatially correlated with each other as well as with G+C content and gene density. Based on these spatial correlations, we have confirmed that recently-described ancestral vs. clade-specific genome territories can be defined by their repeat content. The clade-specific Short Interspersed Nuclear Element correlations were scattered over the genome and appear to have been extensively remodelled. In contrast, territories enriched for ancestral repeats tended to be contiguous domains. To determine if the latter territories were evolutionarily conserved, we compared these results with a similar analysis of the human genome, and observed similar ancestral repeat enriched domains. These results indicate that ancestral, evolutionarily conserved mammalian genome territories can be identified on the basis of repeat content alone. Interspersed repeats of different ages appear to be analogous to geologic strata, allowing identification of ancient vs. newly remodelled regions of mammalian genomes.

MicroRNAs are part of the regulatory network that controls EGF induced apoptosis, including elements of the JAK/STAT pathway, in A431 cells

MiRNAs are known to regulate gene expression and in the context of cancer have been shown to regulate metastasis, cell proliferation and cell death. In this report we describe potential miRNA regulatory roles with respect to induction of cell death by pharmacologic dose of Epidermal Growth Factor (EGF). Our previous work suggested that multiple pathways are involved in the induction of apoptosis, including interferon induced genes, cytokines, cytoskeleton and cell adhesion and TP53 regulated genes. Using miRNA time course expression profiling of EGF treated A431 cells and coupling this to our previous gene expression and proteomic data, we have been able to implicate a number of additional miRNAs in the regulation of apoptosis. Specifically we have linked miR-134, miR-145, miR-146b-5p, miR-432 and miR-494 to the regulation of both apoptotic and anti-apoptotic genes expressed as a function of EGF treatment. Whilst additional miRNAs were differentially expressed, these had the largest number of apoptotic and anti-apoptotic targets. We found 5 miRNAs previously implicated in the regulation of apoptosis and our results indicate that an additional 20 miRNAs are likely to be involved based on their correlated expression with targets. Certain targets were linked to multiple miRNAs, including PEG10, BTG1, ID1, IL32 and NCF2. Some miRNAs that target the interferon pathway were found to be down regulated, consistent with a novel layer of regulation of interferon pathway components downstream of JAK/STAT. We have significantly expanded the repertoire of miRNAs that may regulate apoptosis in cancer cells as a result of this work.

Whole-exome sequencing points to considerable genetic heterogeneity of cerebral palsy

Cerebral palsy (CP) is a common, clinically heterogeneous group of disorders affecting movement and posture. Its prevalence has changed little in 50 years and the causes remain largely unknown. The genetic contribution to CP causation has been predicted to be ~2%. We performed whole-exome sequencing of 183 cases with CP including both parents (98 cases) or one parent (67 cases) and 18 singleton cases (no parental DNA). We identified and validated 61 de novo protein-altering variants in 43 out of 98 (44%) case-parent trios. Initial prioritization of variants for causality was by mutation type, whether they were known or predicted to be deleterious and whether they occurred in known disease genes whose clinical spectrum overlaps CP. Further, prioritization used two multidimensional frameworks-the Residual Variation Intolerance Score and the Combined Annotation-dependent Depletion score. Ten de novo mutations in three previously identified disease genes (TUBA1A (n=2), SCN8A (n=1) and KDM5C (n=1)) and in six novel candidate CP genes (AGAP1, JHDM1D, MAST1, NAA35, RFX2 and WIPI2) were predicted to be potentially pathogenic for CP. In addition, we identified four predicted pathogenic, hemizygous variants on chromosome X in two known disease genes, L1CAM and PAK3, and in two novel candidate CP genes, CD99L2 and TENM1. In total, 14% of CP cases, by strict criteria, had a potentially disease-causing gene variant. Half were in novel genes. The genetic heterogeneity highlights the complexity of the genetic contribution to CP. Function and pathway studies are required to establish the causative role of these putative pathogenic CP genes.

A functional genomics catalogue of activated transcription factors during pathogenesis of pneumococcal disease

Background

Streptococcus pneumoniae (the pneumococcus) is the world’s foremost microbial pathogen, killing more people each year than HIV, TB or malaria. The capacity to penetrate deeper host tissues contributes substantially to the ability of this organism to cause disease. Here we investigated, for the first time, functional genomics modulation of 3 pneumococcal strains (serotype 2 [D39], serotype 4 [WCH43] and serotype 6A [WCH16]) during transition from the nasopharynx to lungs to blood and to brain of mice at both promoter and domain activation levels.

Results

We found 7 highly activated transcription factors (TFs) [argR, codY, hup, rpoD, rr02, scrR and smrC] capable of binding to a large number of up-regulated genes, potentially constituting the regulatory backbone of pneumococcal pathogenesis. Strain D39 showed a distinct profile in employing a large number of TFs during blood infection. Interestingly, the same highly activated TFs used by D39 in blood are also used by WCH16 and WCH43 during brain infection. This indicates that different pneumococcal strains might activate a similar set of TFs and regulatory elements depending on the final site of infection. Hierarchical clustering analysis showed that all the highly activated TFs, except rpoD, clustered together with a high level of similarity in all 3 strains, which might suggest redundancy in the regulatory roles of these TFs during infection. Discriminant function analysis of the TFs in various niches highlights differential regulatory backgrounds of the 3 strains, and pathogenesis data confirms codY as the most significant predictor discriminating between these strains in various niches, particularly in the blood. Moreover, the predicted TF and domain activation profiles of the 3 strains correspond with their distinct pathogenicity characteristics.

Conclusions

Our findings suggest that the pneumococcus changes the short binding sites in the promoter regions of genes in a niche-specific manner to enhance its ability to disseminate from one host niche to another. This study provides a framework for an improved understanding of the dynamics of pneumococcal pathogenesis, and opens a new avenue into similar investigations in other pathogenic bacteria.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-769) contains supplementary material, which is available to authorized users.

ARMC5 mutations are common in familial bilateral macronodular adrenal hyperplasia

CONTEXT

Bilateral macronodular adrenal hyperplasia (BMAH) is a rare form of adrenal Cushing's syndrome. Familial cases have been reported, but at the time we conducted this study, the genetic basis of BMAH was unknown. Recently, germline variants of armadillo repeat containing 5 (ARMC5) in patients with isolated BMAH and somatic, second-hit mutations in tumor nodules, were identified.

OBJECTIVE

Our objective was to identify the genetic basis of familial BMAH.

DESIGN

We performed whole exome capture and sequencing of 2 affected individuals from each of 4 BMAH families (BMAH-01, BMAH-02, BMAH-03, and BMAH-05). Based on clinical evaluation, there were 7, 3, 3, and 4 affected individuals in these families, respectively. Sanger sequencing of ARMC5 was performed in 1 other BMAH kindred, BMAH-06.

RESULTS

Exome sequencing identified novel variants Chr16:g.31477540, c.2139delT, p.(Thr715Leufs*1) (BMAH-02) and Chr16:g.31473811, c.943C→T, p.(Arg315Trp) (BMAH-03) in ARMC5 (GRch37/hg19), validated by Sanger sequencing. BMAH-01 had a recently reported mutation Chr16:g.31476121, c.1777C→T, p.(Arg593Trp). Sanger sequencing of ARMC5 in BMAH-06 identified a previously reported mutation, Chr16:g. 31473688; c.799C→T, p.(Arg267*). The genetic basis of BMAH in BMAH-05 was not identified.

CONCLUSIONS

Our studies have detected ARMC5 mutations in 4 of 5 BMAH families tested, confirming that these mutations are a frequent cause of BMAH. Two of the 4 families had novel mutations, indicating allelic heterogeneity. Preclinical evaluation did not predict mutation status. The ARMC5-negative family had unusual prominent hyperaldosteronism. Further studies are needed to determine the penetrance of BMAH in ARMC5 mutation-positive relatives of affected patients, the practical utility of genetic screening and genotype-phenotype correlations.