Three near-complete genome assemblies reveal substantial centromere dynamics from diploid to tetraploid in Brachypodium genus

BACKGROUND

Centromeres are critical for maintaining genomic stability in eukaryotes, and their turnover shapes genome architectures and drives karyotype evolution. However, the co-evolution of centromeres from different species in allopolyploids over millions of years remains largely unknown.

RESULTS

Here, we generate three near-complete genome assemblies, a tetraploid Brachypodium hybridum and its two diploid ancestors, Brachypodium distachyon and Brachypodium stacei. We detect high degrees of sequence, structural, and epigenetic variations of centromeres at base-pair resolution between closely related Brachypodium genomes, indicating the appearance and accumulation of species-specific centromere repeats from a common origin during evolution. We also find that centromere homogenization is accompanied by local satellite repeats bursting and retrotransposon purging, and the frequency of retrotransposon invasions drives the degree of interspecies centromere diversification. We further investigate the dynamics of centromeres during alloploidization process, and find that dramatic genetics and epigenetics architecture variations are associated with the turnover of centromeres between homologous chromosomal pairs from diploid to tetraploid. Additionally, our pangenomes analysis reveals the ongoing variations of satellite repeats and stable evolutionary homeostasis within centromeres among individuals of each Brachypodium genome with different polyploidy levels.

CONCLUSIONS

Our results provide unprecedented information on the genomic, epigenomic, and functional diversity of highly repetitive DNA between closely related species and their allopolyploid genomes at both coarse and fine scale.

A happy accident: a novel turfgrass reference genome

Poa pratensis, commonly known as Kentucky bluegrass, is a popular cool-season grass species used as turf in lawns and recreation areas globally. Despite its substantial economic value, a reference genome had not previously been assembled due to the genome’s relatively large size and biological complexity that includes apomixis, polyploidy, and interspecific hybridization. We report here a fortuitous de novo assembly and annotation of a P. pratensis genome. Instead of sequencing the genome of a C4 grass, we accidentally sampled and sequenced tissue from a weedy P. pratensis whose stolon was intertwined with that of the C4 grass. The draft assembly consists of 6.09 Gbp with an N50 scaffold length of 65.1 Mbp, and a total of 118 scaffolds, generated using PacBio long reads and Bionano optical map technology. We annotated 256 K gene models and found 58% of the genome to be composed of transposable elements. To demonstrate the applicability of the reference genome, we evaluated population structure and estimated genetic diversity in P. pratensis collected from three North American prairies, two in Manitoba, Canada and one in Colorado, USA. Our results support previous studies that found high genetic diversity and population structure within the species. The reference genome and annotation will be an important resource for turfgrass breeding and study of bluegrasses.

Genomic imbalance modulates transposable element expression in maize

Genomic imbalance refers to the more severe phenotypic consequences of changing part of a chromosome compared with the whole genome set. Previous genome imbalance studies in maize have identified prevalent inverse modulation of genes on the unvaried chromosomes (trans) with both the addition or subtraction of chromosome arms. Transposable elements (TEs) comprise a substantial fraction of the genome, and their reaction to genomic imbalance is therefore of interest. Here, we analyzed TE expression using RNA-seq data of aneuploidy and ploidy series and found that most aneuploidies showed an inverse modulation of TEs, but reductions in monosomy and increases in disomy and trisomy were also common. By contrast, the ploidy series showed little TE modulation. The modulation of TEs and genes in the same experimental group were compared, and TEs showed greater modulation than genes, especially in disomy. Class I and II TEs were differentially modulated in most aneuploidies, and some superfamilies in each TE class also showed differential modulation. Finally, the significantly upregulated TEs in three disomies (TB-7Lb, TB9Lc, and TB-10L19) did not increase the proportion of adjacent gene expression when compared with non-differentially expressed TEs, indicating that modulations of TEs do not compound the effect on genes. These results suggest that the prevalent inverse TE modulation in aneuploidy results from stoichiometric upset of the regulatory machinery used by TEs, similar to the response of core genes to genomic imbalance.

MicroRNAs play regulatory roles in genomic balance

Classic genetics studies found that genomic imbalance caused by changing the dosage of part of the genome (aneuploidy) has more detrimental effects than altering the dosage of the whole genome (ploidy). Previous analysis revealed global modulation of gene expression triggered by aneuploidy across various species, including maize (Zea mays), Arabidopsis, yeast, mammals, etc. Plant microRNAs (miRNAs) are a class of 20- to 24-nt endogenous small noncoding RNAs that carry out post-transcriptional gene expression regulation. That miRNAs and their putative targets are preferentially retained as duplicates after whole-genome duplication, as are many transcription factors and signaling components, indicates miRNAs are likely to be dosage-sensitive and potentially involved in genomic balance networks. This review addresses the following questions regarding the role of miRNAs in genomic imbalance. (1) How do aneuploidy and polyploidy impact the expression of miRNAs? (2) Do miRNAs play a regulatory role in modulating the expression of their targets under genomic imbalance?

Nitrous Oxide-Induced Metaphase Arrest: A Technique for Somatic Chromosome Analysis

Procedures to arrest metaphase chromosomes are used for determining chromosome numbers, chromosomal aberrations, and natural chromosome variation, as well as chromosome sorting. Here is described a technique of nitrous oxide gas treatment of freshly harvested root tips that is highly effective at producing an excellent mitotic index together with well-spread chromosomes. The details of the treatment and equipment used are provided. The metaphase spreads can be used directly for determining chromosome numbers or for in situ hybridization to reveal chromosomal features.

Prospects of Feral Crop De Novo Redomestication

Modern agriculture depends on a narrow variety of crop species, leaving global food and nutritional security highly vulnerable to the adverse effects of climate change and population expansion. Crop improvement using conventional and molecular breeding approaches leveraging plant genetic diversity using crop wild relatives (CWRs) has been one approach to address these issues. However, the rapid pace of the global change requires additional innovative solutions to adapt agriculture to meet global needs. Neodomestication-the rapid and targeted introduction of domestication traits using introgression or genome editing of CWRs-is being explored as a supplementary approach. These methods show promise; however, they have so far been limited in efficiency and applicability. We propose expanding the scope of neodomestication beyond truly wild CWRs to include feral crops as a source of genetic diversity for novel crop development, in this case 'redomestication'. Feral crops are plants that have escaped cultivation and evolved independently, typically adapting to their local environments. Thus, feral crops potentially contain valuable adaptive features while retaining some domestication traits. Due to their genetic proximity to crop species, feral crops may be easier targets for de novo domestication (i.e. neodomestication via genome editing techniques). In this review, we explore the potential of de novo redomestication as an application for novel crop development by genome editing of feral crops. This approach to efficiently exploit plant genetic diversity would access an underutilized reservoir of genetic diversity that could prove important in support of global food insecurity in the face of the climate change.

