Evolutionary Dynamics of Clustered Irregularly Interspaced Short Palindromic Repeat Systems in the Ocean Metagenome

Clustered regularly interspaced short palindromic repeats (CRISPRs) form a recently characterized type of prokaryotic antiphage defense system. The phage-host interactions involving CRISPRs have been studied in experiments with selected bacterial or archaeal species and, computationally, in completely sequenced genomes. However, these studies do not allow one to take prokaryotic population diversity and phage-host interaction dynamics into account. This gap can be filled by using metagenomic data: in particular, the largest existing data set, generated from the Sorcerer II Global Ocean Sampling expedition. The application of three publicly available CRISPR recognition programs to the Global Ocean metagenome produced a large proportion of false-positive results. To address this problem, a filtering procedure was designed. It resulted in about 200 reliable CRISPR cassettes, which were then studied in detail. The repeat consensuses were clustered into several stable classes that differed from the existing classification. Short fragments of DNA similar to the cassette spacers were more frequently present in the same geographical location than in other locations (P, <0.0001). We developed a catalogue of elementary CRISPR-forming events and reconstructed the likely evolutionary history of cassettes that had common spacers. Metagenomic collections allow for relatively unbiased analysis of phage-host interactions and CRISPR evolution. The results of this study demonstrate that CRISPR cassettes retain the memory of the local virus population at a particular ocean location. CRISPR evolution may be described using a limited vocabulary of elementary events that have a natural biological interpretation.Prokaryotes are highly diverse (33). One of the explanations of this diversity is the high extinction rate, due to genetic aggression, which leads to the clearance of ecological niches and, as a result, may allow new prokaryotic species to emerge. In the absence of host defense, viral infection of prokaryotic colonies results in colony extinction or the fixation of a fraction of the invader''s genetic material in the host genome, profoundly affecting the life cycle of the host (32). Thus, bacteria and archaea have developed various kinds of defense mechanisms to resist this pressure; the best studied of these mechanisms is restriction-modification systems (4).Along with well-known prokaryotic defense mechanisms, such as rapid evolution of cell receptors or the use of restriction-modification or toxin-antitoxin systems (see, e.g., references 6, 21, and 25), newly discovered clustered regularly interspaced palindromic repeat (CRISPR) systems seem to play an important role in protecting the cell from archaeal virus or bacteriophage assaults (reviewed in reference 36). A typical CRISPR system is a genetic locus comprising CRISPR-associated (cas) genes coding for proteins of several distinct functional classes (8, 19, 29) and a CRISPR cassette. A CRISPR cassette is formed by almost identical direct repeats with an average length of 32 nucleotides (nt), which are separated by similarly sized, unique spacers. A considerable proportion of spacers is similar to known phage or virus sequences, suggesting that the system is involved in antivirus defense (8, 29, 31). This involvement was experimentally demonstrated when a CRISPR system was shown to be essential for cell survival after invasion by foreign DNA (5). The mechanism is thought to be analogous to eukaryotic RNA interference (29), but it has not been characterized in detail yet.CRISPR cassettes retain information that could be used to reveal the evolutionary history of individual systems. First, it has been shown that CRISPR-associated genes could be divided into eight subtypes according to operon organization and gene phylogeny (19). Second, the repeats of different CRISPR cassettes may be similar, which might indicate a common origin of such cassettes. The first attempt to cluster CRISPR cassettes by the similarity of repeat sequences resulted in 12 clusters (27). In that study, the cassettes were obtained by the application of PILER-CR to completely sequenced genomes. Third, pairwise comparison of spacers could also reveal the specific evolutionary history of individual CRISPR cassettes.So far, most large-scale studies of CRISPR systems have been restricted to well-studied organisms with completely sequenced genomes (5, 9, 20, 28, 30). However, the dynamic interaction between viruses or phages and microorganisms in natural environments is of particular interest (2, 10, 15, 23, 35, 38, 40-42). It may be studied using CRISPRs in a metagenome, that is, sequenced DNA fragments collected in one geographical location and therefore representing one ecological niche with all its inhabitants. This approach is interesting for two reasons. First, metagenomic samples provide a common census of coexisting organisms, i.e., in many cases, both the infecting viruses and phages and their victims. Second, most bacteria and archaea from metagenomic samples cannot be cultivated, and hence little is known about their CRISPR systems.To date, three studies have considered host-virus interactions in metagenomes. One study used two thermophilic Synechococcus isolates from microbial mats in hot springs at Yellowstone National Park to demonstrate fast coevolution of the host and phage genomes (22). Two studies described archaeal and bacterial interactions with viruses and phages, respectively, in acidophilic biofilms (2, 39). All environmental communities analyzed so far are extreme and are dominated by few species. Natural samples containing many diverse coexisting organisms may arguably be more interesting.The largest available metagenome, produced by the Sorcerer II Global Ocean Sampling (GOS) expedition, comprises samples of genetic material collected from more than 50 geographical locations of the Pacific and Atlantic oceans (34). This variety provides an opportunity to study the evolution of phage-host interactions reflected in CRISPRs.Three algorithms, PILER-CR (14), the CRISPR recognition tool (CRT) (7), and CRISPRFinder (18), have been developed as tools for the discovery of new CRISPR cassettes. All these algorithms define candidate CRISPR cassette sequences as short direct repeats separated by short unique spacers; they then use a variety of standard repeat-finding techniques. However, the implementation of specific details is different.PILER-CR constructs local alignments of the input sequence to itself; each hit between two close regions is a candidate for an alignment of a repeat with its neighbor copy. In terms of dynamic programming, taking into account the repeat structure of a CRISPR cassette implies looking for hits only within a relatively narrow band around the main diagonal of the dot plot. This process is followed by several refinement steps.CRT does not use alignments to identify candidate repeats; rather, it derives them directly from the analysis of an input sequence. It is based on finding series of short repeats of a specified length (searching for exact k-mer matches) and then extending these repeats (increasing k-mer length) while allowing for a certain level of mismatches.Finally, CRISPRFinder is based on a suffix-tree-based algorithm for repeat discovery, again with additional refinement.All three algorithms were used for the CRISPR cassette search in this study.     

