A conserved sequence signature is essential for robust plant miRNA biogenesis.

Micro (mi)RNAs are 20-22nt long non-coding RNA molecules involved in post-transcriptional silencing of targets having high base-pair complementarity. Plant miRNAs are processed from long Pol II-transcripts with specific stem-loop structures by Dicer-like (DCL) 1 protein. Although there were reports indicating how a specific region is selected for miRNA biogenesis, molecular details were unclear. Here, we show that the presence of specific GC-rich sequence signature within miRNA/miRNA* region is required for the precise miRNA biogenesis. The involvement of GC-rich signatures in precise processing and abundance of miRNAs was confirmed through detailed molecular and functional analysis. Consistent with the presence of the miRNA-specific GC signature, target RNAs of miRNAs also possess conserved complementary sequence signatures in their miRNA binding motifs. The selection of these GC signatures was dependent on an RNA binding protein partner of DCL1 named HYL1. Finally, we demonstrate a direct application of this discovery for enhancing the abundance and efficiency of artificial miRNAs that are popular in plant functional genomic studies.

The emergence of sequence-dependent structural motifs in stretched, torsionally constrained DNA.

The double-helical structure of DNA results from canonical base pairing and stacking interactions. However, variations from steady-state conformations resulting from mechanical perturbations in cells have physiological relevance but their dependence on sequence remains unclear. Here, we use molecular dynamics simulations showing sequence differences result in markedly different structural motifs upon physiological twisting and stretching. We simulate overextension on different sequences of DNA ((AA)12, (AT)12, (CC)12 and (CG)12) with supercoiling densities at 200 and 50 mM salt concentrations. We find that DNA denatures in the majority of stretching simulations, surprisingly including those with over-twisted DNA. GC-rich sequences are observed to be more stable than AT-rich ones, with the specific response dependent on the base pair order. Furthermore, we find that (AT)12 forms stable periodic structures with non-canonical hydrogen bonds in some regions and non-canonical stacking in others, whereas (CG)12 forms a stacking motif of four base pairs independent of supercoiling density. Our results demonstrate that 20-30% DNA extension is sufficient for breaking B-DNA around and significantly above cellular supercoiling, and that the DNA sequence is crucial for understanding structural changes under mechanical stress. Our findings have important implications for the activities of protein machinery interacting with DNA in all cells.

New Roles for Canonical Transcription Factors in Repeat Expansion Diseases.

The presence of microsatellite repeat expansions within genes is associated with >30 neurological diseases. Of interest, (GGGGCC)>30-repeats within C9orf72 are associated with amyotrophic lateral sclerosis and frontotemporal dementia (ALS/FTD). These expansions can be 100s to 1000s of units long. Thus, it is perplexing how RNA-polymerase II (RNAPII) can successfully transcribe them. Recent investigations focusing on GGGGCC-transcription have identified specific, canonical complexes that may promote RNAPII-transcription at these GC-rich microsatellites: the DSIF complex and PAF1C. These complexes may be important for resolving the unique secondary structures formed by GGGGCC-DNA during transcription. Importantly, this process can produce potentially toxic repeat-containing RNA that can encode potentially toxic peptides, impacting neuron function and health. Understanding how transcription of these repeats occurs has implications for therapeutics in multiple diseases.

Transcription elongation factor AFF2/FMR2 regulates expression of expanded GGGGCC repeat-containing C9ORF72 allele in ALS/FTD.

Expanded GGGGCC (G4C2) repeats in C9ORF72 cause amyotrophic lateral sclerosis (ALS) and frontotemporal dementia (FTD). How RNAs containing expanded G4C2 repeats are transcribed in human neurons is largely unknown. Here we describe a Drosophila model in which poly(GR) expression in adult neurons causes axonal and locomotor defects and premature death without apparent TDP-43 pathology. In an unbiased genetic screen, partial loss of Lilliputian (Lilli) activity strongly suppresses poly(GR) toxicity by specifically downregulating the transcription of GC-rich sequences in Drosophila. Knockout of AFF2/FMR2 (one of four mammalian homologues of Lilli) with CRISPR-Cas9 decreases the expression of the mutant C9ORF72 allele containing expanded G4C2 repeats and the levels of repeat RNA foci and dipeptide repeat proteins in cortical neurons derived from induced pluripotent stem cells of C9ORF72 patients, resulting in rescue of axonal degeneration and TDP-43 pathology. Thus, AFF2/FMR2 regulates the transcription and toxicity of expanded G4C2 repeats in human C9ORF72-ALS/FTD neurons.

