The Chara Genome: Secondary Complexity and Implications for Plant Terrestrialization.

Land plants evolved from charophytic algae, among which Charophyceae possess the most complex body plans. We present the genome of Chara braunii; comparison of the genome to those of land plants identified evolutionary novelties for plant terrestrialization and land plant heritage genes. C. braunii employs unique xylan synthases for cell wall biosynthesis, a phragmoplast (cell separation) mechanism similar to that of land plants, and many phytohormones. C. braunii plastids are controlled via land-plant-like retrograde signaling, and transcriptional regulation is more elaborate than in other algae. The morphological complexity of this organism may result from expanded gene families, with three cases of particular note: genes effecting tolerance to reactive oxygen species (ROS), LysM receptor-like kinases, and transcription factors (TFs). Transcriptomic analysis of sexual reproductive structures reveals intricate control by TFs, activity of the ROS gene network, and the ancestral use of plant-like storage and stress protection proteins in the zygote.

Impact of translational regulation on diel expression revealed by time-series ribosome profiling in Arabidopsis.

Plants have developed the ability to adjust to the day/night cycle through the expression of diel genes, which allow them to effectively respond to environmental changes and optimise their growth and development. Diel oscillations also have substantial implications in many physiological processes, including photosynthesis, floral development, and environmental stress responses. The expression of diel genes is regulated by a combination of the circadian clock and responses to environmental cues, such as light and temperature. A great deal of information is available on the transcriptional regulation of diel gene expression. However, the extent to which translational regulation is involved in controlling diel changes in expression is not yet clear. To investigate the impact of translational regulation on diel expression, we conducted Ribo-seq and RNA-seq analyses on a time-series sample of Arabidopsis shoots cultivated under a 12 h light/dark cycle. Our results showed that translational regulation is involved in about 71% of the genes exhibiting diel changes in mRNA abundance or translational activity, including clock genes, many of which are subject to both translational and transcriptional control. They also revealed that the diel expression of glycosylation and ion-transporter-related genes is mainly established through translational regulation. The expression of several diel genes likely subject to translational regulation through upstream open-reading frames was also determined.

Genome sequence and cell biological toolbox of the highly regenerative, coenocytic green feather alga Bryopsis.

Green feather algae (Bryopsidales) undergo a unique life cycle in which a single cell repeatedly executes nuclear division without cytokinesis, resulting in the development of a thallus (>100 mm) with characteristic morphology called coenocyte. Bryopsis is a representative coenocytic alga that has exceptionally high regeneration ability: extruded cytoplasm aggregates rapidly in seawater, leading to the formation of protoplasts. However, the genetic basis of the unique cell biology of Bryopsis remains poorly understood. Here, we present a high-quality assembly and annotation of the nuclear genome of Bryopsis sp. (90.7 Mbp, 27 contigs, N50 = 6.7 Mbp, 14 034 protein-coding genes). Comparative genomic analyses indicate that the genes encoding BPL-1/Bryohealin, the aggregation-promoting lectin, are heavily duplicated in Bryopsis, whereas homologous genes are absent in other ulvophyceans, suggesting the basis of regeneration capability of Bryopsis. Bryopsis sp. possesses >30 kinesins but only a single myosin, which differs from other green algae that have multiple types of myosin genes. Consistent with this biased motor toolkit, we observed that the bidirectional motility of chloroplasts in the cytoplasm was dependent on microtubules but not actin in Bryopsis sp. Most genes required for cytokinesis in plants are present in Bryopsis, including those in the SNARE or kinesin superfamily. Nevertheless, a kinesin crucial for cytokinesis initiation in plants (NACK/Kinesin-7II) is hardly expressed in the coenocytic part of the thallus, possibly underlying the lack of cytokinesis in this portion. The present genome sequence lays the foundation for experimental biology in coenocytic macroalgae.

Odd-Paired is Involved in Morphological Divergence of Snail-Feeding Beetles.

