Evolution of DNA sequence nonhomologies among maize inbreds

Structural variation discovery in wheat using PacBio high-fidelity sequencing

Structural variations (SVs) pervade plant genomes and contribute substantially to the phenotypic diversity. However, most SVs were ineffectively assayed due to their complex nature and the limitations of early genomic technologies. By applying the PacBio high-fidelity (HiFi) sequencing for wheat genomes, we performed a comprehensive evaluation of mainstream long-read aligners and SV callers in SV detection. The results indicated that the accuracy of deletion discovery is markedly influenced by callers, accounting for 87.73% of the variance, whereas both aligners (38.25%) and callers (49.32%) contributed substantially to the accuracy variance for insertions. Among the aligners, Winnowmap2 and NGMLR excelled in detecting deletions and insertions, respectively. For SV callers, SVIM achieved the best performance. We demonstrated that combining the aligners and callers mentioned above is optimal for SV detection. Furthermore, we evaluated the effect of sequencing depth on the accuracy of SV detection, revealing that low-coverage HiFi sequencing is sufficiently robust for high-quality SV discovery. This study thoroughly evaluated SV discovery approaches and established optimal workflows for investigating structural variations using low-coverage HiFi sequencing in the wheat genome, which will advance SV discovery and decipher the biological functions of SVs in wheat and many other plants.

Duplication and sub-functionalization of flavonoid biosynthesis genes plays important role in Leguminosae root nodule symbiosis evolution

Gene innovation plays an essential role in trait evolution. Rhizobial symbioses, the most important N2-fixing agent in agricultural systems that exists mainly in Leguminosae, is one of the most attractive evolution events. However, the gene innovations underlying Leguminosae root nodule symbiosis (RNS) remain largely unknown. Here, we investigated the gene gain event in Leguminosae RNS evolution through comprehensive phylogenomic analyses. We revealed that Leguminosae-gain genes were acquired by gene duplication and underwent a strong purifying selection. Kyoto Encyclopedia of Genes and Genomes analyses showed that the innovated genes were enriched in flavonoid biosynthesis pathways, particular downstream of chalcone synthase (CHS). Among them, Leguminosae-gain type Ⅱ chalcone isomerase (CHI) could be further divided into CHI1A and CHI1B clades, which resulted from the products of tandem duplication. Furthermore, the duplicated CHI genes exhibited exon-intron structural divergences evolved through exon/intron gain/loss and insertion/deletion. Knocking down CHI1B significantly reduced nodulation in Glycine max (soybean) and Medicago truncatula; whereas, knocking down its duplication gene CHI1A had no effect on nodulation. Therefore, Leguminosae-gain type Ⅱ CHI participated in RNS and the duplicated CHI1A and CHI1B genes exhibited RNS functional divergence. This study provides functional insights into Leguminosae-gain genetic innovation and sub-functionalization after gene duplication that contribute to the evolution and adaptation of RNS in Leguminosae.

Plant pangenomes for crop improvement, biodiversity and evolution

Plant genome sequences catalogue genes and the genetic elements that regulate their expression. Such inventories further research aims as diverse as mapping the molecular basis of trait diversity in domesticated plants or inquiries into the origin of evolutionary innovations in flowering plants millions of years ago. The transformative technological progress of DNA sequencing in the past two decades has enabled researchers to sequence ever more genomes with greater ease. Pangenomes - complete sequences of multiple individuals of a species or higher taxonomic unit - have now entered the geneticists' toolkit. The genomes of crop plants and their wild relatives are being studied with translational applications in breeding in mind. But pangenomes are applicable also in ecological and evolutionary studies, as they help classify and monitor biodiversity across the tree of life, deepen our understanding of how plant species diverged and show how plants adapt to changing environments or new selection pressures exerted by human beings.

Are cereal grasses a single genetic system?

In 1993, a passionate and provocative call to arms urged cereal researchers to consider the taxon they study as a single genetic system and collaborate with each other. Since then, that group of scientists has seen their discipline blossom. In an attempt to understand what unity of genetic systems means and how the notion was borne out by later research, we survey the progress and prospects of cereal genomics: sequence assemblies, population-scale sequencing, resistance gene cloning and domestication genetics. Gene order may not be as extraordinarily well conserved in the grasses as once thought. Still, several recurring themes have emerged. The same ancestral molecular pathways defining plant architecture have been co-opted in the evolution of different cereal crops. Such genetic convergence as much as cross-fertilization of ideas between cereal geneticists has led to a rich harvest of genes that, it is hoped, will lead to improved varieties.

Combined analysis of transposable elements and structural variation in maize genomes reveals genome contraction outpaces expansion

Structural differences between genomes are a major source of genetic variation that contributes to phenotypic differences. Transposable elements, mobile genetic sequences capable of increasing their copy number and propagating themselves within genomes, can generate structural variation. However, their repetitive nature makes it difficult to characterize fine-scale differences in their presence at specific positions, limiting our understanding of their impact on genome variation. Domesticated maize is a particularly good system for exploring the impact of transposable element proliferation as over 70% of the genome is annotated as transposable elements. High-quality transposable element annotations were recently generated for de novo genome assemblies of 26 diverse inbred maize lines. We generated base-pair resolved pairwise alignments between the B73 maize reference genome and the remaining 25 inbred maize line assemblies. From this data, we classified transposable elements as either shared or polymorphic in a given pairwise comparison. Our analysis uncovered substantial structural variation between lines, representing both simple and complex connections between TEs and structural variants. Putative insertions in SNP depleted regions, which represent recently diverged identity by state blocks, suggest some TE families may still be active. However, our analysis reveals that within these recently diverged genomic regions, deletions of transposable elements likely account for more structural variation events and base pairs than insertions. These deletions are often large structural variants containing multiple transposable elements. Combined, our results highlight how transposable elements contribute to structural variation and demonstrate that deletion events are a major contributor to genomic differences.

Whole-genome and dispersed duplication, including transposed duplication, jointly advance the evolution of TLP genes in seven representative Poaceae lineages

BACKGROUND

In the evolutionary study of gene families, exploring the duplication mechanisms of gene families helps researchers understand their evolutionary history. The tubby-like protein (TLP) family is essential for growth and development in plants and animals. Much research has been done on its function; however, limited information is available with regard to the evolution of the TLP gene family. Herein, we systematically investigated the evolution of TLP genes in seven representative Poaceae lineages.

RESULTS

Our research showed that the evolution of TLP genes was influenced not only by whole-genome duplication (WGD) and dispersed duplication (DSD) but also by transposed duplication (TRD), which has been neglected in previous research. For TLP family size, we found an evolutionary pattern of progressive shrinking in the grass family. Furthermore, the evolution of the TLP gene family was at least affected by evolutionary driving forces such as duplication, purifying selection, and base mutations.

