Change of gene structure and function by non-homologous end-joining, homologous recombination, and transposition of DNA

Transcribed enhancer sequences are required for maize p1 paramutation

Paramutation is a transfer of heritable silencing states between interacting endogenous alleles or between endogenous alleles and homologous transgenes. Prior results demonstrated that paramutation occurs at the P1-rr (red pericarp and red cob) allele of the maize p1 (pericarp color 1) gene when exposed to a transgene containing a 1.2-kb enhancer fragment (P1.2) of P1-rr. The paramutable P1-rr allele undergoes transcriptional silencing resulting in a paramutant light-pigmented P1-rr' state. To define more precisely the sequences required to elicit paramutation, the P1.2 fragment was further subdivided, and the fragments transformed into maize plants and crossed with P1-rr. Analysis of the progeny plants showed that the sequences required for paramutation are located within a ∼600-bp segment of P1.2 and that this segment overlaps with a previously identified enhancer that is present in 4 direct repeats in P1-rr. The paramutagenic segment is transcribed in both the expressed P1-rr and the silenced P1-rr'. Transcription is sensitive to α-amanitin, indicating that RNA polymerase II mediates most of the transcription of this sequence. Although transcription within the paramutagenic sequence was similar in all tested genotypes, small RNAs were more abundant in the silenced P1-rr' epiallele relative to the expressed P1-rr allele. In agreement with prior results indicating the association of RNA-mediated DNA methylation in p1 paramutation, DNA blot analyses detected increased cytosine methylation of the paramutant P1-rr' sequences homologous to the transgenic P1.2 subfragments. Together these results demonstrate that the P1-rr enhancer repeats mediate p1 paramutation.

Complex chromosomal rearrangements induced by transposons in maize

Eukaryotic genomes are large and complex, and gene expression can be affected by multiple regulatory elements and their positions within the dynamic chromatin architecture. Transposable elements are known to play important roles in genome evolution, yet questions remain as to how transposable elements alter genome structure and affect gene expression. Previous studies have shown that genome rearrangements can be induced by Reversed Ends Transposition involving termini of Activator and related transposable elements in maize and other plants. Here, we show that complex alleles can be formed by the rapid and progressive accumulation of Activator-induced duplications and rearrangements. The p1 gene enhancer in maize can induce ectopic expression of the nearby p2 gene in pericarp tissue when placed near it via different structural rearrangements. By screening for p2 expression, we identified and studied 5 cases in which multiple sequential transposition events occurred and increased the p1 enhancer copy number. We see active p2 expression due to multiple copies of the p1 enhancer present near p2 in all 5 cases. The p1 enhancer effects are confirmed by the observation that loss of p2 expression is correlated with transposition-induced excision of the p1 enhancers. We also performed a targeted Chromosome Conformation Capture experiment to test the physical interaction between the p1 enhancer and p2 promoter region. Together, our results show that transposon-induced rearrangements can accumulate rapidly and progressively increase genetic variation important for genomic evolution.

Transposon-induced inversions activate gene expression in the maize pericarp

Chromosomal inversions can have considerable biological and agronomic impacts including disrupted gene function, change in gene expression, and inhibited recombination. Here, we describe the molecular structure and functional impact of six inversions caused by Alternative Transpositions between p1 and p2 genes responsible for floral pigmentation in maize. In maize line p1-wwB54, the p1 gene is null and the p2 gene is expressed in anther and silk but not in pericarp, making the kernels white. By screening for kernels with red pericarp, we identified inversions in this region caused by transposition of Ac and fractured Ac (fAc) transposable elements. We hypothesize that these inversions place the p2 gene promoter near a p1 gene enhancer, thereby activating p2 expression in kernel pericarp. To our knowledge, this is the first report of multiple recurrent inversions that change the position of a gene promoter relative to an enhancer to induce ectopic expression in a eukaryote.

Mutagenic repair of double-stranded DNA breaks in vaccinia virus genomes requires cellular DNA ligase IV activity in the cytosol