Body shape and size diversity and their evolutionary rates correlate with species richness at the macroevolutionary scale. However, the molecular genetic mechanisms underlying the morphological diversification across related species are poorly understood. In beetles, which account for one-fourth of the known species, adaptation to different trophic niches through morphological diversification appears to have contributed to species radiation. Here, we explored the key genes for the morphological divergence of the slender to stout body shape related to divergent feeding methods on large to small snails within the genus Carabus. We show that the zinc-finger transcription factor encoded by odd-paired (opa) controls morphological variation in the snail-feeding ground beetle Carabus blaptoides. Specifically, opa was identified as the gene underlying the slender to stout morphological difference between subspecies through genetic mapping and functional analysis via gene knockdown. Further analyses revealed that changes in opa cis-regulatory sequences likely contributed to the differences in body shape and size between C. blaptoides subspecies. Among opa cis-regulatory sequences, single nucleotide polymorphisms on the transcription factor binding sites may be associated with the morphological differences between C. blaptoides subspecies. opa was highly conserved in a wide range of taxa, especially in beetles. Therefore, opa may play an important role in adaptive morphological divergence in beetles.

Synchronized Expansion and Contraction of Olfactory, Vomeronasal, and Taste Receptor Gene Families in Hystricomorph Rodents.

Chemical senses, including olfaction, pheromones, and taste, are crucial for the survival of most animals. There has long been a debate about whether different types of senses might influence each other. For instance, primates with a strong sense of vision are thought to have weakened olfactory abilities, although the oversimplified trade-off theory is now being questioned. It is uncertain whether such interactions between different chemical senses occur during evolution. To address this question, we examined four receptor gene families related to olfaction, pheromones, and taste: olfactory receptor (OR), vomeronasal receptor type 1 and type 2 (V1R and V2R), and bitter taste receptor (T2R) genes in Hystricomorpha, which is morphologically and ecologically the most diverse group of rodents. We also sequenced and assembled the genome of the grasscutter, Thryonomys swinderianus. By examining 16 available genome assemblies alongside the grasscutter genome, we identified orthologous gene groups among hystricomorph rodents for these gene families to separate the gene gain and loss events in each phylogenetic branch of the Hystricomorpha evolutionary tree. Our analysis revealed that the expansion or contraction of the four gene families occurred synchronously, indicating that when one chemical sense develops or deteriorates, the others follow suit. The results also showed that V1R/V2R genes underwent the fastest evolution, followed by OR genes, and T2R genes were the most evolutionarily stable. This variation likely reflects the difference in ligands of V1R/V2Rs, ORs, and T2Rs: species-specific pheromones, environment-based scents, and toxic substances common to many animals, respectively.

Early embryonic mutations reveal dynamics of somatic and germ cell lineages in mice.

De novo mutations accumulate with zygotic cell divisions. However, the occurrence of these mutations and the way they are inherited by somatic cells and germ cells remain unclear. Here, we present a novel method to reconstruct cell lineages. We identified mosaic mutations in mice using deep whole-genome sequencing and reconstructed embryonic cell lineages based on the variant allele frequencies of the mutations. The reconstructed trees were confirmed using nuclear transfer experiments and the genotyping of approximately 50 offspring of each tree. The most detailed tree had 32 terminal nodes and showed cell divisions from the fertilized egg to germ cell- and somatic cell-specific lineages, indicating at least five independent cell lineages that would be selected as founders of the primordial germ cells. The contributions of each lineage to germ cells and offspring varied widely. At the emergence of the germ cell-specific lineages, 10-15 embryonic mutations had accumulated, suggesting that the pregastrulation mutation rate is 1.0 mutation per mitosis. Subsequent mutation rates were 0.7 for germ cells and 13.2 for tail fibroblasts. Our results show a new framework to assess embryonic lineages; further, we suggest an evolutionary strategy for preserving heterogeneity owing to postzygotic mutations in offspring.

Chimera RNA transcribed from integrated HPV18 genome with adjacent host genomic region promotes oncogenic gene expression through condensate formation.

