Genomic Distance Research Articles

The Isolation and Genome Sequencing of Five Novel Bacteriophages From the Rumen Active Against Butyrivibrio fibrisolvens.

Although the prokaryotic communities of the rumen microbiome are being uncovered through genome sequencing, little is known about the resident viral populations. Whilst temperate phages can be predicted as integrated prophages when analyzing bacterial and archaeal genomes, the genetics underpinning lytic phages remain poorly characterized. To the five genomes of bacteriophages isolated from rumen-associated samples sequenced and analyzed previously, this study adds a further five novel genomes and predictions gleaned from them to further the understanding of the rumen phage population. Lytic bacteriophages isolated from fresh ovine and bovine fecal and rumen fluid samples were active against the predominant fibrolytic ruminal bacterium Butyrivibrio fibrisolvens. The double stranded DNA genomes were sequenced and reconstructed into single circular complete contigs. Based on sequence similarity and genome distances, the five phages represent four species from three separate genera, consisting of: (1) Butyrivibrio phages Arian and Bo-Finn; (2) Butyrivibrio phages Idris and Arawn; and (3) Butyrivibrio phage Ceridwen. They were predicted to all belong to the Siphoviridae family, based on evidence in the genomes such as size, the presence of the tail morphogenesis module, genes that share similarity to those in other siphovirus isolates and phylogenetic analysis using phage proteomes. Yet, phylogenomic analysis and sequence similarity of the entire phage genomes revealed that these five phages are unique and novel. These phages have only been observed undergoing the lytic lifecycle, but there is evidence in the genomes of phages Arawn and Idris for the potential to be temperate. However, there is no evidence in the genome of the bacterial host Butyrivibrio fibrisolvens of prophage genes or genes that share similarity with the phage genomes.

The geno-spatio analysis of Mycobacterium tuberculosis complex in hot and cold spots of Guangxi, China

BackgroundAt present, there are few studies on polymorphism of Mycobacterium tuberculosis (Mtb) gene and how it affects the TB epidemic. This study aimed to document the differences of polymorphisms between tuberculosis hot and cold spot areas of Guangxi Zhuang Autonomous Region, China.MethodsThe cold and hot spot areas, each with 3 counties, had been pre-identified by TB incidence for 5 years from the surveillance database. Whole genome sequencing analysis was performed on all sputum Mtb isolates from the detected cases during January and June 2018. Single nucleotide polymorphism (SNP) of each isolate compared to the H37Rv strain were called and used for lineage and sub-lineage identification. Pairwise SNP differences between every pair of isolates were computed. Analyses of Molecular Variance (AMOVA) across counties of the same hot or cold spot area and between the two areas were performed.ResultsAs a whole, 59.8% (57.7% sub-lineage 2.2 and 2.1% sub-lineage 2.1) and 39.8% (17.8% sub-lineage 4.4, 6.5% sub-lineage 4.2 and 15.5% sub-lineage 4.5) of the Mtb strains were Lineage 2 and Lineage 4 respectively. The percentages of sub-lineage 2.2 (Beijing family strains) are significantly higher in hot spots. Through the MDS dimension reduction, the genomic population structure in the three hot spot counties is significantly different from those three cold spot counties (T-test p = 0.05). The median of SNPs distances among Mtb isolates in cold spots was greater than that in hot spots (897 vs 746, Rank-sum test p < 0.001). Three genomic clusters, each with genomic distance ≤12 SNPs, were identified with 2, 3 and 4 consanguineous strains. Two clusters were from hot spots and one was from cold spots.ConclusionNarrower genotype diversity in the hot area may indicate higher transmissibility of the Mtb strains in the area compared to those in the cold spot area.

Identification and Analysis of p53-Regulated Enhancers in Hepatic Carcinoma.

