Diversity of genome size and chromosome number in homothallic and heterothallic strains of the Closterium peracerosum-strigosum-littorale complex (Desmidiales, Zygnematophyceae, Streptophyta).

The evolutionary transitions of mating systems between outcrossing and self-fertilization are often suggested to associate with the cytological and genomic changes, but the empirical reports are limited in multicellular organisms. Here we used the unicellular zygnematophycean algae, the Closterium peracerosum-strigosum-littorale (C. psl.) complex, to address whether genomic properties such as genome sizes and chromosome numbers are associated with mating system transitions between homothallism (self-fertility) and heterothallism (self-sterility). Phylogenetic analysis revealed the polyphyly of homothallic strains, suggesting multiple transitions between homothallism and heterothallism in the C. psl. complex. Flow cytometry analysis identified a more than 2-fold genome size variation, ranging from 0.53 to 1.42 Gbp, which was positively correlated with chromosome number variation between strains. Although we did not find consistent trends in genome size change and mating system transitions, the mean chromosome sizes tend to be smaller in homothallic strains than in their relative heterothallic strains. This result suggests that homothallic strains possibly have more fragmented chromosomes, which is consistent with the argument that self-fertilizing populations may tolerate more chromosomal rearrangements.

Extensive Copy Number Variation Explains Genome Size Variation in the Unicellular Zygnematophycean Alga, Closterium peracerosum-strigosum-littorale Complex.

Genome sizes are known to vary within and among closely related species, but the knowledge about genomic factors contributing to the variation and their impacts on gene functions is limited to only a small number of species. This study identified a more than 2-fold heritable genome size variation among the unicellular Zygnematophycean alga, Closterium peracerosum-strigosum-littorale (C. psl.) complex, based on short-read sequencing analysis of 22 natural strains and F1 segregation analysis. Six de novo assembled genomes revealed that genome size variation is largely attributable to genome-wide copy number variation (CNV) among strains rather than mating type-linked genomic regions or specific repeat sequences such as rDNA. Notably, about 30% of genes showed CNV even between strains that can mate with each other. Transcriptome and gene ontology analysis demonstrated that CNV is distributed nonrandomly in terms of gene functions, such that CNV was more often observed in the gene set with stage-specific expression. Furthermore, in about 30% of these genes with CNV, the expression level does not increase proportionally with the gene copy number, suggesting presence of dosage compensation, which was overrepresented in genes involved in basic biological functions, such as translation. Nonrandom patterns in gene duplications and corresponding expression changes in terms of gene functions may contribute to maintaining the high level of CNV associated with extensive genome size variation in the C. psl. complex, despite its possible detrimental effects.

The draft genome sequence of the Japanese rhinoceros beetle Trypoxylus dichotomus septentrionalis towards an understanding of horn formation.

The Japanese rhinoceros beetle Trypoxylus dichotomus is a giant beetle with distinctive exaggerated horns present on the head and prothoracic regions of the male. T. dichotomus has been used as a research model in various fields such as evolutionary developmental biology, ecology, ethology, biomimetics, and drug discovery. In this study, de novo assembly of 615 Mb, representing 80% of the genome estimated by flow cytometry, was obtained using the 10 × Chromium platform. The scaffold N50 length of the genome assembly was 8.02 Mb, with repetitive elements predicted to comprise 49.5% of the assembly. In total, 23,987 protein-coding genes were predicted in the genome. In addition, de novo assembly of the mitochondrial genome yielded a contig of 20,217 bp. We also analyzed the transcriptome by generating 16 RNA-seq libraries from a variety of tissues of both sexes and developmental stages, which allowed us to identify 13 co-expressed gene modules. We focused on the genes related to horn formation and obtained new insights into the evolution of the gene repertoire and sexual dimorphism as exemplified by the sex-specific splicing pattern of the doublesex gene. This genomic information will be an excellent resource for further functional and evolutionary analyses, including the evolutionary origin and genetic regulation of beetle horns and the molecular mechanisms underlying sexual dimorphism.

