Recombination products suggest the frequent occurrence of aberrant gene replacement in the moss Physcomitrella patens

In gene replacement, a variant of gene targeting, transformed DNA integrates into the genome by homologous recombination (HR) to replace resident sequences. Gene replacement in the moss Physcomitrella patens is extremely efficient, but often large amounts of additional DNA are integrated at the target locus. A detailed analysis of recombination junctions of PpCOL2 gene knockout mutants shows that the integrated DNA can be highly rearranged. Our data suggest that the replaced sequences were excised by HR and became integrated back into the genome by non‐homologous end‐joining (NHEJ). RAD51‐mediated strand‐invasion and subsequent strand‐exchange is central to the two‐end invasion pathway, the major gene replacement pathway in yeast. In this pathway, integration is initiated by the free ends of a single replacement vector‐derived donor molecule which then integrates as an entity. Gene replacement in P. patens is entirely RAD51‐dependent suggesting the existence of a pathway mechanistically similar to two‐end invasion. However, invasion of the two ends does not seem to be stringently coordinated in P. patens. Actually, often only one fragment end became integrated by HR, or one‐sided integration of two independent donor fragments occurred simultaneously leading to a double‐strand break that is subsequently sealed by NHEJ and thus causes the observed rearrangements.     

RECA plays a dual role in the maintenance of chloroplast genome stability in Physcomitrella patens

Chloroplast DNA (cpDNA) encodes essential genes for chloroplast functions, including photosynthesis. Homologous recombination occurs frequently in cpDNA; however, its significance and underlying mechanism remain poorly understood. In this study, we analyzed the role of a nuclear‐encoded chloroplast‐localized homolog of RecA recombinase, which is a key factor in homologous recombination in bacteria, in the moss Physcomitrella patens. Complete knockout (KO) of the P. patens chloroplast RecA homolog RECA2 caused a modest growth defect and conferred sensitivity to methyl methanesulfonate and UV. The KO mutant exhibited low recovery of cpDNA from methyl methanesulfonate damage, suggesting that RECA2 knockout impairs repair of damaged cpDNA. The RECA2 KO mutant also exhibited reduced cpDNA copy number and an elevated level of cpDNA molecule resulting from aberrant recombination between short dispersed repeats (13–63 bp), indicating that the RECA2 KO chloroplast genome was destabilized. Taken together, these data suggest a dual role for RECA2 in the maintenance of chloroplast genome stability: RECA2 suppresses aberrant recombination between short dispersed repeats and promotes repair of damaged DNA.     

Rapid proliferation and nucleolar organizer targeting centromeric retrotransposons in cotton

Centromeric chromatin in most eukaryotes is composed of highly repetitive centromeric retrotransposons and satellite repeats that are highly variable even among closely related species. The evolutionary mechanisms that underlie the rapid evolution of centromeric repeats remain unknown. To obtain insight into the evolution of centromeric repeats following polyploidy, we studied a model diploid progenitor (Gossypium raimondii, D‐genome) of the allopolyploid (AD‐genome) cottons, G. hirsutum and G. barbadense. Sequence analysis of chromatin‐immunoprecipitated DNA showed that the G. raimondii centromeric repeats originated from retrotransposon‐related sequences. Comparative analysis showed that nine of the 10 analyzed centromeric repeats were absent from the centromeres in the A‐genome and related diploid species (B‐, F‐ and G‐genomes), indicating that they colonized the centromeres of D‐genome lineage after the divergence of the A‐ and D‐ ancestral species or that they were ancestrally retained prior to the origin of Gossypium. Notably, six of the nine repeats were present in both the A‐ and D‐subgenomes in tetraploid G. hirsutum, and increased in abundance in both subgenomes. This finding suggests that centromeric repeats may spread and proliferate between genomes subsequent to polyploidization. Two repeats, Gr334 and Gr359 occurred in both the centromeres and nucleolar organizer regions (NORs) in D‐ and AD‐genome species, yet localized to just the NORs in A‐, B‐, F‐, and G‐genome species. Contained within is a story of an established centromeric repeat that is eliminated and allopolyploidization provides an opportunity for reinvasion and reestablishment, which broadens our evolutionary understanding behind the cycles of centromeric repeat establishment and targeting.     

MSH1 maintains organelle genome stability and genetically interacts with RECA and RECG in the moss Physcomitrella patens