Enhancers can act as cis-regulatory elements to control transcriptional regulation by recruiting transcription factors (TFs) in a distance and orientation-independent manner. However, it is still unclear how p53 participates in the enhancer network as TF in hepatic carcinoma under the condition of DNA damage. A total of 14,286 active enhancers were identified through the integration of stable and unstable enhancer RNAs (eRNAs) captured by CAGE and GRO-seq, respectively. Furthermore, 218 p53-bound enhancers (Enhp53) were identified by analyzing p53 ChIP-seq in HepG2 cells after DNA damage. The results showed that the enhancer expression and histone markers of enhancers (H3K4me1, H3K4me2, H3K4me3, H3K9ac, and H3K27ac) revealed significantly higher level on Enhp53 than Enhno−p53 which suggested that p53 participated in regulating enhancer activity and chromatin structure. By analyzing 124 TFs ChIP-seq from ENCODE, 93 TFs were found significantly enriched on Enhp53 such as GATA4, YY1, and CTCF, indicating p53 may co-regulate enhancers with TFs participation. Moreover, significantly differentially expressed 438 miRNAs and 1,264 mRNAs were identified by analyzing small RNA-seq and RNA-seq, and 26 Enhp53-miRNAs and 145 Enhp53-mRNA interactions were identified by the integration of 3D genome data and genomic distance. The functional enrichment analysis showed that these miRNA targets and mRNAs were significantly involved in tumor biological processes and signaling pathways such as DNA replication, p53 signaling pathway, hepatitis B, focal adhesion, etc. The above results indicated that p53 participated in regulating enhancer network in hepatic carcinoma and Enhp53 exhibited significantly different characteristics with Enhno−p53.

Genome-Wide Identification and Analysis of Enhancer-Regulated microRNAs Across 31 Human Cancers.

Enhancers are cis-regulatory DNA elements that positively regulate the transcription of target genes in a tissue-specific manner, and dysregulation of target genes could lead to various diseases, such as cancer. Recent studies have shown that enhancers can regulate microRNAs (miRNAs) and participate in their biological synthesis. However, the network of enhancer-regulated miRNAs across multiple cancers is still unclear. Here, a total of 2,418 proximal enhancer–miRNA interactions and 1,280 distal enhancer–miRNA interactions were identified through the integration of genomic distance, co-expression, and 3D genome data in 31 cancers. The results showed that both proximal and distal interactions exhibited a significant cancer type-specific feature trend at the tissue level rather than at the single-cell level, and there was a noteworthy positive correlation between the expression of the miRNA and the number of enhancers regulating the same miRNA in most cancers. Furthermore, we found that there was a high correlation between the formation of enhancer–miRNA pairs and the expression of enhancer RNAs (eRNAs) whether in distal or proximal regulation. The characteristics analysis showed that miRes (enhancers that regulated miRNAs) and non-miRes presented significant differences in sequence conservation, guanine–cytosine (GC) content, and histone modification signatures. Notably, GC content, H3K4me1, and H3K36me3 were present differently between distal and proximal regulation, suggesting that they might participate in chromosome looping of enhancer–miRNA interactions. Finally, we introduced a case study, enhancer: chr1:1186391–1186507 ∼ miR-200a was highly relevant to the survival of thyroid cancer patients and a cis-eQTL SNP on the enhancer affected the expression of the TNFRSF18 gene as a tumor suppressor.

Genomic nucleotide-based distance analysis for delimiting old world monkey derived herpes simplex virus species

BackgroundHerpes simplex viruses form a genus within the alphaherpesvirus subfamily, with three identified viral species isolated from Old World monkeys (OWM); Macacine alphaherpesvirus 1 (McHV-1; herpes B), Cercopithecine alphaherpesvirus 2 (SA8), and Papiine alphaherpesvirus 2 (PaHV-2; herpes papio). Herpes B is endemic to macaques, while PaHV-2 and SA8 appear endemic to baboons. All three viruses are genetically and antigenically similar, with SA8 and PaHV-2 thought to be avirulent in humans, while herpes B is a biosafety level 4 pathogen. Recently, next-generation sequencing (NGS) has resulted in an increased number of published OWM herpes simplex genomes, allowing an encompassing phylogenetic analysis.ResultsIn this study, phylogenetic networks, in conjunction with a genome-based genetic distance cutoff method were used to examine 27 OWM monkey herpes simplex isolates. Genome-based genetic distances were calculated, resulting in distances between lion and pig-tailed simplex viruses themselves, and versus herpes B core strains that were higher than those between PaHV-2 and SA8 (approximately 14 and 10% respectively). The species distance cutoff was determined to be 8.94%, with the method recovering separate species status for PaHV-2 and SA8 and showed that lion and pig-tailed simplex viruses (vs core herpes B strains) were well over the distance species cutoff.ConclusionsWe propose designating lion and pig-tailed simplex viruses as separate, individual viral species, and that this may be the first identification of viral cryptic species.

