Squalomix: shark and ray genome analysis consortium and its data sharing platform.

The taxon Elasmobranchii (sharks and rays) contains one of the long-established evolutionary lineages of vertebrates with a tantalizing collection of species occupying critical aquatic habitats. To overcome the current limitation in molecular resources, we launched the Squalomix Consortium in 2020 to promote a genome-wide array of molecular approaches, specifically targeting shark and ray species. Among the various bottlenecks in working with elasmobranchs are their elusiveness and low fecundity as well as the large and highly repetitive genomes. Their peculiar body fluid composition has also hindered the establishment of methods to perform routine cell culturing required for their karyotyping. In the Squalomix consortium, these obstacles are expected to be solved through a combination of in-house cytological techniques including karyotyping of cultured cells, chromatin preparation for Hi-C data acquisition, and high fidelity long-read sequencing. The resources and products obtained in this consortium, including genome and transcriptome sequences, a genome browser powered by JBrowse2 to visualize sequence alignments, and comprehensive matrices of gene expression profiles for selected species are accessible through https://github.com/Squalomix/info.

MarpolBase Expression: A Web-Based, Comprehensive Platform for Visualization and Analysis of Transcriptomes in the Liverwort Marchantia polymorpha.

The liverwort Marchantia polymorpha is equipped with a wide range of molecular and genetic tools and resources that have led to its wide use to explore the evo-devo aspects of land plants. Although its diverse transcriptome data are rapidly accumulating, there is no extensive yet user-friendly tool to exploit such a compilation of data and to summarize results with the latest annotations. Here, we have developed a web-based suite of tools, MarpolBase Expression (MBEX, https://marchantia.info/mbex/), where users can visualize gene expression profiles, identify differentially expressed genes, perform co-expression and functional enrichment analyses and summarize their comprehensive output in various portable formats. Using oil body biogenesis as an example, we demonstrated that the results generated by MBEX were consistent with the published experimental evidence and also revealed a novel transcriptional network in this process. MBEX should facilitate the exploration and discovery of the genetic and functional networks behind various biological processes in M. polymorpha and promote our understanding of the evolution of land plants.

Clostridium zeae sp. nov., isolated from corn silage.

A rod-shaped, Gram-stain-negative, strictly anaerobic, catalase-negative and endospore-forming bacterial strain CSC2T was isolated from corn silage preserved in Tochigi, Japan. The strain CSC2T grew at 15-40 °C, at pH 5.0-7.7 and with up to 0.5 % (w/v) NaCl. The main cellular fatty acids were C14 : 0, C16 : 0 and C16 : 0 dimethyl acetal. The cellular polar lipids detected were diphosphatidylglycerol, phosphatidylglycerol, phosphatidylethanolamine, phosphatidic acid, lysophosphatidylethanolamine, phosphatidylserine, lysophosphatidylcholine and two unidentified polar lipids. Phylogenetic analysis of the 16S rRNA gene showed that strain CSC2T was a member of the genus Clostridium and closely related to Clostridium polyendosporum DSM 57272T (95.6 % gene sequence similarity) and Clostridium fallax ATCC 19400T (95.3 %). The genomic DNA G+C content of strain CSC2T was 31.1 mol% (whole genome analysis). The average nucleotide identity based on blast and digital DNA-DNA hybridization values between strain CSC2T and the type strains of phylogenetically related species were below 71 and 24 %, respectively. On the basis of the genotypic, phenotypic and chemotaxonomic characteristics, it is proposed to designate strain CSC2T as representing Clostridium zeae sp. nov. The type strain is CSC2T (=MAFF212476T=JCM 33766T=DSM 111242T).

Identification of the sex-determining factor in the liverwort Marchantia polymorpha reveals unique evolution of sex chromosomes in a haploid system.

Sex determination is a central process for sexual reproduction and is often regulated by a sex determinant encoded on a sex chromosome. Rules that govern the evolution of sex chromosomes via specialization and degeneration following the evolution of a sex determinant have been well studied in diploid organisms. However, distinct predictions apply to sex chromosomes in organisms where sex is determined in the haploid phase of the life cycle: both sex chromosomes, female U and male V, are expected to maintain their gene functions, even though both are non-recombining. This is in contrast to the X-Y (or Z-W) asymmetry and Y (W) chromosome degeneration in XY (ZW) systems of diploids. Here, we provide evidence that sex chromosomes diverged early during the evolution of haploid liverworts and identify the sex determinant on the Marchantia polymorpha U chromosome. This gene, Feminizer, encodes a member of the plant-specific BASIC PENTACYSTEINE transcription factor family. It triggers female differentiation via regulation of the autosomal sex-determining locus of FEMALE GAMETOPHYTE MYB and SUPPRESSOR OF FEMINIZATION. Phylogenetic analyses of Feminizer and other sex chromosome genes indicate dimorphic sex chromosomes had already been established 430 mya in the ancestral liverwort. Feminizer also plays a role in reproductive induction that is shared with its gametolog on the V chromosome, suggesting an ancestral function, distinct from sex determination, was retained by the gametologs. This implies ancestral functions can be preserved after the acquisition of a sex determination mechanism during the evolution of a dominant haploid sex chromosome system.

Chromatin Organization in Early Land Plants Reveals an Ancestral Association between H3K27me3, Transposons, and Constitutive Heterochromatin.

Genome packaging by nucleosomes is a hallmark of eukaryotes. Histones and the pathways that deposit, remove, and read histone modifications are deeply conserved. Yet, we lack information regarding chromatin landscapes in extant representatives of ancestors of the main groups of eukaryotes, and our knowledge of the evolution of chromatin-related processes is limited. We used the bryophyte Marchantia polymorpha, which diverged from vascular plants circa 400 mya, to obtain a whole chromosome genome assembly and explore the chromatin landscape and three-dimensional genome organization in an early diverging land plant lineage. Based on genomic profiles of ten chromatin marks, we conclude that the relationship between active marks and gene expression is conserved across land plants. In contrast, we observed distinctive features of transposons and other repetitive sequences in Marchantia compared with flowering plants. Silenced transposons and repeats did not accumulate around centromeres. Although a large fraction of constitutive heterochromatin was marked by H3K9 methylation as in flowering plants, a significant proportion of transposons were marked by H3K27me3, which is otherwise dedicated to the transcriptional repression of protein-coding genes in flowering plants. Chromatin compartmentalization analyses of Hi-C data revealed that repressed B compartments were densely decorated with H3K27me3 but not H3K9 or DNA methylation as reported in flowering plants. We conclude that, in early plants, H3K27me3 played an essential role in heterochromatin function, suggesting an ancestral role of this mark in transposon silencing.

Structural basis for replication origin unwinding by an initiator primase of plasmid ColE2-P9: duplex DNA unwinding by a single protein.

Duplex DNA is generally unwound by protein oligomers prior to replication. The Rep protein of plasmid ColE2-P9 (34 kDa) is an essential initiator for plasmid DNA replication. This protein binds the replication origin (Ori) in a sequence-specific manner as a monomer and unwinds DNA. Here we present the crystal structure of the DNA-binding domain of Rep (E2Rep-DBD) in complex with Ori DNA. The structure unveils the basis for Ori-specific recognition by the E2Rep-DBD and also reveals that it unwinds DNA by the concerted actions of its three contiguous structural modules. The structure also shows that the functionally unknown PriCT domain, which forms a compact module, plays a central role in DNA unwinding. The conservation of the PriCT domain in the C termini of some archaeo-eukaryotic primases indicates that it probably plays a similar role in these proteins. Thus, this is the first report providing the structural basis for the functional importance of the conserved PriCT domain and also reveals a novel mechanism for DNA unwinding by a single protein.

Concerted interaction between origin recognition complex (ORC), nucleosomes and replication origin DNA ensures stable ORC-origin binding.

Chromosomal replication origins, where DNA replication is initiated, are determined in eukaryotic cells by specific binding of a six-subunit origin recognition complex (ORC). Many biochemical analyses have showed the detailed properties of the ORC-DNA interaction. However, because of the lack of in vitro analysis, the molecular architecture of the ORC-chromatin interaction is unclear. Recently, mainly from in vivo analyses, a role of chromatin in the ORC-origin interaction has been reported, including the existence of a specific pattern of nucleosome positioning around origins and of a specific interaction between chromatin-or core histones-and Orc1, a subunit of ORC. Therefore, to understand how ORC establishes its interaction with origin in vivo, it is essential to know the molecular mechanisms of the ORC-chromatin interaction. Here, we show that ORC purified from yeast binds more stably to origin-containing reconstituted chromatin than to naked DNA and forms a nucleosome-free region at origins. Molecular imaging using atomic force microscopy (AFM) shows that ORC associates with the adjacent nucleosomes and forms a larger complex. Moreover, stable binding of ORC to chromatin requires linker DNA. Thus, ORC establishes its interaction with origin by binding to both nucleosome-free origin DNA and neighboring nucleosomes.

Ctf4 coordinates the progression of helicase and DNA polymerase alpha.

Ctf4 is a protein conserved in eukaryotes and a constituent of the replisome progression complex. It also plays a role in the establishment of sister chromatid cohesion. In our current study, we demonstrate that the replication checkpoint is activated in the absence of Ctf4, and that the interaction between the MCM helicase-go ichi ni san (GINS) complex and DNA polymerase alpha (Pol alpha)-primase is destabilized specifically in a ctf4Delta mutant. An in vitro interaction between GINS and DNA Pol alpha was also found to be mediated by Ctf4. The same interaction was not affected in the absence of the replication checkpoint mediators Tof1 or Mrc1. In ctf4Delta cells, DNA pol alpha became significantly unstable and was barely detectable at the replication forks in HU. In contrast, the quantities of helicase and DNA pol epsilon bound to replication forks were almost unchanged but their localizations were widely and abnormally dispersed in the mutant cells compared with wild type. These results lead us to propose that Ctf4 is a key connector between DNA helicase and Pol alpha and is required for the coordinated progression of the replisome.

The ColE2-P9 Rep protein binds to the origin DNA as a monomer.

The Rep proteins of some plasmid replicons have two functions. Dimers bind to the operator sequences acting as auto-repressors, whereas monomers bind to the iterons to initiate replication of DNA. The ColE2 Rep proteins are present mostly in a dimeric form with some multimers larger than dimers in solution, while the form of Rep binding to Ori is not known. We used an EMSA-based method to determine the molecular weight of Rep in the Rep-Ori complex. The result suggested that Rep binds to Ori as a monomer. In addition, the result of EMSA using the Rep protein fused with the maltose binding protein and the His6-tag also supported this conclusion. We proposed that dimerization of Rep might probably be involved in keeping the copy number of the ColE2 plasmid at the normal low level by limiting the amount of active monomeric forms of Rep in the host cell.

The Rep protein binding elements of the plasmid ColE2-P9 replication origin.

The plasmid ColE2-P9 (ColE2) origin (32bp) is specifically recognized by the plasmid-specified Rep protein that initiates DNA replication. The ColE2 origin is divided into at least three functional subregions (I, II, and III), and three sites (a, b, and c) found in subregions I and II play important roles in Rep protein binding. We performed SELEX experiments of plasmid ColE2 to determine the optimal sequences for specific binding of the Rep protein. From these experiments, we obtained a common 16-bp sequence (5'-TGAGACCANATAAGCC-3'), which corresponds to about one half of the minimal ColE2 origin and contains sites a and b. Gel mobility shift assays using single-point mutant origins and the Rep protein further indicated that high affinity sequence-specific recognition by the Rep protein requires sites a, b, and c, but that mutations in site c were less disruptive to this recognition than those in sites a and b.

Anatomy of the replication origin of plasmid ColE2-P9.

The plasmid ColE2-P9 origin is a 32-bp region which is specifically recognized by the plasmid-specified Rep protein to initiate DNA replication. We analyzed the structural and functional organization of the ColE2 origin by using various derivatives carrying deletions and single-base-pair substitutions. The origin may be divided into three subregions: subregion I, which is important for stable binding of the Rep protein; subregion II, which is important for binding of the Rep protein and for initiation of DNA replication; and subregion III, which is important for DNA replication but apparently not for binding of the Rep protein. The Rep protein might recognize three specific DNA elements in subregions I and II. The relative transformation frequency of the autonomously replicating plasmids carrying deletions in subregion I is lower, and nevertheless the copy numbers of these plasmids in host bacteria are higher than those of the wild-type plasmid. Efficient and stable binding of the Rep protein to the origin might be important for the replication efficiency to be at the normal (low) level. Subregion II might be essential for interaction with the catalytic domain of the Rep protein for primer RNA synthesis. The 8-bp sequence across the border of subregions II and III, including the primer sequence, is conserved in the (putative) origins of many plasmids, the putative Rep proteins of which are related to the ColE2-P9 Rep protein. Subregion III might be required for a step that is necessary after Rep protein binding has taken place.

Whole genome association study of rheumatoid arthritis using 27 039 microsatellites.

A major goal of current human genome-wide studies is to identify the genetic basis of complex disorders. However, the availability of an unbiased, reliable, cost efficient and comprehensive methodology to analyze the entire genome for complex disease association is still largely lacking or problematic. Therefore, we have developed a practical and efficient strategy for whole genome association studies of complex diseases by charting the human genome at 100 kb intervals using a collection of 27,039 microsatellites and the DNA pooling method in three successive genomic screens of independent case-control populations. The final step in our methodology consists of fine mapping of the candidate susceptible DNA regions by single nucleotide polymorphisms (SNPs) analysis. This approach was validated upon application to rheumatoid arthritis, a destructive joint disease affecting up to 1% of the population. A total of 47 candidate regions were identified. The top seven loci, withstanding the most stringent statistical tests, were dissected down to individual genes and/or SNPs on four chromosomes, including the previously known 6p21.3-encoded Major Histocompatibility Complex gene, HLA-DRB1. Hence, microsatellite-based genome-wide association analysis complemented by end stage SNP typing provides a new tool for genetic dissection of multifactorial pathologies including common diseases.

Novel algorithm for automated genotyping of microsatellites.

Microsatellites or short tandem repeats (STRs) are abundant in the human genome with easily assayed polymorphisms, providing powerful genetic tools for mapping both Mendelian and complex traits. Microsatellite genotyping requires detection of the products of polymerase chain reaction (PCR) amplification by electrophoresis, and analysis of the peak data for discrimination of the true allele. A high-throughput genotyping approach requires computer-based automation at both the detection and analysis phases. In order to achieve this, complicated peak patterns from individual alleles must be interpreted in order to assign alleles. Previous methods consider limited types of noise peaks and cannot provide enough accuracy. By pattern recognition of various types of noise peaks, such as stutter peaks and additional peaks, we have achieved an overall average accuracy of 94% for allele calling in our actual data. Our algorithm is crucial for a high-throughput genotyping system for microsatellite markers by reducing manual editing and human errors.