The dinucleotide composition of the Zika virus genome is shaped by conflicting evolutionary pressures in mammalian hosts and mosquito vectors.

Most vertebrate RNA viruses show pervasive suppression of CpG and UpA dinucleotides, closely resembling the dinucleotide composition of host cell transcriptomes. In contrast, CpG suppression is absent in both invertebrate mRNA and RNA viruses that exclusively infect arthropods. Arthropod-borne (arbo) viruses are transmitted between vertebrate hosts by invertebrate vectors and thus encounter potentially conflicting evolutionary pressures in the different cytoplasmic environments. Using a newly developed Zika virus (ZIKV) model, we have investigated how demands for CpG suppression in vertebrate cells can be reconciled with potentially quite different compositional requirements in invertebrates and how this affects ZIKV replication and transmission. Mutant viruses with synonymously elevated CpG or UpA dinucleotide frequencies showed attenuated replication in vertebrate cell lines, which was rescued by knockout of the zinc-finger antiviral protein (ZAP). Conversely, in mosquito cells, ZIKV mutants with elevated CpG dinucleotide frequencies showed substantially enhanced replication compared to wild type. Host-driven effects on virus replication attenuation and enhancement were even more apparent in mouse and mosquito models. Infections with CpG- or UpA-high ZIKV mutants in mice did not cause typical ZIKV-induced tissue damage and completely protected mice during subsequent challenge with wild-type virus, which demonstrates their potential as live-attenuated vaccines. In contrast, the CpG-high mutants displayed enhanced replication in Aedes aegypti mosquitoes and a larger proportion of mosquitoes carried infectious virus in their saliva. These findings show that mosquito cells are also capable of discriminating RNA based on dinucleotide composition. However, the evolutionary pressure on the CpG dinucleotides of viral genomes in arthropod vectors directly opposes the pressure present in vertebrate host cells, which provides evidence that an adaptive compromise is required for arbovirus transmission. This suggests that the genome composition of arbo flaviviruses is crucial to maintain the balance between high-level replication in the vertebrate host and persistent replication in the mosquito vector.

Human exonization through differential nucleosome occupancy.

Nucleosomal modifications have been implicated in fundamental epigenetic regulation, but the roles of nucleosome occupancy in shaping changes through evolution remain to be addressed. Here we present high-resolution nucleosome occupancy profiles for multiple tissues derived from human, macaque, tree shrew, mouse, and pig. Genome-wide comparison reveals conserved nucleosome occupancy profiles across both different species and tissue types. Notably, we found significantly higher levels of nucleosome occupancy in exons than in introns, a pattern correlated with the different exon-intron GC content. We then determined whether this biased occupancy may play roles in the origination of new exons through evolution, rather than being a downstream effect of exonization, through a comparative approach to sequentially trace the order of the exonization and biased nucleosome binding. By identifying recently evolved exons in human but not in macaque using matched RNA sequencing, we found that higher exonic nucleosome occupancy also existed in macaque regions orthologous to these exons. Presumably, such biased nucleosome occupancy facilitates the origination of new exons by increasing the splice strength of the ancestral nonexonic regions through driving a local difference in GC content. These data thus support a model that sites bound by nucleosomes are more likely to evolve into exons, which we term the "nucleosome-first" model.

Lacisediminimonas profundi gen. nov., sp. nov., a member of the family Oxalobacteraceae isolated from freshwater sediment.