Centromeres: From chromosome biology to biotechnology applications and synthetic genomes in plants

Centromeres are the genomic regions that organize and regulate chromosome behaviours during cell cycle, and their variations are associated with genome instability, karyotype evolution and speciation in eukaryotes. The highly repetitive and epigenetic nature of centromeres were documented during the past half century. With the aid of rapid expansion in genomic biotechnology tools, the complete sequence and structural organization of several plant and human centromeres were revealed recently. Here, we systematically summarize the current knowledge of centromere biology with regard to the DNA compositions and the histone H3 variant (CENH3)-dependent centromere establishment and identity. We discuss the roles of centromere to ensure cell division and to maintain the three-dimensional (3D) genomic architecture in different species. We further highlight the potential applications of manipulating centromeres to generate haploids or to induce polyploids offspring in plant for breeding programs, and of targeting centromeres with CRISPR/Cas for chromosome engineering and speciation. Finally, we also assess the challenges and strategies for de novo design and synthesis of centromeres in plant artificial chromosomes. The biotechnology applications of plant centromeres will be of great potential for the genetic improvement of crops and precise synthetic breeding in the future.

The multiple fates of gene duplications: Deletion, hypofunctionalization, subfunctionalization, neofunctionalization, dosage balance constraints, and neutral variation

Gene duplications have long been recognized as a contributor to the evolution of genes with new functions. Multiple copies of genes can result from tandem duplication, from transposition to new chromosomes, or from whole-genome duplication (polyploidy). The most common fate is that one member of the pair is deleted to return the gene to the singleton state. Other paths involve the reduced expression of both copies (hypofunctionalization) that are held in duplicate to maintain sufficient quantity of function. The two copies can split functions (subfunctionalization) or can diverge to generate a new function (neofunctionalization). Retention of duplicates resulting from doubling of the whole genome occurs for genes involved with multicomponent interactions such as transcription factors and signal transduction components. In contrast, these classes of genes are underrepresented in small segmental duplications. This complementary pattern suggests that the balance of interactors affects the fate of the duplicate pair. We discuss the different mechanisms that maintain duplicated genes, which may change over time and intersect.

Dosage-sensitive miRNAs trigger modulation of gene expression during genomic imbalance in maize

The genomic imbalance caused by varying the dosage of individual chromosomes or chromosomal segments (aneuploidy) has more detrimental effects than altering the dosage of complete chromosome sets (ploidy). Previous analysis of maize (Zea mays) aneuploids revealed global modulation of gene expression both on the varied chromosome (cis) and the remainder of the genome (trans). However, little is known regarding the role of microRNAs (miRNAs) under genomic imbalance. Here, we report the impact of aneuploidy and polyploidy on the expression of miRNAs. In general, cis miRNAs in aneuploids present a predominant gene-dosage effect, whereas trans miRNAs trend toward the inverse level, although other types of responses including dosage compensation, increased effect, and decreased effect also occur. By contrast, polyploids show less differential miRNA expression than aneuploids. Significant correlations between expression levels of miRNAs and their targets are identified in aneuploids, indicating the regulatory role of miRNAs on gene expression triggered by genomic imbalance.

The non-Mendelian behavior of plant B chromosomes

B chromosomes, also known as supernumerary chromosomes, are dispensable elements in the genome of many plants, animals, and fungi. Many B chromosomes have evolved one or more drive mechanisms to transmit themselves at a higher frequency than predicted by Mendelian genetics, and these mechanisms counteract the tendency of non-essential genetic elements to be lost over time. The frequency of Bs in a population results from a balance between their effect on host fitness and their transmission rate. Here, we will summarize the findings of the drive process of plant B chromosomes, focusing on maize and rye.

Effect of aneuploidy of a non-essential chromosome on gene expression in maize

The non-essential supernumerary maize (Zea mays) B chromosome (B) has recently been shown to contain active genes and to be capable of impacting gene expression of the A chromosomes. However, the effect of the B chromosome on gene expression is still unclear. In addition, it is unknown whether the accumulation of the B chromosome has a cumulative effect on gene expression. To examine these questions, the global expression of genes, microRNAs (miRNAs), and transposable elements (TEs) of leaf tissue of maize W22 plants with 0-7 copies of the B chromosome was studied. All experimental genotypes with B chromosomes displayed a trend of upregulated gene expression for a subset of A-located genes compared to the control. Over 3000 A-located genes are significantly differentially expressed in all experimental genotypes with the B chromosome relative to the control. Modulations of these genes are largely determined by the presence rather than the copy number of the B chromosome. By contrast, the expression of most B-located genes is positively correlated with B copy number, showing a proportional gene dosage effect. The B chromosome also causes increased expression of A-located miRNAs. Differentially expressed miRNAs potentially regulate their targets in a cascade of effects. Furthermore, the varied copy number of the B chromosome leads to the differential expression of A-located and B-located TEs. The findings provide novel insights into the function and properties of the B chromosome.

Prospects and challenges of epigenomics in crop improvement

BACKGROUND

The advent of high-throughput epigenome mapping techniques has ushered in a new era of multiomics with powerful tools now available to map and record genomic output at different levels. Integrating the different components of the epigenome from these multiomics measures allows investigations of cis-regulatory elements on a genome-scale. Mapping of chromatin state, chromatin accessibility dynamics, and higher-order chromatin structure enables a new level of understanding of cell fate determination, identity and function in normal growth and development, disease resistance, and yield.

OBJECTIVE

In this paper, the recent advances in epigenomics research of rice, maize, and wheat are reviewed, and the development trends of epigenomics of major crops in the coming years are projected.

METHODS

We highlight the role of epigenomics in regulating growth and development and identifying potential distal cis-regulatory elements in three major crops, and discuss the prospects and challenges for new epigenetics-mediated breeding technologies in crop improvement.

