During ribosomal and transfer RNA maturation, external transcribed spacer (ETS) and internal transcribed spacer (ITS) sequences are excised and, as non-functional by-products, are rapidly degraded. However, we report that the 3'ETS of the glyW-cysT-leuZ polycistronic tRNA precursor is highly and specifically enriched by co-purification with at least two different small regulatory RNAs (sRNAs), RyhB and RybB. Both sRNAs are shown to base pair with the same region in the 3'ETS of leuZ (3'ETS(leuZ)). Disrupting the pairing by mutating 3'ETS(leuZ) strongly increased the activity of sRNAs, even under non-inducing conditions. Our results indicate that 3'ETS(leuZ) prevents sRNA-dependent remodeling of tricarboxylic acid (TCA) cycle fluxes and decreases antibiotic sensitivity when sRNAs are transcriptionally repressed. This suggests that 3'ETS(leuZ) functions as a sponge to absorb transcriptional noise from repressed sRNAs. Additional data showing RybB and MicF sRNAs are co-purified with ITS(metZ-metW) and ITS(metW-metV) strongly suggest a wide distribution of this phenomenon.

The Chlorophyceae, an advanced class of chlorophyte green algae, comprises five lineages that form two major clades (Chlamydomonadales + Sphaeropleales and Oedogoniales + Chaetopeltidales + Chaetophorales). The four complete chloroplast DNA (cpDNA) sequences currently available for chlorophyceans uncovered an extraordinarily fluid genome architecture as well as many structural features distinguishing this group from other green algae. We report here the 521,168-bp cpDNA sequence from a member of the Chaetopeltidales (Floydiella terrestris), the sole chlorophycean lineage not previously sampled for chloroplast genome analysis. This genome, which contains 97 conserved genes and 26 introns (19 group I and 7 group II introns), is the largest chloroplast genome ever sequenced. Intergenic regions account for 77.8% of the genome size and are populated by short repeats. Numerous genomic features are shared with the cpDNA of the chaetophoralean Stigeoclonium helveticum, notably the absence of a large inverted repeat and the presence of unique gene clusters and trans-spliced group II introns. Although only one of the Floydiella group I introns encodes a homing endonuclease gene, our finding of five free-standing reading frames having similarity with such genes suggests that chloroplast group I introns endowed with mobility were once more abundant in the Floydiella lineage. Parsimony analysis of structural genomic features and phylogenetic analysis of chloroplast sequence data unambiguously resolved the Oedogoniales as sister to the Chaetopeltidales and Chaetophorales. An evolutionary scenario of the molecular events that shaped the chloroplast genome in the Chlorophyceae is presented.

The Chlorophyceae (sensu Mattox and Stewart) is a morphologically diverse class of the Chlorophyta displaying biflagellate and quadriflagellate motile cells with varying configurations of the flagellar apparatus. Phylogenetic analyses of 18S rDNA data and combined 18S and 26S rDNA data from a broad range of chlorophycean taxa uncovered five major monophyletic groups (Chlamydomonadales, Sphaeropleales, Oedogoniales, Chaetophorales, and Chaetopeltidales) but could not resolve their branching order. To gain insight into the interrelationships of these groups, we analyzed multiple genes encoded by the chloroplast genomes of Chlamydomonas reinhardtii P. A. Dang. and Chlamydomonas moewusii Gerloff (Chlamydomonadales), Scenedesmus obliquus (Turpin) Kütz. (Sphaeropleales), Oedogonium cardiacum Wittr. (Oedogoniales), Stigeoclonium helveticum Vischer (Chaetophorales), and Floydiella terrestris (Groover et Hofstetter) Friedl et O'Kelly (Chaetopeltidales). The C. moewusii, Oedogonium, and Floydiella chloroplast DNAs were partly sequenced using a random strategy. Trees were reconstructed from nucleotide and amino acid data sets derived from 44 protein-coding genes of 11 chlorophytes and nine streptophytes as well as from 57 protein-coding genes of the six chlorophycean taxa. All best trees identified two robustly supported major lineages within the Chlorophyceae: a clade uniting the Chlamydomonadales and Sphaeropleales, and a clade uniting the Oedogoniales, Chaetophorales, and Chaetopeltidales (OCC clade). This dichotomy is independently supported by molecular signatures in chloroplast genes, such as insertions ⁄ deletions and the distribution of transspliced group II introns. Within the OCC clade, the sister relationship observed for the Chaetophorales and Chaetopeltidales is also strengthened by independent data. Character state reconstruction of basal body orientation allowed us to refine hypotheses regarding the evolution of the flagellar apparatus.

The Genome Analysis Toolkit (GATK) is a popular set of programs for discovering and genotyping variants from next-generation sequencing data. The current GATK recommendation for RNA sequencing (RNA-seq) is to perform variant calling from individual samples, with the drawback that only variable positions are reported. Versions 3.0 and above of GATK offer the possibility of calling DNA variants on cohorts of samples using the HaplotypeCaller algorithm in Genomic Variant Call Format (GVCF) mode. Using this approach, variants are called individually on each sample, generating one GVCF file per sample that lists genotype likelihoods and their genome annotations. In a second step, variants are called from the GVCF files through a joint genotyping analysis. This strategy is more flexible and reduces computational challenges in comparison to the traditional joint discovery workflow. Using a GVCF workflow for mining SNP in RNA-seq data provides substantial advantages, including reporting homozygous genotypes for the reference allele as well as missing data. Taking advantage of RNA-seq data derived from primary macrophages isolated from 50 cows, the GATK joint genotyping method for calling variants on RNA-seq data was validated by comparing this approach to a so-called “per-sample” method. In addition, pair-wise comparisons of the two methods were performed to evaluate their respective sensitivity, precision and accuracy using DNA genotypes from a companion study including the same 50 cows genotyped using either genotyping-by-sequencing or with the Bovine SNP50 Beadchip (imputed to the Bovine high density). Results indicate that both approaches are very close in their capacity of detecting reference variants and that the joint genotyping method is more sensitive than the per-sample method. Given that the joint genotyping method is more flexible and technically easier, we recommend this approach for variant calling in RNA-seq experiments. Electronic supplementary material The online version of this article (10.1186/s40104-019-0359-0) contains supplementary material, which is available to authorized users.

The chloroplast genome has experienced many architectural changes during the evolution of chlorophyte green algae, with the class Chlorophyceae displaying the lowest degree of ancestral traits. We have previously shown that the completely sequenced chloroplast DNAs (cpDNAs) of Chamydomonas reinhardtii (Chlamydomonadales) and Scenedesmus obliquus (Sphaeropleales) are highly scrambled in gene order relative to one another. Here, we report the complete cpDNA sequence of Stigeoclonium helveticum (Chaetophorales), a member of a third chlorophycean lineage. This genome, which encodes 97 genes and contains 21 introns (including four putatively trans-spliced group II introns inserted at novel sites), is remarkably rich in derived features and extremely rearranged relative to its chlorophycean counterparts. At 223,902 bp, Stigeoclonium cpDNA is the largest chloroplast genome sequenced thus far, and in contrast to those of Chlamydomonas and Scenedesmus, features no large inverted repeat. Interestingly, the pattern of gene distribution between the DNA strands and the bias in base composition along each strand suggest that the Stigeoclonium genome replicates bidirectionally from a single origin. Unlike most known trans-spliced group II introns, those of Stigeoclonium exhibit breaks in domains I and II. By placing our comparative genome analyses in a phylogenetic framework, we inferred an evolutionary scenario of the mutational events that led to changes in genome architecture in the Chlorophyceae.