The genome of Lyophyllum shimeji provides insight into the initial evolution of ectomycorrhizal fungal genomes.

Mycorrhizae are one of the most fundamental symbioses between plants and fungi, with ectomycorrhizae being the most widespread in boreal forest ecosystems. Ectomycorrhizal fungi are hypothesized to have evolved convergently from saprotrophic ancestors in several fungal clades, especially members of the subdivision Agaricomycotina. Studies on fungal genomes have identified several typical characteristics of mycorrhizal fungi, such as genome size expansion and decreases in plant cell-wall degrading enzymes (PCWDEs). However, genomic changes concerning the evolutionary transition to the ectomycorrhizal lifestyle are largely unknown. In this study, we sequenced the genome of Lyophyllum shimeji, an ectomycorrhizal fungus that is phylogenetically related to saprotrophic species and retains some saprotroph-like traits. We found that the genome of Ly. shimeji strain AT787 lacks both incremental increases in genome size and reduced numbers of PCWDEs. Our findings suggest that the previously reported common genomic traits of mycorrhizal fungi are not essential for the ectomycorrhizal lifestyle, but are a result of abolishing saprotrophic activity. Since Ly. shimeji is commercially consumed as an edible mushroom, the newly available genomic information may also impact research designed to enhance the cultivation of this mushroom.

GRAS transcription factors regulate cell division planes in moss overriding the default rule.

Plant cells are surrounded by a cell wall and do not migrate, which makes the regulation of cell division orientation crucial for development. Regulatory mechanisms controlling cell division orientation may have contributed to the evolution of body organization in land plants. The GRAS family of transcription factors was transferred horizontally from soil bacteria to an algal common ancestor of land plants. SHORTROOT (SHR) and SCARECROW (SCR) genes in this family regulate formative periclinal cell divisions in the roots of flowering plants, but their roles in nonflowering plants and their evolution have not been studied in relation to body organization. Here, we show that SHR cell autonomously inhibits formative periclinal cell divisions indispensable for leaf vein formation in the moss Physcomitrium patens, and SHR expression is positively and negatively regulated by SCR and the GRAS member LATERAL SUPPRESSOR, respectively. While precursor cells of a leaf vein lacking SHR usually follow the geometry rule of dividing along the division plane with the minimum surface area, SHR overrides this rule and forces cells to divide nonpericlinally. Together, these results imply that these bacterially derived GRAS transcription factors were involved in the establishment of the genetic regulatory networks modulating cell division orientation in the common ancestor of land plants and were later adapted to function in flowering plant and moss lineages for their specific body organizations.

A divergent RWP-RK transcription factor determines mating type in heterothallic Closterium.

The Closterium peracerosum-strigosum-littorale complex (Closterium, Zygnematophyceae) has an isogamous mating system. Members of the Zygnematophyceae are the closest relatives to extant land plants and are distantly related to chlorophytic models, for which a genetic basis of mating type (MT) determination has been reported. We thus investigated MT determination in Closterium. We sequenced genomes representing the two MTs, mt+ and mt-, in Closterium and identified CpMinus1, a gene linked to the mt- phenotype. We analyzed its function using reverse genetics methods. CpMinus1 encodes a divergent RWP-RK domain-containing-like transcription factor and is specifically expressed during gamete differentiation. Introduction of CpMinus1 into an mt+ strain was sufficient to convert it to a phenotypically mt- strain, while CpMinus1-knockout mt- strains were phenotypically mt+. We propose that CpMinus1 is the major MT determinant that acts by evoking the mt- phenotype and suppressing the mt+ phenotype in heterothallic Closterium. CpMinus1 likely evolved independently in the Zygnematophyceae lineage, which lost an egg-sperm anisogamous mating system. mt- specific regions possibly constitute an MT locus flanked by common sequences that undergo some recombination.

An ancestral function of strigolactones as symbiotic rhizosphere signals.

