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.

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.

Three genomes in the algal genus Volvox reveal the fate of a haploid sex-determining region after a transition to homothallism.

Transitions between separate sexes (dioecy) and other mating systems are common across eukaryotes. Here, we study a change in a haploid dioecious green algal species with male- and female-determining chromosomes (U and V). The genus Volvox is an oogamous (with large, immotile female gametes and small, motile male gametes) and includes both heterothallic species (with distinct male and female genotypes, associated with a mating-type system that prevents fusion of gametes of the same sex) and homothallic species (bisexual, with the ability to self-fertilize). We date the origin of an expanded sex-determining region (SDR) in Volvox to at least 75 Mya, suggesting that homothallism represents a breakdown of dioecy (heterothallism). We investigated the involvement of the SDR of the U and V chromosomes in this transition. Using de novo whole-genome sequences, we identified a heteromorphic SDR of ca 1 Mbp in male and female genotypes of the heterothallic species Volvox reticuliferus and a homologous region (SDLR) in the closely related homothallic species Volvox africanus, which retained several different hallmark features of an SDR. The V. africanus SDLR includes a large region resembling the female SDR of the presumptive heterothallic ancestor, whereas most genes from the male SDR are absent. However, we found a multicopy array of the male-determining gene, MID, in a different genomic location from the SDLR. Thus, in V. africanus, an ancestrally female genotype may have acquired MID and thereby gained male traits.

Identification of small-sized intrachromosomal segments at the ends of INV-DUP-DEL patterns.

The mechanism of chromosomal rearrangement associated with inverted-duplication-deletion (INV-DUP-DEL) pattern formation has been investigated by many researchers, and several possible mechanisms have been proposed. Currently, fold-back and subsequent dicentric chromosome formation has been established as non-recurrent INV-DUP-DEL pattern formation mechanisms. In the present study, we analyzed the breakpoint junctions of INV-DUP-DEL patterns in five patients using long-read whole-genome sequencing and detected 2.2-6.1 kb copy-neutral regions in all five patients. At the end of the INV-DUP-DEL, two patients exhibited chromosomal translocations, which are recognized as telomere capture, and one patient showed direct telomere healing. The remaining two patients had additional small-sized intrachromosomal segments at the end of the derivative chromosomes. These findings have not been previously reported but they may only be explained by the presence of telomere capture breakage. Further investigations are required to better understand the mechanisms underlying this finding.

Long-read sequence analysis for clustered genomic copy number aberrations revealed architectures of intricately intertwined rearrangements.

Most chromosomal aberrations revealed by chromosomal microarray testing (CMA) are simple; however, very complex chromosomal structural rearrangements can also be found. Although the mechanism of structural rearrangements has been gradually revealed, not all mechanisms have been elucidated. We analyzed the breakpoint-junctions (BJs) of two or more clustered copy number variations (CNVs) in the same chromosome arms to understand their conformation and the mechanism of complex structural rearrangements. Combining CMA with long-read whole-genome sequencing (WGS) analysis, we successfully determined all BJs for the clustered CNVs identified in four patients. Multiple CNVs were intricately intertwined with each other, and clustered CNVs in four patients were involved in global complex chromosomal rearrangements. The BJs of two clustered deletions identified in two patients showed microhomologies, and their characteristics were explained by chromothripsis. In contrast, the BJs in the other two patients, who showed clustered deletions and duplications, consisted of blunt-end and nontemplated insertions. These findings could be explained only by alternative nonhomologous end-joining, a mechanism related to polymerase theta. All the patients had at least one inverted segment. Three patients showed cryptic aberrations involving a disruption and a deletion/duplication, which were not detected by CMA but were first identified by WGS. This result suggested that complex rearrangements should be considered if clustered CNVs are observed in the same chromosome arms. Because CMA has potential limitations in genotype-phenotype correlation analysis, a more detailed analysis by whole genome examination is recommended in cases of suspected complex structural aberrations.

The mechanoreceptor DEG-1 regulates cold tolerance in Caenorhabditis elegans.

Caenorhabditis elegans mechanoreceptors located in ASG sensory neurons have been found to sense ambient temperature, which is a key trait for animal survival. Here, we show that experimental loss of xanthine dehydrogenase (XDH-1) function in AIN and AVJ interneurons results in reduced cold tolerance and atypical neuronal response to changes in temperature. These interneurons connect with upstream neurons such as the mechanoreceptor-expressing ASG. Ca2+ imaging revealed that ASG neurons respond to warm temperature via the mechanoreceptor DEG-1, a degenerin/epithelial Na+ channel (DEG/ENaC), which in turn affects downstream AIN and AVJ circuits. Ectopic expression of DEG-1 in the ASE gustatory neuron results in the acquisition of warm sensitivity, while electrophysiological analysis revealed that DEG-1 and human MDEG1 were involved in warm sensation. Taken together, these results suggest that cold tolerance is regulated by mechanoreceptor-mediated circuit calculation.