Pseudomonas izuensis sp. nov., a novel species isolated from Izu Oshima, Japan.

An aerobic, Gram-stain-negative, rod-shaped bacterial strain, IzPS43_3003T, was isolated from Izu Oshima, an active volcanic island located 22 km east of the Izu Peninsula, Japan. The sequence of its 16S rRNA gene indicated that IzPS43_3003T belongs to the Pseudomonas fluorescens lineage, with its sequence being most similar to that of Pseudomonas vancouverensis DhA-51T (99.79 %). Phylogenetic analysis based on whole genome sequences showed that IzPS43_3003T was a member of the Pseudomonas jessenii subgroup. The average nucleotide identity values and genome-to genome distances between the whole genome sequences of IzPS43_3003T and other type strains showed that the highest correlations were with Pseudomonas moorei DSM 12647T (87.3 and 33.5% respectively). These genotypic and phenotypic analyses indicated that IzPS43_3003T belongs to a novel species, Pseudomonas izuensis sp. nov. Its type strain is IzPS43_3003T (=LMG 31527T,=CECT 9963T).

Uncovering Sub-Structure and Genomic Profiles in Across-Countries Subpopulations of Angus Cattle

Highlighting genomic profiles for geographically distinct subpopulations of the same breed may provide insights into adaptation mechanisms to different environments, reveal genomic regions divergently selected, and offer initial guidance to joint genomic analysis. Here, we characterized similarities and differences between the genomic patterns of Angus subpopulations, born and raised in Canada (N = 382) and Brazil (N = 566). Furthermore, we systematically scanned for selection signatures based on the detection of autozygosity islands common between the two subpopulations, and signals of divergent selection, via FST and varLD tests. The principal component analysis revealed a sub-structure with a close connection between the two subpopulations. The averages of genomic relationships, inbreeding coefficients, and linkage disequilibrium at varying genomic distances were rather similar across them, suggesting non-accentuated differences in overall genomic diversity. Autozygosity islands revealed selection signatures common to both subpopulations at chromosomes 13 (63.77–65.25 Mb) and 14 (22.81–23.57 Mb), which are notably known regions affecting growth traits. Nevertheless, further autozygosity islands along with FST and varLD tests unravel particular sites with accentuated population subdivision at BTAs 7 and 18 overlapping with known QTL and candidate genes of reproductive performance, thermoregulation, and resistance to infectious diseases. Our findings indicate overall genomic similarity between Angus subpopulations, with noticeable signals of divergent selection in genomic regions associated with the adaptation in different environments.

The corrected gene proximity map for analyzing the 3D genome organization using Hi-C data

BackgroundGenome-wide ligation-based assays such as Hi-C provide us with an unprecedented opportunity to investigate the spatial organization of the genome. Results of a typical Hi-C experiment are often summarized in a chromosomal contact map, a matrix whose elements reflect the co-location frequencies of genomic loci. To elucidate the complex structural and functional interactions between those genomic loci, networks offer a natural and powerful framework.ResultsWe propose a novel graph-theoretical framework, the Corrected Gene Proximity (CGP) map to study the effect of the 3D spatial organization of genes in transcriptional regulation. The starting point of the CGP map is a weighted network, the gene proximity map, whose weights are based on the contact frequencies between genes extracted from genome-wide Hi-C data. We derive a null model for the network based on the signal contributed by the 1D genomic distance and use it to “correct” the gene proximity for cell type 3D specific arrangements. The CGP map, therefore, provides a network framework for the 3D structure of the genome on a global scale. On human cell lines, we show that the CGP map can detect and quantify gene co-regulation and co-localization more effectively than the map obtained by raw contact frequencies. Analyzing the expression pattern of metabolic pathways of two hematopoietic cell lines, we find that the relative positioning of the genes, as captured and quantified by the CGP, is highly correlated with their expression change. We further show that the CGP map can be used to form an inter-chromosomal proximity map that allows large-scale abnormalities, such as chromosomal translocations, to be identified.ConclusionsThe Corrected Gene Proximity map is a map of the 3D structure of the genome on a global scale. It allows the simultaneous analysis of intra- and inter- chromosomal interactions and of gene co-regulation and co-localization more effectively than the map obtained by raw contact frequencies, thus revealing hidden associations between global spatial positioning and gene expression. The flexible graph-based formalism of the CGP map can be easily generalized to study any existing Hi-C datasets.