Most cervical cancers are caused by human papillomavirus (HPV) infection. In HeLa cells, the HPV18 viral genome is integrated at chromosome 8q24.21 and activates transcription of the proto-oncogene c-Myc. However, the mechanism of how the integrated HPV genome and its transcribed RNAs exhibit transcription activation function has not been fully elucidated. In this study, we found that HPV18 transcripts contain an enhancer RNA-like function to activate proximal genes including CCAT1-5L and c-Myc. We showed that the human genome-integrated HPV18 genes are activated by transcription coregulators including BRD4 and Mediator. The transcribed HPV18 RNAs form a liquid-like condensate at chromosome 8q24.21 locus, which in turn accumulates RNA polymerase II. Moreover, we focused on a relatively uncharacterized transcript from the upstream region of CCAT1, named URC. The URC RNA is transcribed as a chimera RNA with HPV18 and is composed of the 3'-untranslated region of the HPV18 transcript. We experimentally showed that the URC contributes to stabilization of HPV18 RNAs by supplying a polyadenylation site for the HPV18 transcript. Our findings suggest that integrated HPV18 at 8q24.21 locus produces HPV18-URC chimera RNA and promotes tumorigenesis through RNA-based condensate formation.

Arms race between anti-silencing and RdDM in noncoding regions of transposable elements.

Transposable elements (TEs) are among the most dynamic parts of genomes. Since TEs are potentially deleterious, eukaryotes silence them through epigenetic mechanisms such as repressive histone modifications and DNA methylation. We previously reported that Arabidopsis TEs, called VANDALs, counteract epigenetic silencing through a group of sequence-specific anti-silencing proteins, VANCs. VANC proteins bind to noncoding regions of specific VANDAL copies and induce loss of silent chromatin marks. The VANC-target regions form tandem repeats, which diverge rapidly. Sequence-specific anti-silencing allows these TEs to proliferate with minimum host damage. Here, we show that RNA-directed DNA methylation (RdDM) efficiently targets noncoding regions of VANDAL TEs to silence them de novo. Thus, escape from RdDM could be a primary event leading to the rapid evolution and diversification of sequence-specific anti-silencing systems. We propose that this selfish behavior of TEs paradoxically could make them diverse and less harmful to the host.

The RNA helicase DDX6 controls early mouse embryogenesis by repressing aberrant inhibition of BMP signaling through miRNA-mediated gene silencing.

The evolutionarily conserved RNA helicase DDX6 is a central player in post-transcriptional regulation, but its role during embryogenesis remains elusive. We here show that DDX6 enables proper cell lineage specification from pluripotent cells by analyzing Ddx6 knockout (KO) mouse embryos and employing an in vitro epiblast-like cell (EpiLC) induction system. Our study unveils that DDX6 is an important BMP signaling regulator. Deletion of Ddx6 causes the aberrant upregulation of the negative regulators of BMP signaling, which is accompanied by enhanced expression of Nodal and related genes. Ddx6 KO pluripotent cells acquire higher pluripotency with a strong inclination toward neural lineage commitment. During gastrulation, abnormally expanded Nodal and Eomes expression in the primitive streak likely promotes endoderm cell fate specification while inhibiting mesoderm differentiation. We also genetically dissected major DDX6 pathways by generating Dgcr8, Dcp2, and Eif4enif1 KO models in addition to Ddx6 KO. We found that the miRNA pathway mutant Dgcr8 KO phenocopies Ddx6 KO, indicating that DDX6 mostly works along with the miRNA pathway during early development, whereas its P-body-related functions are dispensable. Therefore, we conclude that DDX6 prevents aberrant upregulation of BMP signaling inhibitors by participating in miRNA-mediated gene silencing processes. Overall, this study delineates how DDX6 affects the development of the three primary germ layers during early mouse embryogenesis and the underlying mechanism of DDX6 function.

Repeated translocation of a supergene underlying rapid sex chromosome turnover in Takifugu pufferfish.