A novel Gram-stain-negative bacterial strain, CHu64-6-4T, was isolated from a 67-cm-long sediment core collected from the Daechung Reservoir at a water depth of 17 m, Daejeon, Republic of Korea. The cells of strain CHu64-6-4T were aerobic nonmotile and formed colorless colonies on R2A agar. The phylogenetic analysis based on 16S rRNA gene sequencing indicated that the strain formed a separate lineage within the family Oxalobacteraceae, exhibiting 97.2% and 97.1% 16S rRNA gene sequence similarities to Glaciimonas singularis and Paraherbaspirillum soli, respectively. Strain CHu64-6-4T showed less than 74.4% average nucleotide identity compared to the type strains of related genera within the family Oxalobacteraceae. In the UPGMA dendrogram based on the ANI values of genomic sequences, strain CHu64-6-4T formed an evolutionary lineage independent of the genera Glaciimonas and some other taxa. The chemotaxonomic results showed Q-8 as the predominant respiratory ubiquinone, phosphatidylglycerol, diphosphatidylglycerol, and phosphatidylethnolamine as the major polar lipids, Summed Feature 3 (C16:1ω7c and/or iso-C15:0 2-OH), C16:0, and C18:1ω7c as the major fatty acids, and a DNA G+C content of 62.1 mol%. The combined genotypic and phenotypic data showed that strain CHu64-6-4T could be distinguished from all genera within the family Oxalobacteraceae and represents a novel genus, Lacisediminimonas profundi gen. nov., with the name Lacisediminimonas profundi sp. nov., in the family Oxalobacteraceae. The type strain is CHu64-6-4T (=KCTC 62287T=JCM 32676T).

Paenibacillus pinistramenti sp. nov., isolated from pine litter.

A Gram-stain-positive bacterium, designated strain ASL46T, was isolated from litter layer of a pine forest located in Anmyondo, Korea. Strain ASL46T was found to be an aerobic, motile, endospore-forming rod which can grow at 20-45 °C (optimum, 37 °C), at pH 6.0-11.0 (optimum, pH 7.0) and at salinities of 0-2% (w/v) NaCl (optimum, 1% NaCl). Phylogenetic analyses based on 16S rRNA gene sequences revealed that strain ASL46T belongs to the genus Paenibacillus, showing highest sequence similarity to P. yonginensis DCY84T (98.3%), P. physcomitrella XBT (97.4%) and P. faecis CIP 101062T (96.6%). The average nucleotide identity (ANI) and DNA-DNA relatedness between the strain ASL46T and P. physcomitrella XBT and P. yonginensis DCY84T yielded ANI values of 84.6 and 84.5% and DNA-DNA relatedness of 11.7 ± 0.7 and 10.9 ± 0.2%, respectively. The DNA G+C content of the genomic DNA of strain ASL46T was 52.1 mol%. The predominant isoprenoid quinone was identified as menaquinone-7 and the major cellular fatty acids were determined to be anteiso-C15:0, C16:0 and iso-C16:0. The major polar lipids were diphosphatidylglycerol, phosphatidylglycerol, phosphatidylethanolamine, five unidentified aminophospholipids, an unidentified phospholipid and an unidentified glycolipid. The whole-cell sugar was found to be ribose and cell wall peptidoglycan contained meso-diaminopimelic acid. On the basis of phylogenetic analyses, and phenotypic and chemotaxonomic characteristics, strain ASL46T represents a novel species of the genus Paenibacillus, for which the name Paenibacillus pinistramenti sp. nov. is proposed. The type strain is ASL46T (= KACC 18701T = NBRC 111876T).

Nocardia panacis sp. nov., a novel actinomycete with antiphytopathogen activity isolated from the rhizosphere of Panax notoginseng.

Strain YIM PH21724T was isolated from the rhizosphere of Panax notoginseng. Phylogenetic analyses based on 16S rRNA gene sequences showed that the strain exhibits close phylogenetic relatedness to Nocardia kroppenstedtii N1286T (97.70%), Nocardia farcinica NCTC 11134T (97.67%) and Nocardia puris DSM 44599T (97.40%). The menaquinones were identified as MK-9 (H4), MK-8 (H4, ω-cyclo) and MK-8 (H4), and the major fatty acids (> 10%) were identified as C16:0, C18:1 ω9c and C18:0 10-methyl. The polar lipids were found to be composed of diphosphatidylglycerol, phosphatidylethanolamine, phosphatidylinositol, phosphatidylinositol mannosides and an unidentified lipid. The G + C content of the genomic DNA was determined to be 67.01 mol%. The phenotypic, chemotaxonomic, phylogenetic and genomic results clearly show strain YIM PH21724T should be classified in the genus Nocardia and represents a novel species, for which the name Nocardia panacis sp. nov. is proposed. The type strain is YIM PH21724T (= DSM 105904T = KCTC 49030T = CCTCC AA 2017043T).

