Chromosome-Length Assembly of the Baikal Seal (Pusa sibirica) Genome Reveals a Historically Large Population Prior to Isolation in Lake Baikal

Pusa sibirica, the Baikal seal, is the only extant, exclusively freshwater, pinniped species. The pending issue is, how and when they reached their current habitat—the rift lake Baikal, more than three thousand kilometers away from the Arctic Ocean. To explore the demographic history and genetic diversity of this species, we generated a de novo chromosome-length assembly, and compared it with three closely related marine pinniped species. Multiple whole genome alignment of the four species compared with their karyotypes showed high conservation of chromosomal features, except for three large inversions on chromosome VI. We found the mean heterozygosity of the studied Baikal seal individuals was relatively low (0.61 SNPs/kbp), but comparable to other analyzed pinniped samples. Demographic reconstruction of seals revealed differing trajectories, yet remarkable variations in Ne occurred during approximately the same time periods. The Baikal seal showed a significantly more severe decline relative to other species. This could be due to the difference in environmental conditions encountered by the earlier populations of Baikal seals, as ice sheets changed during glacial–interglacial cycles. We connect this period to the time of migration to Lake Baikal, which occurred ~3–0.3 Mya, after which the population stabilized, indicating balanced habitat conditions.

Transcriptomic Analysis of Flowering Time Genes in Cultivated Chickpea and Wild Cicer

Chickpea (Cicer arietinum L.) is a major grain legume and a good source of plant-based protein. However, comprehensive knowledge of flowering time control in Cicer is lacking. In this study, we acquire high-throughput transcriptome sequencing data and analyze changes in gene expression during floral transition in the early flowering cultivar ICCV 96029, later flowering C. arietinum accessions, and two wild species, C. reticulatum and C. echinospermum. We identify Cicer orthologs of A. thaliana flowering time genes and analyze differential expression of 278 genes between four species/accessions, three tissue types, and two conditions. Our results show that the differences in gene expression between ICCV 96029 and other cultivated chickpea accessions are vernalization-dependent. In addition, we highlight the role of FTa3, an ortholog of FLOWERING LOCUS T in Arabidopsis, in the vernalization response of cultivated chickpea. A common set of differentially expressed genes was found for all comparisons between wild species and cultivars. The direction of expression change for different copies of the FT-INTERACTING PROTEIN 1 gene was variable in different comparisons, which suggests complex mechanisms of FT protein transport. Our study makes a contribution to the understanding of flowering time control in Cicer, and can provide genetic strategies to further improve this important agronomic trait.

Aging Induces Profound Changes in sncRNA in Rat Sperm and These Changes Are Modified by Perinatal Exposure to Environmental Flame Retardant

Advanced paternal age at fertilization is a risk factor for multiple disorders in offspring and may be linked to age-related epigenetic changes in the father's sperm. An understanding of aging-related epigenetic changes in sperm and environmental factors that modify such changes is needed. Here, we characterize changes in sperm small non-coding RNA (sncRNA) between young pubertal and mature rats. We also analyze the modification of these changes by exposure to environmental xenobiotic 2,2',4,4'-tetrabromodiphenyl ether (BDE-47). sncRNA libraries prepared from epididymal spermatozoa were sequenced and analyzed using DESeq 2. The distribution of small RNA fractions changed with age, with fractions mapping to rRNA and lncRNA decreasing and fractions mapping to tRNA and miRNA increasing. In total, 249 miRNA, 908 piRNA and 227 tRNA-derived RNA were differentially expressed (twofold change, false discovery rate (FDR) p ≤ 0.05) between age groups in control animals. Differentially expressed miRNA and piRNA were enriched for protein-coding targets involved in development and metabolism, while piRNA were enriched for long terminal repeat (LTR) targets. BDE-47 accelerated age-dependent changes in sncRNA in younger animals, decelerated these changes in older animals and increased the variance in expression of all sncRNA. Our results indicate that the natural aging process has profound effects on sperm sncRNA profiles and this effect may be modified by environmental exposure.

Rat liver epigenome programing by perinatal exposure to 2,2',4'4'-tetrabromodiphenyl ether

Perinatal exposures to polybrominated diphenyl ethers permanently reprogram liver metabolism and induce a nonalcoholic fatty liver disease-like phenotype and insulin resistance in rodents. Aim: To test if these changes are associated with altered liver epigenome. Materials & methods: Expression of small RNA and changes in DNA methylation in livers of adult rats were analyzed following perinatal exposure to 2,2',4,4'-tetrabromodiphenyl ether, the polybrominated diphenyl ether congener most prevalent in human tissues. Results: We identified 33 differentially methylated DNA regions and 15 differentially expressed miRNAs. These changes were enriched for terms related to lipid and carbohydrate metabolism, insulin signaling, Type-2 diabetes and nonalcoholic fatty liver disease. Conclusion: Changes in the liver epigenome are a likely candidate mechanism of long-term maintenance of an aberrant metabolic phenotype.