Poxviruses comprise a group of large dsDNA viruses that include members relevant to human and animal health, such as variola virus, monkeypox virus, cowpox virus and vaccinia virus (VACV). Poxviruses are remarkable for their unique replication cycle, which is restricted to the cytoplasm of infected cells. The independence from the host nucleus requires poxviruses to encode most of the enzymes involved in DNA replication, transcription and processing. Here, we use the CRISPR/Cas9 genome engineering system to induce DNA damage to VACV (strain Western Reserve) genomes. We show that targeting CRISPR/Cas9 to essential viral genes limits virus replication efficiently. Although VACV is a strictly cytoplasmic pathogen, we observed extensive viral genome editing at the target site; this is reminiscent of a non-homologous end-joining DNA repair mechanism. This pathway was not dependent on the viral DNA ligase, but critically involved the cellular DNA ligase IV. Our data show that DNA ligase IV can act outside of the nucleus to allow repair of dsDNA breaks in poxvirus genomes. This pathway might contribute to the introduction of mutations within the genome of poxviruses and may thereby promote the evolution of these viruses.

Analysis of tandem gene copies in maize chromosomal regions reconstructed from long sequence reads

Haplotype variation not only involves SNPs but also insertions and deletions, in particular gene copy number variations. However, comparisons of individual genomes have been difficult because traditional sequencing methods give too short reads to unambiguously reconstruct chromosomal regions containing repetitive DNA sequences. An example of such a case is the protein gene family in maize that acts as a sink for reduced nitrogen in the seed. Previously, 41-48 gene copies of the alpha zein gene family that spread over six loci spanning between 30- and 500-kb chromosomal regions have been described in two Iowa Stiff Stalk (SS) inbreds. Analyses of those regions were possible because of overlapping BAC clones, generated by an expensive and labor-intensive approach. Here we used single-molecule real-time (Pacific Biosciences) shotgun sequencing to assemble the six chromosomal regions from the Non-Stiff Stalk maize inbred W22 from a single DNA sequence dataset. To validate the reconstructed regions, we developed an optical map (BioNano genome map; BioNano Genomics) of W22 and found agreement between the two datasets. Using the sequences of full-length cDNAs from W22, we found that the error rate of PacBio sequencing seemed to be less than 0.1% after autocorrection and assembly. Expressed genes, some with premature stop codons, are interspersed with nonexpressed genes, giving rise to genotype-specific expression differences. Alignment of these regions with those from the previous analyzed regions of SS lines exhibits in part dramatic differences between these two heterotic groups.

Comparative genomic de-convolution of the cotton genome revealed a decaploid ancestor and widespread chromosomal fractionation

The 'apparently' simple genomes of many angiosperms mask complex evolutionary histories. The reference genome sequence for cotton (Gossypium spp.) revealed a ploidy change of a complexity unprecedented to date, indeed that could not be distinguished as to its exact dosage. Herein, by developing several comparative, computational and statistical approaches, we revealed a 5× multiplication in the cotton lineage of an ancestral genome common to cotton and cacao, and proposed evolutionary models to show how such a decaploid ancestor formed. The c. 70% gene loss necessary to bring the ancestral decaploid to its current gene count appears to fit an approximate geometrical model; that is, although many genes may be lost by single-gene deletion events, some may be lost in groups of consecutive genes. Gene loss following cotton decaploidy has largely just reduced gene copy numbers of some homologous groups. We designed a novel approach to deconvolute layers of chromosome homology, providing definitive information on gene orthology and paralogy across broad evolutionary distances, both of fundamental value and serving as an important platform to support further studies in and beyond cotton and genomics communities.

P1 Epigenetic Regulation in Leaves of High Altitude Maize Landraces: Effect of UV-B Radiation