How bacterial xenogeneic silencer rok distinguishes foreign from self DNA in its resident genome.

Bacterial xenogeneic silencers play important roles in bacterial evolution by recognizing and inhibiting expression from foreign genes acquired through horizontal gene transfer, thereby buffering against potential fitness consequences of their misregulated expression. Here, the detailed DNA binding properties of Rok, a xenogeneic silencer in Bacillus subtilis, was studied using protein binding microarray, and the solution structure of its C-terminal DNA binding domain was determined in complex with DNA. The C-terminal domain of Rok adopts a typical winged helix fold, with a novel DNA recognition mechanism different from other winged helix proteins or xenogeneic silencers. Rok binds the DNA minor groove by forming hydrogen bonds to bases through N154, T156 at the N-terminal of α3 helix and R174 of wing W1, assisted by four lysine residues interacting electrostatically with DNA backbone phosphate groups. These structural features endow Rok with preference towards DNA sequences harboring AACTA, TACTA, and flexible multiple TpA steps, while rigid A-tracts are disfavored. Correspondingly, the Bacillus genomes containing Rok are rich in A-tracts and show a dramatic underrepresentation of AACTA and TACTA, which are significantly enriched in Rok binding regions. These observations suggest that the xenogeneic silencing protein and its resident genome may have evolved cooperatively.

HIV-1 tolerates changes in A-count in a small segment of the pol gene.

BACKGROUND: The HIV-1 RNA genome has a biased nucleotide composition with a surplus of As. Several hypotheses have been put forward to explain this striking phenomenon, but the A-count of the HIV-1 genome has thus far not been systematically manipulated. The reason for this reservation is the likelihood that known and unknown sequence motifs will be affected by such a massive mutational approach, thus resulting in replication-impaired virus mutants. We present the first attempt to increase and decrease the A-count in a relatively small polymerase (pol) gene segment of HIV-1 RNA. RESULTS: To minimize the mutational impact, a new mutational approach was developed that is inspired by natural sequence variation as present in HIV-1 isolates. This phylogeny-instructed mutagenesis allowed us to create replication-competent HIV-1 mutants with a significantly increased or decreased local A-count. The local A-count of the wild-type (wt) virus (40.2%) was further increased to 46.9% or reduced to 31.7 and 26.3%. These HIV-1 variants replicate efficiently in vitro, despite the fact that the pol changes cause a quite profound move in HIV-SIV sequence space. CONCLUSIONS: Extrapolating these results to the complete 9 kb RNA genome, we may cautiously suggest that the A-rich signature does not have to be maintained. This survey also provided clues that silent codon changes, in particular from G-to-A, determine the subtype-specific sequence signatures.

The mutational spectrum of non-CpG DNA varies with CpG content.

The accumulation of base substitutions (mutations) not subject to natural selection is the neutral mutation rate. Because this rate reflects the in vivo processes involved in maintaining the integrity of genetic information, the factors that affect the neutral mutation rate are of considerable interest. Mammals exhibit two dramatically different neutral mutation rates: the CpG mutation rate, wherein the C of most CpGs (i.e., methyl-CpG) mutate at 10-50 times that of C in any other context or of any other base. The latter mutations constitute the non-CpG rate. The high CpG rate results from the spontaneous deamination of methyl-C to T and incomplete restoration of the ensuing T:G mismatches to C:Gs. Here, we determined the neutral non-CpG mutation rate as a function of CpG content by comparing sequence divergence of thousands of pairs of neutrally evolving chimpanzee and human orthologs that differ primarily in CpG content. Both the mutation rate and the mutational spectrum (transition/transversion ratio) of non-CpG residues change in parallel as sigmoidal (logistic) functions of CpG content. As different mechanisms generate transitions and transversions, these results indicate that both mutation rate and mutational processes are contingent on the local CpG content. We consider several possible mechanisms that might explain how CpG exerts these effects.

Sequence features that drive human promoter function and tissue specificity.