CONCLUSION

In this review, we summarize and analyze recent epigenomic advances in three major crops epigenomics and discuss possibilities and challenges for future research in the field.

One Hundred Years of Gene Balance: How Stoichiometric Issues Affect Gene Expression, Genome Evolution, and Quantitative Traits

A century ago experiments with the flowering plant Datura stramonium and the fruit fly Drosophila melanogaster revealed that adding an extra chromosome to a karyotype was much more detrimental than adding a whole set of chromosomes. This phenomenon was referred to as gene balance and has been recapitulated across eukaryotic species. Here, we retrace some developments in this field. Molecular studies suggest that the basis of balance involves stoichiometric relationships of multi-component interactions. This concept has implication for the mechanisms controlling gene expression, genome evolution, sex chromosome evolution/dosage compensation, speciation mechanisms, and the underlying genetics of quantitative traits.

The supernumerary B chromosome of maize: drive and genomic conflict

The supernumerary B chromosome of maize is dispensable, containing no vital genes, and thus is variable in number and presence in lines of maize. In order to be maintained in populations, it has a drive mechanism consisting of nondisjunction at the pollen mitosis that produces the two sperm cells, and then the sperm with the two B chromosomes has a preference for fertilizing the egg as opposed to the central cell in the process of double fertilization. The sequence of the B chromosome coupled with B chromosomal aberrations has localized features involved with nondisjunction and preferential fertilization, which are present at the centromeric region. The predicted genes from the sequence have paralogues dispersed across all A chromosomes and have widely different divergence times suggesting that they have transposed to the B chromosome over evolutionary time followed by degradation or have been co-opted for the selfish functions of the supernumerary chromosome.

De novo centromere formation on chromosome fragments with an inactive centromere in maize (Zea mays)

The B chromosome of maize undergoes nondisjunction at the second pollen mitosis as part of its accumulation mechanism. Previous work identified 9-Bic-1 (9-B inactivated centromere-1), which comprises an epigenetically silenced B chromosome centromere that was translocated to the short arm of chromosome 9(9S). This chromosome is stable in isolation, but when normal B chromosomes are added to the genotype, it will attempt to undergo nondisjunction during the second pollen mitosis and usually fractures the chromosome in 9S. These broken chromosomes allow a test of whether the inactive centromere is reactivated or whether a de novo centromere is formed elsewhere on the chromosome to allow recovery of fragments. Breakpoint determination on the B chromosome and chromosome 9 showed that mini chromosome B1104 has the same breakpoint as 9-Bic-1 in the B centromere region and includes a portion of 9S. CENH3 binding was found on the B centromere region and on 9S, suggesting both centromere reactivation and de novo centromere formation. Another mini chromosome, B496, showed evidence of rearrangement, but it also only showed evidence for a de novo centromere. Other mini chromosome fragments recovered were directly derived from the B chromosome with breakpoints concentrated near the centromeric knob region, which suggests that the B chromosome is broken at a low frequency due to the failure of the sister chromatids to separate at the second pollen mitosis. Our results indicate that both reactivation and de novo centromere formation could occur on fragments derived from the progenitor possessing an inactive centromere.

GNET2: an R package for constructing gene regulatory networks from transcriptomic data

MOTIVATION

The Gene Network Estimation Tool (GNET) is designed to build gene regulatory networks (GRNs) from transcriptomic gene expression data with a probabilistic graphical model. The data preprocessing, model construction and visualization modules of the original GNET software were developed on different programming platforms, which were inconvenient for users to deploy and use.

RESULTS

Here, we present GNET2, an improved implementation of GNET as an integrated R package. GNET2 provides more flexibility for parameter initialization and regulatory module construction based on the core iterative modeling process of the original algorithm. The data exchange interface of GNET2 is handled within an R session automatically. Given the growing demand for regulatory network reconstruction from transcriptomic data, GNET2 offers a convenient option for GRN inference on large datasets.

AVAILABILITY AND IMPLEMENTATION

The source code of GNET2 is available at https://github.com/jianlin-cheng/GNET2.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Genome-wide mapping reveals R-loops associated with centromeric repeats in maize

R-loops are stable chromatin structures comprising a DNA:RNA hybrid and a displaced single-stranded DNA. R-loops have been implicated in gene expression and chromatin structure, as well as in replication blocks and genome instability. Here, we conducted a genome-wide identification of R-loops and identified more than 700,000 R-loop peaks in the maize (Zea mays) genome. We found that sense R-loops were mainly enriched in promoters and transcription termination sites and relatively less enriched in gene bodies, which is different from the main gene-body localization of sense R-loops in Arabidopsis and Oryza sativa. At the chromosome scale, maize R-loops were enriched in pericentromeric heterochromatin regions, and a significant portion of R-loops were derived from transposable elements. In centromeres, R-loops preferentially formed within the binding regions of the centromere-specific histone CENH3, and centromeric retrotransposons were strongly associated with R-loop formation. Furthermore, centromeric retrotransposon R-loops were observed by applying the single-molecule imaging technique of atomic force microscopy. These findings elucidate the fundamental character of R-loops in the maize genome and reveal the potential role of R-loops in centromeres.

Gene-dosage issues: a recurrent theme in whole genome duplication events

Two recent studies have addressed the long-term consequences of whole genome duplications (WGD). Specifically, they analyzed transcriptomes of the plant Arabidopsis thaliana and of four salmonids to assess the impact of WGD on gene expression. These studies point to commonalities in gene expression adjustments after polyploidization that we outline and discuss below.

Genomic imbalance determines positive and negative modulation of gene expression in diploid maize

Genomic imbalance caused by changing the dosage of individual chromosomes (aneuploidy) has a more detrimental effect than varying the dosage of complete sets of chromosomes (ploidy). We examined the impact of both increased and decreased dosage of 15 distal and 1 interstitial chromosomal regions via RNA-seq of maize (Zea mays) mature leaf tissue to reveal new aspects of genomic imbalance. The results indicate that significant changes in gene expression in aneuploids occur both on the varied chromosome (cis) and the remainder of the genome (trans), with a wider spread of modulation compared with the whole-ploidy series of haploid to tetraploid. In general, cis genes in aneuploids range from a gene-dosage effect to dosage compensation, whereas for trans genes the most common effect is an inverse correlation in that expression is modulated toward the opposite direction of the varied chromosomal dosage, although positive modulations also occur. Furthermore, this analysis revealed the existence of increased and decreased effects in which the expression of many genes under genome imbalance are modulated toward the same direction regardless of increased or decreased chromosomal dosage, which is predicted from kinetic considerations of multicomponent molecular interactions. The findings provide novel insights into understanding mechanistic aspects of gene regulation.

