Long-read assembly and comparative evidence-based reanalysis of Cryptosporidium genome sequences reveals expanded transporter repertoire and duplication of entire chromosome ends including subtelomeric regions

Cryptosporidiosis is a leading cause of waterborne diarrheal disease globally and an important contributor to mortality in infants and the immunosuppressed. Despite its importance, the Cryptosporidium community has only had access to a good, but incomplete, Cryptosporidium parvum IOWA reference genome sequence. Incomplete reference sequences hamper annotation, experimental design and interpretation. We have generated a new C. parvum IOWA genome assembly supported by PacBio and Oxford Nanopore long-read technologies and a new comparative and consistent genome annotation for three closely related species C. parvum, Cryptosporidium hominis and Cryptosporidium tyzzeri. We made 1,926 C. parvum annotation updates based on experimental evidence. They include new transporters, ncRNAs, introns and altered gene structures. The new assembly and annotation revealed a complete Dnmt2 methylase ortholog. Comparative annotation between C. parvum, C. hominis and C. tyzzeri revealed that most "missing" orthologs are found suggesting that the biological differences between the species must result from gene copy number variation, differences in gene regulation and single nucleotide variants (SNVs). Using the new assembly and annotation as reference, 190 genes are identified as evolving under positive selection, including many not detected previously. The new C. parvum IOWA reference genome assembly is larger, gap free and lacks ambiguous bases. This chromosomal assembly recovers all 16 chromosome ends, 13 of which are contiguously assembled. The three remaining chromosome ends are provisionally placed. These ends represent duplication of entire chromosome ends including subtelomeric regions revealing a new level of genome plasticity that will both inform and impact future research.

RTEL1 influences the abundance and localization of TERRA RNA

Telomere repeat containing RNAs (TERRAs) are a family of long non-coding RNAs transcribed from the subtelomeric regions of eukaryotic chromosomes. TERRA transcripts can form R-loops at chromosome ends; however the importance of these structures or the regulation of TERRA expression and retention in telomeric R-loops remain unclear. Here, we show that the RTEL1 (Regulator of Telomere Length 1) helicase influences the abundance and localization of TERRA in human cells. Depletion of RTEL1 leads to increased levels of TERRA RNA while reducing TERRA-containing R loops at telomeres. In vitro, RTEL1 shows a strong preference for binding G-quadruplex structures which form in TERRA. This binding is mediated by the C-terminal region of RTEL1, and is independent of the RTEL1 helicase domain. RTEL1 binding to TERRA appears to be essential for cell viability, underscoring the importance of this function. Degradation of TERRA-containing R-loops by overexpression of RNAse H1 partially recapitulates the increased TERRA levels and telomeric instability associated with RTEL1 deficiency. Collectively, these data suggest that regulation of TERRA is a key function of the RTEL1 helicase, and that loss of that function may contribute to the disease phenotypes of patients with RTEL1 mutations.

Dynamics of SAS-I mediated H4 K16 acetylation during DNA replication in yeast

The acetylation of H4 lysine 16 (H4 K16Ac) in Saccharomyces cerevisiae counteracts the binding of the heterochromatin complex SIR to chromatin and inhibits gene silencing. Contrary to other histone acetylation marks, the H4 K16Ac level is high on genes with low transcription, whereas highly transcribed genes show low H4 K16Ac. Approximately 60% of cellular H4 K16Ac in S. cerevisiae is provided by the SAS-I complex, which consists of the MYST-family acetyltransferase Sas2, Sas4 and Sas5. The absence of SAS-I causes inappropriate spreading of the SIR complex and gene silencing in subtelomeric regions. Here, we investigated the genome-wide dynamics of SAS-I dependent H4 K16Ac during DNA replication. Replication is highly disruptive to chromatin and histone marks, since histones are removed to allow progression of the replication fork, and chromatin is reformed with old and new histones after fork passage. We found that H4 K16Ac appears in chromatin immediately upon replication. Importantly, this increase depends on the presence of functional SAS-I complex. Moreover, the appearance of H4 K16Ac is delayed in genes that are strongly transcribed. This indicates that transcription counteracts SAS-I-mediated H4 K16 acetylation, thus “sculpting” histone modification marks at the time of replication. We furthermore investigated which acetyltransferase acts redundantly with SAS-I to acetylate H4 K16Ac. esa1Δ sds3Δ cells, which were also sas2Δ sir3Δ in order to maintain viability, contained no detectable H4 K16Ac, showing that Esa1 and Sas2 are redundant for cellular H4 K16 acetylation. Furthermore, esa1Δ sds3Δ sas2Δ sir3Δ showed a more pronounced growth defect compared to the already defective esa1Δ sds3Δ sir3Δ. This indicates that SAS-I has cellular functions beyond preventing the spreading of heterochromatin.

Full-length Genome of the Ogataea polymorpha strain HU-11 reveals large duplicated segments in subtelomeric regions