Promoters are important regulatory elements that contain the necessary sequence features for cells to initiate transcription. To functionally characterize a large set of human promoters, we measured the transcriptional activities of 4575 putative promoters across eight cell lines using transient transfection reporter assays. In parallel, we measured gene expression in the same cell lines and observed a significant correlation between promoter activity and endogenous gene expression (r = 0.43). As transient transfection assays directly measure the promoting effect of a defined fragment of DNA sequence, decoupled from epigenetic, chromatin, or long-range regulatory effects, we sought to predict whether a promoter was active using sequence features alone. CG dinucleotide content was highly predictive of ubiquitous promoter activity, necessitating the separation of promoters into two groups: high CG promoters, mostly ubiquitously active, and low CG promoters, mostly cell line-specific. Computational models trained on the binding potential of transcriptional factor (TF) binding motifs could predict promoter activities in both high and low CG groups: average area under the receiver operating characteristic curve (AUC) of the models was 91% and exceeded the AUC of CG content by an average of 23%. Known relationships, for example, between HNF4A and hepatocytes, were recapitulated in the corresponding cell lines, in this case the liver-derived cell line HepG2. Half of the associations between tissue-specific TFs and cell line-specific promoters were new. Our study underscores the importance of collecting functional information from complementary assays and conditions to understand biology in a systematic framework.

The influence of the local sequence environment on RNA loop structures.

RNA folding is assumed to be a hierarchical process. The secondary structure of an RNA molecule, signified by base-pairing and stacking interactions between the paired bases, is formed first. Subsequently, the RNA molecule adopts an energetically favorable three-dimensional conformation in the structural space determined mainly by the rotational degrees of freedom associated with the backbone of regions of unpaired nucleotides (loops). To what extent the backbone conformation of RNA loops also results from interactions within the local sequence context or rather follows global optimization constraints alone has not been addressed yet. Because the majority of base stacking interactions are exerted locally, a critical influence of local sequence on local structure appears plausible. Thus, local loop structure ought to be predictable, at least in part, from the local sequence context alone. To test this hypothesis, we used Random Forests on a nonredundant data set of unpaired nucleotides extracted from 97 X-ray structures from the Protein Data Bank (PDB) to predict discrete backbone angle conformations given by the discretized η/θ-pseudo-torsional space. Predictions on balanced sets with four to six conformational classes using local sequence information yielded average accuracies of up to 55%, thus significantly better than expected by chance (17%-25%). Bases close to the central nucleotide appear to be most tightly linked to its conformation. Our results suggest that RNA loop structure does not only depend on long-range base-pairing interactions; instead, it appears that local sequence context exerts a significant influence on the formation of the local loop structure.

Polymerase ribozyme efficiency increased by G/T-rich DNA oligonucleotides.

The RNA world hypothesis states that the early evolution of life went through a stage where RNA served as genome and as catalyst. The replication of RNA world organisms would have been facilitated by ribozymes that catalyze RNA polymerization. To recapitulate an RNA world in the laboratory, a series of RNA polymerase ribozymes was developed previously. However, these ribozymes have a polymerization efficiency that is too low for self-replication, and the most efficient ribozymes prefer one specific template sequence. The limiting factor for polymerization efficiency is the weak sequence-independent binding to its primer/template substrate. Most of the known polymerase ribozymes bind an RNA heptanucleotide to form the P2 duplex on the ribozyme. By modifying this heptanucleotide, we were able to significantly increase polymerization efficiency. Truncations at the 3'-terminus of this heptanucleotide increased full-length primer extension by 10-fold, on a specific template sequence. In contrast, polymerization on several different template sequences was improved dramatically by replacing the RNA heptanucleotide with DNA oligomers containing randomized sequences of 15 nt. The presence of G and T in the random sequences was sufficient for this effect, with an optimal composition of 60% G and 40% T. Our results indicate that these DNA sequences function by establishing many weak and nonspecific base-pairing interactions to the single-stranded portion of the template. Such low-specificity interactions could have had important functions in an RNA world.

Next-generation tag sequencing for cancer gene expression profiling.