Predominantly inverse modulation of gene expression in genomically unbalanced disomic haploid maize

The phenotypic consequences of the addition or subtraction of part of a chromosome is more severe than changing the dosage of the whole genome. By crossing diploid trisomies to a haploid inducer, we identified 17 distal segmental haploid disomies that cover ∼80% of the maize genome. Disomic haploids provide a level of genomic imbalance that is not ordinarily achievable in multicellular eukaryotes, allowing the impact to be stronger and more easily studied. Transcriptome size estimates revealed that a few disomies inversely modulate most of the transcriptome. Based on RNA sequencing, the expression levels of genes located on the varied chromosome arms (cis) in disomies ranged from being proportional to chromosomal dosage (dosage effect) to showing dosage compensation with no expression change with dosage. For genes not located on the varied chromosome arm (trans), an obvious trans-acting effect can be observed, with the majority showing a decreased modulation (inverse effect). The extent of dosage compensation of varied cis genes correlates with the extent of trans inverse effects across the 17 genomic regions studied. The results also have implications for the role of stoichiometry in gene expression, the control of quantitative traits, and the evolution of dosage-sensitive genes.

Preferential meiotic chromosome pairing among homologous chromosomes with cryptic sequence variation in tetraploid maize

Meiotic chromosome pairing between homoeologous chromosomes was reported in many nascent allopolyploids. Homoeologous pairing is gradually eliminated and replaced by exclusive homologous pairing in well-established allopolyploids, an evolutionary process referred to as the diploidization of allopolyploids. A fundamental question of the diploidization of allopolyploids is whether and to what extent the DNA sequence variation among homoeologous chromosomes contribute to the establishment of exclusive homologous chromosome pairing. We developed aneuploid tetraploid maize lines that contain three copies of chromosome 10 derived from inbred lines B73 and H99. We were able to identify the parental origin of each copy of chromosome 10 in the materials using oligonucleotide-based haplotype-specific chromosome painting. We demonstrate that the two identical copies of chromosome 10 from H99 pair preferentially over chromosome 10 from B73 in different stages of prophase I and metaphase I during meiosis. Thus, homologous chromosome pairing is favored to partners with the most similar DNA sequences and can be discriminated based on cryptic sequence variation. We propose that innate preference of homologous chromosome pairing exists in nascent allopolyploids and serves as the first layer that would eventually block all homoeologous chromosome pairing in allopolyploids.

Phosphorylation of histone H3 by Haspin regulates chromosome alignment and segregation during mitosis in maize

In human cells, Haspin-mediated histone H3 threonine 3 (H3T3) phosphorylation promotes centromeric localization of the chromosomal passenger complex, thereby ensuring proper kinetochore-microtubule attachment. Haspin also binds to PDS5 cohesin-associated factor B (Pds5B), antagonizing the Wings apart-like protein homolog (Wapl)-Pds5B interaction and thus preventing Wapl from releasing centromeric cohesion during mitosis. However, the role of Haspin in plant chromosome segregation is not well understood. Here, we show that in maize (Zea mays) mitotic cells, ZmHaspin localized to the centromere during metaphase and anaphase, whereas it localized to the telomeres during meiosis. These results suggest that ZmHaspin plays different roles during mitosis and meiosis. Knockout of ZmHaspin led to decreased H3T3 phosphorylation and histone H3 serine 10 phosphorylation, and defects in chromosome alignment and segregation in mitosis. These lines of evidence suggest that Haspin regulates chromosome segregation in plants via the mechanism described for humans, namely, H3T3 phosphorylation. Plant Haspin proteins lack the RTYGA and PxVxL motifs needed to bind Pds5B and heterochromatin protein 1, and no obvious cohesion defects were detected in ZmHaspin knockout plants. Taken together, these results highlight the conserved but slightly different roles of Haspin proteins in cell division in plants and in animals.

An empirical bayesian approach for testing gene expression fold change and its application in detecting global dosage effects

We are motivated by biological studies intended to understand global gene expression fold change. Biologists have generally adopted a fixed cutoff to determine the significance of fold changes in gene expression studies (e.g. by using an observed fold change equal to two as a fixed threshold). Scientists can also use a t-test or a modified differential expression test to assess the significance of fold changes. However, these methods either fail to take advantage of the high dimensionality of gene expression data or fail to test fold change directly. Our research develops a new empirical Bayesian approach to substantially improve the power and accuracy of fold-change detection. Specifically, we more accurately estimate gene-wise error variation in the log of fold change. We then adopt a t-test with adjusted degrees of freedom for significance assessment. We apply our method to a dosage study in Arabidopsis and a Down syndrome study in humans to illustrate the utility of our approach. We also present a simulation study based on real datasets to demonstrate the accuracy of our method relative to error variance estimation and power in fold-change detection. Our developed R package with a detailed user manual is publicly available on GitHub at https://github.com/cuiyingbeicheng/Foldseq.

DeepGRN: prediction of transcription factor binding site across cell-types using attention-based deep neural networks

Background

Due to the complexity of the biological systems, the prediction of the potential DNA binding sites for transcription factors remains a difficult problem in computational biology. Genomic DNA sequences and experimental results from parallel sequencing provide available information about the affinity and accessibility of genome and are commonly used features in binding sites prediction. The attention mechanism in deep learning has shown its capability to learn long-range dependencies from sequential data, such as sentences and voices. Until now, no study has applied this approach in binding site inference from massively parallel sequencing data. The successful applications of attention mechanism in similar input contexts motivate us to build and test new methods that can accurately determine the binding sites of transcription factors.

Results