The first mitochondrial genome from Scopuridae (Insecta: Plecoptera) reveals structural features and phylogenetic implications.

This study determined the first complete mitochondrial genome (mitogenome) of a stonefly, Scopura longa (Plecoptera: Scopuridae), and reconstructed a phylogeny based on two datasets of mitogenomes in eighteen available stoneflies to examine the relationships among Plecoptera. The complete mitogenome of S. longa is a circular molecule of 15,798bp in size. It contains 13 protein-coding genes (PCGs), 22 transfer RNA genes (tRNAs), two ribosomal RNA genes (rRNAs) and a control region (CR). Most PCGs used standard ATN start codons and TAN stop codons. All tRNAs could be folded as typical cloverleaf secondary structures except tRNASer(AGN), with the dihydrouridine (DHU) arm missing. Tandem repeat regions, two potential stem-loop (SL) structures, Poly T structures and G+C-rich regions are detected in the control region. Finally, the phylogenetic relationships among the families within the Arctoperlaria were reconstructed. The topological structures of the two trees were almost identical. The present phylogenetic analysis shows that S. longa belongs to the infraorder Euholognatha. The monophyly of each family is generally well supported based on nucleotide sequences. The Pteronarcyidae is sister to ((Peltoperlidae+Styloperlidae)+(Perlidae+(Perlodidae+Chloroperlidae))).

Repeat elements organise 3D genome structure and mediate transcription in the filamentous fungus Epichloë festucae.

Structural features of genomes, including the three-dimensional arrangement of DNA in the nucleus, are increasingly seen as key contributors to the regulation of gene expression. However, studies on how genome structure and nuclear organisation influence transcription have so far been limited to a handful of model species. This narrow focus limits our ability to draw general conclusions about the ways in which three-dimensional structures are encoded, and to integrate information from three-dimensional data to address a broader gamut of biological questions. Here, we generate a complete and gapless genome sequence for the filamentous fungus, Epichloë festucae. We use Hi-C data to examine the three-dimensional organisation of the genome, and RNA-seq data to investigate how Epichloë genome structure contributes to the suite of transcriptional changes needed to maintain symbiotic relationships with the grass host. Our results reveal a genome in which very repeat-rich blocks of DNA with discrete boundaries are interspersed by gene-rich sequences that are almost repeat-free. In contrast to other species reported to date, the three-dimensional structure of the genome is anchored by these repeat blocks, which act to isolate transcription in neighbouring gene-rich regions. Genes that are differentially expressed in planta are enriched near the boundaries of these repeat-rich blocks, suggesting that their three-dimensional orientation partly encodes and regulates the symbiotic relationship formed by this organism.

Developmental and cancer-associated plasticity of DNA replication preferentially targets GC-poor, lowly expressed and late-replicating regions.

The spatiotemporal program of metazoan DNA replication is regulated during development and altered in cancers. We have generated novel OK-seq, Repli-seq and RNA-seq data to compare the DNA replication and gene expression programs of twelve cancer and non-cancer human cell types. Changes in replication fork directionality (RFD) determined by OK-seq are widespread but more frequent within GC-poor isochores and largely disconnected from transcription changes. Cancer cell RFD profiles cluster with non-cancer cells of similar developmental origin but not with different cancer types. Importantly, recurrent RFD changes are detected in specific tumour progression pathways. Using a model for establishment and early progression of chronic myeloid leukemia (CML), we identify 1027 replication initiation zones (IZs) that progressively change efficiency during long-term expression of the BCR-ABL1 oncogene, being twice more often downregulated than upregulated. Prolonged expression of BCR-ABL1 results in targeting of new IZs and accentuation of previous efficiency changes. Targeted IZs are predominantly located in GC-poor, late replicating gene deserts and frequently silenced in late CML. Prolonged expression of BCR-ABL1 results in massive deletion of GC-poor, late replicating DNA sequences enriched in origin silencing events. We conclude that BCR-ABL1 expression progressively affects replication and stability of GC-poor, late-replicating regions during CML progression.

