Search

OBJECTIVES: Mint oil is used in various commercial applications world-wide. Mint oil is typically harvested from commercial clones of peppermint or spearmints. Spearmints are the product of a cross between two diploid species: Mentha longifolia (horse mint) and Mentha suaveolens (apple mint). Peppermints are the product of an additional hybridization step between spearmint and an octoploid Mentha aquatica (water mint). Here, we present a chromosome-scale assembly of the genome of a clone of M. suaveolens. Together with the previously assembled genome of M. longifolia, these assemblies are instrumental in addressing questions regarding the origins of spearmint and peppermint oil and the genomic composition of commercial spearmints, and to start elaborating strategies for mint cultivar improvement. DATA DESCRIPTION: A Falcon assembly of the genome of M. suaveolens was generated from 103X coverage of PacBio long reads. Additional scaffolding was conducted by Dovetail Genomics, using a Chicago library, and a HiC library. The resulting assembly had an N50 of 44.7 Mb, and 98.45% of the 536 Mb of the assembly were contained within 12 large superscaffolds. Finally, a genetic map was applied to correct persistent misjoins. Illumina RNA-Seq libraries from a variety of tissues were used to annotate the genome.

The genetic control of many plant traits can be highly complex. Both allelic variation (sequence change) and dosage variation (copy number change) contribute to a plants phenotype. While numerous studies have investigated the effect of allelic or dosage variation, very few have documented both within the same system, leaving their relative contribution to phenotypic effects unclear. The Populus genome is highly polymorphic, and poplars are fairly tolerant of gene dosage variation. Here, using a previously established Populus hybrid F1 population, we assessed and compared the effect of natural allelic variation and induced dosage variation on biomass, phenology, and leaf morphology traits. We identified QTLs for many of these traits, but our results indicate limited overlap between the QTLs associated with natural allelic variation and induced dosage variation. Additionally, the integration of data from both allelic and dosage variation identifies a larger set of QTLs that together explain a larger percentage of the phenotypic variance. Finally, our results suggest that the effect of the large indels might mask that of allelic QTLs. Our study helps clarify the relationship between allelic and dosage variation and their effects on quantitative traits.

High-throughput sequencing-based methods for bulked segregant analysis (BSA) allow for the rapid identification of genetic markers associated with traits of interest. BSA studies have successfully identified qualitative (binary) and quantitative trait loci (QTLs) using QTL mapping. However, most require population structures that fit the models available and a reference genome. Instead, high-throughput short-read sequencing can be combined with BSA of k-mers (BSA-k-mer) to map traits that appear refractory to standard approaches. This method can be applied to any organism and is particularly useful for species with genomes diverged from the closest sequenced genome. It is also instrumental when dealing with highly heterozygous and potentially polyploid genomes without phased haplotype assemblies and for which a single haplotype can control a trait. Finally, it is flexible in terms of population structure. Here, we apply the BSA-k-mer method for the rapid identification of candidate regions related to seed spot and seed size in diploid potato. Using a mixture of F1 and F2 individuals from a cross between 2 highly heterozygous parents, candidate sequences were identified for each trait using the BSA-k-mer approach. Using parental reads, we were able to determine the parental origin of the loci. Finally, we mapped the identified k-mers to a closely related potato genome to validate the method and determine the genomic loci underlying these sequences. The location identified for the seed spot matches with previously identified loci associated with pigmentation in potato. The loci associated with seed size are novel. Both loci are relevant in future breeding toward true seeds in potato.

INTRODUCTION: Owing to advances in high-throughput genome sequencing, QTL-Seq mapping of salt tolerance traits is a major platform for identifying soil-salinity tolerance QTLs to accelerate marker-assisted selection for salt-tolerant rice varieties. We performed QTL-BSA-Seq in the seedling stage of rice from a genetic cross of the extreme salt-sensitive variety, IR29, and Jao Khao (JK), a Thai salt-tolerant variety. METHODS: A total of 462 F2 progeny grown in soil and treated with 160 mM NaCl were used as the QTL mapping population. Two high- and low-bulk sets, based on cell membrane stability (CMS) and tiller number at the recovery stage (TN), were equally sampled. The genomes of each pool were sequenced, and statistical significance of QTL was calculated using QTLseq and G prime (G) analysis, which is based on calculating the allele frequency differences or Δ(SNP index). RESULTS: Both methods detected the overlapping interval region, wherein CMS-bulk was mapped at two loci in the 38.41-38.85 Mb region with 336 SNPs on chromosome 1 (qCMS1) and the 26.13-26.80 Mb region with 1,011 SNPs on chromosome 3 (qCMS3); the Δ(SNP index) peaks were -0.2709 and 0.3127, respectively. TN-bulk was mapped at only one locus in the overlapping 38.26-38.95 Mb region on chromosome 1 with 575 SNPs (qTN1) and a Δ(SNP index) peak of -0.3544. These identified QTLs in two different genetic backgrounds of segregating populations derived from JK were validated. The results confirmed the colocalization of the qCMS1 and qTN1 traits on chromosome 1. Based on the CMS trait, qCMS1/qTN1 stably expressed 6%-18% of the phenotypic variance in the two validation populations, while qCMS1/qTN1 accounted for 16%-20% of the phenotypic variance in one validation population based on the TN trait. CONCLUSION: The findings confirm that the CMS and TN traits are tightly linked to the long arm of chromosome 1 rather than to chromosome 3. The validated qCMS-TN1 QTL can be used for gene/QTL pyramiding in marker-assisted selection to expedite breeding for salt resistance in rice at the seedling stage.

Minichromosomes are small, sometimes circular, rearranged chromosomes consisting of one centromere and short chromosomal arms formed by treatments that break DNA, including plant transformation. Minichromosomes have the potential to serve as vectors to quickly move valuable genes across a wide range of germplasm, including into adapted crop varieties. To realize this potential, minichromosomes must be reliably generated, easily manipulated, and stably inherited. Here we show a reliable method for minichromosome formation in haploids resulting from CENH3-mediated genome elimination, a process that generates genome instability and karyotypic novelty specifically on one parental genome. First, we identified 2 out of 260 haploids, each containing a single-copy minichromosome originating from centromeric regions of chromosomes 1 and 3, respectively. The chromosome 1 minichromosome we characterized did not pair at meiosis but displayed consistent transmission over nine selfing generations. Next, we demonstrated that CENH3-based haploid induction can produce minichromosomes in a targeted manner. Haploid inducers carrying a selectable pericentromeric marker were used to isolate additional chromosome-specific minichromosomes, which occurred in 3 out of 163 haploids. Our findings document the formation of heritable, rearranged chromosomes, and we provide a method for convenient minichromosome production.

