The Entire Mitochondrial Genome of Macrophthalmus abbreviatus Reveals Insights into the Phylogeny and Gene Rearrangements of Brachyura

Brachyuran crabs comprise the most species-rich clades among extant Decapoda and are divided into several major superfamilies. However, the phylogeny of Brachyuran remains controversial, comprehensive analysis of the overall phylogeny is still lacking. Complete mitochondrial genome (mitogenome) can indicate phylogenetic relationships, as well as useful information for gene rearrangement mechanisms and molecular evolution. In this study, we firstly sequenced and annotated the complete mitogenome of Macrophthalmus abbreviatus (Brachyura; Macrophthalmidae). The mitogenome length of M. abbreviatus is 16,322 bp, containing the entire set of 37 genes and a control region typically observed in Brachyuran mitogenomes. The genome composition of M. abbreviatus was highly A+T biased 76.3% showing positive AT-skew (0.033) and negative GC-skew (??0.351). In M. abbreviatus mitogenome, most tRNA genes were folded into the clover-leaf secondary structure except trnH, trnS1 and trnC, which was similar to the other species in Macrophthalmidae. Phylogenetic analysis showed that all families form a monophyletic, and Varunidae and Macrophthalmidae clustered into a monophyletic clade as sister groups. Comparative analyses of rearrangement among Brachyura revealed that Varunidae (Grapsoidea) and Macrophthalmidae (Ocypodoidea) had the same gene order, which reinforced the result of phylogeny. The combined results of two aspects revealed that the polyphyly of Ocypodoidea and Grapsoidea were well supported. In general, the results obtained in this research will contribute to further studies on molecular based for the classification and gene rearrangements of Macrophthalmidae or even Brachyura.

     

Targeting of >1.5 Mb of Human DNA into the Mouse X Chromosome Reveals Presence of cis-Acting Regulators of Epigenetic Silencing

Regulatory sequences can influence the expression of flanking genes over long distances, and X chromosome inactivation is a classic example of cis-acting epigenetic gene regulation. Knock-ins directed to the Mus musculus Hprt locus offer a unique opportunity to analyze the spread of silencing into different human DNA sequences in the identical genomic environment. X chromosome inactivation of four knock-in constructs, including bacterial artificial chromosome (BAC) integrations of over 195 kb, was demonstrated by both the lack of expression from the inactive X chromosome in females with nonrandom X chromosome inactivation and promoter DNA methylation of the human transgene in females. We further utilized promoter DNA methylation to assess the inactivation status of 74 human reporter constructs comprising >1.5 Mb of DNA. Of the 47 genes examined, only the PHB gene showed female DNA hypomethylation approaching the level seen in males, and escape from X chromosome inactivation was verified by demonstration of expression from the inactive X chromosome. Integration of PHB resulted in lower DNA methylation of the flanking HPRT promoter in females, suggesting the action of a dominant cis-acting escape element. Female-specific DNA hypermethylation of CpG islands not associated with promoters implies a widespread imposition of DNA methylation during X chromosome inactivation; yet transgenes demonstrated differential capacities to accumulate DNA methylation when integrated into the identical location on the inactive X chromosome, suggesting additional cis-acting sequence effects. As only one of the human transgenes analyzed escaped X chromosome inactivation, we conclude that elements permitting ongoing expression from the inactive X are rare in the human genome.     

Whole Genome Sequencing Reveals Local Transmission Patterns of Mycobacterium bovis in Sympatric Cattle and Badger Populations

Whole genome sequencing (WGS) technology holds great promise as a tool for the forensic epidemiology of bacterial pathogens. It is likely to be particularly useful for studying the transmission dynamics of an observed epidemic involving a largely unsampled ‘reservoir’ host, as for bovine tuberculosis (bTB) in British and Irish cattle and badgers. BTB is caused by Mycobacterium bovis, a member of the M. tuberculosis complex that also includes the aetiological agent for human TB. In this study, we identified a spatio-temporally linked group of 26 cattle and 4 badgers infected with the same Variable Number Tandem Repeat (VNTR) type of M. bovis. Single-nucleotide polymorphisms (SNPs) between sequences identified differences that were consistent with bacterial lineages being persistent on or near farms for several years, despite multiple clear whole herd tests in the interim. Comparing WGS data to mathematical models showed good correlations between genetic divergence and spatial distance, but poor correspondence to the network of cattle movements or within-herd contacts. Badger isolates showed between zero and four SNP differences from the nearest cattle isolate, providing evidence for recent transmissions between the two hosts. This is the first direct genetic evidence of M. bovis persistence on farms over multiple outbreaks with a continued, ongoing interaction with local badgers. However, despite unprecedented resolution, directionality of transmission cannot be inferred at this stage. Despite the often notoriously long timescales between time of infection and time of sampling for TB, our results suggest that WGS data alone can provide insights into TB epidemiology even where detailed contact data are not available, and that more extensive sampling and analysis will allow for quantification of the extent and direction of transmission between cattle and badgers.     