CONCLUSIONS

This study presents the first comprehensive evolutionary analysis of the TLP gene family in grasses. We demonstrated that the TLP gene family is also influenced by a transposed duplication mechanism. Several new insights into the evolution of the TLP gene family are presented. This work provides a good reference for studying gene evolution and the origin of duplication.

Large haplotypes highlight a complex age structure within the maize pan-genome

The genomes of maize and other eukaryotes contain stable haplotypes in regions of low recombination. These regions, including centromeres, long heterochromatic blocks, and rDNA arrays, have been difficult to analyze with respect to their diversity and origin. Greatly improved genome assemblies are now available that enable comparative genomics over these and other nongenic spaces. Using 26 complete maize genomes, we developed methods to align intergenic sequences while excluding genes and regulatory regions. The centromere haplotypes (cenhaps) extend for megabases on either side of the functional centromere regions and appear as evolutionary strata, with haplotype divergence/coalescence times dating as far back as 450 thousand years ago (kya). Application of the same methods to other low recombination regions (heterochromatic knobs and rDNA) and all intergenic spaces revealed that deep coalescence times are ubiquitous across the maize pan-genome. Divergence estimates vary over a broad timescale with peaks at ∼16 and 300 kya, reflecting a complex history of gene flow among diverging populations and changes in population size associated with domestication. Cenhaps and other long haplotypes provide vivid displays of this ancient diversity.

Plant pan-genomics and its applications

Plant genomes are so highly diverse that a substantial proportion of genomic sequences are not shared among individuals. The variable DNA sequences, along with the conserved core sequences, compose the more sophisticated pan-genome that represents the collection of all non-redundant DNA in a species. With rapid progress in genome sequencing technologies, pan-genome research in plants is now accelerating. Here we review recent advances in plant pan-genomics, including major driving forces of structural variations that constitute the variable sequences, methodological innovations for representing the pan-genome, and major successes in constructing plant pan-genomes. We also summarize recent efforts toward decoding the remaining dark matter in telomere-to-telomere or gapless plant genomes. These new genome resources, which have remarkable advantages over numerous previously assembled less-than-perfect genomes, are expected to become new references for genetic studies and plant breeding.

Annotation of Siberian Larch (Larix sibirica Ledeb.) Nuclear Genome—One of the Most Cold-Resistant Tree Species in the Only Deciduous GENUS in Pinaceae

The recent release of the nuclear, chloroplast and mitochondrial genome assemblies of Siberian larch (Larix sibirica Ledeb.), one of the most cold-resistant tree species in the only deciduous genus of Pinaceae, with seasonal senescence and a rot-resistant valuable timber widely used in construction, greatly contributed to the development of genomic resources for the larch genus. Here, we present an extensive repeatome analysis and the first annotation of the draft nuclear Siberian larch genome assembly. About 66% of the larch genome consists of highly repetitive elements (REs), with the likely wave of retrotransposons insertions into the larch genome estimated to occur 4–5 MYA. In total, 39,370 gene models were predicted, with 87% of them having homology to the Arabidopsis-annotated proteins and 78% having at least one GO term assignment. The current state of the genome annotations allows for the exploration of the gymnosperm and angiosperm species for relative gene abundance in different functional categories. Comparative analysis of functional gene categories across different angiosperm and gymnosperm species finds that the Siberian larch genome has an overabundance of genes associated with programmed cell death (PCD), autophagy, stress hormone biosynthesis and regulatory pathways; genes that may play important roles in seasonal senescence and stress response to extreme cold in larch. Despite being incomplete, the draft assemblies and annotations of the conifer genomes are at a point of development where they now represent a valuable source for further genomic, genetic and population studies.

The Practical Haplotype Graph, a platform for storing and using pangenomes for imputation

MOTIVATION

Pangenomes provide novel insights for population and quantitative genetics, genomics, and breeding not available from studying a single reference genome. Instead, a species is better represented by a pangenome or collection of genomes. Unfortunately, managing and using pangenomes for genomically diverse species is computationally and practically challenging. We developed a trellis graph representation anchored to the reference genome that represents most pangenomes well and can be used to impute complete genomes from low density sequence or variant data.

RESULTS

The Practical Haplotype Graph (PHG) is a pangenome pipeline, database (PostGRES & SQLite), data model (Java, Kotlin, or R), and Breeding API (BrAPI) web service. The PHG has already been able to accurately represent diversity in four major crops including maize, one of the most genomically diverse species, with up to 1000-fold data compression. Using simulated data, we show that, at even 0.1X coverage, with appropriate reads and sequence alignment, imputation results in extremely accurate haplotype reconstruction. The PHG is a platform and environment for the understanding and application of genomic diversity.

AVAILABILITY

All resources listed here are freely available. The PHG Docker used to generate the simulation results is https://hub.docker.com/ as maizegenetics/phg:0.0.27. PHG source code is at https://bitbucket.org/bucklerlab/practicalhaplotypegraph/src/master/. The code used for the analysis of simulated data is at https://bitbucket.org/bucklerlab/phg-manuscript/src/master/. The PHG database of NAM parent haplotypes is in the CyVerse data store (https://de.cyverse.org/de/) and named /iplant/home/shared/panzea/panGenome/PHG_db_maize/phg_v5Assemblies_20200608.db.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

An SGSGeneloss-Based Method for Constructing a Gene Presence-Absence Table Using Mosdepth

Presence-absence variants (PAV) are genomic regions present in some individuals of a species, but not others. PAVs have been shown to contribute to genomic diversity, especially in bacteria and plants. These structural variations have been linked to traits and can be used to track a species' evolutionary history. PAVs are usually called by aligning short read sequence data from one or more individuals to a reference genome or pangenome assembly, and then comparing coverage. Regions where reads do not align define absence in that individual, and the regions are classified as PAVs. The method below details how to align sequence reads to a reference and how to use the sequencing-coverage calculator Mosdepth to identify PAVs and construct a PAV table for use in downstream comparative genome analysis.

M. S, Beyene Y, Das B, Prasanna BM, Gowda M. Discovery and Validation of a Recessively Inherited Major-Effect QTL Conferring Resistance to Maize Lethal Necrosis (MLN) Disease