Recent studies have revealed a surprising diversity of sex chromosomes in vertebrates. However, the detailed mechanism of their turnover is still elusive. To understand this process, it is necessary to compare closely related species in terms of sex-determining genes and the chromosomes harboring them. Here, we explored the genus Takifugu, in which one strong candidate sex-determining gene, Amhr2, has been identified. To trace the processes involved in transitions in the sex-determination system in this genus, we studied 12 species and found that while the Amhr2 locus likely determines sex in the majority of Takifugu species, three species have acquired sex-determining loci at different chromosomal locations. Nevertheless, the generation of genome assemblies for the three species revealed that they share a portion of the male-specific supergene that contains a candidate sex-determining gene, GsdfY, along with genes that potentially play a role in male fitness. The shared supergene spans ∼100 kb and is flanked by two duplicated regions characterized by CACTA transposable elements. These results suggest that the shared supergene has taken over the role of sex-determining locus from Amhr2 in lineages leading to the three species, and repeated translocations of the supergene underlie the turnover of sex chromosomes in these lineages. These findings highlight the underestimated role of a mobile supergene in the turnover of sex chromosomes in vertebrates.

Turnover of mammal sex chromosomes in the Sry-deficient Amami spiny rat is due to male-specific upregulation of Sox9.

Mammalian sex chromosomes are highly conserved, and sex is determined by SRY on the Y chromosome. Two exceptional rodent groups in which some species lack a Y chromosome and Sry offer insights into how novel sex genes can arise and replace Sry, leading to sex chromosome turnover. However, intensive study over three decades has failed to reveal the identity of novel sex genes in either of these lineages. We here report our discovery of a male-specific duplication of an enhancer of Sox9 in the Amami spiny rat Tokudaia osimensis, in which males and females have only a single X chromosome (XO/XO) and the Y chromosome and Sry are completely lost. We performed a comprehensive survey to detect sex-specific genomic regions in the spiny rat. Sex-related genomic differences were limited to a male-specific duplication of a 17-kb unit located 430 kb upstream of Sox9 on an autosome. Hi-C analysis using male spiny rat cells showed the duplicated region has potential chromatin interaction with Sox9. The duplicated unit harbored a 1,262-bp element homologous to mouse enhancer 14 (Enh14), a candidate Sox9 enhancer that is functionally redundant in mice. Transgenic reporter mice showed that the spiny rat Enh14 can function as an embryonic testis enhancer in mice. Embryonic gonads of XX mice in which Enh14 was replaced by the duplicated spiny rat Enh14 showed increased Sox9 expression and decreased Foxl2 expression. We propose that male-specific duplication of this Sox9 enhancer substituted for Sry function, defining a novel Y chromosome in the spiny rat.

Genomic and transcriptomic analyses of the subterranean termite Reticulitermes speratus: Gene duplication facilitates social evolution.

Termites are model social organisms characterized by a polyphenic caste system. Subterranean termites (Rhinotermitidae) are ecologically and economically important species, including acting as destructive pests. Rhinotermitidae occupies an important evolutionary position within the clade representing a transitional taxon between the higher (Termitidae) and lower (other families) termites. Here, we report the genome, transcriptome, and methylome of the Japanese subterranean termite Reticulitermes speratus Our analyses highlight the significance of gene duplication in social evolution in this termite. Gene duplication associated with caste-biased gene expression was prevalent in the R. speratus genome. The duplicated genes comprised diverse categories related to social functions, including lipocalins (chemical communication), cellulases (wood digestion and social interaction), lysozymes (social immunity), geranylgeranyl diphosphate synthase (social defense), and a novel class of termite lineage-specific genes with unknown functions. Paralogous genes were often observed in tandem in the genome, but their expression patterns were highly variable, exhibiting caste biases. Some of the assayed duplicated genes were expressed in caste-specific organs, such as the accessory glands of the queen ovary and the frontal glands of soldier heads. We propose that gene duplication facilitates social evolution through regulatory diversification, leading to caste-biased expression and subfunctionalization and/or neofunctionalization conferring caste-specialized functions.

WUSCHEL-related homeobox family genes in rice control lateral root primordium size.