A Global Genome Segmentation Method for Exploration of Epigenetic Patterns

Current genome-wide ChIP-seq experiments on different epigenetic marks aim at unraveling the interplay between their regulation mechanisms. Published evaluation tools, however, allow testing for predefined hypotheses only. Here, we present a novel method for annotation-independent exploration of epigenetic data and their inter-correlation with other genome-wide features. Our method is based on a combinatorial genome segmentation solely using information on combinations of epigenetic marks. It does not require prior knowledge about the data (e.g. gene positions), but allows integrating the data in a straightforward manner. Thereby, it combines compression, clustering and visualization of the data in a single tool. Our method provides intuitive maps of epigenetic patterns across multiple levels of organization, e.g. of the co-occurrence of different epigenetic marks in different cell types. Thus, it facilitates the formulation of new hypotheses on the principles of epigenetic regulation. We apply our method to histone modification data on trimethylation of histone H3 at lysine 4, 9 and 27 in multi-potent and lineage-primed mouse cells, analyzing their combinatorial modification pattern as well as differentiation-related changes of single modifications. We demonstrate that our method is capable of reproducing recent findings of gene centered approaches, e.g. correlations between CpG-density and the analyzed histone modifications. Moreover, combining the clustered epigenetic data with information on the expression status of associated genes we classify differences in epigenetic status of e.g. house-keeping genes versus differentiation-related genes. Visualizing the distribution of modification states on the chromosomes, we discover strong patterns for chromosome X. For example, exclusively H3K9me3 marked segments are enriched, while poised and active states are rare. Hence, our method also provides new insights into chromosome-specific epigenetic patterns, opening up new questions how “epigenetic computation” is distributed over the genome in space and time.     

Integrative Analysis Reveals Relationships of Genetic and Epigenetic Alterations in Osteosarcoma

Background

Osteosarcomas are the most common non-haematological primary malignant tumours of bone, and all conventional osteosarcomas are high-grade tumours showing complex genomic aberrations. We have integrated genome-wide genetic and epigenetic profiles from the EuroBoNeT panel of 19 human osteosarcoma cell lines based on microarray technologies.

Principal Findings

The cell lines showed complex patterns of DNA copy number changes, where genomic copy number gains were significantly associated with gene-rich regions and losses with gene-poor regions. By integrating the datasets, 350 genes were identified as having two types of aberrations (gain/over-expression, hypo-methylation/over-expression, loss/under-expression or hyper-methylation/under-expression) using a recurrence threshold of 6/19 (>30%) cell lines. The genes showed in general alterations in either DNA copy number or DNA methylation, both within individual samples and across the sample panel. These 350 genes are involved in embryonic skeletal system development and morphogenesis, as well as remodelling of extracellular matrix. The aberrations of three selected genes, CXCL5, DLX5 and RUNX2, were validated in five cell lines and five tumour samples using PCR techniques. Several genes were hyper-methylated and under-expressed compared to normal osteoblasts, and expression could be reactivated by demethylation using 5-Aza-2′-deoxycytidine treatment for four genes tested; AKAP12, CXCL5, EFEMP1 and IL11RA. Globally, there was as expected a significant positive association between gain and over-expression, loss and under-expression as well as hyper-methylation and under-expression, but gain was also associated with hyper-methylation and under-expression, suggesting that hyper-methylation may oppose the effects of increased copy number for detrimental genes.

Conclusions

Integrative analysis of genome-wide genetic and epigenetic alterations identified dependencies and relationships between DNA copy number, DNA methylation and mRNA expression in osteosarcomas, contributing to better understanding of osteosarcoma biology.     