P1 is a R2R3-MYB transcription factor that regulates the accumulation of a specific group of flavonoids in maize floral tissues, such as flavones and phlobaphenes. P1 is also highly expressed in leaves of maize landraces adapted to high altitudes and higher levels of UV-B radiation. In this work, we analyzed the epigenetic regulation of the P1 gene by UV-B in leaves of different maize landraces. Our results demonstrate that DNA methylation in the P1 proximal promoter, intron1 and intron2 is decreased by UV-B in all lines analyzed; however, the basal DNA methylation levels are lower in the landraces than in B73, a low altitude inbred line. DNA demethylation by UV-B is accompanied by a decrease in H3 methylation at Lys 9 and 27, and by an increase in H3 acetylation. smRNAs complementary to specific regions of the proximal promoter and of intron 2 3' end are also decreased by UV-B; interestingly, P1 smRNA levels are lower in the landraces than in B73 both under control conditions and after UV-B exposure, suggesting that smRNAs regulate P1 expression by UV-B in maize leaves. Finally, we investigated if different P1 targets in flower tissues are also regulated by this transcription factor in response to UV-B. Some targets analyzed show an induction in maize landraces in response to UV-B, with higher basal expression levels in the landraces than in B73; however, not all the transcripts analyzed were found to be regulated by UV-B in leaves.

European flint landraces grown in situ reveal adaptive introgression from modern maize

We have investigated the role of selection in the determination of the detected levels of introgression from modern maize hybrid varieties into maize landraces still cultivated in situ in Italy. We exploited the availability of a historical collection of landraces undertaken before the introduction and widespread use of modern maize, to analyse genomic changes that have occurred in these maize landraces over 50 years of co-existence with hybrid varieties. We have combined a previously published SSR dataset (n=21) with an AFLP loci dataset (n=168) to provide higher resolution power and to obtain a more detailed picture. We show that selection pressures for adaptation have favoured new alleles introduced by migration from hybrids. This shows the potential for analysis of historical introgression even over this short period of 50 years, for an understanding of the evolution of the genome and for the identification of its functionally important regions. Moreover, this demonstrates that landraces grown in situ represent almost unique populations for use for such studies when the focus is on the domesticated plant. This is due to their adaptation, which has arisen from their dynamic evolution under a continuously changing agro-ecological environment, and their capture of new alleles from hybridisation. We have also identified loci for which selection has inhibited introgression from modern germplasm and has enhanced the distinction between landraces and modern maize. These loci indicate that selection acted in the past, during the formation of the flint and dent gene pools. In particular, the locus showing the strongest signals of selection is a Misfit transposable element. Finally, molecular characterisation of the same samples with two different molecular markers has allowed us to compare their performances. Although the genetic-diversity and population-structure analyses provide the same global qualitative pattern, which thus provides the same inferences, there are differences related to their natures and characteristics.

DNA damage and repair in plants under ultraviolet and ionizing radiations

Being sessile, plants are continuously exposed to DNA-damaging agents present in the environment such as ultraviolet (UV) and ionizing radiations (IR). Sunlight acts as an energy source for photosynthetic plants; hence, avoidance of UV radiations (namely, UV-A, 315–400 nm; UV-B, 280–315 nm; and UV-C, <280 nm) is unpreventable. DNA in particular strongly absorbs UV-B; therefore, it is the most important target for UV-B induced damage. On the other hand, IR causes water radiolysis, which generates highly reactive hydroxyl radicals (OH^•) and causes radiogenic damage to important cellular components. However, to maintain genomic integrity under UV/IR exposure, plants make use of several DNA repair mechanisms. In the light of recent breakthrough, the current minireview (a) introduces UV/IR and overviews UV/IR-mediated DNA damage products and (b) critically discusses the biochemistry and genetics of major pathways responsible for the repair of UV/IR-accrued DNA damage. The outcome of the discussion may be helpful in devising future research in the current context.

Birth of three stowaway-like MITE families via microhomology-mediated miniaturization of a Tc1/Mariner element in the yellow fever mosquito