The chromosome axis plays a crucial role in meiotic recombination. Here, we study the function of ASY1, the Arabidopsis homolog of the yeast chromosome axis-associated component Hop1. Specifically, we characterized cross-over (CO) distribution in female and male meiosis by deep sequencing of the progeny of an allelic series of asy1 mutants. Combining data from nearly 1,000 individual plants, we find that reduced ASY1 functionality leads to genomic instability and sometimes drastic genomic rearrangements. We further observed that COs are less frequent and appear in more distal chromosomal regions in plants with no or reduced ASY1 functionality, consistent with previous analyses. However, our sequencing approach revealed that the reduction in CO number is not as dramatic as suggested by cytological analyses. Analysis of double mutants of asy1 with mutants with three other CO factors, MUS81, MSH4, and MSH5, as well as the determination of foci number of the CO regulator MLH1 demonstrates that the majority of the COs in asy1, similar to the situation in the wildtype (WT), largely belong to the class I, which are subject to interference. However, these COs are redistributed in asy1 mutants and typically appear much closer than in the WT. Hence, ASY1 plays a key role in CO interference that spaces COs along a chromosome. Conversely, since a large proportion of chromosomes do not receive any CO, we conclude that CO assurance, the process that ensures the obligatory assignment of one CO per chromosome, is also affected in asy1 mutants.

Chromoanagenesis is a catastrophic event that involves localized chromosomal shattering and reorganization. In this study, we report a case of chromoanagenesis resulting from defective meiosis in the MEIOTIC ASYNAPTIC MUTANT 1 (asy1) background in Arabidopsis thaliana. We provide a detailed characterization of the genomic structure of this individual with a severely shattered segment of chromosome 1. We identified 260 novel DNA junctions in the affected region, most of which affect gene sequence on 1 or both sides of the junction. Our results confirm that asy1-related defective meiosis is a potential trigger for chromoanagenesis. This is the first example of chromoanagenesis associated with female meiosis and indicates the potential for genome evolution during oogenesis. Chromoanagenesis is a complex and catastrophic event that results in severely restructured chromosomes. It has been identified in cancer cells and in some plant samples, after specific triggering events. Here, we identified this kind of genome restructuring in a mutant that exhibits defective meiosis in the model plant system Arabidopsis thaliana.

Diverse morphological, cellular and physiological changes occur during seed maturation in Bixa orellana when the seed tissues form specialized cell glands that produce reddish latex with high bixin amounts. Transcriptomic profiling during seed development in three B. orellana accessions (P12, N4 and N5) with contrasting morphologic characteristics showed enrichment in pathways of triterpenes, sesquiterpenes, and cuticular wax biosynthesis. WGCNA allows groups of all identified genes in six modules the module turquoise, the largest and highly correlated with the bixin content. The high number of genes in this module suggests a diversification of regulatory mechanisms for bixin accumulation with the genes belonging to isoprene, triterpenes and carotene pathways, being more highly correlated with the bixin content. Analysis of key genes of the mevalonate (MVA) and the 2C-methyl-D-erythritol-4-phosphate (MEP) pathways revealed specific activities of orthologs of BoHMGR, BoFFP, BoDXS, and BoHDR. This suggests that isoprenoid production is necessary for compounds included in the reddish latex of developing seeds. The carotenoid-related genes BoPSY2, BoPDS1 and BoZDS displayed a high correlation with bixin production, consistent with the requirement for carotene precursors for apocarotenoid biosynthesis. The BoCCD gene member (BoCCD4-4) and some BoALDH (ALDH2B7.2 and ALDH3I1) and BoMET (BoSABATH1 and BoSABATH8) gene members were highly correlated to bixin in the final seed development stage. This suggested a contributing role for several genes in apocarotenoid production. The results revealed high genetic complexity in the biosynthesis of reddish latex and bixin in specialized seed cell glands in different accessions of B. orellana suggesting gene expression coordination between both metabolite biosynthesis processes.

Wide crosses result in postzygotic elimination of one parental chromosome set, but the mechanisms that result in such differential fate are poorly understood. Here, we show that alterations of centromeric histone H3 (CENH3) lead to its selective removal from centromeres of mature Arabidopsis eggs and early zygotes, while wild-type CENH3 persists. In the hybrid zygotes and embryos, CENH3 and essential centromere proteins load preferentially on the CENH3-rich centromeres of the wild-type parent, while CENH3-depleted centromeres fail to reconstitute new CENH3-chromatin and the kinetochore and are frequently lost. Genome elimination is opposed by E3 ubiquitin ligase VIM1. We propose a model based on cooperative binding of CENH3 to chromatin to explain the differential CENH3 loading rates. Thus, parental CENH3 polymorphisms result in epigenetically distinct centromeres that instantiate a strong mating barrier and produce haploids.

Large-scale structural variations, such as chromosomal translocations, can have profound effects on fitness and phenotype, but are difficult to identify and characterize. Here, we describe a simple and effective method aimed at identifying translocations using only the dosage of sequence reads mapped on the reference genome. We binned reads on genomic segments sized according to sequencing coverage and identified instances when copy number segregated in populations. For each dosage-polymorphic 1 Mb bin, we tested independence, effectively an apparent linkage disequilibrium (LD), with other variable bins. In nine potato (Solanum tuberosum) dihaploid families translocations affecting pericentromeric regions were common and in two cases were due to genomic misassembly. In two populations, we found evidence for translocation affecting euchromatic arms. In cv. PI 310467, a nonreciprocal translocation between chromosomes (chr.) 7 and 8 resulted in a 5-3 copy number change affecting several Mb at the respective chromosome tips. In cv. "Alca Tarma," the terminal arm of chr. 4 translocated to the tip of chr. 1. Using oligonucleotide-based fluorescent in situ hybridization painting probes (oligo-FISH), we tested and confirmed the predicted arrangement in PI 310467. In 192 natural accessions of Arabidopsis thaliana, dosage haplotypes tended to vary continuously and resulted in higher noise, while apparent LD between pericentromeric regions suggested the effect of repeats. This method, LD-CNV, should be useful in species where translocations are suspected because it tests linkage without the need for genotyping.

Chromoanagenesis is a genomic catastrophe that results in chromosomal shattering and reassembly. These extreme single chromosome events were first identified in cancer, and have since been observed in other systems, but have so far only been formally documented in plants in the context of haploid induction crosses. The frequency, origins, consequences, and evolutionary impact of such major chromosomal remodeling in other situations remain obscure. Here, we demonstrate the occurrence of chromoanagenesis in poplar (Populus sp.) trees produced from gamma-irradiated pollen. Specifically, in this population of siblings carrying indel mutations, two individuals exhibited highly frequent copy number variation (CNV) clustered on a single chromosome, one of the hallmarks of chromoanagenesis. Using short-read sequencing, we confirmed the presence of clustered segmental rearrangement. Independently, we identified and validated novel DNA junctions and confirmed that they were clustered and corresponded to these rearrangements. Our reconstruction of the novel sequences suggests that the chromosomal segments have reorganized randomly to produce a novel rearranged chromosome but that two different mechanisms might be at play. Our results indicate that gamma irradiation can trigger chromoanagenesis, suggesting that this may also occur when natural or induced mutagens cause DNA breaks. We further demonstrate that such events can be tolerated in poplar, and even replicated clonally, providing an attractive system for more in-depth investigations of their consequences.