The development of a plastic root system is essential for stable crop production under variable environments. Rice plants have two types of lateral roots (LRs): S-type (short and thin) and L-type (long, thick, and capable of further branching). LR types are determined at the primordium stage, with a larger primordium size in L-types than S-types. Despite the importance of LR types for rice adaptability to variable water conditions, molecular mechanisms underlying the primordium size control of LRs are unknown. Here, we show that two WUSCHEL-related homeobox (WOX) genes have opposing roles in controlling LR primordium (LRP) size in rice. Root tip excision on seminal roots induced L-type LR formation with wider primordia formed from an early developmental stage. QHB/OsWOX5 was isolated as a causative gene of a mutant that is defective in S-type LR formation but produces more L-type LRs than wild-type (WT) plants following root tip excision. A transcriptome analysis revealed that OsWOX10 is highly up-regulated in L-type LRPs. OsWOX10 overexpression in LRPs increased the LR diameter in an expression-dependent manner. Conversely, the mutation in OsWOX10 decreased the L-type LR diameter under mild drought conditions. The qhb mutants had higher OsWOX10 expression than WT after root tip excision. A yeast one-hybrid assay revealed that the transcriptional repressive activity of QHB was lost in qhb mutants. An electrophoresis mobility shift assay revealed that OsWOX10 is a potential target of QHB. These data suggest that QHB represses LR diameter increase, repressing OsWOX10 Our findings could help improve root system plasticity under variable environments.

A leaf-emanated signal orchestrates grain size and number in response to maternal resources.

In plants, variations in seed size and number are outcomes of different reproductive strategies. Both traits are often environmentally influenced, suggesting that a mechanism exists to coordinate these phenotypes in response to available maternal resources. Yet, how maternal resources are sensed and influence seed size and number is largely unknown. Here, we report a mechanism that senses maternal resources and coordinates grain size and number in the wild rice Oryza rufipogon, a wild progenitor of Asian cultivated rice. We showed that FT-like 9 (FTL9) regulates both grain size and number and that maternal photosynthetic assimilates induce FTL9 expression in leaves to act as a long-range signal that increases grain number and reduces size. Our findings highlight a strategy that benefits wild plants to survive in a fluctuating environment. In this strategy, when maternal resources are sufficient, wild plants increase their offspring number while preventing an increase in offspring size by the action of FTL9, which helps expand their habitats. In addition, we found that a loss-of-function allele (ftl9) is prevalent among wild and cultivated populations, offering a new scenario in the history of rice domestication.

MiRiQ Database: A Platform for In Silico Rice Mutant Screening.