Tug of war in a double-nanopore system.

We simulate a tug-of-war (TOW) scenario for a model double-stranded DNA threading through a double nanopore (DNP) system. The DNA, simultaneously captured at both pores, is subject to two equal and opposite forces -f[over ⃗]_{L}=f[over ⃗]_{R} (TOW), where f[over ⃗]_{L} and f[over ⃗]_{R} are the forces applied to the left and the right pore, respectively. Even though the net force on the DNA polymer Δf[over ⃗]_{LR}=f[over ⃗]_{L}+f[over ⃗]_{R}=0, the mean first passage time (MFPT) 〈τ〉 depends on the magnitude of the TOW forces |f_{L}|=|f_{R}|=f_{LR}. We qualitatively explain this dependence of 〈τ〉 on f_{LR} from the known results for the single-pore translocation of a triblock copolymer A-B-A with ℓ_{pB}>ℓ_{pA}, where ℓ_{pA} and ℓ_{pB} are the persistence length of the A and B segments, respectively. We demonstrate that the time of flight of a monomer with index m[〈τ_{LR}(m)〉] from one pore to the other exhibits quasiperiodic structure commensurate with the distance between the pores d_{LR}. Finally, we study the situation where we offset the TOW biases so that Δf[over ⃗]_{LR}=f[over ⃗]_{L}+f[over ⃗]_{R}≠0, and qualitatively reproduce the experimental result of the dependence of the MFPT on Δf[over ⃗]_{LR}. We demonstrate that, for a moderate bias, the MFPT for the DNP system for a chain length N follows the same scaling ansatz as that for the single nanopore, 〈τ〉=(AN^{1+ν}+η_{pore}N)(Δf_{LR})^{-1}, where η_{pore} is the pore friction, which enables us to estimate 〈τ〉 for a long chain. Our Brownian dynamics simulation studies provide fundamental insights and valuable information about the details of the translocation speed obtained from 〈τ_{LR}(m)〉, and accuracy of the translation of the data obtained in the time domain to units of genomic distances.

Determinants of transcription factor regulatory range

Characterization of the genomic distances over which transcription factor (TF) binding influences gene expression is important for inferring target genes from TF chromatin immunoprecipitation followed by sequencing (ChIP-seq) data. Here we systematically examine the relationship between thousands of TF and histone modification ChIP-seq data sets with thousands of gene expression profiles. We develop a model for integrating these data, which reveals two classes of TFs with distinct ranges of regulatory influence, chromatin-binding preferences, and auto-regulatory properties. We find that the regulatory range of the same TF bound within different topologically associating domains (TADs) depend on intrinsic TAD properties such as local gene density and G/C content, but also on the TAD chromatin states. Our results suggest that considering TF type, binding distance to gene locus, as well as chromatin context is important in identifying implicated TFs from GWAS SNPs.

Phylogenetic Reconstruction Based on Synteny Block and Gene Adjacencies