Maize lethal necrosis (MLN) is a viral disease with a devastating effect on maize production. Developing and deploying improved varieties with resistance to the disease is important to effectively control MLN; however, little is known about the causal genes and molecular mechanism(s) underlying MLN resistance. Screening thousands of maize inbred lines revealed KS23-5 and KS23-6 as two of the most promising donors of MLN resistance alleles. KS23-5 and KS23-6 lines were earlier developed at the University of Hawaii, United States, on the basis of a source population constituted using germplasm from Kasetsart University, Thailand. Both linkage mapping and association mapping approaches were used to discover and validate genomic regions associated with MLN resistance. Selective genotyping of resistant and susceptible individuals within large F2 populations coupled with genome-wide association study identified a major-effect QTL (qMLN06_157) on chromosome 6 for MLN disease severity score and area under the disease progress curve values in all three F2 populations involving one of the KS23 lines as a parent. The major-effect QTL (qMLN06_157) is recessively inherited and explained 55%-70% of the phenotypic variation with an approximately 6 Mb confidence interval. Linkage mapping in three F3 populations and three F2 populations involving KS23-5 or KS23-6 as one of the parents confirmed the presence of this major-effect QTL on chromosome 6, demonstrating the efficacy of the KS23 allele at qMLN06.157 in varying populations. This QTL could not be identified in population that was not derived using KS23 lines. Validation of this QTL in six F2 populations with 20 SNPs closely linked with qMLN06.157 was further confirmed its consistent expression across populations and its recessive nature of inheritance. On the basis of the consistent and effective resistance afforded by the KS23 allele at qMLN06.157, the QTL can be used in both marker-assisted forward breeding and marker-assisted backcrossing schemes to improve MLN resistance of breeding populations and key lines for eastern Africa.

High-Throughput and Accurate Determination of Transgene Copy Number and Zygosity in Transgenic Maize: From DNA Extraction to Data Analysis

It is vital to develop high-throughput methods to determine transgene copy numbers initially and zygosity during subsequent breeding. In this study, the target sequence of the previously reported endogenous reference gene hmg was analyzed using 633 maize inbred lines, and two SNPs were observed. These SNPs significantly increased the PCR efficiency, while the newly developed hmg gene assay (hmg-taq-F2/R2) excluding these SNPs reduced the efficiency into normal ranges. The TaqMan amplification efficiency of bar and hmg with newly developed primers was calculated as 0.993 and 1.000, respectively. The inter-assay coefficient of variation (CV) values for the bar and hmg genes varied from 1.18 to 2.94%. The copy numbers of the transgene bar using new TaqMan assays were identical to those using dPCR. Significantly, the precision of one repetition reached 96.7% of that of three repetitions of single-copy plants analyzed by simple random sampling, and the actual accuracy reached 95.8%, confirmed by T₁ and T₂ progeny. With the high-throughput DNA extraction and automated data analysis procedures developed in this study, nearly 2700 samples could be analyzed within eight hours by two persons. The combined results suggested that the new hmg gene assay developed here could be a universal maize reference gene system, and the new assay has high throughput and high accuracy for large-scale screening of maize varieties around the world.

CRISPR-Cas9-mediated 75.5-Mb inversion in maize

CRISPR-Cas is a powerful double-strand-break technology with wide-ranging applications from gene discovery to commercial product development. Thus far, this tool has been almost exclusively used for gene knockouts and deletions, with a few examples of gene edits and targeted gene insertions. Here, we demonstrate the application of CRISPR-Cas9 technology to mediate targeted 75.5-Mb pericentric inversion in chromosome 2 in one of the elite maize inbred lines from Corteva Agriscience. This inversion unlocks a large chromosomal region containing substantial genetic variance for recombination, thus providing opportunities for the development of new maize varieties with improved phenotypes.

The barley pan-genome reveals the hidden legacy of mutation breeding

Genetic diversity is key to crop improvement. Owing to pervasive genomic structural variation, a single reference genome assembly cannot capture the full complement of sequence diversity of a crop species (known as the 'pan-genome'1). Multiple high-quality sequence assemblies are an indispensable component of a pan-genome infrastructure. Barley (Hordeum vulgare L.) is an important cereal crop with a long history of cultivation that is adapted to a wide range of agro-climatic conditions2. Here we report the construction of chromosome-scale sequence assemblies for the genotypes of 20 varieties of barley-comprising landraces, cultivars and a wild barley-that were selected as representatives of global barley diversity. We catalogued genomic presence/absence variants and explored the use of structural variants for quantitative genetic analysis through whole-genome shotgun sequencing of 300 gene bank accessions. We discovered abundant large inversion polymorphisms and analysed in detail two inversions that are frequently found in current elite barley germplasm; one is probably the product of mutation breeding and the other is tightly linked to a locus that is involved in the expansion of geographical range. This first-generation barley pan-genome makes previously hidden genetic variation accessible to genetic studies and breeding.

Pan-genome: A promising resource for noncoding RNA discovery in plants

Plant genomes contain both protein-coding and noncoding sequences including transposable elements (TEs) and noncoding RNAs (ncRNAs). The ncRNAs are recognized as important elements that play fundamental roles in the structural organization and function of plant genomes. Despite various hypotheses, TEs are believed to be a major precursor of ncRNAs. Transposable elements are also prime factors that cause genomic variation among members of a species. Hence, TEs pose a major challenge in the discovery and analysis of ncRNAs. With the increase in the number of sequenced plant genomes, it is now accepted that a single reference genome is insufficient to represent the complete genomic diversity and contents of a species, and exploring the pan-genome of a species is critical. In this review, we summarize the recent progress in the field of plant pan-genomes. We also discuss TEs and their roles in ncRNA biogenesis and present our perspectives on the application of pan-genomes for the discovery of ncRNAs to fully explore and exploit their biological roles in plants.

Genetic differentiation of Mutator insertion polymorphisms and association with agronomic traits in waxy and common maize

BACKGROUND

As waxy maize is considered a key economic crop in Korea, an understanding of its genetic variation and differentiation is fundamental for the selective plant breeding. The maize genome is primarily composed of transposable elements, for which large and stable insertions generate variations that reflect selection during evolution.

OBJECTIVES

This study was to elucidate the genetic diversity based on the contribution of TEs and to investigate the effect of Mu transposition on the genetic divergence of waxy and common maize. We also performed an association analysis on these inbred lines to determine the Mu insertions associated with agronomic traits.

METHODS

In this study, we utilized a Mutator-based transposon display method to study the genetic diversity and population structure of 40 waxy and 40 common inbred lines of maize in the Gangwon Agricultural Research and Extension Services collection at the Maize Research Institute.

RESULTS

We detected polymorphisms in 86.33% of 278 Mutator (Mu) anchored loci, reflecting the activity of the Mu element and its contribution to genetic variation. Common maize showed a substantial amount of genetic diversity, which was greater than that observed in waxy maize. Principal-coordinate and neighbor-joining cluster analyzes consistently supported the presence of two genetically distinct groups. However, the distribution of genetic variation within the populations was much higher than the genetic differentiation among the populations. To explore the contribution of the Mu element to phenotypic variation, we analyzed the associations with ten important agronomical traits. On the basis of the combined results from two models (QGLM and Q + KLM), we found significant associations between seven Mu loci and four different traits.