Eukaryotic genomes contain numerous DNA transposons that move by a cut-and-paste mechanism. The majority of these elements are self-insufficient and dependent on their autonomous relatives to transpose. Miniature inverted repeat transposable elements (MITEs) are often the most numerous nonautonomous DNA elements in a higher eukaryotic genome. Little is known about the origin of these MITE families as few of them are accompanied by their direct ancestral elements in a genome. Analyses of MITEs in the yellow fever mosquito identified its youngest MITE family, designated as Gnome, that contains at least 116 identical copies. Genome-wide search for direct ancestral autonomous elements of Gnome revealed an elusive single copy Tc1/Mariner-like element, named as Ozma, that encodes a transposase with a DD37E triad motif. Strikingly, Ozma also gave rise to two additional MITE families, designated as Elf and Goblin. These three MITE families were derived at different times during evolution and bear internal sequences originated from different regions of Ozma. Upon close inspection of the sequence junctions, the internal deletions during the formation of these three MITE families always occurred between two microhomologous sites (6–8 bp). These results suggest that multiple MITE families may originate from a single ancestral autonomous element, and formation of MITEs can be mediated by sequence microhomology. Ozma and its related MITEs are exceptional candidates for the long sought-after endogenous active transposon tool in genetic control of mosquitoes.

Paramutagenicity of a p1 epiallele in maize

Complex silencing mechanisms in plants and other kingdoms target transposons, repeat sequences, invasive viral nucleic acids and transgenes, but also endogenous genes and genes involved in paramutation. Paramutation occurs in a heterozygote when a transcriptionally active allele heritably adopts the epigenetic state of a transcriptionally and/or post-transcriptionally repressed allele. P1-rr and its silenced epiallele P1-pr, which encode a Myb-like transcription factor mediating pigmentation in floral organs of Zea mays, differ in their cytosine methylation pattern and chromatin structure at a complex enhancer site. Here, we tested whether P1-pr is able to heritably silence its transcriptionally active P1-rr allele in a heterozygote and whether DNA methylation is associated with the establishment and maintenance of P1-rr silencing. We found that P1-pr participates in paramutation as the repressing allele and P1-rr as the sensitive allele. Silencing of P1-rr is highly variable compared to the inducing P1-pr resulting in a wide range of gene expression. Whereas cytosine methylation at P1-rr is negatively correlated with transcription and pigment levels after segregation of P1-pr, methylation lags behind the establishment of the repressed p1 gene expression. We propose a model in which P1-pr paramutation is triggered by changing epigenetic states of transposons immediately adjacent to a P1-rr enhancer sequence. Considering the vast amount of transposable elements in the maize genome close to regulatory elements of genes, numerous loci could undergo paramutation-induced allele silencing, which could also have a significant impact on breeding agronomically important traits.

Epiallele biogenesis in maize

We have correlated cytosine methylation of two epialleles, P1-rr and P1-pr, with variation in gene expression and therefore phenotype. The p1 gene in maize encodes a transcription factor that controls phlobaphene pigment accumulation in floral tissues. While cytosine methylation was assayed in various regions spanning 17 kb, the only difference in DNA methylation pattern between the expressed P1-rr allele and the silenced P1-pr allele was detected in a region that consists of a complex arrangement of transposons and adjacent repeats. This region, which comprises the distal enhancer element of P1-rr, is hypermethylated in P1-pr compared to P1-rr. Based on other precedents, we hypothesize that DNA methylation spreads from the transposable elements into the flanking P1-rr enhancer, thereby transcriptionally silencing the gene. Interestingly, P1-pr is reactivated in mutants of the dominant epigenetic modifier Ufo1. DNA methylation in the distal enhancer sequence is significantly reduced, which inversely correlates with increased transcript levels and pigmentation in P1-pr Ufo1 plants. If in general DNA methylation spreads from transposons into adjacent sequences containing regulatory elements for neighboring genes, the corresponding genes could be silenced by chance. Given the large amount of transposable elements in the maize genome, epialleles may be far more frequent than previously estimated.

Analysis of the P1 promoter in response to UV-B radiation in allelic variants of high-altitude maize

BACKGROUND

Plants living at high altitudes are typically exposed to elevated UV-B radiation, and harbor mechanisms to prevent the induced damage, such as the accumulation of UV-absorbing compounds. The maize R2R3-MYB transcription factor P1 controls the accumulation of several UV-B absorbing phenolics by activating a subset of flavonoid biosynthetic genes in leaves of maize landraces adapted to high altitudes.

RESULTS