Abstract Background The availability of low cost sequencing has spurred its application to discovery and typing of variation, including variation induced by mutagenesis. Mutation discovery is challenging as it requires a substantial amount of sequencing and analysis to detect very rare changes and distinguish them from noise. Also challenging are the cases when the organism of interest has not been sequenced or is highly divergent from the reference. Results We describe the development of a simple method for reduced representation sequencing. Input DNA was digested with a single restriction enzyme and ligated to Y adapters modified to contain a sequence barcode and to provide a compatible overhang for ligation. We demonstrated the efficiency of this method at SNP discovery using rice and arabidopsis. To test its suitability for the discovery of very rare SNP, one control and three mutagenized rice individuals (1, 5 and 10 mM sodium azide) were used to prepare genomic libraries for Illumina sequencers by ligating barcoded adapters to NlaIII restriction sites. For genome-dependent discovery 15-30 million of 80 base reads per individual were aligned to the reference sequence achieving individual sequencing coverage from 7 to 15×. We identified high-confidence base changes by comparing sequences across individuals and identified instances consistent with mutations, i.e. changes that were found in a single treated individual and were solely GC to AT transitions. For genome-independent discovery 70-mers were extracted from the sequence of the control individual and single-copy sequence was identified by comparing the 70-mers across samples to evaluate copy number and variation. This de novo "genome" was used to align the reads and identify mutations as above. Covering approximately 1/5 of the 380 Mb genome of rice we detected mutation densities ranging from 0.6 to 4 per Mb of diploid DNA depending on the mutagenic treatment. Conclusions The combination of a simple and cost-effective library construction method, with Illumina sequencing, and the use of a bioinformatic pipeline allows practical SNP discovery regardless of whether a genomic reference is available.

Abstract Background TILLING (Targeting induced local lesions IN genomes) is an efficient reverse genetics approach for detecting induced mutations in pools of individuals. Combined with the high-throughput of next-generation sequencing technologies, and the resolving power of overlapping pool design, TILLING provides an efficient and economical platform for functional genomics across thousands of organisms. Results We propose a probabilistic method for calling TILLING-induced mutations, and their carriers, from high throughput sequencing data of overlapping population pools, where each individual occurs in two pools. We assign a probability score to each sequence position by applying Bayes' Theorem to a simplified binomial model of sequencing error and expected mutations, taking into account the coverage level. We test the performance of our method on variable quality, high-throughput sequences from wheat and rice mutagenized populations. Conclusions We show that our method effectively discovers mutations in large populations with sensitivity of 92.5% and specificity of 99.8%. It also outperforms existing SNP detection methods in detecting real mutations, especially at higher levels of coverage variability across sequenced pools, and in lower quality short reads sequence data. The implementation of our method is available from: http://www.cs.ucdavis.edu/filkov/CAMBa/.

The Targeting Induced Local Lesions in Genomes (TILLING) technology is a reverse genetic strategy broadly applicable to every kind of genome and represents an attractive tool for functional genomic and agronomic applications. It consists of chemical random mutagenesis followed by high-throughput screening of point mutations in targeted genomic regions. Although multiple methods for mutation discovery in amplicons have been described, next-generation sequencing (NGS) is the tool of choice for mutation detection because it quickly allows for the analysis of a large number of amplicons. The aim of the present work was to screen a previously generated sunflower TILLING population and identify alterations in genes involved in several important and complex physiological processes. Twenty-one candidate sunflower genes were chosen as targets for the screening. The TILLING by sequencing strategy allowed us to identify multiple mutations in selected genes and we subsequently validated 16 mutations in 11 different genes through Sanger sequencing. In addition to addressing challenges posed by outcrossing, our detection and validation of mutations in multiple regulatory loci highlights the importance of this sunflower population as a genetic resource.

Abstract Background Camelina sativa, an oilseed crop in the Brassicaceae family, has inspired renewed interest due to its potential for biofuels applications. Little is understood of the nature of the C. sativa genome, however. A study was undertaken to characterize two genes in the fatty acid biosynthesis pathway, fatty acid desaturase (FAD) 2 and fatty acid elongase (FAE) 1, which revealed unexpected complexity in the C. sativa genome. Results In C. sativa, Southern analysis indicates the presence of three copies of both FAD2 and FAE1 as well as LFY, a known single copy gene in other species. All three copies of both CsFAD2 and CsFAE1 are expressed in developing seeds, and sequence alignments show that previously described conserved sites are present, suggesting that all three copies of both genes could be functional. The regions downstream of CsFAD2 and upstream of CsFAE1 demonstrate co-linearity with the Arabidopsis genome. In addition, three expressed haplotypes were observed for six predicted single-copy genes in 454 sequencing analysis and results from flow cytometry indicate that the DNA content of C. sativa is approximately three-fold that of diploid Camelina relatives. Phylogenetic analyses further support a history of duplication and indicate that C. sativa and C. microcarpa might share a parental genome. Conclusions There is compelling evidence for triplication of the C. sativa genome, including a larger chromosome number and three-fold larger measured genome size than other Camelina relatives, three isolated copies of FAD2, FAE1, and the KCS17-FAE1 intergenic region, and three expressed haplotypes observed for six predicted single-copy genes. Based on these results, we propose that C. sativa be considered an allohexaploid. The characterization of fatty acid synthesis pathway genes will allow for the future manipulation of oil composition of this emerging biofuel crop; however, targeted manipulations of oil composition and general development of C. sativa should consider and, when possible take advantage of, the implications of polyploidy.