In flowering plants, strigolactones (SLs) have dual functions as hormones that regulate growth and development, and as rhizosphere signaling molecules that induce symbiosis with arbuscular mycorrhizal (AM) fungi. Here, we report the identification of bryosymbiol (BSB), an SL from the bryophyte Marchantia paleacea. BSB is also found in vascular plants, indicating its origin in the common ancestor of land plants. BSB synthesis is enhanced at AM symbiosis permissive conditions and BSB deficient mutants are impaired in AM symbiosis. In contrast, the absence of BSB synthesis has little effect on the growth and gene expression. We show that the introduction of the SL receptor of Arabidopsis renders M. paleacea cells BSB-responsive. These results suggest that BSB is not perceived by M. paleacea cells due to the lack of cognate SL receptors. We propose that SLs originated as AM symbiosis-inducing rhizosphere signaling molecules and were later recruited as plant hormone.

Phylogenetic distribution and expression pattern analyses identified a divergent basal body assembly protein involved in land plant spermatogenesis.

Land plant spermatozoids commonly possess characteristic structures such as the spline, which consists of a microtubule array, the multilayered structure (MLS) in which the uppermost layer is a continuum of the spline, and multiple flagella. However, the molecular mechanisms underpinning spermatogenesis remain to be elucidated. We successfully identified candidate genes involved in spermatogenesis, deeply divergent BLD10s, by computational analyses combining multiple methods and omics data. We then examined the functions of BLD10s in the liverwort Marchantia polymorpha and the moss Physcomitrium patens. MpBLD10 and PpBLD10 are required for normal basal body (BB) and flagella formation. Mpbld10 mutants exhibited defects in remodeling of the cytoplasm and nucleus during spermatozoid formation, and thus MpBLD10 should be involved in chromatin reorganization and elimination of the cytoplasm during spermiogenesis. We identified orthologs of MpBLD10 and PpBLD10 in diverse Streptophyta and found that MpBLD10 and PpBLD10 are orthologous to BLD10/CEP135 family proteins, which function in BB assembly. However, BLD10s evolved especially quickly in land plants and MpBLD10 might have acquired additional functions in spermatozoid formation through rapid molecular evolution.

Discovery of ultrafast myosin, its amino acid sequence, and structural features.

Cytoplasmic streaming with extremely high velocity (∼70 µm s-1) occurs in cells of the characean algae (Chara). Because cytoplasmic streaming is caused by myosin XI, it has been suggested that a myosin XI with a velocity of 70 µm s-1, the fastest myosin measured so far, exists in Chara cells. However, the velocity of the previously cloned Chara corallina myosin XI (CcXI) was about 20 µm s-1, one-third of the cytoplasmic streaming velocity in Chara Recently, the genome sequence of Chara braunii has been published, revealing that this alga has four myosin XI genes. We cloned these four myosin XI (CbXI-1, 2, 3, and 4) and measured their velocities. While the velocities of CbXI-3 and CbXI-4 motor domains (MDs) were similar to that of CcXI MD, the velocities of CbXI-1 and CbXI-2 MDs were 3.2 times and 2.8 times faster than that of CcXI MD, respectively. The velocity of chimeric CbXI-1, a functional, full-length CbXI-1 construct, was 60 µm s-1 These results suggest that CbXI-1 and CbXI-2 would be the main contributors to cytoplasmic streaming in Chara cells and show that these myosins are ultrafast myosins with a velocity 10 times faster than fast skeletal muscle myosins in animals. We also report an atomic structure (2.8-Å resolution) of myosin XI using X-ray crystallography. Based on this crystal structure and the recently published cryo-electron microscopy structure of acto-myosin XI at low resolution (4.3-Å), it appears that the actin-binding region contributes to the fast movement of Chara myosin XI. Mutation experiments of actin-binding surface loops support this hypothesis.

Elevated mutation rates underlie the evolution of the aquatic plant family Podostemaceae.