Here, we studied the UV-B regulation of P1 in maize leaves of high altitude landraces, and we investigated how UV-B regulates P1 binding to the CHS promoter in both low and high altitude lines. In addition, we analyzed whether the expansion in the P1 expression domain between these maize landraces and inbred lines is associated to changes in the molecular structure of the proximal promoter, distal enhancer and first intron of P1. Finally, using transient expression experiments in protoplasts from various maize genotypes, we investigated whether the different expression patterns of P1 in the high altitude landraces could be attributed to trans- or cis-acting elements.

CONCLUSIONS

Together, our results demonstrate that, although differences in cis-acting elements exist between the different lines under study, the different patterns of P1 expression are largely a consequence of effects in trans.

Divergence of gene regulation through chromosomal rearrangements

BACKGROUND

The molecular mechanisms that modify genome structures to give birth and death to alleles are still not well understood. To investigate the causative chromosomal rearrangements, we took advantage of the allelic diversity of the duplicated p1 and p2 genes in maize. Both genes encode a transcription factor involved in maysin synthesis, which confers resistance to corn earworm. However, p1 also controls accumulation of reddish pigments in floral tissues and has therefore acquired a new function after gene duplication. p1 alleles vary in their tissue-specific expression, which is indicated in their allele designation: the first suffix refers to red or white pericarp pigmentation and the second to red or white glume pigmentation.

RESULTS

Comparing chromosomal regions comprising p1-ww[4Co63], P1-rw1077 and P1-rr4B2 alleles with that of the reference genome, P1-wr[B73], enabled us to reconstruct additive events of transposition, chromosome breaks and repairs, and recombination that resulted in phenotypic variation and chimeric regulatory signals. The p1-ww[4Co63] null allele is probably derived from P1-wr[B73] by unequal crossover between large flanking sequences. A transposon insertion in a P1-wr-like allele and NHEJ (non-homologous end-joining) could have resulted in the formation of the P1-rw1077 allele. A second NHEJ event, followed by unequal crossover, probably led to the duplication of an enhancer region, creating the P1-rr4B2 allele. Moreover, a rather dynamic picture emerged in the use of polyadenylation signals by different p1 alleles. Interestingly, p1 alleles can be placed on both sides of a large retrotransposon cluster through recombination, while functional p2 alleles have only been found proximal to the cluster.

CONCLUSIONS

Allelic diversity of the p locus exemplifies how gene duplications promote phenotypic variability through composite regulatory signals. Transposition events increase the level of genomic complexity based not only on insertions but also on excisions that cause DNA double-strand breaks and trigger illegitimate recombination.

Ancestral grass karyotype reconstruction unravels new mechanisms of genome shuffling as a source of plant evolution

The comparison of the chromosome numbers of today's species with common reconstructed paleo-ancestors has led to intense speculation of how chromosomes have been rearranged over time in mammals. However, similar studies in plants with respect to genome evolution as well as molecular mechanisms leading to mosaic synteny blocks have been lacking due to relevant examples of evolutionary zooms from genomic sequences. Such studies require genomes of species that belong to the same family but are diverged to fall into different subfamilies. Our most important crops belong to the family of the grasses, where a number of genomes have now been sequenced. Based on detailed paleogenomics, using inference from n = 5-12 grass ancestral karyotypes (AGKs) in terms of gene content and order, we delineated sequence intervals comprising a complete set of junction break points of orthologous regions from rice, maize, sorghum, and Brachypodium genomes, representing three different subfamilies and different polyploidization events. By focusing on these sequence intervals, we could show that the chromosome number variation/reduction from the n = 12 common paleo-ancestor was driven by nonrandom centric double-strand break repair events. It appeared that the centromeric/telomeric illegitimate recombination between nonhomologous chromosomes led to nested chromosome fusions (NCFs) and synteny break points (SBPs). When intervals comprising NCFs were compared in their structure, we concluded that SBPs (1) were meiotic recombination hotspots, (2) corresponded to high sequence turnover loci through repeat invasion, and (3) might be considered as hotspots of evolutionary novelty that could act as a reservoir for producing adaptive phenotypes.

Pervasive gene content variation and copy number variation in maize and its undomesticated progenitor