Functional Identification and Analysis of cis-Acting Sequences Which Mediate Genome Cleavage and Packaging in Human Herpesvirus 6

Sequences present at the genomic termini of herpesviruses become linked during lytic-phase replication and provide the substrate for cleavage and packaging of unit length viral genomes. We have previously shown that homologs of the consensus herpesvirus cleavage-packaging signals, pac1 and pac2, are located at the left and right genomic termini of human herpesvirus 6 (HHV-6), respectively. Immediately adjacent to these elements are two distinct arrays of human telomeric repeat sequences (TRS). We now show that the unique sequence element formed at the junction of HHV-6B genome concatemers (pac2-pac1) is necessary and sufficient for virally mediated cleavage of plasmid DNAs containing the HHV-6B lytic-phase origin of DNA replication (oriLyt). The concatemeric junction sequence also allowed for the packaging of these plasmid molecules into intracellular nucleocapsids as well as mature, infectious viral particles. In addition, this element significantly enhanced the replication efficiency of oriLyt-containing plasmids in virally infected cells. Experiments revealed that the concatemeric junction sequence possesses an unusual, S1 nuclease-sensitive conformation (anisomorphic DNA), which might play a role in this apparent enhancement of DNA replication—although additional studies will be required to test this hypothesis. Finally, we also analyzed whether the presence of flanking viral TRS had any effect on the functional activity of the minimal concatemeric junction (pac2-pac1). These experiments revealed that the TRS motifs, either alone or in combination, had no effect on the efficiency of virally mediated DNA replication or DNA cleavage. Taken together, these data show that the cleavage and packaging of HHV-6 DNA are mediated by cis-acting consensus sequences similar to those found in other herpesviruses, and that these sequences also influence the efficiency of HHV-6 DNA replication. Since the adjacent TRS do not influence either viral cleavage and packaging or viral DNA replication, their function remains uncertain.     

Reactivation of an Inactive Centromere Reveals Epigenetic and Structural Components for Centromere Specification in Maize

Stable maize (Zea mays) chromosomes were recovered from an unstable dicentric containing large and small versions of the B chromosome centromere. In the stable chromosome, the smaller centromere had become inactivated. This inactive centromere can be inherited from one generation to the next attached to the active version and loses all known cytological and molecular properties of active centromeres. When separated from the active centromere by intrachromosomal recombination, the inactive centromere can be reactivated. The reactivated centromere regains the molecular attributes of activity in anaphase I of meiosis. When two copies of the dicentric chromosome with one active and one inactive centromere are present, homologous chromosome pairing reduces the frequency of intrachromosomal recombination and thus decreases, but does not eliminate, the reactivation of inactive centromeres. These findings indicate an epigenetic component to centromere specification in that centromere inactivation can be directed by joining two centromeres in opposition. These findings also indicate a structural aspect to centromere specification revealed by the gain of activity at the site of the previously inactive sequences.     

Regulation of per and cry Genes Reveals a Central Role for the D-Box Enhancer in Light-Dependent Gene Expression

Light serves as a key environmental signal for synchronizing the circadian clock with the day night cycle. The zebrafish represents an attractive model for exploring how light influences the vertebrate clock mechanism. Direct illumination of most fish tissues and cell lines induces expression of a broad range of genes including DNA repair, stress response and key clock genes. We have previously identified D- and E-box elements within the promoter of the zebrafish per2 gene that together direct light-induced gene expression. However, is the combined regulation by E- and D-boxes a general feature for all light-induced gene expression? We have tackled this question by examining the regulation of additional light-inducible genes. Our results demonstrate that with the exception of per2, all other genes tested are not induced by light upon blocking of de novo protein synthesis. We reveal that a single D-box serves as the principal light responsive element within the cry1a promoter. Furthermore, upon inhibition of protein synthesis D-box mediated gene expression is abolished while the E-box confers light driven activation as observed in the per2 gene. Given the existence of different photoreceptors in fish cells, our results implicate the D-box enhancer as a general convergence point for light driven signaling.     

The Mitochondrial Genome of Soybean Reveals Complex Genome Structures and Gene Evolution at Intercellular and Phylogenetic Levels