Characterization of Chicken Tumor Necrosis Factor-α, a Long Missed Cytokine in Birds.

Tumor necrosis factor-α (TNF-α) is a pleiotropic cytokine playing critical roles in host defense and acute and chronic inflammation. It has been described in fish, amphibians, and mammals but was considered to be absent in the avian genomes. Here, we report on the identification and functional characterization of the avian ortholog. The chicken TNF-α (chTNF-α) is encoded by a highly GC-rich gene, whose product shares with its mammalian counterpart 45% homology in the extracellular part displaying the characteristic TNF homology domain. Orthologs of chTNF-α were identified in the genomes of 12 additional avian species including Palaeognathae and Neognathae, and the synteny of the closely adjacent loci with mammalian TNF-α orthologs was demonstrated in the crow (Corvus cornix) genome. In addition to chTNF-α, we obtained full sequences for homologs of TNF-α receptors 1 and 2 (TNFR1, TNFR2). chTNF-α mRNA is strongly induced by lipopolysaccharide (LPS) stimulation of monocyte derived, splenic and bone marrow macrophages, and significantly upregulated in splenic tissue in response to i.v. LPS treatment. Activation of T-lymphocytes by TCR crosslinking induces chTNF-α expression in CD4+ but not in CD8+ cells. To gain insights into its biological activity, we generated recombinant chTNF-α in eukaryotic and prokaryotic expression systems. Both, the full-length cytokine and the extracellular domain rapidly induced an NFκB-luciferase reporter in stably transfected CEC-32 reporter cells. Collectively, these data provide strong evidence for the existence of a fully functional TNF-α/TNF-α receptor system in birds thus filling a gap in our understanding of the evolution of cytokine systems.

Refining the Ambush Hypothesis: Evidence That GC- and AT-Rich Bacteria Employ Different Frameshift Defence Strategies.

Stop codons are frequently selected for beyond their regular termination function for error control. The "ambush hypothesis" proposes out-of-frame stop codons (OSCs) terminating frameshifted translations are selected for. Although early indirect evidence was partially supportive, recent evidence suggests OSC frequencies are not exceptional when considering underlying nucleotide content. However, prior null tests fail to control amino acid/codon usages or possible local mutational biases. We therefore return to the issue using bacterial genomes, considering several tests defining and testing against a null. We employ simulation approaches preserving amino acid order but shuffling synonymous codons or preserving codons while shuffling amino acid order. Additionally, we compare codon usage in amino acid pairs, where one codon can but the next, otherwise identical codon, cannot encode an OSC. OSC frequencies exceed expectations typically in AT-rich genomes, the +1 frame and for TGA/TAA but not TAG. With this complex evidence, simply rejecting or accepting the ambush hypothesis is not warranted. We propose a refined post hoc model, whereby AT-rich genomes have more accidental frameshifts, handled by RF2-RF3 complexes (associated with TGA/TAA) and are mostly +1 (or -2) slips. Supporting this, excesses positively correlate with in silico predicted frameshift probabilities. Thus, we propose a more viable framework, whereby genomes broadly adopt one of the two strategies to combat frameshifts: preventing frameshifting (GC-rich) or permitting frameshifts but minimizing impacts when most are caught early (AT-rich). Our refined framework holds promise yet some features, such as the bias of out-of-frame sense codons, remain unexplained.

CG dinucleotide suppression enables antiviral defence targeting non-self RNA.