Thirteen New Plastid Genomes from Mixotrophic and Autotrophic Species Provide Insights into Heterotrophy Evolution in Neottieae Orchids

Mixotrophic species use both organic and mineral carbon sources. Some mixotrophic plants combine photosynthesis and a nutrition called mycoheterotrophy, where carbon is obtained from fungi forming mycorrhizal symbiosis with their roots. These species can lose photosynthetic abilities and evolve full mycoheterotrophy. Besides morphological changes, the latter transition is associated with a deep alteration of the plastid genome. Photosynthesis-related genes are lost first, followed by housekeeping genes, eventually resulting in a highly reduced genome. Whether relaxation of selective constraints already occurs for the plastid genome of mixotrophic species, which remain photosynthetic, is unclear. This is partly due to the difficulty of comparing plastid genomes of autotrophic, mixotrophic, and mycoheterotrophic species in a narrow phylogenetic framework. We address this question in the orchid tribe Neottieae, where this large assortment of nutrition types occurs. We sequenced 13 new plastid genomes, including 9 mixotrophic species and covering all 6 Neottieae genera. We investigated selective pressure on plastid genes in each nutrition type and conducted a phylogenetic inference of the group. Surprisingly, photosynthesis-related genes did not experience selection relaxation in mixotrophic species compared with autotrophic relatives. Conversely, we observed evidence for selection intensification for some plastid genes. Photosynthesis is thus still under purifying selection, maybe because of its role in fruit formation and thus reproductive success. Phylogenetic analysis resolved most relationships, but short branches at the base of the tree suggest an evolutionary radiation at the beginning of Neottieae history, which, we hypothesize, may be linked to mixotrophy emergence.

Single-nucleotide-resolution mapping of DNA gyrase cleavage sites across the Escherichia coli genome

An important antibiotic target, DNA gyrase is an essential bacterial enzyme that introduces negative supercoils into DNA and relaxes positive supercoils accumulating in front of moving DNA and RNA polymerases. By altering the superhelical density, gyrase may regulate expression of bacterial genes. The information about how gyrase is distributed along genomic DNA and whether its distribution is affected by drugs is scarce. During catalysis, gyrase cleaves both DNA strands forming a covalently bound intermediate. By exploiting the ability of several topoisomerase poisons to stabilize this intermediate we developed a ChIP-Seq-based approach to locate, with single nucleotide resolution, DNA gyrase cleavage sites (GCSs) throughout the Escherichia coli genome. We identified an extended gyrase binding motif with phased 10-bp G/C content variation, indicating that bending ability of DNA contributes to gyrase binding. We also found that GCSs are enriched in extended regions located downstream of highly transcribed operons. Transcription inhibition leads to redistribution of gyrase suggesting that the enrichment is functionally significant. Our method can be applied for precise mapping of prokaryotic and eukaryotic type II topoisomerases cleavage sites in a variety of organisms and paves the way for future studies of various topoisomerase inhibitors.

Peripubertal serum dioxin concentrations and subsequent sperm methylome profiles of young Russian adults

BACKGROUND

The association of exposure to endocrine disrupting chemicals in the peripubertal period with subsequent sperm DNA methylation is unknown.

OBJECTIVE

We examined the association of peripubertal serum 2,3,7,8-tetrachlorodibenzo-p-dioxin (TCDD) concentrations with whole-genome bisulfite sequencing (WGBS) of sperm collected in young adulthood.

METHODS

The Russian Children's Study is a prospective cohort of 516 boys who were enrolled at 8-9 years of age and provided semen samples at 18-19 years of age. WGBS of sperm was conducted to identify differentially methylated regions (DMR) between highest (n = 4) and lowest (n = 4) peripubertal TCDD groups.

RESULTS

We found 52 DMRs that distinguished lowest and highest peripubertal serum TCDD concentrations. One of the top scoring networks, "Cellular Assembly and Organization, Cellular Function and Maintenance, Carbohydrate Metabolism", identified estrogen receptor alpha as its central regulator.

CONCLUSION

Findings from our limited sample size suggest that peripubertal environmental exposures are associated with sperm DNA methylation in young adults.

Spacer-length DNA intermediates are associated with Cas1 in cells undergoing primed CRISPR adaptation