Comparative analysis using the draft genome sequence of California poppy (Eschscholzia californica) for exploring the candidate genes involved in benzylisoquinoline alkaloid biosynthesis.

Genome characterization of California poppy (Eschscholzia californica cv. "Hitoezaki"), which produces pharmaceutically important benzylisoquinoline alkaloids (BIAs), was carried out using the draft genome sequence. The numbers of tRNA and rRNA genes were close to those of the other plant species tested, whereas the frequency of repetitive sequences was distinct from those species. Comparison of the predicted genes with those of Amborella trichopoda, Nelumbo nucifera, Solanum lycopersicum, and Arabidopsis thaliana, and analyses of gene ontology and Kyoto Encyclopedia of Genes and Genomes pathway indicated that the enzyme genes involved in BIA biosynthesis were highly enriched in the California poppy genome. Further comparative analysis using the genome information of Papaver somniferum and Aquilegia coerulea, both BIA-producing plants, revealed that many genes encoding BIA biosynthetic enzymes, transcription factors, transporters, and candidate proteins, possibly related to BIA biosynthesis, were specifically distributed in these plant species.

Endoribonuclease ENDU-2 regulates multiple traits including cold tolerance via cell autonomous and nonautonomous controls in Caenorhabditis elegans.

Environmental temperature acclimation is essential to animal survival, yet thermoregulation mechanisms remain poorly understood. We demonstrate cold tolerance in Caenorhabditis elegans as regulated by paired ADL chemosensory neurons via Ca2+-dependent endoribonuclease (EndoU) ENDU-2. Loss of ENDU-2 function results in life span, brood size, and synaptic remodeling abnormalities in addition to enhanced cold tolerance. Enzymatic ENDU-2 defects localized in the ADL and certain muscle cells led to increased cold tolerance in endu-2 mutants. Ca2+ imaging revealed ADL neurons were responsive to temperature stimuli through transient receptor potential (TRP) channels, concluding that ADL function requires ENDU-2 action in both cell-autonomous and cell-nonautonomous mechanisms. ENDU-2 is involved in caspase expression, which is central to cold tolerance and synaptic remodeling in dorsal nerve cord. We therefore conclude that ENDU-2 regulates cell type-dependent, cell-autonomous, and cell-nonautonomous cold tolerance.

A genome database for a Japanese population of the larvacean Oikopleura dioica.