Gene order can be used as an informative character to reconstruct phylogenetic relationships between species independently from the local information present in gene/protein sequences. PhyChro is a reconstruction method based on chromosomal rearrangements, applicable to a wide range of eukaryotic genomes with different gene contents and levels of synteny conservation. For each synteny breakpoint issued from pairwise genome comparisons, the algorithm defines two disjoint sets of genomes, named partial splits, respectively, supporting the two block adjacencies defining the breakpoint. Considering all partial splits issued from all pairwise comparisons, a distance between two genomes is computed from the number of partial splits separating them. Tree reconstruction is achieved through a bottom-up approach by iteratively grouping sister genomes minimizing genome distances. PhyChro estimates branch lengths based on the number of synteny breakpoints and provides confidence scores for the branches. PhyChro performance is evaluated on two data sets of 13 vertebrates and 21 yeast genomes by using up to 130,000 and 179,000 breakpoints, respectively, a scale of genomic markers that has been out of reach until now. PhyChro reconstructs very accurate tree topologies even at known problematic branching positions. Its robustness has been benchmarked for different synteny block reconstruction methods. On simulated data PhyChro reconstructs phylogenies perfectly in almost all cases, and shows the highest accuracy compared with other existing tools. PhyChro is very fast, reconstructing the vertebrate and yeast phylogenies in <15 min.

HiC-Hiker: a probabilistic model to determine contig orientation in chromosome-length scaffolds with Hi-C.

Motivation De novo assembly of reference-quality genomes used to require enormously laborious tasks. In particular, it is extremely time-consuming to build genome markers for ordering assembled contigs along chromosomes; thus, they are only available for well-established model organisms. To resolve this issue, recent studies demonstrated that Hi-C could be a powerful and cost-effective means to output chromosome-length scaffolds for non-model species with no genome marker resources, because the Hi-C contact frequency between a pair of two loci can be a good estimator of their genomic distance, even if there is a large gap between them. Indeed, state-of-the-art methods such as 3D-DNA are now widely used for locating contigs in chromosomes. However, it remains challenging to reduce errors in contig orientation because shorter contigs have fewer contacts with their neighboring contigs. These orientation errors lower the accuracy of gene prediction, read alignment, and synteny block estimation in comparative genomics.ResultsTo reduce these contig orientation errors, we propose a new algorithm, named HiC-Hiker, which has a firm grounding in probabilistic theory, rigorously models Hi-C contacts across contigs, and effectively infers the most probable orientations via the Viterbi algorithm. We compared HiC-Hiker and 3D-DNA using human and worm genome contigs generated from short reads, evaluated their performances, and observed a remarkable reduction in the contig orientation error rate from 4.3% (3D-DNA) to 1.7% (HiC-Hiker). Our algorithm can consider long-range information between distal contigs and precisely estimates Hi-C read contact probabilities among contigs, which may also be useful for determining the ordering of contigs.Availability and implementationHiC-Hiker is freely available at: https://github.com/ryought/hic_hiker.

Serpentine: a flexible 2D binning method for differential Hi-C analysis.

MotivationHi-C contact maps reflect the relative contact frequencies between pairs of genomic loci, quantified through deep sequencing. Differential analyses of these maps enable downstream biological interpretations. However, the multi-fractal nature of the chromatin polymer inside the cellular envelope results in contact frequency values spanning several orders of magnitude: contacts between loci pairs separated by large genomic distances are much sparser than closer pairs. The same is true for poorly covered regions, such as repeated sequences. Both distant and poorly covered regions translate into low signal-to-noise ratios. There is no clear consensus to address this limitation.ResultsWe present Serpentine, a fast, flexible procedure operating on raw data, which considers the contacts in each region of a contact map. Binning is performed only when necessary on noisy regions, preserving informative ones. This results in high-quality, low-noise contact maps that can be conveniently visualized for rigorous comparative analyses.Availability and implementationSerpentine is available on the PyPI repository and https://github.com/koszullab/serpentine; documentation and tutorials are provided at https://serpentine.readthedocs.io/en/latest/.Supplementary information Supplementary data are available at Bioinformatics online.

T-Gene: improved target gene prediction.