Abstract Background Wheat (Triticum ssp.) is an important food source for humans in many regions around the world. However, the ability to understand and modify gene function for crop improvement is hindered by the lack of available genomic resources. TILLING is a powerful reverse genetics approach that combines chemical mutagenesis with a high-throughput screen for mutations. Wheat is specially well-suited for TILLING due to the high mutation densities tolerated by polyploids, which allow for very efficient screens. Despite this, few TILLING populations are currently available. In addition, current TILLING screening protocols require high-throughput genotyping platforms, limiting their use. Results We developed mutant populations of pasta and common wheat and organized them for TILLING. To simplify and decrease costs, we developed a non-denaturing polyacrylamide gel set-up that uses ethidium bromide to detect fragments generated by crude celery juice extract digestion of heteroduplexes. This detection method had similar sensitivity as traditional LI-COR screens, suggesting that it represents a valid alternative. We developed genome-specific primers to circumvent the presence of multiple homoeologous copies of our target genes. Each mutant library was characterized by TILLING multiple genes, revealing high mutation densities in both the hexaploid (~1/38 kb) and tetraploid (~1/51 kb) populations for 50% GC targets. These mutation frequencies predict that screening 1,536 lines for an effective target region of 1.3 kb with 50% GC content will result in ~52 hexaploid and ~39 tetraploid mutant alleles. This implies a high probability of obtaining knock-out alleles (P = 0.91 for hexaploid, P = 0.84 for tetraploid), in addition to multiple missense mutations. In total, we identified over 275 novel alleles in eleven targeted gene/genome combinations in hexaploid and tetraploid wheat and have validated the presence of a subset of them in our seed stock. Conclusion We have generated reverse genetics TILLING resources for pasta and bread wheat and achieved a high mutation density in both populations. We also developed a modified screening method that will lower barriers to adopt this promising technology. We hope that the use of this reverse genetics resource will enable more researchers to pursue wheat functional genomics and provide novel allelic diversity for wheat improvement.

In angiosperms, endosperm development comprises a series of developmental transitions controlled by genetic and epigenetic mechanisms that are initiated after double fertilization. Polycomb repressive complex 2 (PRC2) is a key component of these mechanisms that mediate histone H3 lysine 27 trimethylation (H3K27me3); the action of PRC2 is well described in Arabidopsis thaliana but remains uncertain in cereals. In this study, we demonstrate that mutation of the rice (Oryza sativa) gene EMBRYONIC FLOWER2a (OsEMF2a), encoding a zinc-finger containing component of PRC2, causes an autonomous endosperm phenotype involving proliferation of the central cell nuclei with separate cytoplasmic domains, even in the absence of fertilization. Detailed cytological and transcriptomic analyses revealed that the autonomous endosperm can produce storage compounds, starch granules, and protein bodies specific to the endosperm. These events have not been reported in Arabidopsis. After fertilization, we observed an abnormally delayed developmental transition in the endosperm. Transcriptome and H3K27me3 ChIP-seq analyses using endosperm from the emf2a mutant identified downstream targets of PRC2. These included >100 transcription factor genes such as type-I MADS-box genes, which are likely required for endosperm development. Our results demonstrate that OsEMF2a-containing PRC2 controls endosperm developmental programs before and after fertilization.

Mint oil is a key source of natural flavors with wide industrial applications. Two unbalanced polyploid cultivars named Native (Mentha Spicata L) and Scotch (M. × gracilis Sole) are the main producers of spearmint type oil, which is characterized by high levels of the monoterpenes (-)-carvone and (-)-limonene. These cultivars have been the backbone of spearmint oil production for decades, while breeding and improvement remained largely unexplored, in part, due to sterility in cultivated lines. Here we show that sexual breeding at the diploid level can be leveraged to develop new varieties that produce spearmint type oil, along with the improvement of other important traits. Using field trials and GC-FID oil analysis we characterized plant materials from a public germplasm repository and identified a diploid accession that exhibited 89.5% increase in oil yield, compared to the industry standard, and another that produces spearmint type oil. Spearmint-type oil was present at high frequency in a segregating F2 population (32/160) produced from these two accessions. Field-testing of ten of these F2 lines showed segregation for oil yield and confirmed the production of spearmint-type oil profiles. Two of these lines combined high yield and spearmint-type oil with acceptable analytic and sensory profiles. These results demonstrate that spearmint-type oil can be produced in a diploid background with high yield potential, providing a simpler genetic system for the development of improved spearmint varieties.

Abstract Background Soybean (Glycine max L. Merr.) is an important nitrogen-fixing crop that provides much of the world's protein and oil. However, the available tools for investigation of soybean gene function are limited. Nevertheless, chemical mutagenesis can be applied to soybean followed by screening for mutations in a target of interest using a strategy known as Targeting Induced Local Lesions IN Genomes (TILLING). We have applied TILLING to four mutagenized soybean populations, three of which were treated with ethyl methanesulfonate (EMS) and one with N-nitroso-N-methylurea (NMU). Results We screened seven targets in each population and discovered a total of 116 induced mutations. The NMU-treated population and one EMS mutagenized population had similar mutation density (~1/140 kb), while another EMS population had a mutation density of ~1/250 kb. The remaining population had a mutation density of ~1/550 kb. Because of soybean's polyploid history, PCR amplification of multiple targets could impede mutation discovery. Indeed, one set of primers tested in this study amplified more than a single target and produced low quality data. To address this problem, we removed an extraneous target by pretreating genomic DNA with a restriction enzyme. Digestion of the template eliminated amplification of the extraneous target and allowed the identification of four additional mutant alleles compared to untreated template. Conclusion The development of four independent populations with considerable mutation density, together with an additional method for screening closely related targets, indicates that soybean is a suitable organism for high-throughput mutation discovery even with its extensively duplicated genome.

Abstract Background The composition of the individual eukaryote's genome and its variation within a species remain poorly defined. Even for a sequenced genome such as that of the model plant Arabidopsis thaliana accession Col-0, the large arrays of heterochromatic repeats are incompletely sequenced, with gaps of uncertain size persisting in them. Results Using geographically separate populations of A. thaliana, we assayed variation in the heterochromatic repeat arrays using two independent methods and identified significant polymorphism among them, with variation by as much as a factor of two in the centromeric 180 bp repeat, in the 45S rDNA arrays and in the Athila retroelements. In the accession with highest genome size as measured by flow cytometry, Loh-0, we found more than a two-fold increase in 5S RNA gene copies relative to Col-0; results from fluorescence in situ hybridization with 5S probes were consistent with the existence of size polymorphism between Loh-0 and Col-0 at the 5S loci. Comparative genomic hybridization results of Loh-0 and Col-0 did not support contiguous variation in copy number of protein-coding genes on the scale needed to explain their observed genome size difference. We developed a computational data model to test whether the variation we measured in the repeat fractions could account for the different genome sizes determined with flow cytometry, and found that this proposed relationship could account for about 50% of the variance in genome size among the accessions. Conclusion Our analyses are consistent with substantial repeat number polymorphism for 5S and 45S ribosomal genes among accession of A. thaliana. Differences are also suggested for centromeric and pericentromeric repeats. Our analysis also points to the difficulties in measuring the repeated fraction of the genome and suggests that independent validation of genome size should be sought in addition to flow cytometric measurements.