Background: Currently, methylotrophic yeasts (e.g., Pichia pastoris, Hansenula polymorpha, and Candida boindii) are subjects of intense genomics studies in basic research and industrial applications. In the genus Ogataea, most research is focused on three basic O. polymorpha strains—CBS4732, NCYC495, and DL-1. However, these three strains are of independent origin and unclear relationship. As a high-yield engineered O. polymorpha strain, HU-11 can be regarded as identical to CBS4732, because the only difference between them is a 5-bp insertion. Results: In the present study, we have assembled the full-length genome of O. polymorpha HU-11 using high-depth PacBio and Illumina data. Long terminal repeat (LTR) retrotransposons, rDNA, 5' and 3' telomeric, subtelomeric, low complexity and other repeat regions were curated to improve the genome quality. We took advantage of the full-length HU-11 genome sequence for the genome annotation and comparison. Particularly, we determined the exact location of the rDNA genes and LTR retrotransposons in seven chromosomes and detected large duplicated segments in the subtelomeric regions. Three novel findings are: (1) O. polymorpha NCYC495 is so phylogenetically close to CBS4732/HU-11 that the syntenic regions covers nearly 100% of their genomes with a nucleotide identity of 99.5%, while NCYC495 is significantly distinct from DL-1; (2) large segment duplication in subtelomeric regions is the main reason for genome expansion in yeasts; and (3) the duplicated segments in subtelomeric regions may be integrated at telomeric tandem repeats (TRs) through a molecular mechanism, which can be used to develop a simple and highly efficient genome editing system to integrate or cleave large segments into yeast genomes. Conclusions: Our findings provide new opportunities for in-depth understanding of genome evolution in methylotrophic yeasts and lay the foundations for the industrial applications of O. polymorpha HU-11 and CBS4732. The full-length genome of the O. polymorpha strain HU-11 should be included into the NCBI RefSeq database for future studies of O. polymorpha CBS4732, NCYC495, and their derivative strains.

Redistribution of Meiotic Crossovers Along Wheat Chromosomes by Virus-Induced Gene Silencing

Meiotic recombination is the main driver of genetic diversity in wheat breeding. The rate and location of crossover (CO) events are regulated by genetic and epigenetic factors. In wheat, most COs occur in subtelomeric regions but are rare in centromeric and pericentric areas. The aim of this work was to increase COs in both “hot” and “cold” chromosomal locations. We used Virus-Induced gene Silencing (VIGS) to downregulate the expression of recombination-suppressing genes XRCC2 and FANCM and of epigenetic maintenance genes MET1 and DDM1 during meiosis. VIGS suppresses genes in a dominant, transient and non-transgenic manner, which is convenient in wheat, a hard-to-transform polyploid. F1 hybrids of a cross between two tetraploid lines whose genome was fully sequenced (wild emmer and durum wheat), were infected with a VIGS vector ∼ 2 weeks before meiosis. Recombination was measured in F2 seedlings derived from F1-infected plants and non-infected controls. We found significant up and down-regulation of CO rates along subtelomeric regions as a result of silencing either MET1, DDM1 or XRCC2 during meiosis. In addition, we found up to 93% increase in COs in XRCC2-VIGS treatment in the pericentric regions of some chromosomes. Silencing FANCM showed no effect on CO. Overall, we show that CO distribution was affected by VIGS treatments rather than the total number of COs which did not change. We conclude that transient silencing of specific genes during meiosis can be used as a simple, fast and non-transgenic strategy to improve breeding abilities in specific chromosomal regions.

Application of an optimized annotation pipeline to the Cryptococcus deuterogattii genome reveals dynamic primary metabolic gene clusters and genomic impact of RNAi loss

Evaluating the quality of a de novo annotation of a complex fungal genome based on RNA-seq data remains a challenge. In this study, we sequentially optimized a Cufflinks-CodingQuary-based bioinformatics pipeline fed with RNA-seq data using the manually annotated model pathogenic yeasts Cryptococcus neoformans and Cryptococcus deneoformans as test cases. Our results show that the quality of the annotation is sensitive to the quantity of RNA-seq data used and that the best quality is obtained with 5–10 million reads per RNA-seq replicate. We also showed that the number of introns predicted is an excellent a priori indicator of the quality of the final de novo annotation. We then used this pipeline to annotate the genome of the RNAi-deficient species Cryptococcus deuterogattii strain R265 using RNA-seq data. Dynamic transcriptome analysis revealed that intron retention is more prominent in C. deuterogattii than in the other RNAi-proficient species C. neoformans and C. deneoformans. In contrast, we observed that antisense transcription was not higher in C. deuterogattii than in the two other Cryptococcus species. Comparative gene content analysis identified 21 clusters enriched in transcription factors and transporters that have been lost. Interestingly, analysis of the subtelomeric regions in these three annotated species identified a similar gene enrichment, reminiscent of the structure of primary metabolic clusters. Our data suggest that there is active exchange between subtelomeric regions, and that other chromosomal regions might participate in adaptive diversification of Cryptococcus metabolite assimilation potential.