Chloroplast and mitochondrial DNA encodes genes that are essential for photosynthesis and respiration, respectively. Thus, loss of integrity of the genomic DNA of organelles leads to a decline in organelle function and alteration of organelle genetic information. RECA (RECA1 and RECA2) and RECG, which are homologs of bacterial homologous recombination repair (HRR) factors RecA and RecG, respectively, play an important role in the maintenance of integrity of the organelle genome by suppressing aberrant recombination between short dispersed repeats (SDRs) in the moss Physcomitrella patens. On the other hand, MutS homolog 1 (MSH1), a plant‐specific MSH with a C‐terminal GIY‐YIG endonuclease domain, is involved in the maintenance of integrity of the organelle genome in the angiosperm Arabidopsis thaliana. Here, we address the role of the duplicated MSH1 genes, MSH1A and MSH1B, in P. patens, in which MSH1A lacks the C‐terminal endonuclease domain. MSH1A and MSH1B localized to both chloroplast and mitochondrial nucleoids in protoplast cells. Single and double knockout (KO) mutants of MSH1A and MSH1B showed no obvious morphological defects; however, MSH1B KO and double KO mutants, as well as MSH1B GIY‐YIG deletion mutants, exhibited genomic instability due to recombination between SDRs in chloroplasts and mitochondria. Creating double KO mutations of each combination of MSH1B, RECA2 and RECG synergistically increased recombination between SDRs in chloroplasts and mitochondria. These results show the role of MSH1 in the maintenance of integrity of the organelle genome and the genetic interaction between MSH1 and homologs of HRR factors in the basal land plant P. patens.     

The bryophytes Physcomitrium patens and Marchantia polymorpha as model systems for studying evolutionary cell and developmental biology in plants

Bryophytes are nonvascular spore-forming plants. Unlike in flowering plants, the gametophyte (haploid) generation of bryophytes dominates the sporophyte (diploid) generation. A comparison of bryophytes with flowering plants allows us to answer some fundamental questions raised in evolutionary cell and developmental biology. The moss Physcomitrium patens was the first bryophyte with a sequenced genome. Many cell and developmental studies have been conducted in this species using gene targeting by homologous recombination. The liverwort Marchantia polymorpha has recently emerged as an excellent model system with low genomic redundancy in most of its regulatory pathways. With the development of molecular genetic tools such as efficient genome editing, both P. patens and M. polymorpha have provided many valuable insights. Here, we review these advances with a special focus on polarity formation at the cell and tissue levels. We examine current knowledge regarding the cellular mechanisms of polarized cell elongation and cell division, including symmetric and asymmetric cell division. We also examine the role of polar auxin transport in mosses and liverworts. Finally, we discuss the future of evolutionary cell and developmental biological studies in plants.

A review of the cell biological and developmental mechanisms of bryophytes, including Physcomitrium patens and Marchantia polymorpha.     

A lineage‐specific centromere retrotransposon in Oryza brachyantha

Most eukaryotic centromeres contain large quantities of repetitive DNA, such as satellite repeats and retrotransposons. Unlike most transposons in plant genomes, the centromeric retrotransposon (CR) family is conserved over long evolutionary periods among a majority of the grass species. CR elements are highly concentrated in centromeres, and are likely to play a role in centromere function. In order to study centromere evolution in the Oryza (rice) genus, we sequenced the orthologous region to centromere 8 of Oryza sativa from a related species, Oryza brachyantha. We found that O. brachyantha does not have the canonical CRR (CR of rice) found in the centromeres of all other Oryza species. Instead, a new Ty3‐gypsy (Metaviridae) retroelement (FRetro3) was found to colonize the centromeres of this species. This retroelement is found in high copy numbers in the O. brachyantha genome, but not in other Oryza genomes, and based on the dating of long terminal repeats (LTRs) of FRetro3 it was amplified in the genome in the last few million years. Interestingly, there is a high level of removal of FRetro3 based on solo‐LTRs to full‐length elements, and this rapid turnover may have played a role in the replacement of the canonical CRR with the new element by active deletion. Comparison with previously described ChIP cloning data revealed that FRetro3 is found in CENH3‐associated chromatin sequences. Thus, within a single lineage of the Oryza genus, the canonical component of grass centromeres has been replaced with a new retrotransposon that has all the hallmarks of a centromeric retroelement.     

Shifts in the evolutionary rate and intensity of purifying selection between two Brassica genomes revealed by analyses of orthologous transposons and relics of a whole genome triplication