We describe a new method, Tag-seq, which employs ultra high-throughput sequencing of 21 base pair cDNA tags for sensitive and cost-effective gene expression profiling. We compared Tag-seq data to LongSAGE data and observed improved representation of several classes of rare transcripts, including transcription factors, antisense transcripts, and intronic sequences, the latter possibly representing novel exons or genes. We observed increases in the diversity, abundance, and dynamic range of such rare transcripts and took advantage of the greater dynamic range of expression to identify, in cancers and normal libraries, altered expression ratios of alternative transcript isoforms. The strand-specific information of Tag-seq reads further allowed us to detect altered expression ratios of sense and antisense (S-AS) transcripts between cancer and normal libraries. S-AS transcripts were enriched in known cancer genes, while transcript isoforms were enriched in miRNA targeting sites. We found that transcript abundance had a stronger GC-bias in LongSAGE than Tag-seq, such that AT-rich tags were less abundant than GC-rich tags in LongSAGE. Tag-seq also performed better in gene discovery, identifying >98% of genes detected by LongSAGE and profiling a distinct subset of the transcriptome characterized by AT-rich genes, which was expressed at levels below those detectable by LongSAGE. Overall, Tag-seq is sensitive to rare transcripts, has less sequence composition bias relative to LongSAGE, and allows differential expression analysis for a greater range of transcripts, including transcripts encoding important regulatory molecules.

Unusual composition of a yeast chromosome arm is associated with its delayed replication.

The 11.3-Mb genome of the yeast Lachancea (Saccharomyces) kluyveri displays an intriguing compositional heterogeneity: a region of approximately 1 Mb, covering almost the whole left arm of chromosome C (C-left), has an average GC content of 52.9%, which is significantly higher than the 40.4% global GC content of the rest of the genome. This region contains the MAT locus, which remains normal in composition. The excess of GC base pairs affects both coding and noncoding sequences, and thus is not due to selective pressure acting on protein sequences. It leads to a strong codon usage bias and alters the amino acid composition of the 457 proteins encoded on C-left that do not show obvious bias for functional categories, or the presence of paralogs or orthologs of essential genes of Saccharomyces cerevisiae. They share significant synteny conservation with other species of the Saccharomycetaceae, and phylogenetic analysis indicates that C-left originates from a Lachancea species. In contrast, there is a complete absence of transposable elements in C-left, whereas 18 elements per megabase are distributed across the rest of the genome. Comparative hybridization of synchronized cells using high-density genome arrays reveals that C-left is replicated later during S phase than the rest of the genome. Two possible primary causes of this major compositional heterogeneity are discussed: an ancient hybridization of two related species with very distinct GC composition, or an intrinsic mechanism, possibly associated with the loss of the silent cassettes from C-left that progressively increased the GC content and generated the delayed replication of this chromosomal arm.

[Oxidation and deamination of nucleobases as an epigenetic tool]. / Oksydacyjnie modyfikowane i deaminowane zasady DNA jako czynnik epigenetyczny.

 Recent discoveries have demonstrated that 5-methylcytosine (5mC) may be hydroxymethylated to 5-hydroxymethylcytosine (5hmC) in mammals and that genomic DNA may contain about 0.02-0.7% of 5hmC. The aforementioned modification is the key intermediate of active DNA demethylation and has been named "the sixth base in DNA". Although active DNA demethylation in mammals is still controversial, the most plausible mechanism/s of active 5mC demethylation include involvement of three families of enzymes; i) Tet, which is involved in hydroxylation of 5mC to form 5hmC, which can be further oxidized to 5-formylcytosine (5fC) and 5-carboxylcytosine (5caC); ii) deamination of 5mC (or 5hmC) by AID/APOBEC to form thymine or 5-hydroxymethyluracil (5hmU) mispaired with guanine; iii) the BER pathway induced by involvement of TDG glycosylase to replace the above described base modification (5fC, 5caC, 5hmU) with cytosine to demethylate DNA. A plausible scenario for engagement of TDG glycosylase (or some other G-T glycosylase) is through prior deamination of 5-mC to thymine, which generates a G: T substrate for the enzyme. Here cytidine deaminase of the AID/APOBEC family was implicated in the deamination step. It is possible that TDG may act in concert with these deaminases. It seems that mutations are not the only effect of oxidatively modified DNA bases. These, as yet, understudied aspects of the damage suggest a potential for 8-oxoguanine (8-oxoGua) to affect gene expression via chromatin relaxation. It is possible that 8-oxoGua presence in specific DNA sequences may be widely used for transcription regulation, which suggests the epigenetic nature of 8-oxoGua presence in DNA.