In this study, we propose a novel tool (named DeepGRN) for transcription factors binding site prediction based on the combination of two components: single attention module and pairwise attention module. The performance of our methods is evaluated on the ENCODE-DREAM in vivo Transcription Factor Binding Site Prediction Challenge datasets. The results show that DeepGRN achieves higher unified scores in 6 of 13 targets than any of the top four methods in the DREAM challenge. We also demonstrate that the attention weights learned by the model are correlated with potential informative inputs, such as DNase-Seq coverage and motifs, which provide possible explanations for the predictive improvements in DeepGRN.

Conclusions

DeepGRN can automatically and effectively predict transcription factor binding sites from DNA sequences and DNase-Seq coverage. Furthermore, the visualization techniques we developed for the attention modules help to interpret how critical patterns from different types of input features are recognized by our model.

Development of a Transformable Fast-Flowering Mini-Maize as a Tool for Maize Gene Editing

Maize (Zea mays ssp. mays) is a popular genetic model due to its ease of crossing, well-established toolkits, and its status as a major global food crop. Recent technology developments for precise manipulation of the genome are further impacting both basic biological research and biotechnological application in agriculture. Crop gene editing often requires a process of genetic transformation in which the editing reagents are introduced into plant cells. In maize, this procedure is well-established for a limited number of public lines that are amenable for genetic transformation. Fast-Flowering Mini-Maize (FFMM) lines A and B were recently developed as an open-source tool for maize research by reducing the space requirements and the generation time. Neither line of FFMM were competent for genetic transformation using traditional protocols, a necessity to its status as a complete toolkit for public maize genetic research. Here we report the development of new lines of FFMM that have been bred for amenability to genetic transformation. By hybridizing a transformable maize genotype high Type-II callus parent A (Hi-II A) with line A of FFMM, we introgressed the ability to form embryogenic callus from Hi-II A into the FFMM-A genetic background. Through multiple generations of iterative self-hybridization or doubled-haploid method, we established maize lines that have a strong ability to produce embryogenic callus from immature embryos and maintain resemblance to FFMM-A in flowering time and stature. Using an Agrobacterium-mediated standard transformation method, we successfully introduced the CRISPR-Cas9 reagents into immature embryos and generated transgenic and mutant lines displaying the expected mutant phenotypes and genotypes. The transformation frequencies of the tested genotypes, defined as the numbers of transgenic event producing T1 seeds per 100 infected embryos, ranged from 0 to 17.1%. Approximately 80% of transgenic plants analyzed in this study showed various mutation patterns at the target site. The transformable FFMM line, FFMM-AT, can serve as a useful genetic and genomic resource for the maize community.

A transposon surveillance mechanism that safeguards plant male fertility during stress

Although plants are able to withstand a range of environmental conditions, spikes in ambient temperature can impact plant fertility causing reductions in seed yield and notable economic losses^1,2. Therefore, understanding the precise molecular mechanisms that underpin plant fertility under environmental constraints is critical to safeguarding future food production³. Here, we identified two Argonaute-like proteins whose activities are required to sustain male fertility in maize plants under high temperatures. We found that MALE-ASSOCIATED ARGONAUTE-1 and -2 associate with temperature-induced phased secondary small RNAs in pre-meiotic anthers and are essential to controlling the activity of retrotransposons in male meiocyte initials. Biochemical and structural analyses revealed how male-associated Argonaute activity and its interaction with retrotransposon RNA targets is modulated through the dynamic phosphorylation of a set of highly conserved, surface-located serine residues. Our results demonstrate that an Argonaute-dependent, RNA-guided surveillance mechanism is critical in plants to sustain male fertility under environmentally constrained conditions, by controlling the mutagenic activity of transposons in male germ cells.

The Gene Balance Hypothesis: Epigenetics and Dosage Effects in Plants

Dosage effects in plants are caused by changes in the copy number of chromosomes, segments of chromosomes, or multiples of individual genes. Genes often exhibit a dosage effect in which the amount of product is closely correlated with the number of copies present. However, when larger segments of chromosomes are varied, there are trans-acting effects across the genome that are unleashed that modulate gene expression in cascading effects. These appear to be mediated by the stoichiometric relationship of gene regulatory machineries. There are both positive and negative modulations of target gene expression, but the latter is the plurality effect. When this inverse effect is combined with a dosage effect, compensation for a gene can occur in which its expression is similar to the normal diploid regardless of the change in chromosomal dosage. In contrast, changing the whole genome in a polyploidy series has fewer relative effects as the stoichiometric relationship is not disrupted. Together, these observations suggest that the stoichiometry of gene regulation is important as a reflection of the mode of assembly of the individual subunits involved in the effective regulatory macromolecular complexes. This principle has implications for gene expression mechanisms, quantitative trait genetics, and the evolution of genes depending on the mode of duplication, either segmentally or via whole-genome duplication.

Inbreeding Depression in Genotypically Matched Diploid and Tetraploid Maize

The genetic and molecular basis of heterosis has long been studied but without a consensus about mechanism. The opposite effect, inbreeding depression, results from repeated self-pollination and leads to a reduction in vigor. A popular explanation for this reaction is the homozygosis of recessive, slightly deleterious alleles upon inbreeding. However, extensive studies in alfalfa indicated that inbreeding between diploids and autotetraploids was similar despite the fact that homozygosis of alleles would be dramatically different. The availability of tetraploid lines of maize generated directly from various inbred lines provided the opportunity to examine this issue in detail in perfectly matched diploid and tetraploid hybrids and their parallel inbreeding regimes. Identical hybrids at the diploid and tetraploid levels were inbred in triplicate for seven generations. At the conclusion of this regime, F1 hybrids and selected representative generations (S1, S3, S5, S7) were characterized phenotypically in randomized blocks during the same field conditions. Quantitative measures of the multiple generations of inbreeding provided little evidence for a distinction in the decline of vigor between the diploids and the tetraploids. The results suggest that the homozygosis of completely recessive, slightly deleterious alleles is an inadequate hypothesis to explain inbreeding depression in general.

Rapid Birth or Death of Centromeres on Fragmented Chromosomes in Maize