Molecular evolutionary rates vary among lineages and influence the evolutionary process. Here, we report elevated genome-wide mutation rates in Podostemaceae, a family of aquatic plants with a unique body plan that allows members to live on submerged rocks in fast-flowing rivers. Molecular evolutionary analyses using 1640 orthologous gene groups revealed two historical increases in evolutionary rates: the first at the emergence of the family and the second at the emergence of Podostemoideae, which is the most diversified subfamily. In both branches, synonymous substitution rates were elevated, indicating higher mutation rates. On early branches, mutations were biased in favour of AT content, which is consistent with a role for ultraviolet light-induced mutation and habitat shift. In ancestors of Podostemoideae, DNA-repair genes were enriched in genes under positive selection, which may have responded to the meristem architectural changes.

Upregulation of the Long Noncoding RNA HULC by Hepatitis C Virus and Its Regulation of Viral Replication.

BACKGROUND: Many long noncoding RNAs (lncRNAs) have important roles in biological processes. The lncRNA HULC was found to be upregulated in human hepatoma tissues. HULC is thought to be involved in multiple steps of hepatoma development and progression; however, the relationship between HULC and hepatitis C virus (HCV) infection, which is a leading cause of hepatoma, remains unclear. METHODS: We examined the effect of HCV replication on HULC expression and the underlying mechanism using cell culture systems. Subsequently, we tested the effect of HULC suppression and overexpression on HCV replication. Finally, we examined the impact of HCV eradication on HULC expression using human liver tissue and blood samples. RESULTS: HCV replication increased HULC expression in cell cultures. A promoter assay showed that an HCV nonstructural protein, NS5A, increased HULC transcription. HULC suppression inhibited HCV replication; conversely, its overexpression enhanced HCV replication. These effects on HCV replication seemed to occur by the modification of HCV translation. Measurements from human liver and blood samples showed that HCV eradication significantly reduced HULC levels in the liver and blood. CONCLUSIONS: HCV infection increases HULC expression in vitro and in vivo. HULC modulates HCV replication through an HCV internal ribosome entry site-directed translation step.

IFN-γ‒Induced APOBEC3B Contributes to Merkel Cell Polyomavirus Genome Mutagenesis in Merkel Cell Carcinoma.

Merkel cell polyomavirus (MCPyV) is the causative agent of an aggressive skin tumor, Merkel cell carcinoma. The viral genome is integrated into the tumor genome and harbors nonsense mutations in the helicase domain of large T antigen. However, the molecular mechanisms by which the viral genome gains the tumor-specific mutations remain to be elucidated. Focusing on host cytosine deaminases APOBEC3s, we find that A3A, A3B, or A3G introduces A3-specific mutations into episomal MCPyV genomes in MCPyV-replicating 293-derivative cells. Sequence analysis of MCPyV genomes retrieved from the NCBI database revealed a decrease of TpC dinucleotide, a preferred target for A3A and A3B, in the 3'-region of the large T antigen‒coding sequence. The viral DNA isolated from tumors contained mutated cytosines, with a remarkable bias toward TpC dinucleotide. Analysis of publicly available microarray data showed that expression of IFN-γ and cytotoxic T lymphocyte markers was positively correlated with the A3A, A3B, and A3G levels in MCPyV-positive but not in MCPyV-negative tumors. Finally, IFN-γ treatment induced A3B and A3G expression in the MCPyV-positive Merkel cell carcinoma cell line MS-1. These results suggest that the IFN-γ-A3B axis plays pivotal roles in evolutionally shaping MCPyV genomic sequences and in generating tumor-specific large T antigen mutations during development of Merkel cell carcinoma.

The complete plastid genome sequence of the enigmatic moss, Takakia lepidozioides (Takakiopsida, Bryophyta): evolutionary perspectives on the largest collection of genes in mosses and the intensive RNA editing.

