Long-read assembly of a Great Dane genome highlights the contribution of GC-rich sequence and mobile elements to canine genomes.

Technological advances have allowed improvements in genome reference sequence assemblies. Here, we combined long- and short-read sequence resources to assemble the genome of a female Great Dane dog. This assembly has improved continuity compared to the existing Boxer-derived (CanFam3.1) reference genome. Annotation of the Great Dane assembly identified 22,182 protein-coding gene models and 7,049 long noncoding RNAs, including 49 protein-coding genes not present in the CanFam3.1 reference. The Great Dane assembly spans the majority of sequence gaps in the CanFam3.1 reference and illustrates that 2,151 gaps overlap the transcription start site of a predicted protein-coding gene. Moreover, a subset of the resolved gaps, which have an 80.95% median GC content, localize to transcription start sites and recombination hotspots more often than expected by chance, suggesting the stable canine recombinational landscape has shaped genome architecture. Alignment of the Great Dane and CanFam3.1 assemblies identified 16,834 deletions and 15,621 insertions, as well as 2,665 deletions and 3,493 insertions located on secondary contigs. These structural variants are dominated by retrotransposon insertion/deletion polymorphisms and include 16,221 dimorphic canine short interspersed elements (SINECs) and 1,121 dimorphic long interspersed element-1 sequences (LINE-1_Cfs). Analysis of sequences flanking the 3' end of LINE-1_Cfs (i.e., LINE-1_Cf 3'-transductions) suggests multiple retrotransposition-competent LINE-1_Cfs segregate among dog populations. Consistent with this conclusion, we demonstrate that a canine LINE-1_Cf element with intact open reading frames can retrotranspose its own RNA and that of a SINEC_Cf consensus sequence in cultured human cells, implicating ongoing retrotransposon activity as a driver of canine genetic variation.

Cell-line specific RNA editing patterns in Trypanosoma brucei suggest a unique mechanism to generate protein variation in a system intolerant to genetic mutations.

Trypanosoma brucei possesses a highly complex RNA editing system that uses guide RNAs to direct the insertion and deletion of uridines in mitochondrial mRNAs. These changes extensively alter the target mRNAs and can more than double them in length. Recently, analyses showed that several of the edited genes possess the capacity to encode two different protein products. The overlapped reading frames can be accessed through alternative RNA editing that shifts the translated reading frame. In this study, we analyzed the editing patterns of three putative dual-coding genes, ribosomal protein S12 (RPS12), the 5' editing domain of NADH dehydrogenase subunit 7 (ND7 5'), and C-rich region 3 (CR3). We found evidence that alternatively 5'-edited ND7 5' and CR3 transcripts are present in the transcriptome, providing evidence for the use of dual ORFs in these transcripts. Moreover, we found that CR3 has a complex set of editing pathways that vary substantially between cell lines. These findings suggest that alternative editing can work to introduce genetic variation in a system that selects against nucleotide mutations.

Intrinsic and regulated properties of minimally edited trypanosome mRNAs.

Most mitochondrial mRNAs in kinetoplastids require extensive uridine insertion/deletion editing to generate translatable open reading frames. Editing is specified by trans-acting gRNAs and involves a complex machinery including basal and accessory factors. Here, we utilize high-throughput sequencing to analyze editing progression in two minimally edited mRNAs that provide a simplified system due their requiring only two gRNAs each for complete editing. We show that CYb and MURF2 mRNAs exhibit barriers to editing progression that differ from those previously identified for pan-edited mRNAs, primarily at initial gRNA usage and gRNA exchange. We demonstrate that mis-edited junctions arise through multiple pathways including mis-alignment of cognate gRNA, incorrect and sometimes promiscuous gRNA utilization and inefficient gRNA anchoring. We then examined the roles of accessory factors RBP16 and MRP1/2 in maintaining edited CYb and MURF2 populations. RBP16 is essential for initiation of CYb and MURF2 editing, as well as MURF2 editing progression. In contrast, MRP1/2 stabilizes both edited mRNA populations, while further promoting progression of MURF2 mRNA editing. We also analyzed the effects of RNA Editing Substrate Binding Complex components, TbRGG2 and GAP1, and show that both proteins modestly impact progression of editing on minimally edited mRNAs, suggesting a novel function for GAP1.

Mitochondrial dual-coding genes in Trypanosoma brucei.

Trypanosoma brucei is transmitted between mammalian hosts by the tsetse fly. In the mammal, they are exclusively extracellular, continuously replicating within the bloodstream. During this stage, the mitochondrion lacks a functional electron transport chain (ETC). Successful transition to the fly, requires activation of the ETC and ATP synthesis via oxidative phosphorylation. This life cycle leads to a major problem: in the bloodstream, the mitochondrial genes are not under selection and are subject to genetic drift that endangers their integrity. Exacerbating this, T. brucei undergoes repeated population bottlenecks as they evade the host immune system that would create additional forces of genetic drift. These parasites possess several unique genetic features, including RNA editing of mitochondrial transcripts. RNA editing creates open reading frames by the guided insertion and deletion of U-residues within the mRNA. A major question in the field has been why this metabolically expensive system of RNA editing would evolve and persist. Here, we show that many of the edited mRNAs can alter the choice of start codon and the open reading frame by alternative editing of the 5' end. Analyses of mutational bias indicate that six of the mitochondrial genes may be dual-coding and that RNA editing allows access to both reading frames. We hypothesize that dual-coding genes can protect genetic information by essentially hiding a non-selected gene within one that remains under selection. Thus, the complex RNA editing system found in the mitochondria of trypanosomes provides a unique molecular strategy to combat genetic drift in non-selective conditions.

Analysis of the Trypanosoma brucei EATRO 164 Bloodstream Guide RNA Transcriptome.

The mitochondrial genome of Trypanosoma brucei contains many cryptogenes that must be extensively edited following transcription. The RNA editing process is directed by guide RNAs (gRNAs) that encode the information for the specific insertion and deletion of uridylates required to generate translatable mRNAs. We have deep sequenced the gRNA transcriptome from the bloodstream form of the EATRO 164 cell line. Using conventionally accepted fully edited mRNA sequences, ~1 million gRNAs were identified. In contrast, over 3 million reads were identified in our insect stage gRNA transcriptome. A comparison of the two life cycle transcriptomes show an overall ratio of procyclic to bloodstream gRNA reads of 3.5:1. This ratio varies significantly by gene and by gRNA populations within genes. The variation in the abundance of the initiating gRNAs for each gene, however, displays a trend that correlates with the developmental pattern of edited gene expression. A comparison of related major classes from each transcriptome revealed a median value of ten single nucleotide variations per gRNA. Nucleotide variations were much less likely to occur in the consecutive Watson-Crick anchor region, indicating a very strong bias against G:U base pairs in this region. This work indicates that gRNAs are expressed during both life cycle stages, and that differential editing patterns observed for the different mitochondrial mRNA transcripts are not due to the presence or absence of gRNAs. However, the abundance of certain gRNAs may be important in the developmental regulation of RNA editing.