Comparative genomics has revealed common occurrences in karyotype evolution such as chromosomal end-to-end fusions and insertions of one chromosome into another near the centromere, as well as many cases of de novo centromeres that generate positional polymorphisms. However, how rearrangements such as dicentrics and acentrics persist without being destroyed or lost remains unclear. Here, we sought experimental evidence for the frequency and timeframe for inactivation and de novo formation of centromeres in maize (Zea mays). The pollen from plants with supernumerary B chromosomes was gamma-irradiated and then applied to normal maize silks of a line without B chromosomes. In ∼8,000 first-generation seedlings, we found many B-A translocations, centromere expansions, and ring chromosomes. We also found many dicentric chromosomes, but a fraction of these show only a single primary constriction, which suggests inactivation of one centromere. Chromosomal fragments were found without canonical centromere sequences, revealing de novo centromere formation over unique sequences; these were validated by immunolocalization with Thr133-phosphorylated histone H2A, a marker of active centromeres, and chromatin immunoprecipitation-sequencing with the CENH3 antibody. These results illustrate the regular occurrence of centromere birth and death after chromosomal rearrangement during a narrow window of one to potentially only a few cell cycles for the rearranged chromosomes to be recognized in this experimental regime.

The deposition of CENH3 in maize is stringently regulated

The centromere, as an essential element to mediate chromosome segregation, is epigenetically determined by CENH3-containing nucleosomes as a functional marker; therefore the accurate deposition of CENH3 is crucial for chromosome transmission. We characterized the deposition of CENH3 in maize by over-expression and mutational analysis. Our results revealed that over-expressing CENH3 in callus is lethal while over-expressing GFP-CENH3 and CENH3-YFP in callus and plants is not and can be partly deposited normally. Different mutations of GFP-CENH3 demonstrated that CENH3-Thr4 in the N-terminus was needed for the deposition as a positive phosphorylation site and the last five amino acids in the C-terminus are necessary for deposition. The C-terminal tail of CENH3 is confirmed to be responsible for the interaction of CENH3 and histone H4, which indicates that CENH3 maintains deposition in centromeres via interacting with H4 to form stable nucleosomes. For GFP-CENH3 and CENH3-YFP, the fused tags at the termini probably affect the structure of CENH3 and reduce its interaction with other proteins, which in turn could decrease proper deposition. Taken together, multiple amino acids or motifs were shown to play essential roles in CENH3 deposition, which is suggested to be affected by numerous factors in maize.

Site-specific recombinase genome engineering toolkit in maize

Site-specific recombinase enzymes function in heterologous cellular environments to initiate strand-switching reactions between unique DNA sequences termed recombinase binding sites. Depending on binding site position and orientation, reactions result in integrations, excisions, or inversions of targeted DNA sequences in a precise and predictable manner. Here, we established five different stable recombinase expression lines in maize through Agrobacterium-mediated transformation of T-DNA molecules that contain coding sequences for Cre, R, FLPe, phiC31 Integrase, and phiC31 excisionase. Through the bombardment of recombinase activated DsRed transient expression constructs, we have determined that all five recombinases are functional in maize plants. These recombinase expression lines could be utilized for a variety of genetic engineering applications, including selectable marker removal, targeted transgene integration into predetermined locations, and gene stacking.

Magnitude of modulation of gene expression in aneuploid maize depends on the extent of genomic imbalance

Aneuploidy has profound effects on an organism, typically more so than polyploidy, and the basis of this contrast is not fully understood. A dosage series of the maize long arm of chromosome 1 (1L) was used to compare relative global gene expression in different types and degrees of aneuploidy to gain insights into how the magnitude of genomic imbalance as well as hypoploidy affects global gene expression. While previously available methods require a selective examination of specific genes, RNA sequencing provides a whole-genome view of gene expression in aneuploids. Most studies of global aneuploidy effects have concentrated on individual types of aneuploids because multiple dose aneuploidies of the same genomic region are difficult to produce in most model genetic organisms. The genetic toolkit of maize allows the examination of multiple ploidies and 1-4 doses of chromosome arms. Thus, a detailed examination of expression changes both on the varied chromosome arms and elsewhere in the genome is possible, in both hypoploids and hyperploids, compared with euploid controls. Previous studies observed the inverse trans effect, in which genes not varied in DNA dosage were expressed in a negative relationship to the varied chromosomal region. This response was also the major type of changes found globally in this study. Many genes varied in dosage showed proportional expression changes, though some were seen to be partly or fully dosage compensated. It was also found that the effects of aneuploidy were progressive, with more severe aneuploids producing effects of greater magnitude.

Engineered minichromosomes in plants

Artificial chromosome platforms are described in plants. Because the function of centromeres is largely epigenetic, attempts to produce artificial chromosomes with plant centromere DNA have failed. The removal of the centromeric sequences from the cell strips off the centromeric histone that is the apparent biochemical marker of centromere activity. Thus, engineered minichromosomes have been produced by telomere mediated chromosomal truncation. The introduction of telomere repeats will cleave the chromosome at the site of insertion and attach the accompanying transgenes in the process. Such truncation events have been documented in maize, Arabidopsis, barley, rice, Brassica and wheat. Truncation of the nonvital supernumerary B chromosome of maize is a favorite target but engineered minichromosomes derived from the normal A chromosomes have also been recovered. Transmission through mitosis of small chromosomes is apparently normal but there is loss during meiosis. Potential solutions to address this issue are discussed. With procedures now well established to produce the foundation for artificial chromosomes in plants, current efforts are directed at building them up to specification using gene stacking methods and editing techniques.

Meiotic crossovers characterized by haplotype-specific chromosome painting in maize

Meiotic crossovers (COs) play a critical role in generating genetic variation and maintaining faithful segregation of homologous chromosomes during meiosis. We develop a haplotype-specific fluorescence in situ hybridization (FISH) technique that allows visualization of COs directly on metaphase chromosomes. Oligonucleotides (oligos) specific to chromosome 10 of maize inbreds B73 and Mo17, respectively, are synthesized and labeled as FISH probes. The parental and recombinant chromosome 10 in B73 x Mo17 F₁ hybrids and F₂ progenies can be unambiguously identified by haplotype-specific FISH. Analysis of 58 F₂ plants reveals lack of COs in the entire proximal half of chromosome 10. However, we detect COs located in regions very close to the centromere in recombinant inbred lines from an intermated B73 x Mo17 population, suggesting effective accumulation of COs in recombination-suppressed chromosomal regions through intermating and the potential to generate favorable allelic combinations of genes residing in these regions.