During primed CRISPR adaptation spacers are preferentially selected from DNA recognized by CRISPR interference machinery, which in the case of Type I CRISPR-Cas systems consists of CRISPR RNA (crRNA) bound effector Cascade complex that locates complementary targets, and Cas3 executor nuclease/helicase. A complex of Cas1 and Cas2 proteins is capable of inserting new spacers in the CRISPR array. Here, we show that in Escherichia coli cells undergoing primed adaptation, spacer-sized fragments of foreign DNA are associated with Cas1. Based on sensitivity to digestion with nucleases, the associated DNA is not in a standard double-stranded state. Spacer-sized fragments are cut from one strand of foreign DNA in Cas1- and Cas3-dependent manner. These fragments are generated from much longer S1-nuclease sensitive fragments of foreign DNA that require Cas3 for their production. We propose that in the course of CRISPR interference Cas3 generates fragments of foreign DNA that are recognized by the Cas1-Cas2 adaptation complex, which excises spacer-sized fragments and channels them for insertion into CRISPR array.

Dynamics of Escherichia coli type I-E CRISPR spacers over 42 000 years

CRISPR-Cas are nucleic acid-based prokaryotic immune systems. CRISPR arrays accumulate spacers from foreign DNA and provide resistance to mobile genetic elements containing identical or similar sequences. Thus, the set of spacers present in a given bacterium can be regarded as a record of encounters of its ancestors with genetic invaders. Such records should be specific for different lineages and change with time, as earlier acquired spacers get obsolete and are lost. Here, we studied type I-E CRISPR spacers of Escherichia coli from extinct pachyderm. We find that many spacers recovered from intestines of a 42 000-year-old mammoth match spacers of present-day E. coli. Present-day CRISPR arrays can be reconstructed from palaeo sequences, indicating that the order of spacers has also been preserved. The results suggest that E. coli CRISPR arrays were not subject to intensive change through adaptive acquisition during this time.

Effect of method of deduplication on estimation of differential gene expression using RNA-seq

Background

RNA-seq is a useful tool for analysis of gene expression. However, its robustness is greatly affected by a number of artifacts. One of them is the presence of duplicated reads.

Results

To infer the influence of different methods of removal of duplicated reads on estimation of gene expression in cancer genomics, we analyzed paired samples of hepatocellular carcinoma (HCC) and non-tumor liver tissue. Four protocols of data analysis were applied to each sample: processing without deduplication, deduplication using a method implemented in SAMtools, and deduplication based on one or two molecular indices (MI). We also analyzed the influence of sequencing layout (single read or paired end) and read length. We found that deduplication without MI greatly affects estimated expression values; this effect is the most pronounced for highly expressed genes.

Conclusion

The use of unique molecular identifiers greatly improves accuracy of RNA-seq analysis, especially for highly expressed genes. We developed a set of scripts that enable handling of MI and their incorporation into RNA-seq analysis pipelines. Deduplication without MI affects results of differential gene expression analysis, producing a high proportion of false negative results. The absence of duplicate read removal is biased towards false positives. In those cases where using MI is not possible, we recommend using paired-end sequencing layout.

Foreign DNA acquisition by the I-F CRISPR-Cas system requires all components of the interference machinery

CRISPR immunity depends on acquisition of fragments of foreign DNA into CRISPR arrays. For type I-E CRISPR–Cas systems two modes of spacer acquisition, naïve and primed adaptation, were described. Naïve adaptation requires just two most conserved Cas1 and Cas2 proteins; it leads to spacer acquisition from both foreign and bacterial DNA and results in multiple spacers incapable of immune response. Primed adaptation requires all Cas proteins and a CRISPR RNA recognizing a partially matching target. It leads to selective acquisition of spacers from DNA molecules recognized by priming CRISPR RNA, with most spacers capable of protecting the host. Here, we studied spacer acquisition by a type I-F CRISPR–Cas system. We observe both naïve and primed adaptation. Both processes require not just Cas1 and Cas2, but also intact Csy complex and CRISPR RNA. Primed adaptation shows a gradient of acquisition efficiency as a function of distance from the priming site and a strand bias that is consistent with existence of single-stranded adaption intermediates. The results provide new insights into the mechanism of spacer acquisition and illustrate surprising mechanistic diversity of related CRISPR–Cas systems.

The Cas6e ribonuclease is not required for interference and adaptation by the E. coli type I-E CRISPR-Cas system

CRISPR-Cas are small RNA-based adaptive prokaryotic immunity systems protecting cells from foreign DNA or RNA. Type I CRISPR-Cas systems are composed of a multiprotein complex (Cascade) that, when bound to CRISPR RNA (crRNA), can recognize double-stranded DNA targets and recruit the Cas3 nuclease to destroy target-containing DNA. In the Escherichia coli type I-E CRISPR-Cas system, crRNAs are generated upon transcription of CRISPR arrays consisting of multiple palindromic repeats and intervening spacers through the function of Cas6e endoribonuclease, which cleaves at specific positions of repeat sequences of the CRISPR array transcript. Cas6e is also a component of Cascade. Here, we show that when mature unit-sized crRNAs are provided in a Cas6e-independent manner by transcription termination, the CRISPR-Cas system can function without Cas6e. The results should allow facile interrogation of various targets by type I-E CRISPR-Cas system in E. coli using unit-sized crRNAs generated by transcription.