Recent sequencing of the Brassica rapa and Brassica oleracea genomes revealed extremely contrasting genomic features such as the abundance and distribution of transposable elements between the two genomes. However, whether and how these structural differentiations may have influenced the evolutionary rates of the two genomes since their split from a common ancestor are unknown. Here, we investigated and compared the rates of nucleotide substitution between two long terminal repeats (LTRs) of individual orthologous LTR‐retrotransposons, the rates of synonymous and non‐synonymous substitution among triplicated genes retained in both genomes from a shared whole genome triplication event, and the rates of genetic recombination estimated/deduced by the comparison of physical and genetic distances along chromosomes and ratios of solo LTRs to intact elements. Overall, LTR sequences and genic sequences showed more rapid nucleotide substitution in B. rapa than in B. oleracea. Synonymous substitution of triplicated genes retained from a shared whole genome triplication was detected at higher rates in B. rapa than in B. oleracea. Interestingly, non‐synonymous substitution was observed at lower rates in the former than in the latter, indicating shifted densities of purifying selection between the two genomes. In addition to evolutionary asymmetry, orthologous genes differentially regulated and/or disrupted by transposable elements between the two genomes were also characterized. Our analyses suggest that local genomic and epigenomic features, such as recombination rates and chromatin dynamics reshaped by independent proliferation of transposable elements and elimination between the two genomes, are perhaps partially the causes and partially the outcomes of the observed inter‐specific asymmetric evolution.     

DNA methylation of retrotransposons,DNA transposons and genes in sugar beet (Beta vulgaris L.)

The methylation of cytosines shapes the epigenetic landscape of plant genomes, coordinates transgenerational epigenetic inheritance, represses the activity of transposable elements (TEs), affects gene expression and, hence, can influence the phenotype. Sugar beet (Beta vulgaris ssp. vulgaris), an important crop that accounts for 30% of worldwide sugar needs, has a relatively small genome size (758 Mbp) consisting of approximately 485 Mbp repetitive DNA (64%), in particular satellite DNA, retrotransposons and DNA transposons. Genome‐wide cytosine methylation in the sugar beet genome was studied in leaves and leaf‐derived callus with a focus on repetitive sequences, including retrotransposons and DNA transposons, the major groups of repetitive DNA sequences, and compared with gene methylation. Genes showed a specific methylation pattern for CG, CHG (H = A, C, and T) and CHH sites, whereas the TE pattern differed, depending on the TE class (class 1, retrotransposons and class 2, DNA transposons). Along genes and TEs, CG and CHG methylation was higher than that of adjacent genomic regions. In contrast to the relatively low CHH methylation in retrotransposons and genes, the level of CHH methylation in DNA transposons was strongly increased, pointing to a functional role of asymmetric methylation in DNA transposon silencing. Comparison of genome‐wide DNA methylation between sugar beet leaves and callus revealed a differential methylation upon tissue culture. Potential epialleles were hypomethylated (lower methylation) at CG and CHG sites in retrotransposons and genes and hypermethylated (higher methylation) at CHH sites in DNA transposons of callus when compared with leaves.     

Analysis of transposable elements and organellar DNA in male and female genomes of a species with a huge Y chromosome reveals distinct Y centromeres

Few angiosperms have distinct Y chromosomes. Among those that do are Silene latifolia (Caryophyllaceae), Rumex acetosa (Polygonaceae) and Coccinia grandis (Cucurbitaceae), the latter having a male/female difference of 10% of the total genome (female individuals have a 0.85 pg genome, male individuals 0.94 pg), due to a Y chromosome that arose about 3 million years ago. We compared the sequence composition of male and female C. grandis plants and determined the chromosomal distribution of repetitive and organellar DNA with probes developed from 21 types of repetitive DNA, including 16 mobile elements. The size of the Y chromosome is largely due to the accumulation of certain repeats, such as members of the Ty1/copia and Ty3/gypsy superfamilies, an unclassified element and a satellite, but also plastome‐ and chondriome‐derived sequences. An abundant tandem repeat with a unit size of 144 bp stains the centromeres of the X chromosome and the autosomes, but is absent from the Y centromere. Immunostaining with pericentromere‐specific markers for anti‐histone H3Ser10ph and H2AThr120ph revealed a Y‐specific extension of these histone marks. That the Y centromere has a different make‐up from all the remaining centromeres raises questions about its spindle attachment, and suggests that centromeric or pericentromeric chromatin might be involved in the suppression of recombination.     

