Hierarchical length and sequence preferences establish a single major piRNA 3'-end

PIWI-interacting RNAs (piRNAs) guard germline genomes against the deleterious action of retroviruses and other mobile genetic elements. How piRNAs faithfully discriminate between self and non-self to restrict all mobile elements while sparing essential genes remains a key outstanding question in genome biology. PiRNAs use extensive base-pairing to recognize their targets and variable 3'ends could change the specificity and efficacy of piRNA silencing. Here, we identify conserved rules that ensure the generation of a single major piRNA 3'end in flies and mice. Our data suggest that the PIWI proteins initially define a short interval on pre-piRNAs that grants access to the ZUC-processor complex. Within this Goldilocks zone, the preference of the ZUC-processor to cut in front of a Uridine determines the ultimate processing site. We observe a mouse-specific roadblock that relocates the Goldilocks zone and generates an opportunity for consecutive trimming by PNLDC1. Our data reveal a conserved hierarchy between length and sequence preferences that controls the piRNA sequence space. The unanticipated precision of 3'end formation bolsters the emerging understanding that the functional piRNA sequence space is tightly controlled to ensure effective defense.

Characterising genome architectures using Genome Decomposition Analysis

Genome architecture describes how genes and other features are arranged in genomes. These arrangements reflect the evolutionary pressures on genomes and underlie biological processes such as chromosomal segregation and the regulation of gene expression. We present a new tool called Genome Decomposition Analysis (GDA) that characterises genome architectures and acts as an accessible approach for discovering hidden features of a genome assembly. With the imminent deluge of high quality genome assemblies from projects such as the Darwin Tree of Life and the Earth BioGenome Project, GDA has been designed to facilitate their exploration and the discovery of novel genome biology. We highlight the effectiveness of our approach in characterising the genome architectures of single-celled eukaryotic parasites from the phylum Apicomplexa and show that it scales well to large genomes.SignificanceGenome sequencing has revealed that there are functionally important arrangements of genes, repetitive elements and regulatory sequences within chromosomes. Identifying these arrangements requires extensive computation and analysis. Furthermore, improvements in genome sequencing technology and the establishment of consortia aiming to sequence all species of eukaryotes mean that there is a need for high throughput methods for discovering new genome biology. Here we present a software pipeline, named GDA, which determines the patterns of genomic features across chromosomes and uses these to characterise genome architecture. We show that it recapitulates the known genome architecture of several Apicomplexan parasites and use it to identify features in a recently sequenced, less well-characterised genome. GDA scales well to large genomes and is freely available.

Long reads and Hi-C sequencing illuminate the two-compartment genome of the model arbuscular mycorrhizal symbiont Rhizophagus irregularis

Chromosome folding links genome structure with gene function by generating distinct nuclear compartments and topologically associating domains (TADs). In mammals, these domains undergo preferential interactions and regulate gene expression, however in fungi the role of chromosome folding in genome biology is unclear. Here, we combine Nanopore (ONT) sequencing with chromatin conformation capture sequencing (Hi-C) to reveal chromosome diversity in a group of obligate plant symbionts with a multinucleate mycelium; the arbuscular mycorrhizal fungi (AMF). We find that phylogenetically distinct strains of the model AMF Rhizophagus irregularis all carry 33 chromosomes. Homologous chromosomes show within species variability in size, as well as in gene and repeat content. Strain-specific Hi-C sequencing reveals that all strains have a 3D genome organization that resembles a checkerboard structure with two distinct (A/B) chromatin compartments. Each compartment differs in the level of gene transcription, regulation of candidate effectors and methylation rate. The A-compartment is more gene-dense and contains most core genes, while the B-compartment is more repeat-rich and has higher rates of chromosomal rearrangement. While the B-compartment is transcriptionally repressed, it has significantly more secreted proteins and in planta up-regulated candidate effectors suggesting a possible host-induced change in chromosome conformation. Overall, this study provides a fine-scale view into the genome biology and evolution of prominent plant symbionts, and opens avenues to study the mechanisms that generate and modify chromosome folding during host-microbe interactions.

Re-assembly of nineteenth-century smallpox vaccine genomes reveals the contemporaneous use of horsepox and horsepox-related viruses in the USA

According to a recent article published in Genome Biology, Duggan and coworkers sequenced and partially assembled five genomes of smallpox vaccines from the nineteenth century. No information regarding the ends of genomes was presented, and they are important to understand the evolutionary relationship of the different smallpox vaccine genomes during the centuries. We re-assembled the genomes, which include the largest genomes in the vaccinia lineage and one true horsepox strain. Moreover, the assemblies reveal a diverse genetic structure in the genome ends. Our data emphasize the concurrent use of horsepox and horsepox-related viruses as the smallpox vaccine in the nineteenth century.

IDR2D identifies reproducible genomic interactions

Chromatin interaction data from protocols such as ChIA-PET, HiChIP and Hi-C provide valuable insights into genome organization and gene regulation, but can include spurious interactions that do not reflect underlying genome biology. We introduce an extension of the Irreproducible Discovery Rate (IDR) method called IDR2D that identifies replicable interactions shared by chromatin interaction experiments. IDR2D provides a principled set of interactions and eliminates artifacts from single experiments. The method is available as a Bioconductor package for the R community, as well as an online service at https://idr2d.mit.edu.

Longest protein, longest transcript or most expression, for accurate gene reconstruction of transcriptomes?

Methods of transcript assembly and reduction filters are compared for recovery of reference gene sets of human, pig and plant, including longest coding sequence with EvidentialGene, longest transcript with CD-HIT, and most RNA-seq with TransRate. EvidentialGene methods are the most accurate in recovering reference genes, and maintain accuracy for alternate transcripts and paralogs. In comparison, filtering large over-assemblies by longest RNA measures, and most RNA-seq expression measures, discards a large portion of accurate models, especially alternates and paralogs. Accuracy of protein calculations is compared, with errors found in popular methods, as is accuracy of transcript assemblers. Gene reconstruction accuracy depends upon the underlying measurements, where protein criteria, including homology among species, have the strength of evolutionary biology that other criteria lack. EvidentialGene provides a gene reconstruction algorithm that is consistent with genome biology.