Vertebrate genomes exhibit marked CG suppression-that is, lower than expected numbers of 5'-CG-3' dinucleotides. This feature is likely to be due to C-to-T mutations that have accumulated over hundreds of millions of years, driven by CG-specific DNA methyl transferases and spontaneous methyl-cytosine deamination. Many RNA viruses of vertebrates that are not substrates for DNA methyl transferases mimic the CG suppression of their hosts. This property of viral genomes is unexplained. Here we show, using synonymous mutagenesis, that CG suppression is essential for HIV-1 replication. The deleterious effect of CG dinucleotides on HIV-1 replication was cumulative, associated with cytoplasmic RNA depletion, and was exerted by CG dinucleotides in both translated and non-translated exonic RNA sequences. A focused screen using small inhibitory RNAs revealed that zinc-finger antiviral protein (ZAP) inhibited virion production by cells infected with CG-enriched HIV-1. Crucially, HIV-1 mutants containing segments whose CG content mimicked random nucleotide sequence were defective in unmanipulated cells, but replicated normally in ZAP-deficient cells. Crosslinking-immunoprecipitation-sequencing assays demonstrated that ZAP binds directly and selectively to RNA sequences containing CG dinucleotides. These findings suggest that ZAP exploits host CG suppression to identify non-self RNA. The dinucleotide composition of HIV-1, and perhaps other RNA viruses, appears to have adapted to evade this host defence.

Avian Genomes Revisited: Hidden Genes Uncovered and the Rates versus Traits Paradox in Birds.

According to current assemblies, avian genomes differ from those of the other lineages of amniotes in 1) containing a lower number of genes; 2) displaying a high stability of karyotype and recombination map; and 3) lacking any correlation between evolutionary rates (dN/dS) and life-history traits, unlike mammals and nonavian reptiles. We question the reality of the bird missing genes and investigate whether insufficient representation of bird gene content might have biased previous evolutionary analyses. Mining RNAseq data, we show that the vast majority of the genes missing from avian genome assemblies are actually present in most species of birds. These mainly correspond to the GC-rich fraction of the bird genome, which is the most difficult to sequence, assemble and annotate. With the inclusion of these genes in a phylogenomic analysis of high-quality alignments, we uncover a positive and significant correlation between the ratio of nonsynonymous to synonymous substitution rate (dN/dS) and life-history traits in Neoaves. We report a strong effect of GC-biased gene conversion on the dN/dS ratio in birds and a peculiar behavior of Palaeognathae (ostrich and allies) and Galloanserae (chickens, ducks and allies). Avian genomes do not contain fewer genes than mammals or nonavian reptiles. Previous analyses have overlooked ∼15% of the bird gene complement. GC-rich regions, which are the most difficult to access, are a key component of amniote genomes. They experience peculiar molecular processes and must be included for unbiased functional and comparative genomic analyses in birds.

DNA helix: the importance of being AT-rich.

The AT-rich DNA is mostly associated with condensed chromatin, whereas the GC-rich sequence is preferably located in the dispersed chromatin. The AT-rich genes are prone to be tissue-specific (silenced in most tissues), while the GC-rich genes tend to be housekeeping (expressed in many tissues). This paper reports another important property of DNA base composition, which can affect repertoire of genes with high AT content. The GC-rich sequence is more liable to mutation. We found that Spearman correlation between human gene GC content and mutation probability is above 0.9. The change of base composition even in synonymous sites affects mutation probability of nonsynonymous sites and thus of encoded proteins. There is a unique type of housekeeping genes, which are especially unsafe when prone to mutation. Natural selection which usually removes deleterious mutations, in the case of these genes only increases the hazard because it can descend to suborganismal (cellular) level. These are cell cycle-related genes. In accordance with the proposed concept, they have low GC content of synonymous sites (despite them being housekeeping). The gene-centred protein interaction enrichment analysis (PIEA) showed the core clusters of genes whose interactants are modularly enriched in genes with AT-rich synonymous codons. This interconnected network is involved in double-strand break repair, DNA integrity checkpoints and chromosome pairing at mitosis. The damage of these genes results in genome and chromosome instability leading to cancer and other 'error catastrophes'. Reducing the nonsynonymous mutations, the usage of AT-rich synonymous codons can decrease probability of cancer by above 20-fold.

Short DNA sequence patterns accurately identify broadly active human enhancers.