Construction of Agropyron Gaertn. genetic linkage maps using a wheat 660K SNP array reveals a homoeologous relationship with the wheat genome

Agropyron Gaertn. (P genome) is a wild relative of wheat that harbours many genetic variations that could be used to increase the genetic diversity of wheat. To agronomically transfer important genes from the P genome to a wheat chromosome by induced homoeologous pairing and recombination, it is necessary to determine the chromosomal relationships between Agropyron and wheat. Here, we report using the wheat 660K single nucleotide polymorphism (SNP) array to genotype a segregating Agropyron F₁ population derived from an interspecific cross between two cross‐pollinated diploid collections ‘Z1842’ [A. cristatum (L.) Beauv.] (male parent) and ‘Z2098’ [A. mongolicum Keng] (female parent) and 35 wheat–A. cristatum addition/substitution lines. Genetic linkage maps were constructed using 913 SNP markers distributed among seven linkage groups spanning 839.7 cM. The average distance between adjacent markers was 1.8 cM. The maps identified the homoeologous relationship between the P genome and wheat and revealed that the P and wheat genomes are collinear and relatively conserved. In addition, obvious rearrangements and introgression spread were observed throughout the P genome compared with the wheat genome. Combined with genotyping data, the complete set of wheat–A. cristatum addition/substitution lines was characterized according to their homoeologous relationships. In this study, the homoeologous relationship between the P genome and wheat was identified using genetic linkage maps, and the detection mean for wheat–A. cristatum introgressions might significantly accelerate the introgression of genetic variation from Agropyron into wheat for exploitation in wheat improvement programmes.     

Effects of sequence diversity and recombination on the accuracy of phylogenetic trees estimated by kSNP

kSNP v2 is a powerful tool for single nucleotide polymorphism (SNP) identification from complete microbial genomes and for estimating phylogenetic trees from the identified SNPs. kSNP can analyse finished genomes, genome assemblies, raw reads or any combination of those and does not require either genome alignment or reference genomes. This study uses sequence evolution simulations to evaluate the topological accuracy of kSNP trees and to assess the effects of diversity and recombination on that accuracy. The accuracies of kSNP trees are strongly affected by increasing diversity, with parsimony accuracy > maximum‐likelihood accuracy > neighbour‐joining accuracy. Accuracy is also strongly influenced by recombination; as recombination increases accuracy decreases. Reliable trees are arbitrarily defined as those that have ≥ 90% topological accuracy. It is determined that the best predictor of topological accuracy is the ratio of r/m, a measure of the effect of recombination, to FCK (the fraction of core kmers), a measure of diversity. Tools are available to allow investigators to determine both r/m and FCK, and the relationship between topological accuracy and the ratio of r/m to FCK is determined. The practical implication of this study is that kSNP is an effective tool for estimating phylogenetic trees from microbial genome sequences provided that both recombination and sequence diversity are within acceptable ranges.     

Biological implications of the occurrence of 32 members of the XTH (xyloglucan endotransglucosylase/hydrolase) family of proteins in the bryophyte Physcomitrella patens

This comprehensive overview of the xyloglucan endotransglucosylase/hydrolase (XTH) family of genes and proteins in bryophytes, based on research using genomic resources that are newly available for the moss Physcomitrella patens, provides new insights into plant evolution. In angiosperms, the XTH genes are found in large multi‐gene families, probably reflecting the diverse roles of individual XTHs in various cell types. As there are fewer cell types in P. patens than in angiosperms such as Arabidopsis and rice, it is tempting to deduce that there are fewer XTH family genes in bryophytes. However, the present study unexpectedly identified as many as 32 genes that potentially encode XTH family proteins in the genome of P. patens, constituting a fairly large multi‐gene family that is comparable in size with those of Arabidopsis and rice. In situ localization of xyloglucan endotransglucosylase activity in this moss indicates that some P. patens XTH proteins exhibit biochemical functions similar to those found in angiosperms, and that their expression profiles are tissue‐dependent. However, comparison of structural features of families of XTH genes between P. patens and angiosperms demonstrated the existence of several bryophyte‐specific XTH genes with distinct structural and functional features that are not found in angiosperms. These bryophyte‐specific XTH genes might have evolved to meet morphological and functional needs specific to the bryophyte. These findings raise interesting questions about the biological implications of the XTH family of proteins in non‐seed plants.     