Genetic studies using mutant resources have significantly contributed to elucidating plant gene function. Massive mutant libraries sequenced by next-generation sequencing technology facilitate mutant identification and functional analysis of genes of interest. Here, we report the creation and release of an open-access database (https://miriq.agr.kyushu-u.ac.jp/index.php), called Mutation-induced Rice in Kyushu University (MiRiQ), designed for in silico mutant screening based on a whole-genome-sequenced mutant library. This database allows any user to easily find mutants of interest without laborious efforts such as large-scale screening by PCR. The initial version of the MiRiQ database (version 1.0) harbors a total of 1.6 million single-nucleotide variants (SNVs) and InDels of 721 M1 plants that were mutagenized by N-methyl-N-nitrosourea treatment of the rice cultivar Nipponbare (Oryza sativa ssp. japonica). The SNVs were distributed among 87% of all 35,630 annotated protein-coding genes of the Nipponbare genome and were predicted to induce missense and nonsense mutations. The MiRiQ database provides built-in tools, such as a search tool by keywords and JBrowse for mutation searches. Users can request mutant seeds in the M2 or M3 generations from a request form linked to this database. We believe that the availability of a wide range of gene mutations in this database will benefit the plant science community and breeders worldwide by accelerating functional genomic research and crop improvement.

Transcriptome Analysis of Rice Root Tips Reveals Auxin, Gibberellin and Ethylene Signaling Underlying Nutritropism.

Nutritropism is a positive tropism toward nutrients in plant roots. An NH4+ gradient is a nutritropic stimulus in rice (Oryza sativa L.). When rice roots are exposed to an NH4+ gradient generated around nutrient sources, root tips bend toward and coil around the sources. The molecular mechanisms are largely unknown. Here, we analyzed the transcriptomes of the inside and outside of bending root tips exhibiting nutritropism to reveal nutritropic signal transduction. Tissues facing the nutrient sources (inside) and away (outside) were separately collected by laser microdissection. Principal component analysis revealed distinct transcriptome patterns between the two tissues. Annotations of 153 differentially expressed genes implied that auxin, gibberellin and ethylene signaling were activated differentially between the sides of the root tips under nutritropism. Exogenous application of transport and/or biosynthesis inhibitors of these phytohormones largely inhibited the nutritropism. Thus, signaling and de novo biosynthesis of the three phytohormones are necessary for nutritropism. Expression patterns of IAA genes implied that auxins accumulated more in the inside tissues, meaning that ammonium stimulus is transduced to auxin signaling in nutritropism similar to gravity stimulus in gravitropism. SAUR and expansin genes, which are known to control cell wall modification and to promote cell elongation in shoot gravitropism, were highly expressed in the inside tissues rather than the outside tissues, and our transcriptome data are unexplainable for differential elongation in root nutritropism.

Shared evolutionary trajectories of three independent neo-sex chromosomes in Drosophila.

Dosage compensation (DC) on the X Chromosome counteracts the deleterious effects of gene loss on the Y Chromosome. However, DC is not efficient if the X Chromosome also degenerates. This indeed occurs in Drosophila miranda, in which both the neo-Y and the neo-X are under accelerated pseudogenization. To examine the generality of this pattern, we investigated the evolution of two additional neo-sex chromosomes that emerged independently in D. albomicans and D. americana and reanalyzed neo-sex chromosome evolution in D. miranda Comparative genomic and transcriptomic analyses revealed that the pseudogenization rate on the neo-X is also accelerated in D. albomicans and D. americana although to a lesser extent than in D. miranda In males, neo-X-linked genes whose neo-Y-linked homologs are pseudogenized tended to be up-regulated more than those whose neo-Y-linked homologs remain functional. Moreover, genes under strong functional constraint and genes highly expressed in the testis tended to remain functional on the neo-X and neo-Y, respectively. Focusing on the D. miranda and D. albomicans neo-sex chromosomes that emerged independently from the same autosome, we further found that the same genes tend to become pseudogenized in parallel on the neo-Y. These genes include Idgf6 and JhI-26, which may be unnecessary or even harmful in males. Our results indicate that neo-sex chromosomes in Drosophila share a common evolutionary trajectory after their emergence, which may prevent sex chromosomes from being an evolutionary dead end.

Loss of One X and the Y Chromosome Changes the Configuration of the X Inactivation Center in the Genus Tokudaia.

INTRODUCTION: X chromosome inactivation (XCI) is an essential mechanism for dosage compensation between females and males in mammals. In females, XCI is controlled by a complex, conserved locus termed the X inactivation center (Xic), in which the lncRNA Xist is the key regulator. However, little is known about the Xic in species with unusual sex chromosomes. The genus Tokudaia includes three rodent species endemic to Japan. Tokudaia osimensis and Tokudaia tokunoshimensis lost the Y chromosome (XO/XO), while Tokudaia muenninki (TMU) acquired a neo-X region by fusion of the X chromosome and an autosome (XX/XY). We compared the gene location and structure in the Xic among Tokudaia species. METHODS: Gene structure of nine genes in Xic was predicted, and the gene location and genome sequences of Xic were compared between mouse and Tokudaia species. The expression level of the gene was confirmed by transcripts per million calculation using RNA-seq data. RESULTS: Compared to mouse, the Xic gene order and location were conserved in Tokudaia species. However, remarkable structure changes were observed in lncRNA genes, Xist and Tsix, in the XO/XO species. In Xist, important functional repeats, B-, C-, D-, and E-repeats, were partially or completely lost due to deletions in these species. RNA-seq data showed that female-specific expression patterns of Xist and Tsix were confirmed in TMU, however, not in the XO/XO species. Additionally, three deletions and one inversion were confirmed in the intergenic region between Jpx and Ftx in the XO/XO species. CONCLUSION: Our findings indicate that even if the Xist and Tsix lncRNAs are expressed, they are incapable of producing a successful and lasting XCI in the XO/XO species. We hypothesized that the significant structure change in the intergenic region of Jpx-Ftx resulted in the inability to perform the XCI, and, as a result, a lack of Xist expression. Our results collectively suggest that structural changes in the Xic occurred in the ancestral lineage of XO/XO species, likely due to the loss of one X chromosome and the Y chromosome as a consequence of the degradation of the XCI system.

Behaviors of nucleosomes with mutant histone H4s in euchromatic domains of living human cells.

Since Robert Feulgen first stained DNA in the cell, visualizing genome chromatin has been a central issue in cell biology to uncover how chromatin is organized and behaves in the cell. To approach this issue, we have developed single-molecule imaging of nucleosomes, a basic unit of chromatin, to unveil local nucleosome behavior in living cells. In this study, we investigated behaviors of nucleosomes with various histone H4 mutants in living HeLa cells to address the role of H4 tail acetylation, including H4K16Ac and others, which are generally associated with more transcriptionally active chromatin regions. We ectopically expressed wild-type (wt) or mutated H4s (H4K16 point; H4K5,8,12,16 quadruple; and H4 tail deletion) fused with HaloTag in HeLa cells. Cells that expressed wtH4-Halo, H4K16-Halo mutants, and multiple H4-Halo mutants had euchromatin-concentrated distribution. Consistently, the genomic regions of the wtH4-Halo nucleosomes corresponded to Hi-C contact domains (or topologically associating domains, TADs) with active chromatin marks (A-compartment). Utilizing single-nucleosome imaging, we found that none of the H4 deacetylation or acetylation mimicked H4 mutants altered the overall local nucleosome motion. This finding suggests that H4 mutant nucleosomes embedded in the condensed euchromatic domains with excess endogenous H4 nucleosomes cannot cause an observable change in the local motion. Interestingly, H4 with four lysine-to-arginine mutations displayed a substantial freely diffusing fraction in the nucleoplasm, whereas H4 with a truncated N-terminal tail was incorporated in heterochromatic regions as well as euchromatin. Our study indicates the power of single-nucleosome imaging to understand individual histone/nucleosome behavior reflecting chromatin environments in living cells.

W chromosome sequences of two bombycid moths provide an insight into the origin of Fem.

Fem is a W-linked gene that encodes a piRNA precursor, and its product, Fem piRNA, is a master factor of female determination in Bombyx mori. Fem has low similarity to any known sequences, and the origin of Fem remains unclear. So far, two hypotheses have been proposed for the origin of Fem: The first hypothesis is that Fem is an allele of Masc, which assumes that the W chromosome was originally a homologous chromosome of the Z chromosome. The second hypothesis is that Fem arose by the transposition of Masc to the W chromosome. To explore the origin of Fem, we determined the W chromosome sequences of B. mori and, as a comparison, a closely relative bombycid species of Trilocha varians with a Fem-independent sex determination system. To our surprise, although the sequences of W and Z chromosomes show no homology to each other, a few pairs of homologues are shared by W and Z chromosomes, indicating the W chromosome of both species originated from Z chromosome. In addition, the W chromosome of T. varians lacks Fem, while the W chromosome of B. mori has over 100 copies of Fem. The high-quality assembly of the W chromosome of B. mori arose the third hypothesis about the origin of Fem: Fem is a chimeric sequence of multiple transposons. More than half of one transcriptional unit of Fem shows a significant homology to RTE-BovB. Moreover, the Fem piRNA-producing region could correspond to the boundary of the two transposons, gypsy and satellite DNA.