Single-parent expression drives dynamic gene expression complementation in maize hybrids.

Single-parent expression (SPE) is defined as gene expression in only one of the two parents. SPE can arise from differential expression between parental alleles, termed non-presence/absence (non-PAV) SPE, or from the physical absence of a gene in one parent, termed PAV SPE. We used transcriptome data of diverse Zea mays (maize) inbreds and hybrids, including 401 samples from five different tissues, to test for differences between these types of SPE genes. Although commonly observed, SPE is highly genotype and tissue specific. A positive correlation was observed between the genetic distance of the two inbred parents and the number of SPE genes identified. Regulatory analysis showed that PAV SPE and non-PAV SPE genes are mainly regulated by cis effects, with a small fraction under trans regulation. Polymorphic transposable element insertions in promoter sequences contributed to the high level of cis regulation for PAV SPE and non-PAV SPE genes. PAV SPE genes were more frequently expressed in hybrids than non-PAV SPE genes. The expression of parentally silent alleles in hybrids of non-PAV SPE genes was relatively rare but occurred in most hybrids. Non-PAV SPE genes with expression of the silent allele in hybrids are more likely to exhibit above high parent expression level than hybrids that do not express the silent allele, leading to non-additive expression. This study provides a comprehensive understanding of the nature of non-PAV SPE and PAV SPE genes and their roles in gene expression complementation in maize hybrids.

Transposable elements contribute to dynamic genome content in maize.

Transposable elements (TEs) are ubiquitous components of eukaryotic genomes and can create variation in genome organization and content. Most maize genomes are composed of TEs. We developed an approach to define shared and variable TE insertions across genome assemblies and applied this method to four maize genomes (B73, W22, Mo17 and PH207) with uniform structural annotations of TEs. Among these genomes we identified approximately 400 000 TEs that are polymorphic, encompassing 1.6 Gb of variable TE sequence. These polymorphic TEs include a combination of recent transposition events as well as deletions of older TEs. There are examples of polymorphic TEs within each of the superfamilies of TEs and they are found distributed across the genome, including in regions of recent shared ancestry among individuals. There are many examples of polymorphic TEs within or near maize genes. In addition, there are 2380 gene annotations in the B73 genome that are located within variable TEs, providing evidence for the role of TEs in contributing to the substantial differences in annotated gene content among these genotypes. TEs are highly variable in our survey of four temperate maize genomes, highlighting the major contribution of TEs in driving variation in genome organization and gene content. 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://github.com/SNAnderson/maizeTE_variation; https://mcstitzer.github.io/maize_TEs.

Using multiple reference genomes to identify and resolve annotation inconsistencies.

BACKGROUND: Advances in sequencing technologies have led to the release of reference genomes and annotations for multiple individuals within more well-studied systems. While each of these new genome assemblies shares significant portions of synteny between each other, the annotated structure of gene models within these regions can differ. Of particular concern are split-gene misannotations, in which a single gene is incorrectly annotated as two distinct genes or two genes are incorrectly annotated as a single gene. These misannotations can have major impacts on functional prediction, estimates of expression, and many downstream analyses. RESULTS: We developed a high-throughput method based on pairwise comparisons of annotations that detect potential split-gene misannotations and quantifies support for whether the genes should be merged into a single gene model. We demonstrated the utility of our method using gene annotations of three reference genomes from maize (B73, PH207, and W22), a difficult system from an annotation perspective due to the size and complexity of the genome. On average, we found several hundred of these potential split-gene misannotations in each pairwise comparison, corresponding to 3-5% of gene models across annotations. To determine which state (i.e. one gene or multiple genes) is biologically supported, we utilized RNAseq data from 10 tissues throughout development along with a novel metric and simulation framework. The methods we have developed require minimal human interaction and can be applied to future assemblies to aid in annotation efforts. CONCLUSIONS: Split-gene misannotations occur at appreciable frequency in maize annotations. We have developed a method to easily identify and correct these misannotations. Importantly, this method is generic in that it can utilize any type of short-read expression data. Failure to account for split-gene misannotations has serious consequences for biological inference, particularly for expression-based analyses.

The limited role of differential fractionation in genome content variation and function in maize (Zea mays L.) inbred lines.