Abstract Background Rice is both a food source for a majority of the world's population and an important model system. Available functional genomics resources include targeted insertion mutagenesis and transgenic tools. While these can be powerful, a non-transgenic, unbiased targeted mutagenesis method that can generate a range of allele types would add considerably to the analysis of the rice genome. TILLING (Targeting Induced Local Lesions in Genomes), a general reverse genetic technique that combines traditional mutagenesis with high throughput methods for mutation discovery, is such a method. Results To apply TILLING to rice, we developed two mutagenized rice populations. One population was developed by treatment with the chemical mutagen ethyl methanesulphonate (EMS), and the other with a combination of sodium azide plus methyl-nitrosourea (Az-MNU). To find induced mutations, target regions of 0.7–1.5 kilobases were PCR amplified using gene specific primers labeled with fluorescent dyes. Heteroduplexes were formed through denaturation and annealing of PCR products, mismatches digested with a crude preparation of CEL I nuclease and cleaved fragments visualized using denaturing polyacrylamide gel electrophoresis. In 10 target genes screened, we identified 27 nucleotide changes in the EMS-treated population and 30 in the Az-MNU population. Conclusion We estimate that the density of induced mutations is two- to threefold higher than previously reported rice populations (about 1/300 kb). By comparison to other plants used in public TILLING services, we conclude that the populations described here would be suitable for use in a large scale TILLING project.

Abstract Background Going from a gene sequence to its function in the context of a whole organism requires a strategy for targeting mutations, referred to as reverse genetics. Reverse genetics is highly desirable in the modern genomics era; however, the most powerful methods are generally restricted to a few model organisms. Previously, we introduced a reverse-genetic strategy with the potential for general applicability to organisms that lack well-developed genetic tools. Our TILLING (Targeting Induced Local Lesions IN Genomes) method uses chemical mutagenesis followed by screening for single-base changes to discover induced mutations that alter protein function. TILLING was shown to be an effective reverse genetic strategy by the establishment of a high-throughput TILLING facility and the delivery of thousands of point mutations in hundreds of Arabidopsis genes to members of the plant biology community. Results We demonstrate that high-throughput TILLING is applicable to maize, an important crop plant with a large genome but with limited reverse-genetic resources currently available. We screened pools of DNA samples for mutations in 1-kb segments from 11 different genes, obtaining 17 independent induced mutations from a population of 750 pollen-mutagenized maize plants. One of the genes targeted was the DMT102 chromomethylase gene, for which we obtained an allelic series of three missense mutations that are predicted to be strongly deleterious. Conclusions Our findings indicate that TILLING is a broadly applicable and efficient reverse-genetic strategy. We are establishing a public TILLING service for maize modeled on the existing Arabidopsis TILLING Project.

Creating true-breeding lines is a critical step in plant breeding. Novel, completely homozygous true-breeding lines can be generated by doubled haploid technology in single generation. Haploid induction through modification of the centromere-specific histone 3 variant (CENH3), including chimeric proteins, expression of non-native CENH3 and single amino acid substitutions, has been shown to induce, on outcrossing to wild type, haploid progeny possessing only the genome of the wild-type parent, in Arabidopsis thaliana. Here, we report the characterization of 31 additional EMS-inducible amino acid substitutions in CENH3 for their ability to complement a knockout in the endogenous CENH3 gene and induce haploid progeny when pollinated by the wild type. We also tested the effect of double amino acid changes, which might be generated through a second round of EMS mutagenesis. Finally, we report on the effects of CRISPR/Cas9-mediated in-frame deletions in the αN helix of the CENH3 histone fold domain. Remarkably, we found that complete deletion of the αN helix, which is conserved throughout angiosperms, results in plants which exhibit normal growth and fertility while acting as excellent haploid inducers when pollinated by wild-type pollen. Both of these technologies, CRISPR mutagenesis and EMS mutagenesis, represent non-transgenic approaches to the generation of haploid inducers.

The challenges of breeding autotetraploid potato (Solanum tuberosum) have motivated the development of alternative breeding strategies. A common approach is to obtain uniparental dihaploids from a tetraploid of interest through pollination with S. tuberosum Andigenum Group (formerly S. phureja) cultivars. The mechanism underlying haploid formation of these crosses is unclear, and questions regarding the frequency of paternal DNA transmission remain. Previous reports have described aneuploid and euploid progeny that, in some cases, displayed genetic markers from the haploid inducer (HI). Here, we surveyed a population of 167 presumed dihaploids for large-scale structural variation that would underlie chromosomal addition from the HI, and for small-scale introgression of genetic markers. In 19 progeny, we detected 10 of the 12 possible trisomies and, in all cases, demonstrated the noninducer parent origin of the additional chromosome. Deep sequencing indicated that occasional, short-tract signals appearing to be of HI origin were better explained as technical artifacts. Leveraging recurring copy number variation patterns, we documented subchromosomal dosage variation indicating segregation of polymorphic maternal haplotypes. Collectively, 52% of the assayed chromosomal loci were classified as dosage variable. Our findings help elucidate the genomic consequences of potato haploid induction and suggest that most potato dihaploids will be free of residual pollinator DNA.

Genomic imprinting regulates parent-specific transcript dosage during seed development and is mainly confined to the endosperm. Elucidation of the function of many imprinted genes has been hampered by the lack of corresponding mutant phenotypes, and the role of imprinting is mainly associated with genome dosage regulation or allocation of resources. Disruption of imprinted genes has also been suggested to mediate endosperm-based post-zygotic hybrid barriers depending on genetic variation and gene dosage. Here, we have analyzed the conservation of a clade from the MADS-box type I class transcription factors in the closely related species Arabidopsis arenosa, A. lyrata, and A. thaliana, and show that AGL36-like genes are imprinted and maternally expressed in seeds of Arabidopsis species and in hybrid seeds between outbreeding species. In hybridizations between outbreeding and inbreeding species the paternally silenced allele of the AGL36-like gene is reactivated in the hybrid, demonstrating that also maternally expressed imprinted genes are perturbed during hybridization and that such effects on imprinted genes are specific to the species combination. Furthermore, we also demonstrate a quantitative effect of genetic diversity and temperature on the strength of the post-zygotic hybridization barrier. Markedly, a small decrease in temperature during seed development increases the survival of hybrid F1 seeds, suggesting that abiotic and genetic parameters play important roles in post-zygotic species barriers, pointing at evolutionary scenarios favoring such effects. OPEN RESEARCH BADGES: This article has earned an Open Data Badge for making publicly available the digitally-shareable data necessary to reproduce the reported results. The data is available at https://www.ncbi.nlm.nih.gov/bioproject/?term=PRJNA562212. All sequences generated in this study have been deposited in the National Center for Biotechnology Information Sequence Read Archive (https://www.ncbi.nlm.nih.gov/sra/) with project number PRJNA562212.