Identifying the genes regulated by a given transcription factor (TF) (its 'target genes') is a key step in developing a comprehensive understanding of gene regulation. Previously, we developed a method (CisMapper) for predicting the target genes of a TF based solely on the correlation between a histone modification at the TF's binding site and the expression of the gene across a set of tissues or cell lines. That approach is limited to organisms for which extensive histone and expression data are available, and does not explicitly incorporate the genomic distance between the TF and the gene. We present the T-Gene algorithm, which overcomes these limitations. It can be used to predict which genes are most likely to be regulated by a TF, and which of the TF's binding sites are most likely involved in regulating particular genes. T-Gene calculates a novel score that combines distance and histone/expression correlation, and we show that this score accurately predicts when a regulatory element bound by a TF is in contact with a gene's promoter, achieving median precision above 60%. T-Gene is easy to use via its web server or as a command-line tool, and can also make accurate predictions (median precision above 40%) based on distance alone when extensive histone/expression data is not available for the organism. T-Gene provides an estimate of the statistical significance of each of its predictions. The T-Gene web server, source code, histone/expression data and genome annotation files are provided at http://meme-suite.org. Supplementary data are available at Bioinformatics online.

Evaluating Enhancer Function and Transcription.

Cell-type- and condition-specific profiles of gene expression require coordination between protein-coding gene promoters and cis-regulatory sequences called enhancers. Enhancers can stimulate gene activity at great genomic distances from their targets, raising questions about how enhancers communicate with specific gene promoters and what molecular mechanisms underlie enhancer function. Characterization of enhancer loci has identified the molecular features of active enhancers that accompany the binding of transcription factors and local opening of chromatin. These characteristics include coactivator recruitment, histone modifications, and noncoding RNA transcription. However, it remains unclear which of these features functionally contribute to enhancer activity. Here, we discuss what is known about how enhancers regulate their target genes and how enhancers and promoters communicate. Further, we describe recent data demonstrating many similarities between enhancers and the gene promoters they control, and we highlight unanswered questions in the field, such as the potential roles of transcription at enhancers.

Sources of Multidrug Resistance in Patients With Previous Isoniazid-Resistant Tuberculosis Identified Using Whole Genome Sequencing: A Longitudinal Cohort Study.

BackgroundMeta-analysis of patients with isoniazid-resistant tuberculosis (TB) given standard first-line anti-TB treatment indicated an increased risk of multidrug-resistant TB (MDR-TB) emerging (8%), compared to drug-sensitive TB (0.3%). Here we use whole genome sequencing (WGS) to investigate whether treatment of patients with preexisting isoniazid-resistant disease with first-line anti-TB therapy risks selecting for rifampicin resistance, and hence MDR-TB.MethodsPatients with isoniazid-resistant pulmonary TB were recruited and followed up for 24 months. Drug susceptibility testing was performed by microscopic observation drug susceptibility assay, mycobacterial growth indicator tube, and by WGS on isolates at first presentation and in the case of re-presentation. Where MDR-TB was diagnosed, WGS was used to determine the genomic relatedness between initial and subsequent isolates. De novo emergence of MDR-TB was assumed where the genomic distance was 5 or fewer single-nucleotide polymorphisms (SNPs), whereas reinfection with a different MDR-TB strain was assumed where the distance was 10 or more SNPs.ResultsTwo hundred thirty-nine patients with isoniazid-resistant pulmonary TB were recruited. Fourteen (14/239 [5.9%]) patients were diagnosed with a second episode of TB that was multidrug resistant. Six (6/239 [2.5%]) were identified as having evolved MDR-TB de novo and 6 as having been reinfected with a different strain. In 2 cases, the genomic distance was between 5 and 10 SNPs and therefore indeterminate.ConclusionsIn isoniazid-resistant TB, de novo emergence and reinfection of MDR-TB strains equally contributed to MDR development. Early diagnosis and optimal treatment of isoniazid-resistant TB are urgently needed to avert the de novo emergence of MDR-TB during treatment.

The DLO Hi-C Tool for Digestion-Ligation-Only Hi-C Chromosome Conformation Capture Data Analysis.

It is becoming increasingly important to understand the mechanism of regulatory elements on target genes in long-range genomic distance. 3C (chromosome conformation capture) and its derived methods are now widely applied to investigate three-dimensional (3D) genome organizations and gene regulation. Digestion-ligation-only Hi-C (DLO Hi-C) is a new technology with high efficiency and cost-effectiveness for whole-genome chromosome conformation capture. Here, we introduce the DLO Hi-C tool, a flexible and versatile pipeline for processing DLO Hi-C data from raw sequencing reads to normalized contact maps and for providing quality controls for different steps. It includes more efficient iterative mapping and linker filtering. We applied the DLO Hi-C tool to different DLO Hi-C datasets and demonstrated its ability in processing large data with multithreading. The DLO Hi-C tool is suitable for processing DLO Hi-C and in situ DLO Hi-C datasets. It is convenient and efficient for DLO Hi-C data processing.