KEY MESSAGE: Complete chloroplast genome sequence of a moss, Takakia lepidozioides (Takakiopsida) is reported. The largest collection of genes in mosses and the intensive RNA editing were discussed from evolutionary perspectives. We assembled the entire plastid genome sequence of Takakia lepidozioides (Takakiopsida), emerging from the first phylogenetic split among extant mosses. The genome sequences were assembled into a circular molecule 149,016 bp in length, with a quadripartite structure comprising a large and a small single-copy region separated by inverted repeats. It contained 88 genes coding for proteins, 32 for tRNA, four for rRNA, two open reading frames, and at least one pseudogene (tufA). This is the largest number of genes of all sequenced plastid genomes in mosses and Takakia is the only moss that retains the seven coding genes ccsA, cysA, cysT, petN rpoA, rps16 and trnPGGG. Parsimonious interpretation of gene loss suggests that the last common ancestor of bryophytes had all seven genes and that mosses lost at least three of them during their diversification. Analyses of the plastid transcriptome identified the extraordinary frequency of RNA editing with more than 1100 sites. We indicated a close correlation between the monoplastidy of vegetative tissue and the intensive RNA editing sites in the plastid genome in land plant lineages. Here, we proposed a hypothesis that the small population size of plastids in each vegetative cell of some early diverging land plants, including Takakia, might cause the frequent fixation of mutations in plastid genome through the intracellular genetic drift and that deleterious mutations might be continuously compensated by RNA editing during or following transcription.

The genome of Shorea leprosula (Dipterocarpaceae) highlights the ecological relevance of drought in aseasonal tropical rainforests.

Hyperdiverse tropical rainforests, such as the aseasonal forests in Southeast Asia, are supported by high annual rainfall. Its canopy is dominated by the species-rich tree family of Dipterocarpaceae (Asian dipterocarps), which has both ecological (e.g., supports flora and fauna) and economical (e.g., timber production) importance. Recent ecological studies suggested that rare irregular drought events may be an environmental stress and signal for the tropical trees. We assembled the genome of a widespread but near threatened dipterocarp, Shorea leprosula, and analyzed the transcriptome sequences of ten dipterocarp species representing seven genera. Comparative genomic and molecular dating analyses suggested a whole-genome duplication close to the Cretaceous-Paleogene extinction event followed by the diversification of major dipterocarp lineages (i.e. Dipterocarpoideae). Interestingly, the retained duplicated genes were enriched for genes upregulated by no-irrigation treatment. These findings provide molecular support for the relevance of drought for tropical trees despite the lack of an annual dry season.

An Agrobacterium-mediated stable transformation technique for the hornwort model Anthoceros agrestis.

Despite their key phylogenetic position and their unique biology, hornworts have been widely overlooked. Until recently there was no hornwort model species amenable to systematic experimental investigation. Anthoceros agrestis has been proposed as the model species to study hornwort biology. We have developed an Agrobacterium-mediated method for the stable transformation of A. agrestis, a hornwort model species for which a genetic manipulation technique was not yet available. High transformation efficiency was achieved by using thallus tissue grown under low light conditions. We generated a total of 274 transgenic A. agrestis lines expressing the ß-glucuronidase (GUS), cyan, green, and yellow fluorescent proteins under control of the CaMV 35S promoter and several endogenous promoters. Nuclear and plasma membrane localization with multiple color fluorescent proteins was also confirmed. The transformation technique described here should pave the way for detailed molecular and genetic studies of hornwort biology, providing much needed insight into the molecular mechanisms underlying symbiosis, carbon-concentrating mechanism, RNA editing and land plant evolution in general.

Chloroplast acquisition without the gene transfer in kleptoplastic sea slugs, Plakobranchus ocellatus.

Some sea slugs sequester chloroplasts from algal food in their intestinal cells and photosynthesize for months. This phenomenon, kleptoplasty, poses a question of how the chloroplast retains its activity without the algal nucleus. There have been debates on the horizontal transfer of algal genes to the animal nucleus. To settle the arguments, this study reported the genome of a kleptoplastic sea slug, Plakobranchus ocellatus, and found no evidence of photosynthetic genes encoded on the nucleus. Nevertheless, it was confirmed that light illumination prolongs the life of mollusk under starvation. These data presented a paradigm that a complex adaptive trait, as typified by photosynthesis, can be transferred between eukaryotic kingdoms by a unique organelle transmission without nuclear gene transfer. Our phylogenomic analysis showed that genes for proteolysis and immunity undergo gene expansion and are up-regulated in chloroplast-enriched tissue, suggesting that these molluskan genes are involved in the phenotype acquisition without horizontal gene transfer.