BackgroundSalt stress, a major plant environmental stress, is a critical constraint for rice productivity. Dissecting the genetic loci controlling salt tolerance in rice for improving productivity, especially at the flowering stage, remains challenging. Here, we conducted a genome-wide association study (GWAS) of salt tolerance based on exome sequencing of the Thai rice accessions.ResultsPhotosynthetic parameters and cell membrane stability under salt stress at the flowering stage; and yield-related traits of 104 Thai rice (Oryza sativa L.) accessions belonging to the indica subspecies were evaluated. The rice accessions were subjected to exome sequencing, resulting in 112,565 single nucleotide polymorphisms (SNPs) called with a minor allele frequency of at least 5%. LD decay analysis of the panel indicates that the average LD for SNPs at 20 kb distance from each other was 0.34 (r2), which decayed to its half value (~ 0.17) at around 80 kb. By GWAS performed using mixed linear model, two hundred loci containing 448 SNPs on exons were identified based on the salt susceptibility index of the net photosynthetic rate at day 6 after salt stress; and the number of panicles, filled grains and unfilled grains per plant. One hundred and forty six genes, which accounted for 73% of the identified loci, co-localized with the previously reported salt quantitative trait loci (QTLs). The top four regions that contained a high number of significant SNPs were found on chromosome 8, 12, 1 and 2. While many are novel, their annotation is consistent with potential involvement in plant salt tolerance and in related agronomic traits. These significant SNPs greatly help narrow down the region within these QTLs where the likely underlying candidate genes can be identified.ConclusionsInsight into the contribution of potential genes controlling salt tolerance from this GWAS provides further understanding of salt tolerance mechanisms of rice at the flowering stage, which can help improve yield productivity under salinity via gene cloning and genomic selection.

Unfavourable environmental conditions, including soil salinity, lead to decreased rice (Oryza sativa L.) productivity, especially at the reproductive stage. In this study, we examined 30 rice varieties, which revealed significant differences in the photosynthetic performance responses under salt stress conditions during the reproductive stage, which ultimately affected yield components after recovery. In rice with a correlation between net photosynthetic rate (PN) and intercellular CO2 concentration (Ci) under salt stress, PN was found to be negatively correlated with filled grain number after recovery. Applying stringent criteria, we identified 130,317 SNPs and 15,396 InDels between two "high-yield rice" varieties and two "low-yield rice" varieties with contrasting photosynthesis and grain yield characteristics. A total of 2,089 genes containing high- and moderate-impact SNPs or InDels were evaluated by gene ontology (GO) enrichment analysis, resulting in over-represented terms in the apoptotic process and kinase activity. Among these genes, 262 were highly expressed in reproductive tissues, and most were annotated as receptor-like protein kinases. These findings highlight the importance of variations in signaling components in the genome and these loci can serve as potential genes in rice breeding to produce a variety with salt avoidance that leads to increased yield in saline soil.

Gene dosage variation and the associated changes in gene expression influence a wide variety of traits, ranging from cancer in humans to yield in plants. It is also expected to affect important traits of ecological and agronomic importance in forest trees, but this variation has not been systematically characterized or exploited. Here we performed a comprehensive scan of the Populus genome for dosage-sensitive loci affecting quantitative trait variation for spring and fall phenology and biomass production. The study population was a large collection of clonally propagated F1 hybrid lines of Populus that saturate the genome 10-fold with deletions and insertions (indels) of known sizes and positions. As a group, the phenotypic means of the indel lines consistently differed from control nonindel lines, with an overall negative effect of both insertions and deletions on all biomass-related traits but more diverse effects and an overall wider phenotypic distribution of the indel lines for the phenology-related traits. We also investigated the correlation between gene dosage at specific chromosomal locations and phenotype, to identify dosage quantitative trait loci (dQTL). Such dQTL were detected for most phenotypes examined, but stronger effect dQTL were identified for the phenology-related traits than for the biomass traits. Our genome-wide screen for dosage sensitivity in a higher eukaryote demonstrates the importance of global genomic balance and the impact of dosage on life history traits.

Rice nucleolin (OsNUC1), consisting of two isoforms, OsNUC1-L and OsNUC1-S, is a multifunctional protein involved in salt-stress tolerance. Here, OsNUC1-S's function was investigated using transgenic rice lines overexpressing OsNUC1-S. Under non-stress conditions, the transgenic lines showed a lower yield, but higher net photosynthesis rates, stomatal conductance, and transpiration rates than wild type only in the second leaves, while in the flag leaves, these parameters were similar among the lines. However, under salt-stress conditions at the booting stage, the higher yields in transgenic lines were detected. Moreover, the gas exchange parameters of the transgenic lines were higher in both flag and second leaves, suggesting a role for OsNUC1-S overexpression in photosynthesis adaptation under salt-stress conditions. Moreover, the overexpression lines could maintain light-saturation points under salt-stress conditions, while a decrease in the light-saturation point owing to salt stress was found in wild type. Based on a transcriptome comparison between wild type and a transgenic line, after 3 and 9 days of salt stress, the significantly differentially expressed genes were enriched in the metabolic process of nucleic acid and macromolecule, photosynthesis, water transport, and cellular homeostasis processes, leading to the better performance of photosynthetic processes under salt-stress conditions at the booting stage.

BackgroundCalmodulin (CaM) is an important calcium sensor protein that transduces Ca2+ signals in plant stress signaling pathways. A previous study has revealed that transgenic rice over-expressing the calmodulin gene OsCam1-1 (LOC_Os03g20370) is more tolerant to salt stress than wild type. To elucidate the role of OsCam1-1 in the salt stress response mechanism, downstream components of the OsCam1-1-mediated response were identified and investigated by transcriptome profiling and target identification.ResultsTranscriptome profiling of transgenic 'Khao Dawk Mali 105' rice over-expressing OsCam1-1 and wild type rice showed that overexpression of OsCam1-1 widely affected the expression of genes involved in several cellular processes under salt stress, including signaling, hormone-mediated regulation, transcription, lipid metabolism, carbohydrate metabolism, secondary metabolism, photosynthesis, glycolysis, tricarboxylic acid (TCA) cycle and glyoxylate cycle. Under salt stress, the photosynthesis rate in the transgenic rice was slightly lower than in wild type, while sucrose and starch contents were higher, suggesting that energy and carbon metabolism were affected by OsCam1-1 overexpression. Additionally, four known and six novel CaM-interacting proteins were identified by cDNA expression library screening with the recombinant OsCaM1. GO terms enriched in their associated proteins that matched those of the differentially expressed genes affected by OsCam1-1 overexpression revealed various downstream cellular processes that could potentially be regulated by OsCaM1 through their actions.ConclusionsThe diverse cellular processes affected by OsCam1-1 overexpression and possessed by the identified CaM1-interacting proteins corroborate the notion that CaM signal transduction pathways compose a complex network of downstream components involved in several cellular processes. These findings suggest that under salt stress, CaM activity elevates metabolic enzymes involved in central energy pathways, which promote or at least maintain the production of energy under the limitation of photosynthesis.