c-Myc and Sp1 contribute to proviral latency by recruiting histone deacetylase 1 to the human immunodeficiency virus type 1 promoter

Histone deacetylase (HDAC) inhibitors such as valproic acid (VPA) induce the expression of quiescent proviral human immunodeficiency virus type 1 (HIV-1) and may deplete proviral infection in vivo. To uncover novel molecular mechanisms that maintain HIV latency, we sought cellular mRNAs whose expression was diminished in resting CD4(+) T cells of HIV-1-infected patients exposed to VPA. c-Myc was prominent among genes markedly downregulated upon exposure to VPA. c-Myc expression repressed HIV-1 expression in chronically infected cell lines. Chromatin immunoprecipitation (ChIP) assays revealed that c-Myc and HDAC1 are coordinately resident at the HIV-1 long terminal repeat (LTR) promoter and absent from the promoter after VPA treatment in concert with histone acetylation, RNA polymerase II recruitment, and LTR expression. Sequential ChIP assays demonstrated that c-Myc, Sp1, and HDAC1 coexist in the same DNA-protein complex at the HIV promoter. Short hairpin RNA inhibition of c-Myc reduces both c-Myc and HDAC1 occupancy, blocks c-Myc repression of Tat activation, and increases LTR expression. These results expand the understanding of mechanisms that recruit HDAC and maintain the latency of HIV-1, suggesting novel therapeutic approaches against latent proviral HIV infection.     

Cajal bodies and the nucleolus are required for a plant virus systemic infection

The nucleolus and Cajal bodies (CBs) are prominent interacting subnuclear domains involved in a number of crucial aspects of cell function. Certain viruses interact with these compartments but the functions of such interactions are largely uncharacterized. Here, we show that the ability of the groundnut rosette virus open reading frame (ORF) 3 protein to move viral RNA long distances through the phloem strictly depends on its interaction with CBs and the nucleolus. The ORF3 protein targets and reorganizes CBs into multiple CB-like structures and then enters the nucleolus by causing fusion of these structures with the nucleolus. The nucleolar localization of the ORF3 protein is essential for subsequent formation of viral ribonucleoprotein (RNP) particles capable of virus long-distance movement and systemic infection. We provide a model whereby the ORF3 protein utilizes trafficking pathways involving CBs to enter the nucleolus and, along with fibrillarin, exit the nucleus to form viral 'transport-competent' RNP particles in the cytoplasm.     

Metamorphosis of cinctoblastula larvae (Homoscleromorpha, porifera)