BACKGROUND: Enhancers are DNA regulatory elements that influence gene expression. There is substantial diversity in enhancers' activity patterns: some enhancers drive expression in a single cellular context, while others are active across many. Sequence characteristics, such as transcription factor (TF) binding motifs, influence the activity patterns of regulatory sequences; however, the regulatory logic through which specific sequences drive enhancer activity patterns is poorly understood. Recent analysis of Drosophila enhancers suggested that short dinucleotide repeat motifs (DRMs) are general enhancer sequence features that drive broad regulatory activity. However, it is not known whether the regulatory role of DRMs is conserved across species. RESULTS: We performed a comprehensive analysis of the relationship between short DNA sequence patterns, including DRMs, and human enhancer activity in 38,538 enhancers across 411 different contexts. In a machine-learning framework, the occurrence patterns of short sequence motifs accurately predicted broadly active human enhancers. However, DRMs alone were weakly predictive of broad enhancer activity in humans and showed different enrichment patterns than in Drosophila. In general, GC-rich sequence motifs were significantly associated with broad enhancer activity, and consistent with this enrichment, broadly active human TFs recognize GC-rich motifs. CONCLUSIONS: Our results reveal the importance of specific sequence motifs in broadly active human enhancers, demonstrate the lack of evolutionary conservation of the role of DRMs, and provide a computational framework for investigating the logic of enhancer sequences.

Synthetic ligand promotes gene expression by affecting GC sequence in promoter.

A naphthyridine carbamate tetramer (NCT8) is a synthetic compound, which selectively binds to nucleic acids containing CGG/CGG sequence. Although NCT8 is a promising compound for a wide range of DNA and RNA based biotechnology such as modulation of specific gene expression, little is known about its behavior in human cells. In the present study, we investigated the changes induced in gene expression by NCT8. Genes differentially expressed in the presence of NCT8 in HeLa cells were identified by whole-transcriptome analysis. The whole-transcriptome analysis showed that NCT8 significantly induced up-regulation of specific genes, whose promoter region has GC-rich sequence.

Retrotransposons are specified as DNA replication origins in the gene-poor regions of Arabidopsis heterochromatin.

Genomic stability depends on faithful genome replication. This is achieved by the concerted activity of thousands of DNA replication origins (ORIs) scattered throughout the genome. The DNA and chromatin features determining ORI specification are not presently known. We have generated a high-resolution genome-wide map of 3230 ORIs in cultured Arabidopsis thaliana cells. Here, we focused on defining the features associated with ORIs in heterochromatin. In pericentromeric gene-poor domains ORIs associate almost exclusively with the retrotransposon class of transposable elements (TEs), in particular of the Gypsy family. ORI activity in retrotransposons occurs independently of TE expression and while maintaining high levels of H3K9me2 and H3K27me1, typical marks of repressed heterochromatin. ORI-TEs largely colocalize with chromatin signatures defining GC-rich heterochromatin. Importantly, TEs with active ORIs contain a local GC content higher than the TEs lacking them. Our results lead us to conclude that ORI colocalization with retrotransposons is determined by their transposition mechanism based on transcription, and a specific chromatin landscape. Our detailed analysis of ORIs responsible for heterochromatin replication has implications on the mechanisms of ORI specification in other multicellular organisms in which retrotransposons are major components of heterochromatin and of the entire genome.

Transcriptome Analysis of Core Dinoflagellates Reveals a Universal Bias towards "GC" Rich Codons.

Although dinoflagellates are a potential source of pharmaceuticals and natural products, the mechanisms for regulating and producing these compounds are largely unknown because of extensive post-transcriptional control of gene expression. One well-documented mechanism for controlling gene expression during translation is codon bias, whereby specific codons slow or even terminate protein synthesis. Approximately 10,000 annotatable genes from fifteen "core" dinoflagellate transcriptomes along a range of overall guanine and cytosine (GC) content were used for codonW analysis to determine the relative synonymous codon usage (RSCU) and the GC content at each codon position. GC bias in the analyzed dataset and at the third codon position varied from 51% and 54% to 66% and 88%, respectively. Codons poor in GC were observed to be universally absent, but bias was most pronounced for codons ending in uracil followed by adenine (UA). GC bias at the third codon position was able to explain low abundance codons as well as the low effective number of codons. Thus, we propose that a bias towards codons rich in GC bases is a universal feature of core dinoflagellates, possibly relating to their unique chromosome structure, and not likely a major mechanism for controlling gene expression.