CONCLUSIONS

These results will assist waxy maize breeders in choosing parental lines and be useful for marker-assisted selection.

Integrated annotations and analyses of small RNA-producing loci from 47 diverse plants

Plant endogenous small RNAs (sRNAs) are important regulators of gene expression. There are two broad categories of plant sRNAs: microRNAs (miRNAs) and endogenous short interfering RNAs (siRNAs). MicroRNA loci are relatively well-annotated but compose only a small minority of the total sRNA pool; siRNA locus annotations have lagged far behind. Here, we used a large data set of published and newly generated sRNA sequencing data (1333 sRNA-seq libraries containing more than 20 billion reads) and a uniform bioinformatic pipeline to produce comprehensive sRNA locus annotations of 47 diverse plants, yielding more than 2.7 million sRNA loci. The two most numerous classes of siRNA loci produced mainly 24- and 21-nucleotide (nt) siRNAs, respectively. Most often, 24-nt-dominated siRNA loci occurred in intergenic regions, especially at the 5′-flanking regions of protein-coding genes. In contrast, 21-nt-dominated siRNA loci were most often derived from double-stranded RNA precursors copied from spliced mRNAs. Genic 21-nt-dominated loci were especially common from disease resistance genes, including from a large number of monocots. Individual siRNA sequences of all types showed very little conservation across species, whereas mature miRNAs were more likely to be conserved. We developed a web server where our data and several search and analysis tools are freely accessible.

Haplotype structure in commercial maize breeding programs in relation to key founder lines

High-density haplotype analysis revealed significant haplotype sharing between ex-PVPs registered from 1976 to 1992 and key maize founders, and uncovered similarities and differences in haplotype sharing patterns by company and heterotic group. Proprietary inbreds developed by the private seed industry have been the major source for driving genetic gain in successful North American maize hybrids for decades. Much of the history of industry germplasm can be traced back to key founder lines, some of which were pivotal in the development of prominent heterotic groups. Previous studies have summarized pedigree-based relationships, genetic diversity and population structure among commercial inbreds with expired Plant Variety Protection (ex-PVP). However, less is known about the extent of haplotype sharing between historical founders and ex-PVPs. A better understanding of the relationships between founders and ex-PVPs provides insight into the haplotype and heterotic group structure among industry germplasm. We performed high-density haplotype analysis with 11.3 million SNPs on 212 maize inbreds, which included 157 ex-PVPs registered 1976-1992 and 55 public lines relevant to PVPs. Among these lines were 12 key founders identified in literature review: 207, A632, B14, B37, B73, LH123HT, LH82, Mo17, Oh43, OH7, PHG39 and Wf9. Our results revealed that, on average, 81.6% of an ex-PVP's genome is shared with at least 1 of these 12 founder lines and more than half when limited to B73, Mo17 and 207. Quantifiable similarities and contrasts among heterotic groups and major US seed industry companies were also observed. The results from this study provide high-resolution haplotype data on ex-PVP germplasm, confirm founder relationship trends observed in previous studies, uncover region-specific haplotype structure differences and demonstrate how haplotype sharing analysis can be used as a tool to explore germplasm diversity.

The Landscape of Copia and Gypsy Retrotransposon During Maize Domestication and Improvement

The release of genomic sequences in the maize HapMap3 population provides an opportunity to study the genetic diversity of maize. In this study, retrotransposon insertion polymorphisms (RIPs) were mapped against the maize genome sequence. In total, 27 retrotransposon families were identified, and more than 170,000 RIPs were discovered in teosinte, landrace, and improved groups. Interestingly, the copy number of transposable elements (TEs) were more abundant in landrace groups than in teosinte or improved groups, suggesting that TEs experienced amplification during domestication and contraction during improvement. Landrace accessions exhibited higher TE insertion frequency compared to the other groups. Furthermore, the position of TE insertions were closer to genes and more abundant in the centromeres of landrace groups compared to the other groups. The three groups could be clearly distinguished by RIPs. These results demonstrate that TEs were amplified and contracted during maize domestication and improvement, respectively.

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.

High throughput genotyping of structural variations in a complex plant genome using an original Affymetrix® axiom® array

Background

Insertions/deletions (InDels) and more specifically presence/absence variations (PAVs) are pervasive in several species and have strong functional and phenotypic effect by removing or drastically modifying genes. Genotyping of such variants on large panels remains poorly addressed, while necessary for approaches such as association mapping or genomic selection.

Results

We have developed, as a proof of concept, a new high-throughput and affordable approach to genotype InDels. We first identified 141,000 InDels by aligning reads from the B73 line against the genome of three temperate maize inbred lines (F2, PH207, and C103) and reciprocally. Next, we designed an Affymetrix® Axiom® array to target these InDels, with a combination of probes selected at breakpoint sites (13%) or within the InDel sequence, either at polymorphic (25%) or non-polymorphic sites (63%) sites. The final array design is composed of 662,772 probes and targets 105,927 InDels, including PAVs ranging from 35 bp to 129kbp. After Affymetrix® quality control, we successfully genotyped 86,648 polymorphic InDels (82% of all InDels interrogated by the array) on 445 maize DNA samples with 422,369 probes. Genotyping InDels using this approach produced a highly reliable dataset, with low genotyping error (~ 3%), high call rate (~ 98%), and high reproducibility (> 95%). This reliability can be further increased by combining genotyping of several probes calling the same InDels (< 0.1% error rate and > 99.9% of call rate for 5 probes). This “proof of concept” tool was used to estimate the kinship matrix between 362 maize lines with 57,824 polymorphic InDels. This InDels kinship matrix was highly correlated with kinship estimated using SNPs from Illumina 50 K SNP arrays.

Conclusions

We efficiently genotyped thousands of small to large InDels on a sizeable number of individuals using a new Affymetrix® Axiom® array. This powerful approach opens the way to studying the contribution of InDels to trait variation and heterosis in maize. The approach is easily extendable to other species and should contribute to decipher the biological impact of InDels at a larger scale.

Interspecific hybridisation and LTR-retrotransposon mobilisation-related structural variation in plants: A case study