The metamorphosis of the cinctoblastula of Homoscleromorpha is studied in five species belonging to three genera. The different steps of metamorphosis are similar in all species. The metamorphosis occurs by the invagination and involution of either the anterior epithelium or the posterior epithelium of the larva. During metamorphosis, morphogenetic polymorphism was observed, which has an individual character and does not depend on either external or species specific factors. In the rhagon, the development of the aquiferous system occurs only by epithelial morphogenesis and subsequent differentiation of cells. Mesohylar cells derive from flagellated cells after ingression. The formation of pinacoderm and choanoderm occurs by the differentiation of the larval flagellated epithelium. This is possibly due to the conservation of cell junctions in the external surface of the larval flagellated cells and of the basement membrane in their internal surface. The main difference in homoscleromorph metamorphosis compared with Demospongiae is the persistence of the flagellated epithelium throughout this process and even in the adult since exo- and endopinacoderm remain flagellated. The antero-posterior axis of the larva corresponds to the baso-apical axis of the adult in Homoscleromorpha.     

A potent and orally active HIV-1 integrase inhibitor

A 1,6-naphthyridine inhibitor of HIV-1 integrase has been discovered with excellent inhibitory activity in cells, good pharmacokinetics, and an excellent ability to inhibit virus with mutant enzyme.     

Long-term data reveal effects of climate,road access,and latitude on mountain goat horn size

Potential negative artificial selection on horn size is a concern for many harvested ungulates. The mountain goat (Oreamnos americanus) has distinct black horns, but targeting animals based on horn size in the field can be challenging. We analyzed over 23,000 horn records that included base circumference and total length, from which we also derived horn volume, from mountain goats harvested in Alaska, British Columbia, and the Northwest Territories from 1980 to 2016. We tested 3 potential drivers of horn size variation: geographical location, environmental conditions, and artificial selection. We found no support for a latitudinal effect with surprisingly little variation across the sampling distribution. The Pacific Decadal Oscillation had the largest effect outside age in the model, suggesting a role of climate in shaping variation. Mountain goats harvested closer to roads had larger horns, indicating that ease of access might allow hunters to be more selective, though the effect size was small. Our findings reinforce the value of accurate and complete record keeping on horn size, age, and sex of harvested animals, and highlight the importance of explicitly considering climate and accessibility when devising management strategies for the mountain goat.     

Feeding Ecology and Regurgitation–Reingestion Behavior of the Critically Endangered Indri indri in the Maromizaha Protected Area,Eastern Madagascar

International Journal of Primatology - Deforestation around the world is a major threat to primates. Understanding primate species’ habitat and dietary requirements is critical in creating...     

Quantitative analysis of mRNA expression levels and DNA methylation profiles of three neighboring genes: FUS1, NPRL2/G21 and RASSF1A in non-small cell lung cancer patients

Background

Tumor suppressor gene (TSG) inactivation plays a crucial role in carcinogenesis. FUS1, NPRL2/G21 and RASSF1A are TSGs from LUCA region at 3p21.3, a critical chromosomal region in lung cancer development. The aim of the study was to analyze and compare the expression levels of these 3 TSGs in NSCLC, as well as in macroscopically unchanged lung tissue surrounding the primary lesion, and to look for the possible epigenetic mechanism of TSG inactivation via gene promoter methylation.

Methods

Expression levels of 3 TSGs and 2 DNA methyltransferases, DNMT1 and DNMT3B, were assessed using real-time PCR method (qPCR) in 59 primary non-small cell lung tumors and the matched macroscopically unchanged lung tissue samples. Promoter methylation status of TSGs was analyzed using methylation-specific PCRs (MSP method) and Methylation Index (MI) value was calculated for each gene.

Results

The expression of all three TSGs were significantly different between NSCLC subtypes: RASSF1A and FUS1 expression levels were significantly lower in squamous cell carcinoma (SCC), and NPRL2/G21 in adenocarcinoma (AC). RASSF1A showed significantly lower expression in tumors vs macroscopically unchanged lung tissues. Methylation frequency was 38–76 %, depending on the gene. The highest MI value was found for RASSF1A (52 %) and the lowest for NPRL2/G21 (5 %). The simultaneous decreased expression and methylation of at least one RASSF1A allele was observed in 71 % tumor samples. Inverse correlation between gene expression and promoter methylation was found for FUS1 (rs = −0.41) in SCC subtype. Expression levels of DNMTs were significantly increased in 75–92 % NSCLCs and were significantly higher in tumors than in normal lung tissue. However, no correlation between mRNA expression levels of DNMTs and DNA methylation status of the studied TSGs was found.

Conclusions

The results indicate the potential role of the studied TSGs in the differentiation of NSCLC histopathological subtypes. The significant differences in RASSF1A expression levels between NSCLC and macroscopically unchanged lung tissue highlight its possible diagnostic role in lung cancer in situ recognition. High percentage of lung tumor samples with simultaneous RASSF1A decreased expression and gene promoter methylation indicates its epigenetic silencing. However, DNMT overexpression doesn’t seem to be a critical determinate of its promoter hypermethylation.