PCR amplification of GC-rich DNA regions using the nucleotide analog N4-methyl-2'-deoxycytidine 5'-triphosphate.

GC-rich DNA regions were PCR-amplified with Taq DNA polymerase using either the canonical set of deoxynucleoside triphosphates or mixtures in which the dCTP had been partially or completely replaced by its N4-methylated analog, N4-methyl-2'-deoxycytidine 5'-triphosphate (N4me-dCTP). In the case of a particularly GC-rich region (78.9% GC), the PCR mixtures containing N4me-dCTP produced the expected amplicon in high yield, while mixtures containing the canonical set of nucleotides produced numerous alternative amplicons. For another GC-rich DNA region (80.6% GC), the target amplicon was only generated by re-amplifying a gel-purified sample of the original amplicon with N4me-dCTP-containing PCR mixtures. In a direct PCR comparison on a highly GC-rich template, mixtures containing N4me-dCTP clearly performed better than did solutions containing the canonical set of nucleotides mixed with various organic additives (DMSO, betaine, or ethylene glycol) that have been reported to resolve or alleviate problems caused by secondary structures in the DNA. This nucleotide analog was also tested in PCR amplification of DNA regions with intermediate GC content, producing the expected amplicon in each case with a melting temperature (Tm) clearly below the Tm of the same amplicon synthesized exclusively with the canonical bases.

Human CST Facilitates Genome-wide RAD51 Recruitment to GC-Rich Repetitive Sequences in Response to Replication Stress.

The telomeric CTC1/STN1/TEN1 (CST) complex has been implicated in promoting replication recovery under replication stress at genomic regions, yet its precise role is unclear. Here, we report that STN1 is enriched at GC-rich repetitive sequences genome-wide in response to hydroxyurea (HU)-induced replication stress. STN1 deficiency exacerbates the fragility of these sequences under replication stress, resulting in chromosome fragmentation. We find that upon fork stalling, CST proteins form distinct nuclear foci that colocalize with RAD51. Furthermore, replication stress induces physical association of CST with RAD51 in an ATR-dependent manner. Strikingly, CST deficiency diminishes HU-induced RAD51 foci formation and reduces RAD51 recruitment to telomeres and non-telomeric GC-rich fragile sequences. Collectively, our findings establish that CST promotes RAD51 recruitment to GC-rich repetitive sequences in response to replication stress to facilitate replication restart, thereby providing insights into the mechanism underlying genome stability maintenance.

Guanine-rich sequences inhibit proofreading DNA polymerases.

DNA polymerases with proofreading activity are important for accurate amplification of target DNA. Despite numerous efforts have been made to improve the proofreading DNA polymerases, they are more susceptible to be failed in PCR than non-proofreading DNA polymerases. Here we showed that proofreading DNA polymerases can be inhibited by certain primers. Further analysis showed that G-rich sequences such as GGGGG and GGGGHGG can cause PCR failure using proofreading DNA polymerases but not Taq DNA polymerase. The inhibitory effect of these G-rich sequences is caused by G-quadruplex and is dose dependent. G-rich inhibitory sequence-containing primers can be used in PCR at a lower concentration to amplify its target DNA fragment.

Chromosome Banding in Amphibia. XXXIII. Demonstration of 5-Methylcytosine-Rich Heterochromatin in Anura.

An experimental approach using monoclonal anti-5-methylcytosine (5-MeC) antibodies and indirect immunofluorescence was elaborated for detecting 5-MeC-rich chromosome regions in anuran chromosomes. This technique was applied to mitotic metaphases of 6 neotropical frog species belonging to 6 genera and 4 families. The hypermethylation patterns were compared with a variety of banding patterns obtained by conventional banding techniques. The hypermethylated DNA sequences are species-specific and located exclusively in constitutive heterochromatin. They are found in centromeric, pericentromeric, telomeric, and interstitial positions of the chromosomes and adjacent to nucleolus organizer regions. 5-MeC-rich DNA sequences can be embedded both in AT- and GC-rich repetitive DNA. The experimental parameters that have major influence on the reproducibility and quality of the anti-5-MeC antibody labeling are discussed.