The rice chromosome segment substitution lines (CSSLs) of 'Khao Dawk Mali 105' ('KDML105') genetic background were developed by backcrossing with 'KDML105' rice and transferring the region from chromosome 1 of DH212 which was expected to contain the full putative salt tolerance genetic region. Line of CSSL11, CSSL12, and CSSL16 contained the full putative salt tolerance genetic region were evaluated with the parental lines, 'KDML105' and DH212 at seedling stage of rice. The physiological responses in rice plants were grown under normal condition and 75 mM of NaCl, and then comparative photosynthetic parameters, chlorophyll fluorescence parameters, PhiPS2, ETR, NPQ, as well as growth analysis. In this article, the data of physiological response evaluation in rice at seedling stage after salt stress treatment can be found. This can be useful as the information of the photosynthesis response to salt stress to other rice cultivars and related species.

The CRISPR/Cas9 nuclease system is a powerful and flexible tool for genome editing, and novel applications of this system are being developed rapidly. Here, we used CRISPR/Cas9 to target the FAD2 gene in Arabidopsis thaliana and in the closely related emerging oil seed plant, Camelina sativa, with the goal of improving seed oil composition. We successfully obtained Camelina seeds in which oleic acid content was increased from 16% to over 50% of the fatty acid composition. These increases were associated with significant decreases in the less desirable polyunsaturated fatty acids, linoleic acid (i.e. a decrease from ~16% to

During cell division, spindle fibers attach to chromosomes at centromeres. The DNA sequence at regional centromeres is fast evolving with no conserved genetic signature for centromere identity. Instead CENH3, a centromere-specific histone H3 variant, is the epigenetic signature that specifies centromere location across both plant and animal kingdoms. Paradoxically, CENH3 is also adaptively evolving. An ongoing question is whether CENH3 evolution is driven by a functional relationship with the underlying DNA sequence. Here, we demonstrate that despite extensive protein sequence divergence, CENH3 histones from distant species assemble centromeres on the same underlying DNA sequence. We first characterized the organization and diversity of centromere repeats in wild-type Arabidopsis thaliana We show that A. thaliana CENH3-containing nucleosomes exhibit a strong preference for a unique subset of centromeric repeats. These sequences are largely missing from the genome assemblies and represent the youngest and most homogeneous class of repeats. Next, we tested the evolutionary specificity of this interaction in a background in which the native A. thaliana CENH3 is replaced with CENH3s from distant species. Strikingly, we find that CENH3 from Lepidium oleraceum and Zea mays, although specifying epigenetically weaker centromeres that result in genome elimination upon outcrossing, show a binding pattern on A. thaliana centromere repeats that is indistinguishable from the native CENH3. Our results demonstrate positional stability of a highly diverged CENH3 on independently evolved repeats, suggesting that the sequence specificity of centromeres is determined by a mechanism independent of CENH3.

Introduction: Owing to advances in high-throughput genome sequencing, QTLSeq mapping of salt tolerance traits is a major platform for identifying soil-salinity tolerance QTLs to accelerate marker-assisted selection for salt-tolerant rice varieties. We performed QTL-BSA-Seq in the seedling stage of rice from a genetic cross of the extreme salt-sensitive variety, IR29, and "Jao Khao" (JK), a Thai salt-tolerant variety. Methods: A total of 462 F2 progeny grown in soil and treated with 160 mM NaCl were used as the QTL mapping population. Two high- and low-bulk sets, based on cell membrane stability (CMS) and tiller number at the recovery stage (TN), were equally sampled. The genomes of each pool were sequenced, and statistical significance of QTL was calculated using QTLseq and G prime (G') analysis, which is based on calculating the allele frequency differences or D(SNP index). Results: Both methods detected the overlapping interval region, wherein CMSbulk was mapped at two loci in the 38.41-38.85 Mb region with 336 SNPs on chromosome 1 (qCMS1) and the 26.13-26.80 Mb region with 1,011 SNPs on chromosome 3 (qCMS3); the D(SNP index) peaks were -0.2709 and 0.3127, respectively. TN-bulk was mapped at only one locus in the overlapping 38.26-38.95 Mb region on chromosome 1 with 575 SNPs (qTN1) and a D(SNP index) peak of -0.3544. These identified QTLs in two different genetic backgrounds of segregating populations derived from JK were validated. The results confirmed the colocalization of the qCMS1 and qTN1 traits on chromosome 1. Based on the CMS trait, qCMS1/qTN1 stably expressed 6%-18% of the phenotypic variance in the two validation populations, while qCMS1/qTN1 accounted for 16%-20% of the phenotypic variance in one validation population based on the TN trait. Conclusion: The findings confirm that the CMS and TN traits are tightly linked to the long arm of chromosome 1 rather than to chromosome 3. The validated qCMS-TN1 QTL can be used for gene/QTL pyramiding in marker-assisted selection to expedite breeding for salt resistance in rice at the seedling stage. [ABSTRACT FROM AUTHOR]

Lettuce (Lactuca sativa) seeds exhibit thermoinhibition, or failure to complete germination when imbibed at warm temperatures. Chemical mutagenesis was employed to develop lettuce lines that exhibit germination thermotolerance. Two independent thermotolerant lettuce seed mutant lines, TG01 and TG10, were generated through ethyl methanesulfonate mutagenesis. Genetic and physiological analyses indicated that these two mutations were allelic and recessive. To identify the causal gene(s), we applied bulked segregant analysis by whole genome sequencing. For each mutant, bulked DNA samples of segregating thermotolerant (mutant) seeds were sequenced and analyzed for homozygous single-nucleotide polymorphisms. Two independent candidate mutations were identified at different physical positions in the zeaxanthin epoxidase gene (ABSCISIC ACID DEFICIENT 1/ZEAXANTHIN EPOXIDASE, or ABA1/ZEP) in TG01 and TG10. The mutation in TG01 caused an amino acid replacement, whereas the mutation in TG10 resulted in alternative mRNA splicing. Endogenous abscisic acid contents were reduced in both mutants, and expression of the ABA1 gene from wild-type lettuce under its own promoter fully complemented the TG01 mutant. Conventional genetic mapping confirmed that the causal mutations were located near the ZEP/ABA1 gene, but the bulked segregant whole genome sequencing approach more efficiently identified the specific gene responsible for the phenotype.

Bixin is an apocarotenoid obtained from the seed aril of Bixa orellana L., a tropical plant known as achiote in Mexico. This compound is the second most commonly used natural colouring for food and pharmaceutical industries. B. orellana is an outcrossing species that displays high genetic variability. Recently, the colour traits of sexual organs were associated with the biosynthesis and accumulation of bixin in mature seeds. Herein, we describe a new approach for genotype–phenotype association by surveying lycopene beta-cyclase (Boβ-LCY1) gene variation in sixteen achiote accessions divided into three groups according to contrasting traits, such as flower colour, fruit colour and bixin production. Using a combination of single-strand conformational polymorphism techniques and the sequencing of polymorphic bands, we identified several single-nucleotide polymorphisms that divided the accessions into three haplotypes. Surprisingly, we observed that these three haplotypes were consistent with the same three groups previously characterized by phenotypic traits. We derived a putative sequence for the Boβ-LCY1 gene and surveyed the variations in this sequence. The heterozygosity of Boβ-LCY1 alleles resulted in a higher bixin content, likely associated with heterosis for this metabolite. These findings augment the toolbox available for the selection and genetic improvement of B. orellana and provide a reliable phenotype–genotype association method for commercial varietal selection, contributing to the development of laboratory techniques to identify desirable traits of commercial plant species.

BackgroundBixin or annatto is a commercially important natural orange-red pigment derived from lycopene that is produced and stored in seeds of Bixa orellana L. An enzymatic pathway for bixin biosynthesis was inferred from homology of putative proteins encoded by differentially expressed seed cDNAs. Some activities were later validated in a heterologous system. Nevertheless, much of the pathway remains to be clarified. For example, it is essential to identify the methylerythritol phosphate (MEP) and carotenoid pathways genes.ResultsIn order to investigate the MEP, carotenoid, and bixin pathways genes, total RNA from young leaves and two different developmental stages of seeds from B. orellana were used for the construction of indexed mRNA libraries, sequenced on the Illumina HiSeq 2500 platform and assembled de novo using Velvet, CLC Genomics Workbench and CAP3 software. A total of 52,549 contigs were obtained with average length of 1,924 bp. Two phylogenetic analyses of inferred proteins, in one case encoded by thirteen general, single-copy cDNAs, in the other from carotenoid and MEP cDNAs, indicated that B. orellana is closely related to sister Malvales species cacao and cotton. Using homology, we identified 7 and 14 core gene products from the MEP and carotenoid pathways, respectively. Surprisingly, previously defined bixin pathway cDNAs were not present in our transcriptome. Here we propose a new set of gene products involved in bixin pathway.ConclusionThe identification and qRT-PCR quantification of cDNAs involved in annatto production suggest a hypothetical model for bixin biosynthesis that involve coordinated activation of some MEP, carotenoid and bixin pathway genes. These findings provide a better understanding of the mechanisms regulating these pathways and will facilitate the genetic improvement of B. orellana.

The centromeric histone 3 variant (CENH3, aka CENP-A) is essential for the segregation of sister chromatids during mitosis and meiosis. To better define CENH3 functional constraints, we complemented a null allele in Arabidopsis with a variety of mutant alleles, each inducing a single amino acid change in conserved residues of the histone fold domain. Many of these transgenic missense lines displayed wild-type growth and fertility on self-pollination, but exhibited frequent post-zygotic death and uniparental inheritance when crossed with wild-type plants. The failure of centromeres marked by these missense mutation in the histone fold domain of CENH3 reproduces the genome elimination syndromes described with chimeric CENH3 and CENH3 from diverged species. Additionally, evidence that a single point mutation is sufficient to generate a haploid inducer provide a simple one-step method for the identification of non-transgenic haploid inducers in existing mutagenized collections of crop species. As proof of the extreme simplicity of this approach to create haploid-inducing lines, we performed an in silico search for previously identified point mutations in CENH3 and identified an Arabidopsis line carrying the A86V substitution within the histone fold domain. This A87V non-transgenic line, while fully fertile on self-pollination, produced postzygotic death and uniparental haploids when crossed to wild type.

Genome instability is associated with mitotic errors and cancer. This phenomenon can lead to deleterious rearrangements, but also genetic novelty, and many questions regarding its genesis, fate and evolutionary role remain unanswered. Here, we describe extreme chromosomal restructuring during genome elimination, a process resulting from hybridization of Arabidopsis plants expressing different centromere histones H3. Shattered chromosomes are formed from the genome of the haploid inducer, consistent with genomic catastrophes affecting a single, laggard chromosome compartmentalized within a micronucleus. Analysis of breakpoint junctions implicates breaks followed by repair through non-homologous end joining (NHEJ) or stalled fork repair. Furthermore, mutation of required NHEJ factor DNA Ligase 4 results in enhanced haploid recovery. Lastly, heritability and stability of a rearranged chromosome suggest a potential for enduring genomic novelty. These findings provide a tractable, natural system towards investigating the causes and mechanisms of complex genomic rearrangements similar to those associated with several human disorders.

The point of attachment of spindle microtubules to metaphase chromosomes is known as the centromere. Plant and animal centromeres are epigenetically specified by a centromere-specific variant of Histone H3, CENH3 (a.k.a. CENP-A). Unlike canonical histones that are invariant, CENH3 proteins are accumulating substitutions at an accelerated rate. This diversification of CENH3 is a conundrum since its role as the key determinant of centromere identity remains a constant across species. Here, we ask whether naturally occurring divergence in CENH3 has functional consequences. We performed functional complementation assays on cenh3-1, a null mutation in Arabidopsis thaliana, using untagged CENH3s from increasingly distant relatives. Contrary to previous results using GFP-tagged CENH3, we find that the essential functions of CENH3 are conserved across a broad evolutionary landscape. CENH3 from a species as distant as the monocot Zea mays can functionally replace A. thaliana CENH3. Plants expressing variant CENH3s that are fertile when selfed show dramatic segregation errors when crossed to a wild-type individual. The progeny of this cross include hybrid diploids, aneuploids with novel genetic rearrangements and haploids that inherit only the genome of the wild-type parent. Importantly, it is always chromosomes from the plant expressing the divergent CENH3 that missegregate. Using chimeras, we show that it is divergence in the fast-evolving N-terminal tail of CENH3 that is causing segregation errors and genome elimination. Furthermore, we analyzed N-terminal tail sequences from plant CENH3s and discovered a modular pattern of sequence conservation. From this we hypothesize that while the essential functions of CENH3 are largely conserved, the N-terminal tail is evolving to adapt to lineage-specific centromeric constraints. Our results demonstrate that this lineage-specific evolution of CENH3 causes inviability and sterility of progeny in crosses, at the same time producing karyotypic variation. Thus, CENH3 evolution can contribute to postzygotic reproductive barriers.

Interspecific hybridization often induces epigenetic remodeling that leads to transposon activation, gene expression changes, and loss of imprinting. These genomic changes can be deleterious and contribute to postzygotic hybrid incompatibility. In Arabidopsis, loss of genomic imprinting of PHERES1 and presumed failure of Polycomb Repressive Complex contributes to seed inviability observed in A. thaliana X A. arenosa interspecific hybrids. We used this species pair to further analyze the relationship between parentally biased gene expression and postzygotic hybrid incompatibility using two A. thaliana accessions, Col-0 and C24, with differential seed survival. We found that parentally biased expression was perturbed to a similar degree in both A. thaliana hybrids for PHERES1, HDG3, and six other normally paternally expressed genes. We propose that early genome remodeling and loss of imprinting of seed development genes induces lethality in both compatible and incompatible hybrids.