Determining mitochondrial genomes is important for elucidating vital activities of seed plants. Mitochondrial genomes are specific to each plant species because of their variable size, complex structures and patterns of gene losses and gains during evolution. This complexity has made research on the soybean mitochondrial genome difficult compared with its nuclear and chloroplast genomes. The present study helps to solve a 30-year mystery regarding the most complex mitochondrial genome structure, showing that pairwise rearrangements among the many large repeats may produce an enriched molecular pool of 760 circles in seed plants. The soybean mitochondrial genome harbors 58 genes of known function in addition to 52 predicted open reading frames of unknown function. The genome contains sequences of multiple identifiable origins, including 6.8 kb and 7.1 kb DNA fragments that have been transferred from the nuclear and chloroplast genomes, respectively, and some horizontal DNA transfers. The soybean mitochondrial genome has lost 16 genes, including nine protein-coding genes and seven tRNA genes; however, it has acquired five chloroplast-derived genes during evolution. Four tRNA genes, common among the three genomes, are derived from the chloroplast. Sizeable DNA transfers to the nucleus, with pericentromeric regions as hotspots, are observed, including DNA transfers of 125.0 kb and 151.6 kb identified unambiguously from the soybean mitochondrial and chloroplast genomes, respectively. The soybean nuclear genome has acquired five genes from its mitochondrial genome. These results provide biological insights into the mitochondrial genome of seed plants, and are especially helpful for deciphering vital activities in soybean.     

Homologous Series by Chromosome Number and the Genome Rearrangements in the Phylogeny of Salmonoidei

Published data on chromosome numbers of Salmonoidei are summarized. The existence of homologous variation of chromosome number in different phyletic lines of this suborder is substantiated. It is suggested that the origin of homologous series is related to major genome rearrangements (simultaneous fusion of several chromosomes).     

Transcriptome Analysis of Honeybee (Apis Mellifera) Haploid and Diploid Embryos Reveals Early Zygotic Transcription during Cleavage

In honeybees, the haplodiploid sex determination system promotes a unique embryogenesis process wherein females develop from fertilized eggs and males develop from unfertilized eggs. However, the developmental strategies of honeybees during early embryogenesis are virtually unknown. Similar to most animals, the honeybee oocytes are supplied with proteins and regulatory elements that support early embryogenesis. As the embryo develops, the zygotic genome is activated and zygotic products gradually replace the preloaded maternal material. The analysis of small RNA and mRNA libraries of mature oocytes and embryos originated from fertilized and unfertilized eggs has allowed us to explore the gene expression dynamics in the first steps of development and during the maternal-to-zygotic transition (MZT). We localized a short sequence motif identified as TAGteam motif and hypothesized to play a similar role in honeybees as in fruit flies, which includes the timing of early zygotic expression (MZT), a function sustained by the presence of the zelda ortholog, which is the main regulator of genome activation. Predicted microRNA (miRNA)-target interactions indicated that there were specific regulators of haploid and diploid embryonic development and an overlap of maternal and zygotic gene expression during the early steps of embryogenesis. Although a number of functions are highly conserved during the early steps of honeybee embryogenesis, the results showed that zygotic genome activation occurs earlier in honeybees than in Drosophila based on the presence of three primary miRNAs (pri-miRNAs) (ame-mir-375, ame-mir-34 and ame-mir-263b) during the cleavage stage in haploid and diploid embryonic development.     

Machine Learning Methods for Exploring Sequence Determinants of 3D Genome Organization

In higher eukaryotic cells, chromosomes are folded inside the nucleus. Recent advances in whole-genome mapping technologies have revealed the multiscale features of 3D genome organization that are intertwined with fundamental genome functions. However, DNA sequence determinants that modulate the formation of 3D genome organization remain poorly characterized. In the past few years, predicting 3D genome organization based on DNA sequence features has become an active area of research. Here, we review the recent progress in computational approaches to unraveling important sequence elements for 3D genome organization. In particular, we discuss the rapid development of machine learning-based methods that facilitate the connections between DNA sequence features and 3D genome architectures at different scales. While much progress has been made in developing predictive models for revealing important sequence features for 3D genome organization, new research is urgently needed to incorporate multi-omic data and enhance model interpretability, further advancing our understanding of gene regulation mechanisms through the lens of 3D genome organization.