Meiotic crossovers (COs) are essential for proper chromosome segregation and generating novel combinations of alleles. Here, the authors develop haplotype-specific oligos on maize chromosome 10 for fluorescence in situ hybridization and analyze CO patterns in an intermated recombinant population derived from B73 and Mo17.

Hybrid Decay: A Transgenerational Epigenetic Decline in Vigor and Viability Triggered in Backcross Populations of Teosinte with Maize

In the course of generating populations of maize with teosinte chromosomal introgressions, an unusual sickly plant phenotype was noted in individuals from crosses with two teosinte accessions collected near Valle de Bravo, Mexico. The plants of these Bravo teosinte accessions appear phenotypically normal themselves and the F1 plants appear similar to typical maize × teosinte F1s. However, upon backcrossing to maize, the BC1 and subsequent generations display a number of detrimental characteristics including shorter stature, reduced seed set, and abnormal floral structures. This phenomenon is observed in all BC individuals and there is no chromosomal segment linked to the sickly plant phenotype in advanced backcross generations. Once the sickly phenotype appears in a lineage, normal plants are never again recovered by continued backcrossing to the normal maize parent. Whole-genome shotgun sequencing reveals a small number of genomic sequences, some with homology to transposable elements, that have increased in copy number in the backcross populations. Transcriptome analysis of seedlings, which do not have striking phenotypic abnormalities, identified segments of 18 maize genes that exhibit increased expression in sickly plants. A de novo assembly of transcripts present in plants exhibiting the sickly phenotype identified a set of 59 upregulated novel transcripts. These transcripts include some examples with sequence similarity to transposable elements and other sequences present in the recurrent maize parent (W22) genome as well as novel sequences not present in the W22 genome. Genome-wide profiles of gene expression, DNA methylation, and small RNAs are similar between sickly plants and normal controls, although a few upregulated transcripts and transposable elements are associated with altered small RNA or methylation profiles. This study documents hybrid incompatibility and genome instability triggered by the backcrossing of Bravo teosinte with maize. We name this phenomenon "hybrid decay" and present ideas on the mechanism that may underlie it.

Progressive heterosis in genetically defined tetraploid maize

Progressive heterosis, i.e., the additional hybrid vigor in double-cross tetraploid hybrids not found in their single-cross tetraploid parents, has been documented in a number of species including alfalfa, potato, and maize. In this study, four artificially induced maize tetraploids, directly derived from standard inbred lines, were crossed in pairs to create two single-cross hybrids. These hybrids were then crossed to create double-cross hybrids containing genetic material from all four original lines. Replicated field-based phenotyping of the materials over four years indicated a strong progressive heterosis phenotype in tetraploids but not in their diploid counterparts. In particular, the above ground dry weight phenotype of double-cross tetraploid hybrids was on average 34% and 56% heavier than that of the single-cross tetraploid hybrids and the double-cross diploid counterparts, respectively. Additionally, whole-genome resequencing of the original inbred lines and further analysis of these data did not show the expected spectrum of alleles to explain tetraploid progressive heterosis under the complementation of complete recessive model. These results underscore the reality of the progressive heterosis phenotype, its potential utility for increasing crop biomass production, and the need for exploring alternative hypothesis to explain it at a molecular level.

Genomic Balance and Speciation

The role of genomic balance in accumulating species hybrid incompatibilities is discussed. Aneuploidy has been shown to produce more global modulations than polyploidy with the responsible genes being transcription factors and signaling components involved in molecular complexes, illustrating a stoichiometric component to gene expression. Genomic imbalance is usually detrimental to the organism and in many cases results in lethality. Here, it is proposed that once gene flow is prevented between or within populations by various speciation initiating processes, the stoichiometric relationship of members of macromolecular complexes can change via compensatory drift with the eventual result of newly established functional balances. However, when these new relationships are brought together in interspecific hybrids, detrimental consequences will occur. We suggest that these detrimental interactions contribute to hybrid incompatibilities.

Global impacts of chromosomal imbalance on gene expression in Arabidopsis and other taxa

Changes in dosage of part of the genome (aneuploidy) have long been known to produce much more severe phenotypic consequences than changes in the number of whole genomes (ploidy). To examine the basis of these differences, global gene expression in mature leaf tissue for all five trisomies and in diploids, triploids, and tetraploids of Arabidopsis thaliana was studied. The trisomies displayed a greater spread of expression modulation than the ploidy series. In general, expression of genes on the varied chromosome ranged from compensation to dosage effect, whereas genes from the remainder of the genome ranged from no effect to reduced expression approaching the inverse level of chromosomal imbalance (2/3). Genome-wide DNA methylation was examined in each genotype and found to shift most prominently with trisomy 4 but otherwise exhibited little change, indicating that genetic imbalance is generally mechanistically unrelated to DNA methylation. Independent analysis of gene functional classes demonstrated that ribosomal, proteasomal, and gene body methylated genes were less modulated compared with all classes of genes, whereas transcription factors, signal transduction components, and organelle-targeted protein genes were more tightly inversely affected. Comparing transcription factors and their targets in the trisomies and in expression networks revealed considerable discordance, illustrating that altered regulatory stoichiometry is a major contributor to genetic imbalance. Reanalysis of published data on gene expression in disomic yeast and trisomic mouse cells detected similar stoichiometric effects across broad phylogenetic taxa, and indicated that these effects reflect normal gene regulatory processes.

High-efficiency genome editing using a dmc1 promoter-controlled CRISPR/Cas9 system in maize

Previous studies revealed that the promoters for driving both Cas9 and sgRNAs are quite important for efficient genome editing by CRISPR/Cas9 in plants. Here, we report our results of targeted genome editing using the maize dmc1 gene promoter combined with the U3 promoter for Cas9 and sgRNA, respectively. Three loci in the maize genome were selected for targeting. The T0 plants regenerated were highly efficiently edited at the target sites with homozygous or bi-allelic mutants accounting for about 66%. The mutations in T0 plants could be stably transmitted to the T1 generation, and new mutations could be generated in gametes or zygotes. Whole-genome resequencing indicated that no off-target mutations could be detected in the predicted loci with sequence similarity to the targeted site. Our results show that the dmc1 promoter-controlled (DPC) CRISPR/Cas9 system is highly efficient in maize and provide further evidence that the optimization of the promoters used for the CRISPR/Cas9 system is important for enhancing the efficiency of targeted genome editing in plants. The evolutionary conservation of the dmc1 gene suggests its potential for use in other plant species.