Individuals of the same species are generally thought to have very similar genomes. However, there is growing evidence that structural variation in the form of copy number variation (CNV) and presence-absence variation (PAV) can lead to variation in the genome content of individuals within a species. Array comparative genomic hybridization (CGH) was used to compare gene content and copy number variation among 19 diverse maize inbreds and 14 genotypes of the wild ancestor of maize, teosinte. We identified 479 genes exhibiting higher copy number in some genotypes (UpCNV) and 3410 genes that have either fewer copies or are missing in the genome of at least one genotype relative to B73 (DownCNV/PAV). Many of these DownCNV/PAV are examples of genes present in B73, but missing from other genotypes. Over 70% of the CNV/PAV examples are identified in multiple genotypes, and the majority of events are observed in both maize and teosinte, suggesting that these variants predate domestication and that there is not strong selection acting against them. Many of the genes affected by CNV/PAV are either maize specific (thus possible annotation artifacts) or members of large gene families, suggesting that the gene loss can be tolerated through buffering by redundant functions encoded elsewhere in the genome. While this structural variation may not result in major qualitative variation due to genetic buffering, it may significantly contribute to quantitative variation.

Tissue culture-induced novel epialleles of a Myb transcription factor encoded by pericarp color1 in maize

Plants regenerated from tissue culture often display somaclonal variation, that is, somatic and often meiotically heritable phenotypic variation that can result from both genetic and epigenetic modifications. To better understand the molecular basis of somaclonal variation, we have characterized four unique tissue culture-derived epialleles of the pericarp color1 (p1) gene of maize (Zea mays L.). The progenitor p1 allele, P1-wr, is composed of multiple head-to-tail tandemly arranged copies of the complete gene unit and specifies brick-red phlobaphene pigmentation in the cob glumes. The novel epialleles identified in progeny plants regenerated from tissue culture showed partial to complete loss of p1 function indicated by pink or colorless cob glumes. Loss of pigmentation was correlated with nearly complete loss of p1 steady-state transcripts. DNA gel-blot analysis and genomic bisulfite sequencing showed that silencing of the epialleles was associated with hypermethylation of a region in the second intron of P1-wr. Presence of Unstable factor for orange1 (Ufo1), an unlinked epigenetic modifier of p1, restored the cob glume pigmentation in the silenced alleles, and such reactivation was accompanied by hypomethylation of the p1 sequence. This observation confirmed that silencing of the epialleles is indeed due to epigenetic modifications and that the p1 epialleles were capable of functioning in the presence of the correct trans-acting factors. While the low-copy regions of the genome generally undergo hypomethylation during tissue culture, our study shows that the tandemly repeated genes are also prone to hypermethylation and epigenetic silencing.

Gene structure induced epigenetic modifications of pericarp color1 alleles of maize result in tissue-specific mosaicism

Background

The pericarp color1 (p1) gene encodes for a myb-homologous protein that regulates the biosynthesis of brick-red flavonoid pigments called phlobahpenes. The pattern of pigmentation on the pericarp and cob glumes depends upon the allelic constitution at the p1 locus. p1 alleles have unique gene structure and copy number which have been proposed to influence the epigenetic regulation of tissue-specific gene expression. For example, the presence of tandem-repeats has been correlated with the suppression of pericarp pigmentation though a mechanism associated with increased DNA methylation.

Methodology/Principal Findings

Herein, we extensively characterize a p1 allele called P1-mosaic (P1-mm) that has mosaic pericarp and light pink or colorless cob glumes pigmentation. Relative to the P1-wr (white pericarp and red cob glumes), we show that the tandem repeats of P1-mm have a modified gene structure containing a reduced number of repeats. The P1-mm has reduced DNA methylation at a distal enhancer and elevated DNA methylation downstream of the transcription start site.

Conclusions/Significance

Mosaic gene expression occurs in many eukaryotes. Herein we use maize p1 gene as model system to provide further insight about the mechanisms that govern expression mosaicism. We suggest that the gene structure of P1-mm is modified in some of its tandem gene repeats. It is known that repeated genes are susceptible to chromatin-mediated regulation of gene expression. We discuss how the modification to the tandem repeats of P1-mm may have disrupted the epigenetic mechanisms that stably confer tissue-specific expression.