Ehrlichia chaffeensis TRP120 binds a G+C-rich motif in host cell DNA and exhibits eukaryotic transcriptional activator function.

Ehrlichia chaffeensis is an obligately intracellular bacterium that modulates host cell gene transcription in the mononuclear phagocyte, but the host gene targets and mechanisms involved in transcriptional modulation are not well-defined. In this study, we identified a novel tandem repeat DNA-binding domain in the E. chaffeensis 120-kDa tandem repeat protein (TRP120) that directly binds host cell DNA. TRP120 was observed by immunofluorescent microscopy in the nucleus of E. chaffeensis-infected host cells and was detected in nuclear extracts by Western immunoblotting with TRP120-specific antibody. The TRP120 binding sites and associated host cell target genes were identified using high-throughput deep sequencing (Illumina) of immunoprecipitated DNA (chromatin immunoprecipitation and high-throughput DNA sequencing). Multiple em motif elicitation (MEME) analysis of the most highly enriched TRP120-bound sequences revealed a G+C-rich DNA motif, and recombinant TRP120 specifically bound synthetic oligonucleotides containing the motif. TRP120 target gene binding sites were mapped most frequently to intersecting regions (intron/exon; 49%) but were also identified in upstream regulatory regions (25%) and downstream locations (26%). Genes targeted by TRP120 were most frequently associated with transcriptional regulation, signal transduction, and apoptosis. TRP120 targeted inflammatory chemokine genes, CCL2, CCL20, and CXCL11, which were strongly upregulated during E. chaffeensis infection and were also upregulated by direct transfection with recombinant TRP120. This study reveals that TRP120 is a novel DNA-binding protein that is involved in a host gene transcriptional regulation strategy.

Low-complexity sequences and single amino acid repeats: not just "junk" peptide sequences.

For decades proteins were thought to interact in a "lock and key" system, which led to the definition of a paradigm linking stable three-dimensional structure to biological function. As a consequence, any non-structured peptide was considered to be nonfunctional and to evolve neutrally. Surprisingly, the most commonly shared peptides between eukaryotic proteomes are low-complexity sequences that in most conditions do not present a stable three-dimensional structure. However, because these sequences evolve rapidly and because the size variation of a few of them can have deleterious effects, low-complexity sequences have been suggested to be the target of selection. Here we review evidence that supports the idea that these simple sequences should not be considered just "junk" peptides and that selection drives the evolution of many of them.

Acidithiobacillus ferrivorans, sp. nov.; facultatively anaerobic, psychrotolerant iron-, and sulfur-oxidizing acidophiles isolated from metal mine-impacted environments.

Phenotypic and genotypic analysis was carried out on four iron- and sulfur-oxidizing acidophilic bacteria (the "NO-37 group") isolated from different parts of the world. 16S rRNA phylogeny showed that they are highly related to each other, but are less related to the type strain of Acidithiobacillus ferrooxidans. The NO-37 group isolates are obligate chemolithoautotrophs, facultative anaerobes, diazotrophic, and psychrotolerant. They are less tolerant of extremely low pH, and in contrast to At. ferrooxidans (T), all of the NO-37 group isolates are motile. The GC contents of genomic DNA of the NO-37 group isolates were around 56 mol% and the DNA-DNA hybridization value between genomic DNA of isolate NO-37 and At. ferrooxidans (T) was 37%. It also appears that the bacteria of the NO-37 group have a different biochemical mechanism for oxidizing ferrous iron than At. ferrooxidans (T); the gene coding for the archetypal rusticyanin (RusA) was not detected in any of the NO-37 group isolates, rather a gene coding for a homologous protein (RusB) was amplified from three of the four novel isolates. Isolates of the NO-37 group clearly belong to a species that is different to those already recognized in the genus Acidithiobacillus, for which the name Acidithiobacillus ferrivorans is proposed.