A CRISPR/LbCas12a‐based method for highly efficient multiplex gene editing in Physcomitrella patens

Due to their high efficiency, specificity, and flexibility, programmable nucleases, such as those of the clustered regularly interspaced short palindromic repeats (CRISPR)/Cas12a (Cpf1) system, have greatly expanded the applicability of editing the genomes of various organisms. Genes from different gene families or genes with redundant functions in the same gene family can be examined by assembling multiple CRISPR RNAs (crRNAs) in a single vector. However, the activity and efficiency of CRISPR/Cas12a in the non‐vascular plant Physcomitrella patens are largely unknown. Here, we demonstrate that LbCas12a together with its mature crRNA can target multiple loci simultaneously in P. patens with high efficiency via co‐delivery of LbCas12a and a crRNA expression cassette in vivo. The mutation frequencies induced by CRISPR/LbCas12a at a single locus ranged from 26.5 to 100%, with diverse deletions being the most common type of mutation. Our method expands the repertoire of genome editing tools available for P. patens and facilitates the creation of loss‐of‐function mutants of multiple genes from different gene families.     

When a clonal genome finds its way back to a sexual species: evidence from ongoing but rare introgression in the hybridogenetic water frog complex

Besides several exceptions, asexual metazoans are usually viewed as ephemeral sinks for genomes, which become ‘frozen’ in clonal lineages after their emergence from ancestral sexual species. Here, we investigated whether and at what rate the asexuals are able to introgress their genomes back into the parental sexual population, thus more or less importantly affecting the gene pools of sexual species. We focused on hybridogenetic hybrids of western Palaearctic water frogs (Pelophylax esculentus), which originate through hybridization between P. ridibundus and P. lessonae, but transmit only clonal ridibundus genome into their gametes. Although usually mating with P. lessonae, P. esculentus may upon mating with P. ridibundus or another hybrid produce sexually reproducing P. ridibundus offspring with the introgressed ex‐clonal genome. We compared the rate of nuclear amplified fragment length polymorphism (AFLP) and mitochondrial introgression in two types of populations, that is, those where P. ridibundus occurs in isolation and those where it lives with the hybridogens. Although significant differentiation (Φpt) between sexual and clonal ridibundus genomes suggested limited gene flow between sexuals and hybridogens, a non‐negligible (~5%) proportion of P. ridibundus bore introgressed mtDNA and AFLP markers. Whereas transfer of mtDNA was exclusively unidirectional, introgression of nuclear markers was bidirectional. The proportion of introgressed P. ridibundus was highest in syntopic populations with P. esculentus, proving an ongoing and site‐specific interspecific genetic transfer mediated by hybridogenetic hybrids. It turns out that asexual hybrids are not just a sink for genes of sexual species, but may significantly influence the genetic architecture of their sexual counterparts.     

Centromeric DNA characterization in the model grass Brachypodium distachyon provides insights on the evolution of the genus

Brachypodium distachyon is a well‐established model monocot plant, and its small and compact genome has been used as an accurate reference for the much larger and often polyploid genomes of cereals such as Avena sativa (oats), Hordeum vulgare (barley) and Triticum aestivum (wheat). Centromeres are indispensable functional units of chromosomes and they play a core role in genome polyploidization events during evolution. As the Brachypodium genus contains about 20 species that differ significantly in terms of their basic chromosome numbers, genome size, ploidy levels and life strategies, studying their centromeres may provide important insight into the structure and evolution of the genome in this interesting and important genus. In this study, we isolated the centromeric DNA of the B. distachyon reference line Bd21 and characterized its composition via the chromatin immunoprecipitation of the nucleosomes that contain the centromere‐specific histone CENH3. We revealed that the centromeres of Bd21 have the features of typical multicellular eukaryotic centromeres. Strikingly, these centromeres contain relatively few centromeric satellite DNAs; in particular, the centromere of chromosome 5 (Bd5) consists of only ~40 kb. Moreover, the centromeric retrotransposons in B. distachyon (CRBds) are evolutionarily young. These transposable elements are located both within and adjacent to the CENH3 binding domains, and have similar compositions. Moreover, based on the presence of CRBds in the centromeres, the species in this study can be grouped into two distinct lineages. This may provide new evidence regarding the phylogenetic relationships within the Brachypodium genus.