Maize is a diverse paleotetraploid species with considerable presence/absence variation and copy number variation. One mechanism through which presence/absence variation can arise is differential fractionation. Fractionation refers to the loss of duplicate gene pairs from one of the maize subgenomes during diploidization. Differential fractionation refers to non-shared gene loss events between individuals following a whole-genome duplication event. We investigated the prevalence of presence/absence variation resulting from differential fractionation in the syntenic portion of the genome using two whole-genome de novo assemblies of the inbred lines B73 and PH207. Between these two genomes, syntenic genes were highly conserved with less than 1% of syntenic genes being subject to differential fractionation. The few variably fractionated syntenic genes that were identified are unlikely to contribute to functional phenotypic variation, as there is a significant depletion of these genes in annotated gene sets. In further comparisons of 60 diverse inbred lines, non-syntenic genes were six times more likely to be variable than syntenic genes, suggesting that comparisons among additional genome assemblies are not likely to result in the discovery of large-scale presence/absence variation among syntenic genes.

Draft Assembly of Elite Inbred Line PH207 Provides Insights into Genomic and Transcriptome Diversity in Maize.

Intense artificial selection over the last 100 years has produced elite maize (Zea mays) inbred lines that combine to produce high-yielding hybrids. To further our understanding of how genome and transcriptome variation contribute to the production of high-yielding hybrids, we generated a draft genome assembly of the inbred line PH207 to complement and compare with the existing B73 reference sequence. B73 is a founder of the Stiff Stalk germplasm pool, while PH207 is a founder of Iodent germplasm, both of which have contributed substantially to the production of temperate commercial maize and are combined to make heterotic hybrids. Comparison of these two assemblies revealed over 2500 genes present in only one of the two genotypes and 136 gene families that have undergone extensive expansion or contraction. Transcriptome profiling revealed extensive expression variation, with as many as 10,564 differentially expressed transcripts and 7128 transcripts expressed in only one of the two genotypes in a single tissue. Genotype-specific genes were more likely to have tissue/condition-specific expression and lower transcript abundance. The availability of a high-quality genome assembly for the elite maize inbred PH207 expands our knowledge of the breadth of natural genome and transcriptome variation in elite maize inbred lines across heterotic pools.

Genome-Wide Association and Gene Co-expression Network Analyses Reveal Complex Genetics of Resistance to Goss's Wilt of Maize.

Goss's bacterial wilt and leaf blight is a disease of maize caused by the gram positive bacterium Clavibacter michiganensis subsp. nebraskensis (Cmn). First discovered in Nebraska, Goss's wilt has now spread to major maize growing states in the United States and three provinces in Canada. Previous studies conducted using elite maize inbred lines and their hybrids have shown that resistance to Goss's wilt is a quantitative trait. The objective of this study was to further our understanding of the genetic basis of resistance to Goss's wilt by using a combined approach of genome-wide association mapping and gene co-expression network analysis. Genome-wide association analysis was accomplished using a diversity panel consisting of 555 maize inbred lines and a set of 450 recombinant inbred lines (RILs) from three bi-parental mapping populations, providing the most comprehensive screening of Goss's wilt resistance to date. Three SNPs in the diversity panel and 10 SNPs in the combined dataset, including the diversity panel and RILs, were found to be significantly associated with Goss's wilt resistance. Each significant SNP explained 1-5% of the phenotypic variation for Goss's wilt (total of 8-11%). To augment the results of genome-wide association mapping and help identify candidate genes, a time course RNA sequencing experiment was conducted using resistant (N551) and susceptible (B14A) maize inbred lines. Gene co-expression network analysis of this time course experiment identified one module of 141 correlated genes that showed differential regulation in response to Cmn inoculations in both resistant and susceptible lines. SNPs inside and flanking these genes explained 13.3% of the phenotypic variation. Among 1,000 random samples of genes, only 8% of samples explained more phenotypic variance for Goss's wilt resistance than those implicated by the co-expression network analysis. While a statistically significant enrichment was not observed (P < 0.05), these results suggest a possible role for these genes in quantitative resistance at the field level and warrant more research on combining gene co-expression network analysis with quantitative genetic analyses to dissect complex disease resistance traits. The results of the GWAS and co-expression analysis both support the complex nature of resistance to this important disease of maize.

Tandem Duplicate Genes in Maize Are Abundant and Date to Two Distinct Periods of Time.

Tandem duplicate genes are proximally duplicated and as such occur in similar genomic neighborhoods. Using the maize B73 and PH207 de novo genome assemblies, we identified thousands of tandem gene duplicates that account for ∼10% of the annotated genes. These tandem duplicates have a bimodal distribution of ages, which coincide with ancient allopolyploidization and more recent domestication. Tandem duplicates are smaller on average and have a higher probability of containing LTR elements than other genes, suggesting origins in nonhomologous recombination. Within relatively recent tandem duplicate genes, ∼26% appear to be undergoing degeneration or divergence in function from the ancestral copy. Our results show that tandem duplicates are abundant in maize, arose in bursts throughout maize evolutionary history under multiple potential mechanisms, and may provide a substrate for novel phenotypic variation.