Sites of strong Rec12/Spo11 binding in the fission yeast genome are associated with meiotic recombination and with centromeres.

Meiotic recombination arises from Rec12/Spo11-dependent formation of DNA double-strand breaks (DSBs) and their subsequent repair. We identified Rec12-binding peaks across the Schizosaccharomyces pombe genome using chromatin immunoprecipitation after reversible formaldehyde cross-linking combined with whole-genome DNA microarrays. Strong Rec12 binding coincided with previously identified DSBs at the recombination hotspots ura4A, mbs1, and mbs2 and correlated with DSB formation at a new site. In addition, Rec12 binding corresponded to eight novel conversion hotspots and correlated with crossover density in segments of chromosome I. Notably, Rec12 binding inversely correlated with guanine-cytosine (GC) content, contrary to findings in Saccharomyces cerevisiae. Although both replication origins and Rec12-binding sites preferred AT-rich gene-free regions, they seemed to exclude each other. We also uncovered a connection between binding sites of Rec12 and meiotic cohesin Rec8. Rec12-binding peaks lay often within 2.5 kb of a Rec8-binding peak. Rec12 binding showed preference for large intergenic regions and was found to bind preferentially near to genes expressed strongly in meiosis. Surprisingly, Rec12 binding was also detected in centromeric core regions, which raises the intriguing possibility that Rec12 plays additional roles in meiotic chromosome dynamics.

High guanine and cytosine content increases mRNA levels in mammalian cells.

Mammalian genes are highly heterogeneous with respect to their nucleotide composition, but the functional consequences of this heterogeneity are not clear. In the previous studies, weak positive or negative correlations have been found between the silent-site guanine and cytosine (GC) content and expression of mammalian genes. However, previous studies disregarded differences in the genomic context of genes, which could potentially obscure any correlation between GC content and expression. In the present work, we directly compared the expression of GC-rich and GC-poor genes placed in the context of identical promoters and UTR sequences. We performed transient and stable transfections of mammalian cells with GC-rich and GC-poor versions of Hsp70, green fluorescent protein, and IL2 genes. The GC-rich genes were expressed several-fold to over a 100-fold more efficiently than their GC-poor counterparts. This effect was not due to different translation rates of GC-rich and GC-poor mRNA. On the contrary, the efficient expression of GC-rich genes resulted from their increased steady-state mRNA levels. mRNA degradation rates were not correlated with GC content, suggesting that efficient transcription or mRNA processing is responsible for the high expression of GC-rich genes. We conclude that silent-site GC content correlates with gene expression efficiency in mammalian cells.

A systematic approach re-analyzing the effects of temperature disturbance on the microbial community of mesophilic anaerobic digestion.

Microbial communities are key drivers of ecosystem processes, but their behavior in disturbed environments is difficult to measure. How microbial community composition and function respond disturbances is a common challenge in biomedical, environmental, agricultural, and bioenergy research. A novel way to solve this problem is to use a systems-level perspective and describe microbial communities as networks. Based on a mesophilic anaerobic digestion system of swine manure as a tool, we propose a simple framework to investigate changes in microbial communities via compositions, metabolic pathways, genomic properties and interspecies relationships in response to a long-term temperature disturbance. After temperature disturbance, microbial communities tend towards a competitive interaction network with higher GC content and larger genome size. Based on microbial interaction networks, communities responded to the disturbance by showing a transition from acetotrophic (Methanotrichaceae and Methanosarcinaceae) to methylotrophic methanogens (Methanomassiliicoccaceae and Methanobacteriaceae) and a fluctuation in rare biosphere taxa. To conclude, this study may be important for exploring the dynamic relationships between disturbance and microbial communities as a whole, as well as for providing researchers with a better understanding of how changes in microbial communities relate to ecological processes.