Meiotic Studies on Combinations of Chromosomes With Different Sized Centromeres in Maize

Multiple centromere misdivision derivatives of a translocation between the supernumerary B chromosome and the short arm of chromosome 9 (TB-9Sb) permit investigation of how centromeres of different sizes behave in meiosis in opposition or in competition with each other. In the first analysis, heterozygotes were produced between the normal TB-9Sb and derivatives of it that resulted from centromere misdivision that reduced the amounts of centromeric DNA. These heterozygotes could test whether these drastic differences would result in meiotic drive of the larger chromosome in female meiosis. Cytological determinations of the segregation of large and small centromeres among thousands of progeny of four combinations were made. The recovery of the larger centromere was at a few percent higher frequency in two of four combinations. However, examination of phosphorylated histone H2A-Thr133, a characteristic of active centromeres, showed a lack of correlation with the size of the centromeric DNA, suggesting an expansion of the basal protein features of the kinetochore in two of the three cases despite the reduction in the size of the underlying DNA. In the second analysis, plants containing different sizes of the B chromosome centromere were crossed to plants with TB-9Sb with a foldback duplication of 9S (TB-9Sb-Dp9). In the progeny, plants containing large and small versions of the B chromosome centromere were selected by FISH. A meiotic "tug of war" occurred in hybrid combinations by recombination between the normal 9S and the foldback duplication in those cases in which pairing occurred. Such pairing and recombination produce anaphase I bridges but in some cases the large and small centromeres progressed to the same pole. In one combination, new dicentric chromosomes were found in the progeny. Collectively, the results indicate that the size of the underlying DNA of a centromere does not dramatically affect its segregation properties or its ability to progress to the poles in meiosis potentially because the biochemical features of centromeres adjust to the cellular conditions.

Parallel altitudinal clines reveal trends in adaptive evolution of genome size in Zea mays

While the vast majority of genome size variation in plants is due to differences in repetitive sequence, we know little about how selection acts on repeat content in natural populations. Here we investigate parallel changes in intraspecific genome size and repeat content of domesticated maize (Zea mays) landraces and their wild relative teosinte across altitudinal gradients in Mesoamerica and South America. We combine genotyping, low coverage whole-genome sequence data, and flow cytometry to test for evidence of selection on genome size and individual repeat abundance. We find that population structure alone cannot explain the observed variation, implying that clinal patterns of genome size are maintained by natural selection. Our modeling additionally provides evidence of selection on individual heterochromatic knob repeats, likely due to their large individual contribution to genome size. To better understand the phenotypes driving selection on genome size, we conducted a growth chamber experiment using a population of highland teosinte exhibiting extensive variation in genome size. We find weak support for a positive correlation between genome size and cell size, but stronger support for a negative correlation between genome size and the rate of cell production. Reanalyzing published data of cell counts in maize shoot apical meristems, we then identify a negative correlation between cell production rate and flowering time. Together, our data suggest a model in which variation in genome size is driven by natural selection on flowering time across altitudinal clines, connecting intraspecific variation in repetitive sequence to important differences in adaptive phenotypes.

Genome size in plants can vary by orders of magnitude, but this variation has long been considered to be of little functional consequence. Studying three independent adaptations to high altitude in Zea mays, we find that genome size experiences parallel pressures from natural selection, causing a reduction in genome size with increasing altitude. Though reductions in overall repetitive content are responsible for the genome size change, we find that only those individual loci contributing most to the variation in genome size are individually targeted by selection. To identify the phenotype influenced by genome size, we study how variation in genome size within a single wild population impacts leaf growth and cell division. We find that genome size variation correlates negatively with the rate of cell division, suggesting that individuals with larger genomes require longer to complete a mitotic cycle. Finally, we reanalyze data from maize inbreds to show that faster cell division is correlated with earlier flowering, connecting observed variation in genome size to an important adaptive phenotype.

Parallel altitudinal clines reveal trends in adaptive evolution of genome size in Zea mays

While the vast majority of genome size variation in plants is due to differences in repetitive sequence, we know little about how selection acts on repeat content in natural populations. Here we investigate parallel changes in intraspecific genome size and repeat content of domesticated maize (Zea mays) landraces and their wild relative teosinte across altitudinal gradients in Mesoamerica and South America. We combine genotyping, low coverage whole-genome sequence data, and flow cytometry to test for evidence of selection on genome size and individual repeat abundance. We find that population structure alone cannot explain the observed variation, implying that clinal patterns of genome size are maintained by natural selection. Our modeling additionally provides evidence of selection on individual heterochromatic knob repeats, likely due to their large individual contribution to genome size. To better understand the phenotypes driving selection on genome size, we conducted a growth chamber experiment using a population of highland teosinte exhibiting extensive variation in genome size. We find weak support for a positive correlation between genome size and cell size, but stronger support for a negative correlation between genome size and the rate of cell production. Reanalyzing published data of cell counts in maize shoot apical meristems, we then identify a negative correlation between cell production rate and flowering time. Together, our data suggest a model in which variation in genome size is driven by natural selection on flowering time across altitudinal clines, connecting intraspecific variation in repetitive sequence to important differences in adaptive phenotypes.

Barbara McClintock's Unsolved Chromosomal Mysteries: Parallels to Common Rearrangements and Karyotype Evolution

Two obscure studies on chromosomal behavior by Barbara McClintock are revisited in light of subsequent studies and evolutionary genomics of chromosome number reduction. The phenomenon of deficiency recovery in which adjacent genetic markers lost in the zygote reappear in later developmental sectors is discussed in light of de novo centromere formation on chromosomal fragments. Second, McClintock described a small chromosome, which she postulated carried an "X component," that fostered specific types of chromosomal rearrangements mainly involving centromere changes and attachments to the termini of chromosomes. These findings are cast in the context of subsequent studies on centromere misdivision, the tendency of broken fragments to join chromosome ends, and the realization from genomic sequences that nested chromosomal insertion and end-to-end chromosomal fusions are common features of karyotype evolution. Together, these results suggest a synthesis that centromere breaks, inactivation, and de novo formation together with telomeres-acting under some circumstances as double-strand DNA breaks that join with others-is the underlying basis of these chromosomal phenomena.