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.