The larvacean Oikopleura dioica is a planktonic chordate and is a tunicate that belongs to the closest relatives to vertebrates. Its simple and transparent body, invariant embryonic cell lineages, and short life cycle of 5 days make it a promising model organism for the study of developmental biology. The genome browser OikoBase was established in 2013 using Norwegian O. dioica. However, genome information for other populations is not available, even though many researchers have studied local populations. In the present study, we sequenced using Illumina and PacBio RSII technologies the genome of O. dioica from a southwestern Japanese population that was cultured in our laboratory for 3 years. The genome of Japanese O. dioica was assembled into 576 scaffold sequences with a total length and N50 length of 56.6 and 1.5 Mb, respectively. A total of 18,743 gene models (transcript models) were predicted in the genome assembly, named OSKA2016. In addition, 19,277 non-redundant transcripts were assembled using RNA-seq data. The OSKA2016 has global sequence similarity of only 86.5% when compared with the OikoBase, highlighting the sequence difference between the two far distant O. dioica populations on the globe. The genome assembly, transcript assembly, and transcript models were incorporated into ANISEED (https://www.aniseed.cnrs.fr/) for genome browsing and BLAST searches. Mapping of reads obtained from male- or female-specific genome libraries yielded male-specific scaffolds in the OSKA2016 and revealed that over 2.6 Mb of sequence were included in the male-specific Y-region. The genome and transcriptome resources from two distinct populations will be useful datasets for developmental biology, evolutionary biology, and molecular ecology using this model organism.

Tengunoibacter tsumagoiensis gen. nov., sp. nov., Dictyobacter kobayashii sp. nov., Dictyobacter alpinus sp. nov., and description of Dictyobacteraceae fam. nov. within the order Ktedonobacterales isolated from Tengu-no-mugimeshi, a soil-like granular mass of micro-organisms, and emended descriptions of the genera Ktedonobacter and Dictyobacter.

Three mesophilic, Gram-stain-positive, aerobic bacterial strains, designated Uno3T, Uno11T and Uno16T, were isolated from a soil-like granular micro-organism mass (termed Tengu-no-mugimeshi) collected from Tsumagoi, Gunma, Japan. They grow at 11-37 °C and pH 4.0-8.0, form branched mycelia, and have a G+C content between 49.4-50.3 mol%. The major menaquinone and fatty acid of Uno3T are MK-9 and iso-C16 : 0, respectively, whereas Uno11T and Uno16T share MK-9 (H2) and C16 : 1-2OH. The major cell-wall sugars are mannose (Uno3T and Uno11T) and glucose (Uno16T). Phylogenetic analysis based on 16S rRNA gene sequences indicated that these three strains belong to the order Ktedonobacterales and are most closely related to Dictyobacter aurantiacus S-27T (sequence similarity of 91.3, 96.4 and 95.5 %). Average nucleotide identity values were <79.9 % among Uno11T, Uno16T and D. aurantiacus S-27T, well below the 95-96 % species circumscription threshold. Based on phenotypic features and phylogenetic positions, we propose that Uno3T represents a novel genus and species, Tengunoibacter tsumagoiensis gen. nov., sp. nov. (type strain Uno3T=NBRC 113152T=LMG 30471T=BCRC 81113T) within the new family Dictyobacteraceae fam. nov. Strains Uno11T and Uno16T are also considered to represent novel species: Dictyobacterkobayashii sp. nov. (type strain Uno11T=NBRC 113153T=LMG 30472T=BCRC 81114T) and Dictyobacteralpinus sp. nov. (type strain Uno16T=NBRC 113154T=BCRC 81115T). We also propose an emended description of the genus Dictyobacter, classifying it within family Dictyobacteraceae, and provide emended descriptions of the genera Dictyobacter and Ktedonobacter.

Mining of the Uncharacterized Cytochrome P450 Genes Involved in Alkaloid Biosynthesis in California Poppy Using a Draft Genome Sequence.

Land plants produce specialized low molecular weight metabolites to adapt to various environmental stressors, such as UV radiation, pathogen infection, wounding and animal feeding damage. Due to the large variety of stresses, plants produce various chemicals, particularly plant species-specific alkaloids, through specialized biosynthetic pathways. In this study, using a draft genome sequence and querying known biosynthetic cytochrome P450 (P450) enzyme-encoding genes, we characterized the P450 genes involved in benzylisoquinoline alkaloid (BIA) biosynthesis in California poppy (Eschscholzia californica), as P450s are key enzymes involved in the diversification of specialized metabolism. Our in silico studies showed that all identified enzyme-encoding genes involved in BIA biosynthesis were found in the draft genome sequence of approximately 489 Mb, which covered approximately 97% of the whole genome (502 Mb). Further analyses showed that some P450 families involved in BIA biosynthesis, i.e. the CYP80, CYP82 and CYP719 families, were more enriched in the genome of E. californica than in the genome of Arabidopsis thaliana, a plant that does not produce BIAs. CYP82 family genes were highly abundant, so we measured the expression of CYP82 genes with respect to alkaloid accumulation in different plant tissues and two cell lines whose BIA production differs to estimate the functions of the genes. Further characterization revealed two highly homologous P450s (CYP82P2 and CYP82P3) that exhibited 10-hydroxylase activities with different substrate specificities. Here, we discuss the evolution of the P450 genes and the potential for further genome mining of the genes encoding the enzymes involved in BIA biosynthesis.

Germline mutation rates and the long-term phenotypic effects of mutation accumulation in wild-type laboratory mice and mutator mice.

The germline mutation rate is an important parameter that affects the amount of genetic variation and the rate of evolution. However, neither the rate of germline mutations in laboratory mice nor the biological significance of the mutation rate in mammalian populations is clear. Here we studied genome-wide mutation rates and the long-term effects of mutation accumulation on phenotype in more than 20 generations of wild-type C57BL/6 mice and mutator mice, which have high DNA replication error rates. We estimated the base-substitution mutation rate to be 5.4 × 10(-9) (95% confidence interval = 4.6 × 10(-9)-6.5 × 10(-9)) per nucleotide per generation in C57BL/6 laboratory mice, about half the rate reported in humans. The mutation rate in mutator mice was 17 times that in wild-type mice. Abnormal phenotypes were 4.1-fold more frequent in the mutator lines than in the wild-type lines. After several generations, the mutator mice reproduced at substantially lower rates than the controls, exhibiting low pregnancy rates, lower survival rates, and smaller litter sizes, and many of the breeding lines died out. These results provide fundamental information about mouse genetics and reveal the impact of germline mutation rates on phenotypes in a mammalian population.

Molecular Characterization of Eye Pigmentation-Related ABC Transporter Genes in the Ladybird Beetle Harmonia axyridis Reveals Striking Gene Duplication of the white Gene.

Many species of ladybird beetles (Coccinellidae) possess vivid body colors. These colors and patterns show diversity between coccinellid species, or even within species. However, the molecular underpinnings of these striking body colors are scarcely understood. One of the candidate pigmentation molecules responsible for ladybird body color is ommochrome pigment, which is well known as the red pigment molecule responsible for the red eyes of Drosophila. Various insects also use ommochrome in body coloration. It is known that ommochrome pigment precursors are imported into appropriate cells by the ATP binding cassette (ABC) transporter proteins White and Scarlet. Thus, these ABC transporter genes are potentially involved in various color and pattern expressions seen in ladybird beetle species. In this study, in order to identify the repertory of ABC transporter genes responsible for such body colors, we performed molecular characterization of pigment-related ABC transporter genes, especially white and scarlet, in the coccinellid Harmonia axyridis. By using whole genome data for H. axyridis and subsequent RACE-PCR, six white orthologs and one scarlet ortholog were successfully identified. According to the results of functional analyses via RNA interference (RNAi), only one of these genes had a major function in eye pigmentation. Specific effects on body color and pattern were not detected by our RNAi experiments of any of these genes. This is the first report of this striking duplication of white genes and their functional analyses in insects.

Genome-wide analysis of intraspecific DNA polymorphism in 'Micro-Tom', a model cultivar of tomato (Solanum lycopersicum).

Tomato (Solanum lycopersicum) is regarded as a model plant of the Solanaceae family. The genome sequencing of the tomato cultivar 'Heinz 1706' was recently completed. To accelerate the progress of tomato genomics studies, systematic bioresources, such as mutagenized lines and full-length cDNA libraries, have been established for the cultivar 'Micro-Tom'. However, these resources cannot be utilized to their full potential without the completion of the genome sequencing of 'Micro-Tom'. We undertook the genome sequencing of 'Micro-Tom' and here report the identification of single nucleotide polymorphisms (SNPs) and insertion/deletions (indels) between 'Micro-Tom' and 'Heinz 1706'. The analysis demonstrated the presence of 1.23 million SNPs and 0.19 million indels between the two cultivars. The density of SNPs and indels was high in chromosomes 2, 5 and 11, but was low in chromosomes 6, 8 and 10. Three known mutations of 'Micro-Tom' were localized on chromosomal regions where the density of SNPs and indels was low, which was consistent with the fact that these mutations were relatively new and introgressed into 'Micro-Tom' during the breeding of this cultivar. We also report SNP analysis for two 'Micro-Tom' varieties that have been maintained independently in Japan and France, both of which have served as standard lines for 'Micro-Tom' mutant collections. Approximately 28,000 SNPs were identified between these two 'Micro-Tom' lines. These results provide high-resolution DNA polymorphic information on 'Micro-Tom' and represent a valuable contribution to the 'Micro-Tom'-based genomics resources.

Expanded male sex-determining region conserved during the evolution of homothallism in the green alga Volvox.

Male and female genotypes in heterothallic (self-incompatible) species of haploid organisms, such as algae and bryophytes, are generally determined by male and female sex-determining regions (SDRs) in the sex chromosomes. To resolve the molecular genetic basis for the evolution of homothallic (bisexual and self-compatible) species from a heterothallic ancestor, we compared whole-genome data from Thai and Japanese genotypes within the homothallic green alga Volvox africanus. The Thai and Japanese algae harbored expanded ancestral male and female SDRs of ∼1 Mbp each, representing a direct heterothallic ancestor. Therefore, the expanded male and female ancestral SDRs may originate from the ancient (∼75 mya) heterothallic ancestor, and either might have been conserved during the evolution of each homothallic genotype. An expanded SDR-like region seems essential for homothallic sexual reproduction in V. africanus, irrespective of male or female origin. Our study stimulates future studies to elucidate the biological significance of such expanded genomic regions.

Erratum: Expanded male sex-determining region conserved during the evolution of homothallism in the green alga Volvox.

[This corrects the article DOI: 10.1016/j.isci.2023.106893.].

Genome sequencing of the multicellular alga Astrephomene provides insights into convergent evolution of germ-soma differentiation.

Germ-soma differentiation evolved independently in many eukaryotic lineages and contributed to complex multicellular organizations. However, the molecular genetic bases of such convergent evolution remain unresolved. Two multicellular volvocine green algae, Volvox and Astrephomene, exhibit convergent evolution of germ-soma differentiation. The complete genome sequence is now available for Volvox, while genome information is scarce for Astrephomene. Here, we generated the de novo whole genome sequence of Astrephomene gubernaculifera and conducted RNA-seq analysis of isolated somatic and reproductive cells. In Volvox, tandem duplication and neofunctionalization of the ancestral transcription factor gene (RLS1/rlsD) might have led to the evolution of regA, the master regulator for Volvox germ-soma differentiation. However, our genome data demonstrated that Astrephomene has not undergone tandem duplication of the RLS1/rlsD homolog or acquisition of a regA-like gene. Our RNA-seq analysis revealed the downregulation of photosynthetic and anabolic gene expression in Astrephomene somatic cells, as in Volvox. Among genes with high expression in somatic cells of Astrephomene, we identified three genes encoding putative transcription factors, which may regulate somatic cell differentiation. Thus, the convergent evolution of germ-soma differentiation in the volvocine algae may have occurred by the acquisition of different regulatory circuits that generate a similar division of labor.

A chromosome-level genome sequence of Chrysanthemum seticuspe, a model species for hexaploid cultivated chrysanthemum.

Chrysanthemums are one of the most industrially important cut flowers worldwide. However, their segmental allopolyploidy and self-incompatibility have prevented the application of genetic analysis and modern breeding strategies. We thus developed a model strain, Gojo-0 (Chrysanthemum seticuspe), which is a diploid and self-compatible pure line. Here, we present the 3.05 Gb chromosome-level reference genome sequence, which covered 97% of the C. seticuspe genome. The genome contained more than 80% interspersed repeats, of which retrotransposons accounted for 72%. We identified recent segmental duplication and retrotransposon expansion in C. seticuspe, contributing to arelatively large genome size. Furthermore, we identified a retrotransposon family, SbdRT, which was enriched in gene-dense genome regions and had experienced a very recent transposition burst. We also demonstrated that the chromosome-level genome sequence facilitates positional cloning in C. seticuspe. The genome sequence obtained here can greatly contribute as a reference for chrysanthemum in front-line breeding including genome editing.

Comparative genomic analyses illuminate the distinct evolution of megabats within Chiroptera.

The revision of the sub-order Microchiroptera is one of the most intriguing outcomes in recent mammalian molecular phylogeny. The unexpected sister-taxon relationship between rhinolophoid microbats and megabats, with the exclusion of other microbats, suggests that megabats arose in a relatively short period of time from a microbat-like ancestor. In order to understand the genetic mechanism underlying adaptive evolution in megabats, we determined the whole-genome sequences of two rousette megabats, Leschenault's rousette (Rousettus leschenaultia) and the Egyptian fruit bat (R. aegyptiacus). The sequences were compared with those of 22 other mammals, including nine bats, available in the database. We identified that megabat genomes are distinct in that they have extremely low activity of SINE retrotranspositions, expansion of two chemosensory gene families, including the trace amine receptor (TAAR) and olfactory receptor (OR), and elevation of the dN/dS ratio in genes for immunity and protein catabolism. The adaptive signatures discovered in the genomes of megabats may provide crucial insight into their distinct evolution, including key processes such as virus resistance, loss of echolocation, and frugivorous feeding.

Characteristics of induced mutations in offspring derived from irradiated mouse spermatogonia and mature oocytes.

The exposure of germ cells to radiation introduces mutations in the genomes of offspring, and a previous whole-genome sequencing study indicated that the irradiation of mouse sperm induces insertions/deletions (indels) and multisite mutations (clustered single nucleotide variants and indels). However, the current knowledge on the mutation spectra is limited, and the effects of radiation exposure on germ cells at stages other than the sperm stage remain unknown. Here, we performed whole-genome sequencing experiments to investigate the exposure of spermatogonia and mature oocytes. We compared de novo mutations in a total of 24 F1 mice conceived before and after the irradiation of their parents. The results indicated that radiation exposure, 4 Gy of gamma rays, induced 9.6 indels and 2.5 multisite mutations in spermatogonia and 4.7 indels and 3.1 multisite mutations in mature oocytes in the autosomal regions of each F1 individual. Notably, we found two types of deletions, namely, small deletions (mainly 1~12 nucleotides) in non-repeat sequences, many of which showed microhomology at the breakpoint junction, and single-nucleotide deletions in mononucleotide repeat sequences. The results suggest that these deletions and multisite mutations could be a typical signature of mutations induced by parental irradiation in mammals.