The genome sequence of Samia ricini, a new model species of lepidopteran insect.

Samia ricini, a gigantic saturniid moth, has the potential to be a novel lepidopteran model species. Samia ricini is far more resistant to diseases than the current model species Bombyx mori, and therefore can be more easily reared. In addition, genetic resources available for S. ricini rival those for B. mori: at least 26 ecoraces of S. ricini are reported and S. ricini can hybridize with wild Samia species, which are distributed throughout Asian countries, and produce fertile progenies. Physiological traits such as food preference, integument colour and larval spot pattern differ among S. ricini strains and wild Samia species so that those traits can be targeted in forward genetic analyses. To facilitate genetic research in S. ricini, we determined its whole genome sequence. The assembled genome of S. ricini was 458 Mb with 155 scaffolds, and the scaffold N50 length of the assembly was ~ 21 Mb. In total, 16,702 protein coding genes were predicted. While the S. ricini genome was mostly collinear with that of B. mori with some rearrangements and few S. ricini-specific genes were discovered, chorion genes and fibroin genes seemed to have expanded in the S. ricini lineage. As the first step of genetic analyses, causal genes for "Blue," "Yellow," "Spot," and "Red cocoon" phenotypes were mapped to chromosomes.

Corrigendum: Networking and Specificity-Changing DNA Methyltransferases in Helicobacter pylori.

[This corrects the article DOI: 10.3389/fmicb.2020.01628.].

Ethylene signaling mediates host invasion by parasitic plants.

Parasitic plants form a specialized organ, a haustorium, to invade host tissues and acquire water and nutrients. To understand the molecular mechanism of haustorium development, we performed a forward genetics screening to isolate mutants exhibiting haustorial defects in the model parasitic plant Phtheirospermum japonicum. We isolated two mutants that show prolonged and sometimes aberrant meristematic activity in the haustorium apex, resulting in severe defects on host invasion. Whole-genome sequencing revealed that the two mutants respectively have point mutations in homologs of ETHYLENE RESPONSE 1 (ETR1) and ETHYLENE INSENSITIVE 2 (EIN2), signaling components in response to the gaseous phytohormone ethylene. Application of the ethylene signaling inhibitors also caused similar haustorial defects, indicating that ethylene signaling regulates cell proliferation and differentiation of parasite cells. Genetic disruption of host ethylene production also perturbs parasite invasion. We propose that parasitic plants use ethylene as a signal to invade host roots.

Networking and Specificity-Changing DNA Methyltransferases in Helicobacter pylori.

Epigenetic DNA base methylation plays important roles in gene expression regulation. We here describe a gene expression regulation network consisting of many DNA methyltransferases each frequently changing its target sequence-specificity. Our object Helicobacter pylori, a bacterium responsible for most incidence of stomach cancer, carries a large and variable repertoire of sequence-specific DNA methyltransferases. By creating a dozen of single-gene knockout strains for the methyltransferases, we revealed that they form a network controlling methylome, transcriptome and adaptive phenotype sets. The methyltransferases interact with each other in a hierarchical way, sometimes regulated positively by one methyltransferase but negatively with another. Motility, oxidative stress tolerance and DNA damage repair are likewise regulated by multiple methyltransferases. Their regulation sometimes involves translation start and stop codons suggesting coupling of methylation, transcription and translation. The methyltransferases frequently change their sequence-specificity through gene conversion of their target recognition domain and switch their target sets to remodel the network. The emerging picture of a metamorphosing gene regulation network, or firework, consisting of epigenetic systems ever-changing their specificity in search for adaptation, provides a new paradigm in understanding global gene regulation and adaptive evolution.