Pituitary-Derived Circular RNAs Expression and Regulatory Network Prediction During the Onset of Puberty in Landrace × Yorkshire Crossbred Pigs.

Being the center of the hypothalamus-pituitary-ovary (HPO) axis, the pituitary plays a key role in the onset of puberty. Recent studies show that circular RNAs (circRNAs) can perform as miRNA sponges to regulate development in animals. However, the function of pituitary-derived circRNAs in first estrus remains unclear in pigs. In this study, we performed a genome-wide identification and characterization of circRNAs using pituitaries from Landrace × Yorkshire crossbred pigs at three stages: pre-, in-, and post-puberty, to describe such pituitary-derived circRNAs in pigs. A total of 5148 circRNAs were found in the gilts' pituitaries, averaging 18 682 bp in genomic distance, which consisted of approximately 91% exonic, 6% intergenic, and 3% intronic circRNAs. Furthermore, 158 novel circRNAs were identified for the first time and classified as putative pituitary-specific circRNAs. Their expression levels during the onset of puberty, significantly exceeded those of the other circRNAs, and the parental genes of these putative pituitary-specific circRNAs were enriched in “ssc04917: prolactin signaling pathway,” “ssc04080: neuroactive ligand-receptor interaction,” and “ssc04728: dopaminergic synapse” pathways, all of which were consistent with pituitary functioning. Additionally, 17 differentially regulated circRNAs were found and investigated for their potential interaction with miRNAs, along with genes, by constructing a circRNA-targeted miRNA-gene network. Taken together, these results provide new insight into the circRNA-mediated timing of puberty in gilts at the pituitary level.

Measuring significant changes in chromatin conformation with ACCOST.

Chromatin conformation assays such as Hi-C cannot directly measure differences in 3D architecture between cell types or cell states. For this purpose, two or more Hi-C experiments must be carried out, but direct comparison of the resulting Hi-C matrices is confounded by several features of Hi-C data. Most notably, the genomic distance effect, whereby contacts between pairs of genomic loci that are proximal along the chromosome exhibit many more Hi-C contacts that distal pairs of loci, dominates every Hi-C matrix. Furthermore, the form that this distance effect takes often varies between different Hi-C experiments, even between replicate experiments. Thus, a statistical confidence measure designed to identify differential Hi-C contacts must accurately account for the genomic distance effect or risk being misled by large-scale but artifactual differences. ACCOST (Altered Chromatin COnformation STatistics) accomplishes this goal by extending the statistical model employed by DEseq, re-purposing the ‘size factors,’ which were originally developed to account for differences in read depth between samples, to instead model the genomic distance effect. We show via analysis of simulated and real data that ACCOST provides unbiased statistical confidence estimates that compare favorably with competing methods such as diffHiC, FIND and HiCcompare. ACCOST is freely available with an Apache license at https://bitbucket.org/noblelab/accost.

A Well-Mixed E. coli Genome: Widespread Contacts Revealed by Tracking Mu Transposition

The three-dimensional structures of chromosomes are increasingly being recognized as playing a major role in cellular regulatory states. The efficiency and promiscuity of phage Mu transposition was exploited to directly measure invivo interactions between genomic loci in E.coli. Two global organizing principles have emerged: first, the chromosome is well-mixed and uncompartmentalized, with transpositions occurring freely between all measured loci; second, several gene families/regions show "clustering": strong three-dimensional co-localization regardless of linear genomic distance. The activities of the SMC/condensin protein MukB and nucleoid-compacting protein subunit HU-α are essential for the well-mixed state; HU-α is also needed for clustering of 6/7 ribosomal RNA-encoding loci. The data are explained by a model in which the chromosomal structure is driven by dynamic competition between DNA replication and chromosomal relaxation, providing a foundation for determining how region-specific properties contribute to both chromosomal structure and gene regulation.