The dynamics of long-terminal-repeat retrotransposons in two poplar species (Populus deltoides and P. nigra) and in an interspecific hybrid, recently synthesized, were investigated by analyzing the genomic abundance and transcription levels of a collection of 828 full-length retroelements identified in the genome sequence of P. trichocarpa, all occurring also in the genomes of P. deltoides and P. nigra. Overall, genomic abundance and transcription levels of many retrotransposons in the hybrid resulted higher or lower than expected by calculating the mean of the parental values. A bioinformatics procedure was established to ascertain the occurrence of the same retrotransposon loci in the three genotypes. The results indicated that retrotransposon abundance variations between the hybrid and the mean value of the parents were due to i) co-segregation of retrotransposon high- or low-abundant haplotypes; ii) new retroelement insertions; iii) retrotransposon loss. Concerning retrotransposon expression, this was generally low, with only 14/828 elements over- or under-expressed in the hybrid than expected by calculating the mean of the parents. It is concluded that interspecific hybridisation between the two poplar species determine quantitative variation and differential expression of some retrotransposons, with possible consequences for the genetic differentiation of the hybrid.

The Second Site Modifier, Sympathy for the ligule, Encodes a Homolog of Arabidopsis ENHANCED DISEASE RESISTANCE4 and Rescues the Liguleless narrow Maize Mutant

Liguleless narrow1 encodes a plasma membrane-localized receptor-like kinase required for normal development of maize (Zea mays) leaves, internodes, and inflorescences. The semidominant Lgn-R mutation lacks kinase activity, and phenotypic severity is dependent on inbred background. We created near isogenic lines and assayed the phenotype in multiple environments. Lgn-R plants that carry the B73 version of Sympathy for the ligule (Sol-B) fail to grow under hot conditions, but those that carry the Mo17 version (Sol-M) survive at hot temperatures and are significantly taller at cool temperatures. To identify Sol, we used recombinant mapping and analyzed the Lgn-R phenotype in additional inbred backgrounds. We identified amino acid sequence variations in GRMZM2G075262 that segregate with severity of the Lgn-R phenotypes. This gene is expressed at high levels in Lgn-R B73, but expression drops to nonmutant levels with one copy of Sol-M An EMS mutation solidified the identity of SOL as a maize homolog of Arabidopsis (Arabidopsis thaliana) ENHANCED DISEASE RESISTANCE4 (EDR4). SOL, like EDR4, is induced in response to pathogen-associated molecular patterns such as flg22. Integrated transcriptomic and phosphoproteomic analyses suggest that Lgn-R plants constitutively activate an immune signaling cascade that induces temperature-sensitive responses in addition to defects in leaf development. We propose that aspects of the severe Lgn-R developmental phenotype result from constitutive defense induction and that SOL potentially functions in repressing this response in Mo17 but not B73. Identification of LGN and its interaction with SOL provides insight into the integration of developmental control and immune responses.

Multiple Maize Reference Genomes Impact the Identification of Variants by Genome-Wide Association Study in a Diverse Inbred Panel

Use of a single reference genome for genome-wide association studies (GWAS) limits the gene space represented to that of a single accession. This limitation can complicate identification and characterization of genes located within presence-absence variations (PAVs). In this study, we present the draft de novo genome assembly of 'PHJ89', an 'Oh43'-type inbred line of maize ( L.). From three separate reference genome assemblies ('B73', 'PH207', and PHJ89) that represent the predominant germplasm groups of maize, we generated three separate whole-seedling gene expression profiles and single nucleotide polymorphism (SNP) matrices from a panel of 942 diverse inbred lines. We identified 34,447 (B73), 39,672 (PH207), and 37,436 (PHJ89) transcripts that are not present in the respective reference genome assemblies. Genome-wide association studies were conducted in the 942 inbred panel with both the SNP and expression data values to map (SCMV) resistance. Highlighting the impact of alternative reference genomes in gene discovery, the GWAS results for SCMV resistance with expression values as a surrogate measure of PAV resulted in robust detection of the physical location of a known resistance gene when the B73 reference that contains the gene was used, but not the PH207 reference. This study provides the valuable resource of the Oh43-type PHJ89 genome assembly as well as SNP and expression data for 942 individuals generated from three different reference genomes.

The HuangZaoSi Maize Genome Provides Insights into Genomic Variation and Improvement History of Maize

Maize is a globally important crop that was a classic model plant for genetic studies. Here, we report a 2.2 Gb draft genome sequence of an elite maize line, HuangZaoSi (HZS). Hybrids bred from HZS-improved lines (HILs) are planted in more than 60% of maize fields in China. Proteome clustering of six completed sequenced maize genomes show that 638 proteins fall into 264 HZS-specific gene families with the majority of contributions from tandem duplication events. Resequencing and comparative analysis of 40 HZS-related lines reveals the breeding history of HILs. More than 60% of identified selective sweeps were clustered in identity-by-descent conserved regions, and yield-related genes/QTLs were enriched in HZS characteristic selected regions. Furthermore, we demonstrated that HZS-specific family genes were not uniformly distributed in the genome but enriched in improvement/function-related genomic regions. This study provides an important and novel resource for maize genome research and expands our knowledge on the breadth of genomic variation and improvement history of maize.

Prospects of pan-genomics in barley

The concept of a pan-genome refers to intraspecific diversity in genome content and structure, encompassing both genes and intergenic space. Pan-genomic studies employ a combination of de novo sequence assembly and reference-based alignment to discover and genotype structural variants. The large size and complex structure of Triticeae genomes were for a long time an obstacle for genomic research in barley and its relatives. Now that a reference genome is available, computational pipelines for high-quality sequence assembly are in place, and sequence costs continue to drop, investigations into the structural diversity of the barley genome seem within reach. Here, we review the recent progress on pan-genomics in the model grass Brachypodium distachyon, and the cereal crops rice and maize, and devise a multi-tiered strategy for a pan-genome project in barley. Our design involves: (1) the construction of high-quality de novo sequence assemblies for a small core set of representative genotypes, (2) short-read sequencing of a large diversity panel of genebank accessions to medium coverage and (3) the use of complementary methods such as chromosome-conformation capture sequencing and k-mer-based association genetics. The in silico representation of the barley pan-genome may inform about the mechanisms of structural genome evolution in the Triticeae and supplement quantitative genetics models of crop performance for better accuracy and predictive ability.

Exploring and Exploiting Pan-genomics for Crop Improvement

Genetic variation ranging from single-nucleotide polymorphisms to large structural variants (SVs) can cause variation of gene content among individuals within the same species. There is an increasing appreciation that a single reference genome is insufficient to capture the full landscape of genetic diversity of a species. Pan-genome analysis offers a platform to evaluate the genetic diversity of a species via investigation of its entire genome repertoire. Although a recent wave of pan-genomic studies has shed new light on crop diversity and improvement using advanced sequencing technology, the potential applications of crop pan-genomics in crop improvement are yet to be fully exploited. In this review, we highlight the progress achieved in understanding crop pan-genomics, discuss biological activities that cause SVs, review important agronomical traits affected by SVs, and present our perspective on the application of pan-genomics in crop improvement.

Systematic survey of plant LTR-retrotransposons elucidates phylogenetic relationships of their polyprotein domains and provides a reference for element classification