Chromatin modifiers and recombination factors promote a telomere fold-back structure, that is lost during replicative senescence

Telomeres have the ability to adopt a lariat conformation and hence, engage in long and short distance intra-chromosome interactions. Budding yeast telomeres were proposed to fold back into subtelomeric regions, but a robust assay to quantitatively characterize this structure has been lacking. Therefore, it is not well understood how the interactions between telomeres and non-telomeric regions are established and regulated. We employ a telomere chromosome conformation capture (Telo-3C) approach to directly analyze telomere folding and its maintenance in S. cerevisiae. We identify the histone modifiers Sir2, Sin3 and Set2 as critical regulators for telomere folding, which suggests that a distinct telomeric chromatin environment is a major requirement for the folding of yeast telomeres. We demonstrate that telomeres are not folded when cells enter replicative senescence, which occurs independently of short telomere length. Indeed, Sir2, Sin3 and Set2 protein levels are decreased during senescence and their absence may thereby prevent telomere folding. Additionally, we show that the homologous recombination machinery, including the Rad51 and Rad52 proteins, as well as the checkpoint component Rad53 are essential for establishing the telomere fold-back structure. This study outlines a method to interrogate telomere-subtelomere interactions at a single unmodified yeast telomere. Using this method, we provide insights into how the spatial arrangement of the chromosome end structure is established and demonstrate that telomere folding is compromised throughout replicative senescence.

Telomere length dynamics in response to DNA damage in malaria parasites

SUMMARYMalaria remains a major cause of morbidity and mortality within the developing world. Recent work has implicated chromosome end stability and the repair of DNA breaks through telomere healing as potent drivers of variant antigen diversification, thus associating basic mechanisms for maintaining genome integrity with aspects of host-parasite interactions. Here we applied long-read sequencing technology to precisely examine the dynamics of telomere addition and chromosome end stabilization in response to double strand breaks within subtelomeric regions. We observed that the process of telomere healing induces the initial synthesis of telomere repeats well in excess of the minimal number required for end stability. However once stabilized, these newly created telomeres appear to function normally, eventually returning to a length nearing that of intact chromosome ends. These results parallel recent observations in humans, suggesting an evolutionarily conserved mechanism for chromosome end repair.

Application of an optimized annotation pipeline to the Cryptococcus deuterogattii genome reveals dynamic primary metabolic gene clusters and genomic impact of RNAi loss

Evaluating the quality of a de novo annotation of a complex fungal genome based on RNA-seq data remains a challenge. In this study, we sequentially optimized a Cufflinks-CodingQuary based bioinformatics pipeline fed with RNA-seq data using the manually annotated model pathogenic yeasts Cryptococcus neoformans and Cryptococcus deneoformans as test cases. Our results demonstrate that the quality of the annotation is sensitive to the quantity of RNA-seq data used and that the best quality is obtained with 5 to 10 million reads per RNA-seq replicate. We also demonstrated that the number of introns predicted is an excellent a priori indicator of the quality of the final de novo annotation. We then used this pipeline to annotate the genome of the RNAi-deficient species Cryptococcus deuterogattii strain R265 using RNA-seq data. Dynamic transcriptome analysis revealed that intron retention is more prominent in C. deuterogattii than in the other RNAi-proficient species C. neoformans and C. deneoformans. In contrast, we observed that antisense transcription was not higher in C. deuterogattii than in the two other Cryptococcus species. Comparative gene content analysis identified 21 clusters enriched in transcription factors and transporters that have been lost. Interestingly, analysis of the subtelomeric regions in these three annotated species identified a similar gene enrichment, reminiscent of the structure of primary metabolic clusters. Our data suggest that there is active exchange between subtelomeric regions, and that other chromosomal regions might participate in adaptive diversification of Cryptococcus metabolite assimilation potential.

The silencing factor Sir3 is a molecular bridge that sticks together distant loci

ABSTRACTPhysical contacts between distant loci contribute to regulate genome function. However, the molecular mechanisms responsible for settling and maintaining such interactions remain poorly understood. Here we investigate the well conserved interactions between heterochromatin loci. In budding yeast, the 32 telomeres cluster in 3-5 foci in exponentially growing cells. This clustering is functionally linked to the formation of heterochromatin in subtelomeric regions through the recruitment of the silencing complex SIR composed of Sir2/3/4. Combining microscopy and Hi-C on strains expressing different alleles of SIR3, we show that the binding of Sir3 directly promotes long range contacts between distant regions, including the rDNA, telomeres, and internal Sir3 bound sites. Furthermore, we unveil a new property of Sir3 in promoting rDNA compaction. Finally, using a synthetic approach we demonstrate that Sir3 can bond loci belonging to different chromosomes together, when targeted to these loci, independently of its interaction with its known partners (Rap1, and Sir4), Sir2 activity or chromosome context. Altogether these data suggest that Sir3 represents an uncommon example of protein able to bridge directly distant loci.