BACKGROUND

Plant LTR-retrotransposons are classified into two superfamilies, Ty1/copia and Ty3/gypsy. They are further divided into an enormous number of families which are, due to the high diversity of their nucleotide sequences, usually specific to a single or a group of closely related species. Previous attempts to group these families into broader categories reflecting their phylogenetic relationships were limited either to analyzing a narrow range of plant species or to analyzing a small numbers of elements. Furthermore, there is no reference database that allows for similarity based classification of LTR-retrotransposons.

RESULTS

We have assembled a database of retrotransposon encoded polyprotein domains sequences extracted from 5410 Ty1/copia elements and 8453 Ty3/gypsy elements sampled from 80 species representing major groups of green plants (Viridiplantae). Phylogenetic analysis of the three most conserved polyprotein domains (RT, RH and INT) led to dividing Ty1/copia and Ty3/gypsy retrotransposons into 16 and 14 lineages respectively. We also characterized various features of LTR-retrotransposon sequences including additional polyprotein domains, extra open reading frames and primer binding sites, and found that the occurrence and/or type of these features correlates with phylogenies inferred from the three protein domains.

CONCLUSIONS

We have established an improved classification system applicable to LTR-retrotransposons from a wide range of plant species. This system reflects phylogenetic relationships as well as distinct sequence and structural features of the elements. A comprehensive database of retrotransposon protein domains (REXdb) that reflects this classification provides a reference for efficient and unified annotation of LTR-retrotransposons in plant genomes. Access to REXdb related tools is implemented in the RepeatExplorer web server (https://repeatexplorer-elixir.cerit-sc.cz/) or using a standalone version of REXdb that can be downloaded seaparately from RepeatExplorer web page (http://repeatexplorer.org/).

Where to Cross Over? Defining Crossover Sites in Plants

It is believed that recombination in meiosis serves to reshuffle genetic material from both parents to increase genetic variation in the progeny. At the same time, the number of crossovers is usually kept at a very low level. As a consequence, many organisms need to make the best possible use from the one or two crossovers that occur per chromosome in meiosis. From this perspective, the decision of where to allocate rare crossover events becomes an important issue, especially in self-pollinating plant species, which experience limited variation due to inbreeding. However, the freedom in crossover allocation is significantly limited by other, genetic and non-genetic factors, including chromatin structure. Here we summarize recent progress in our understanding of those processes with a special emphasis on plant genomes. First, we focus on factors which influence the distribution of recombination initiation sites and discuss their effects at both, the single hotspot level and at the chromosome scale. We also briefly explain the aspects of hotspot evolution and their regulation. Next, we analyze how recombination initiation sites translate into the development of crossovers and their location. Moreover, we provide an overview of the sequence polymorphism impact on crossover formation and chromosomal distribution.

A Megabase-Scale Deletion is Associated with Phenotypic Variation of Multiple Traits in Maize

Genomic deletions are pervasive in the maize (Zea mays L.) genome, and play important roles in phenotypic variation and adaptive evolution. However, little is known about the biological functions of these genomic deletions. Here, we report the biological function of a megabase-scale deletion, which we identified by position-based cloning of the multi-trait weakened (muw) mutant, which is inherited as a single recessive locus. MUW was mapped to a 5.16-Mb region on chromosome 2. The 5.16-Mb deletion in the muw mutant led to the loss of 48 genes and was responsible for a set of phenotypic abnormities, including wilting leaves, poor yield performance, reduced plant height, increased stomatal density, and rapid water loss. While muw appears to have resulted from double-stranded break repair that was not dependent on intragenomic DNA homology, extensive duplication of maize genes may have mitigated its effects and facilitated its survival.

Shared and genetically distinct Zea mays transcriptome responses to ongoing and past low temperature exposure

BACKGROUND

Cold temperatures and their alleviation affect many plant traits including the abundance of protein coding gene transcripts. Transcript level changes that occur in response to cold temperatures and their alleviation are shared or vary across genotypes. In this study we identify individual transcripts and groups of functionally related transcripts that consistently respond to cold and its alleviation. Genes that respond differently to temperature changes across genotypes may have limited functional importance. We investigate if these genes share functions, and if their genotype-specific gene expression levels change in magnitude or rank across temperatures.

RESULTS

We estimate transcript abundances from over 22,000 genes in two unrelated Zea mays inbred lines during and after cold temperature exposure. Genotype and temperature contribute to many genes' abundances. Past cold exposure affects many fewer genes. Genes up-regulated in cold encode many cytokinin glucoside biosynthesis enzymes, transcription factors, signalling molecules, and proteins involved in diverse environmental responses. After cold exposure, protease inhibitors and cuticular wax genes are newly up-regulated, and environmentally responsive genes continue to be up-regulated. Genes down-regulated in response to cold include many photosynthesis, translation, and DNA replication associated genes. After cold exposure, DNA replication and translation genes are still preferentially downregulated. Lignin and suberin biosynthesis are newly down-regulated. DNA replication, reactive oxygen species response, and anthocyanin biosynthesis genes have strong, genotype-specific temperature responses. The ranks of genotypes' transcript abundances often change across temperatures.

CONCLUSIONS

We report a large, core transcriptome response to cold and the alleviation of cold. In cold, many of the core suite of genes are up or downregulated to control plant growth and photosynthesis and limit cellular damage. In recovery, core responses are in part to prepare for future stress. Functionally related genes are consistently and greatly up-regulated in a single genotype in response to cold or its alleviation, suggesting positive selection has driven genotype-specific temperature responses in maize.

Construction of genetic linkage map and identification of QTLs related to agronomic traits in maize using DNA transposon-based markers

Transposable elements (TEs), are a rich source for molecular marker development as they constitute a significant fraction of the eukaryotic genome and impact the overall genome structure. Here, we utilize Mutator-based transposon display (Mu-TD), and CACTA-derived sequence-characterized amplified regions (SCAR) anchored by simple sequence repeats and single nucleotide polymorphisms to locate quantitative trait loci (QTLs) linked to agriculturally important traits on a genetic map. Specifically, we studied recombinant inbred line populations derived from a cross between dent corn and waxy corn. The resulting linkage map included 259 Mu-anchored fragments, 34 SCARs, and 614 SSR markers distributed throughout the ten maize chromosomes. Linkage analysis revealed three SNP loci associated with kernel starch synthesis genes (sh2, su1, wx1) linked to either Mu-TD loci or SSR markers, which may be useful for maize breeding programs. In addition, we used QTL analysis to determine the chromosomal location of traits related to grain yield and kernel quality. We identified 24 QTLs associated with nine traits located on nine out of ten maize chromosomes. Among these, 13 QTLs involved Mu loci and two involved SCARs. This study demonstrates the potential use of DNA transposon-based markers to construct linkage maps and identify QTLs linked to agronomic traits.

Specific LTR-Retrotransposons Show Copy Number Variations between Wild and Cultivated Sunflowers

The relationship between variation of the repetitive component of the genome and domestication in plant species is not fully understood. In previous work, variations in the abundance and proximity to genes of long terminal repeats (LTR)-retrotransposons of sunflower (Helianthus annuus L.) were investigated by Illumina DNA sequencingtocompare cultivars and wild accessions. In this study, we annotated and characterized 22 specific retrotransposon families whose abundance varies between domesticated and wild genotypes. These families mostly belonged to the Chromovirus lineage of the Gypsy superfamily and were distributed overall chromosomes. They were also analyzed in respect to their proximity to genes. Genes close to retrotransposon were classified according to biochemical pathways, and differences between domesticated and wild genotypes are shown. These data suggest that structural variations related to retrotransposons might have occurred to produce phenotypic variation between wild and domesticated genotypes, possibly by affecting the expression of genes that lie close to inserted or deleted retrotransposons and belong to specific biochemical pathways as those involved in plant stress responses.

Extensive intraspecific gene order and gene structural variations between Mo17 and other maize genomes

Maize is an important crop with a high level of genome diversity and heterosis. The genome sequence of a typical female line, B73, was previously released. Here, we report a de novo genome assembly of a corresponding male representative line, Mo17. More than 96.4% of the 2,183 Mb assembled genome can be accounted for by 362 scaffolds in ten pseudochromosomes with 38,620 annotated protein-coding genes. Comparative analysis revealed large gene-order and gene structural variations: approximately 10% of the annotated genes were mutually nonsyntenic, and more than 20% of the predicted genes had either large-effect mutations or large structural variations, which might cause considerable protein divergence between the two inbred lines. Our study provides a high-quality reference-genome sequence of an important maize germplasm, and the intraspecific gene order and gene structural variations identified should have implications for heterosis and genome evolution.

Diverse Chromosomal Locations of Quantitative Trait Loci for Tolerance to Maize chlorotic mottle virus in Five Maize Populations

The recent rapid emergence of maize lethal necrosis (MLN), caused by coinfection of maize with Maize chlorotic mottle virus (MCMV) and a second virus usually from the family Potyviridae, is causing extensive losses for farmers in East Africa, Southeast Asia, and South America. Although the genetic basis of resistance to potyviruses is well understood in maize, little was known about resistance to MCMV. The responses of five maize inbred lines (KS23-5, KS23-6, N211, DR, and Oh1VI) to inoculation with MCMV, Sugarcane mosaic virus, and MLN were characterized. All five lines developed fewer symptoms than susceptible controls after inoculation with MCMV; however, the virus was detected in systemic leaf tissue from each of the lines similarly to susceptible controls, indicating that the lines were tolerant of MCMV rather than resistant to it. Except for KS23-5, the inbred lines also developed fewer symptoms after inoculation with MLN than susceptible controls. To identify genetic loci associated with MCMV tolerance, large F₂ or recombinant inbred populations were evaluated for their phenotypic responses to MCMV, and the most resistant and susceptible plants were genotyped by sequencing. One to four quantitative trait loci (QTL) were identified in each tolerant population using recombination frequency and positional mapping strategies. In contrast to previous studies of virus resistance in maize, the chromosomal positions and genetic character of the QTL were unique to each population. The results suggest that different, genotype-specific mechanisms are associated with MCMV tolerance in maize. These results will allow for the development of markers for marker-assisted selection of MCMV- and MLN-tolerant maize hybrids for disease control.

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.

Sequence analysis of European maize inbred line F2 provides new insights into molecular and chromosomal characteristics of presence/absence variants

BACKGROUND

Maize is well known for its exceptional structural diversity, including copy number variants (CNVs) and presence/absence variants (PAVs), and there is growing evidence for the role of structural variation in maize adaptation. While PAVs have been described in this important crop species, they have been only scarcely characterized at the sequence level and the extent of presence/absence variation and relative chromosomal landscape of inbred-specific regions remain to be elucidated.

RESULTS

De novo genome sequencing of the French F2 maize inbred line revealed 10,044 novel genomic regions larger than 1 kb, making up 88 Mb of DNA, that are present in F2 but not in B73 (PAV). This set of maize PAV sequences allowed us to annotate PAV content and to analyze sequence breakpoints. Using PAV genotyping on a collection of 25 temperate lines, we also analyzed Linkage Disequilibrium in PAVs and flanking regions, and PAV frequencies within maize genetic groups.

CONCLUSIONS

We highlight the possible role of MMEJ-type double strand break repair in maize PAV formation and discover 395 new genes with transcriptional support. Pattern of linkage disequilibrium within PAVs strikingly differs from this of flanking regions and is in accordance with the intuition that PAVs may recombine less than other genomic regions. We show that most PAVs are ancient, while some are found only in European Flint material, thus pinpointing structural features that may be at the origin of adaptive traits involved in the success of this material. Characterization of such PAVs will provide useful material for further association genetic studies in European and temperate maize.

Location of low copy genes in chromosomes of Brachiaria spp

Repetitive DNA sequences have been widely used in cytogenetic analyses. The use of gene sequences with a low-copy-number, however, is little explored especially in plants. To date, the karyotype details in Brachiaria spp. are limited to the location of rDNA sites. The challenge lies in developing new probes based on incomplete sequencing data for the genus or complete sequencing of related species, since there are no model species with a sequenced genome in Brachiaria spp. The present study aimed at the physical location of conserved genes in chromosomes of Brachiaria ruziziensis, Brachiaria brizantha, and Brachiaria decumbens using RNAseq data, as well as sequences of Setaria italica and Sorghum bicolor through the fluorescent in situ hybridization technique. Five out of approximately 90 selected sequences generated clusters in the chromosomes of the species of Brachiaria studied. We identified genes in synteny with 5S and 45S rDNA sites, which contributed to the identification of chromosome pairs carrying these genes. In some cases, the species of Brachiaria evaluated had syntenic segments conserved across the chromosomes. The use of genomic sequencing data is essential for the enhancement of cytogenetic analyses.

Contributions of Zea mays subspecies mexicana haplotypes to modern maize

Maize was domesticated from lowland teosinte (Zea mays ssp. parviglumis), but the contribution of highland teosinte (Zea mays ssp. mexicana, hereafter mexicana) to modern maize is not clear. Here, two genomes for Mo17 (a modern maize inbred) and mexicana are assembled using a meta-assembly strategy after sequencing of 10 lines derived from a maize-teosinte cross. Comparative analyses reveal a high level of diversity between Mo17, B73, and mexicana, including three Mb-size structural rearrangements. The maize spontaneous mutation rate is estimated to be 2.17 × 10^-8 ~3.87 × 10^-8 per site per generation with a nonrandom distribution across the genome. A higher deleterious mutation rate is observed in the pericentromeric regions, and might be caused by differences in recombination frequency. Over 10% of the maize genome shows evidence of introgression from the mexicana genome, suggesting that mexicana contributed to maize adaptation and improvement. Our data offer a rich resource for constructing the pan-genome of Zea mays and genetic improvement of modern maize varieties.

Cultivar-specific transcriptome prediction and annotation in Ficus carica L

The availability of transcriptomic data sequence is a key step for functional genomics studies. Recently, a repertoire of predicted genes of a Japanese cultivar of fig (Ficus carica L.) was released. Because of the great phenotypic variability that can be found in this species, we decided to study another fig genotype, the Italian cv. Dottato, in order to perform comparative studies between the two cultivars and extend the pan genome of this species. We isolated, sequenced and assembled fig genomic DNA from young fruits of cv. Dottato. Then, putative gene sequences were predicted and annotated. Finally, a comparison was performed between cvs. Dottato and Horaishi predicted transcriptomes. Our data provide a resource (available at the Sequence Read Archive database under SRP109082) to be used for functional genomics of fig, in order to fill the gap of knowledge still existing in this species concerning plant development, defense and adaptation to the environment.

Improved maize reference genome with single-molecule technologies

Complete and accurate reference genomes and annotations provide fundamental tools for characterization of genetic and functional variation. These resources facilitate the determination of biological processes and support translation of research findings into improved and sustainable agricultural technologies. Many reference genomes for crop plants have been generated over the past decade, but these genomes are often fragmented and missing complex repeat regions. Here we report the assembly and annotation of a reference genome of maize, a genetic and agricultural model species, using single-molecule real-time sequencing and high-resolution optical mapping. Relative to the previous reference genome, our assembly features a 52-fold increase in contig length and notable improvements in the assembly of intergenic spaces and centromeres. Characterization of the repetitive portion of the genome revealed more than 130,000 intact transposable elements, allowing us to identify transposable element lineage expansions that are unique to maize. Gene annotations were updated using 111,000 full-length transcripts obtained by single-molecule real-time sequencing. In addition, comparative optical mapping of two other inbred maize lines revealed a prevalence of deletions in regions of low gene density and maize lineage-specific genes.

Improved maize reference genome with single-molecule technologies

An improved reference genome for maize, using single-molecule sequencing and high-resolution optical mapping, enables characterization of structural variation and repetitive regions, and identifies lineage expansions of transposable elements that are unique to maize.

The maize genome was initially reported in 2009 but with some accuracy limitations. Doreen Ware and colleagues report a new reference genome for maize using single-molecule sequencing and high-resolution optical mapping. The technique shows improvements in the gene space including resolution of gaps and misassemblies and correction of order and orientation of genes. The authors characterize structural variation and repetitive regions, and identify transposable element lineage expansions unique to maize.

Complete and accurate reference genomes and annotations provide fundamental tools for characterization of genetic and functional variation¹. These resources facilitate the determination of biological processes and support translation of research findings into improved and sustainable agricultural technologies. Many reference genomes for crop plants have been generated over the past decade, but these genomes are often fragmented and missing complex repeat regions². Here we report the assembly and annotation of a reference genome of maize, a genetic and agricultural model species, using single-molecule real-time sequencing and high-resolution optical mapping. Relative to the previous reference genome³, our assembly features a 52-fold increase in contig length and notable improvements in the assembly of intergenic spaces and centromeres. Characterization of the repetitive portion of the genome revealed more than 130,000 intact transposable elements, allowing us to identify transposable element lineage expansions that are unique to maize. Gene annotations were updated using 111,000 full-length transcripts obtained by single-molecule real-time sequencing⁴. In addition, comparative optical mapping of two other inbred maize lines revealed a prevalence of deletions in regions of low gene density and maize lineage-specific 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.

Characterization of the Poplar Pan-Genome by Genome-Wide Identification of Structural Variation

Many recent studies have emphasized the important role of structural variation (SV) in determining human genetic and phenotypic variation. In plants, studies aimed at elucidating the extent of SV are still in their infancy. Evidence has indicated a high presence and an active role of SV in driving plant genome evolution in different plant species.

With the aim of characterizing the size and the composition of the poplar pan-genome, we performed a genome-wide analysis of structural variation in three intercrossable poplar species: Populus nigra, Populus deltoides, and Populus trichocarpa. We detected a total of 7,889 deletions and 10,586 insertions relative to the P. trichocarpa reference genome, covering respectively 33.2 Mb and 62.9 Mb of genomic sequence, and 3,230 genes affected by copy number variation (CNV). The majority of the detected variants are inter-specific in agreement with a recent origin following separation of species.

Insertions and deletions (INDELs) were preferentially located in low-gene density regions of the poplar genome and were, for the majority, associated with the activity of transposable elements. Genes affected by SV showed lower-than-average expression levels and higher levels of dN/dS, suggesting that they are subject to relaxed selective pressure or correspond to pseudogenes.

Functional annotation of genes affected by INDELs showed over-representation of categories associated with transposable elements activity, while genes affected by genic CNVs showed enrichment in categories related to resistance to stress and pathogens. This study provides a genome-wide catalogue of SV and the first insight on functional and structural properties of the poplar pan-genome.

Genome Reduction Uncovers a Large Dispensable Genome and Adaptive Role for Copy Number Variation in Asexually Propagated Solanum tuberosum

Clonally reproducing plants have the potential to bear a significantly greater mutational load than sexually reproducing species. To investigate this possibility, we examined the breadth of genome-wide structural variation in a panel of monoploid/doubled monoploid clones generated from native populations of diploid potato (Solanum tuberosum), a highly heterozygous asexually propagated plant. As rare instances of purely homozygous clones, they provided an ideal set for determining the degree of structural variation tolerated by this species and deriving its minimal gene complement. Extensive copy number variation (CNV) was uncovered, impacting 219.8 Mb (30.2%) of the potato genome with nearly 30% of genes subject to at least partial duplication or deletion, revealing the highly heterogeneous nature of the potato genome. Dispensable genes (>7000) were associated with limited transcription and/or a recent evolutionary history, with lower deletion frequency observed in genes conserved across angiosperms. Association of CNV with plant adaptation was highlighted by enrichment in gene clusters encoding functions for environmental stress response, with gene duplication playing a part in species-specific expansions of stress-related gene families. This study revealed unique impacts of CNV in a species with asexual reproductive habits and how CNV may drive adaption through evolution of key stress pathways.