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

Editing the 19 kDa alpha-zein gene family generates non-opaque2-based quality protein maize

Maize grain is deficient in lysine. While the opaque2 mutation increases grain lysine, o2 is a transcription factor that regulates a wide network of genes beyond zeins, which leads to pleiotropic and often negative effects. Additionally, the drastic reduction in 19 kDa and 22 kDa alpha-zeins causes a floury kernel, unsuitable for agricultural use. Quality protein maize (QPM) overcame the undesirable kernel texture through the introgression of modifying alleles. However, QPM still lacks a functional o2 transcription factor, which has a penalty on non-lysine amino acids due to the o2 mutation. CRISPR/cas9 gives researchers the ability to directly target genes of interest. In this paper, gene editing was used to specifically target the 19 kDa alpha zein gene family. This allows for proteome rebalancing to occur without an o2 mutation and without a total alpha-zein knockout. The results showed that editing some, but not all, of the 19 kDa zeins resulted in up to 30% more lysine. An edited line displayed an increase of 30% over the wild type. While not quite the 55% lysine increase displayed by QPM, the line had little collateral impact on other amino acid levels compared to QPM. Additionally, the edited line containing a partially reduced 19 kDa showed an advantage in kernel texture that had a complete 19 kDa knockout. These results serve as proof of concept that editing the 19 kDa alpha-zein family alone can enhance lysine while retaining vitreous endosperm and a functional O2 transcription factor.

THP9 enhances seed protein content and nitrogen-use efficiency in maize

Teosinte, the wild ancestor of maize (Zea mays subsp. mays), has three times the seed protein content of most modern inbreds and hybrids, but the mechanisms that are responsible for this trait are unknown^1,2. Here we use trio binning to create a contiguous haplotype DNA sequence of a teosinte (Zea mays subsp. parviglumis) and, through map-based cloning, identify a major high-protein quantitative trait locus, TEOSINTE HIGH PROTEIN 9 (THP9), on chromosome 9. THP9 encodes an asparagine synthetase 4 enzyme that is highly expressed in teosinte, but not in the B73 inbred, in which a deletion in the tenth intron of THP9-B73 causes incorrect splicing of THP9-B73 transcripts. Transgenic expression of THP9-teosinte in B73 significantly increased the seed protein content. Introgression of THP9-teosinte into modern maize inbreds and hybrids greatly enhanced the accumulation of free amino acids, especially asparagine, throughout the plant, and increased seed protein content without affecting yield. THP9-teosinte seems to increase nitrogen-use efficiency, which is important for promoting a high yield under low-nitrogen conditions.

Tandem duplicate expression patterns are conserved between maize haplotypes of the α-zein gene family

Tandem duplication gives rise to copy number variation and subsequent functional novelty among genes as well as diversity between individuals in a species. Functional novelty can result from either divergence in coding sequence or divergence in patterns of gene transcriptional regulation. Here, we investigate conservation and divergence of both gene sequence and gene regulation between the copies of the α-zein gene family in maize inbreds B73 and W22. We used RNA-seq data generated from developing, self-pollinated kernels at three developmental stages timed to coincide with early and peak zein expression. The reference genome annotations for B73 and W22 were modified to ensure accurate inclusion of their respective α-zein gene models to accurately assess copy-specific expression. Expression analysis indicated that although the total expression of α-zeins is higher in W22, the pattern of expression in both lines is conserved. Additional analysis of publicly available RNA-seq data from a diverse population of maize inbreds also demonstrates variation in absolute expression, but conservation of expression patterns across a wide range of maize genotypes and α-zein haplotypes.

Full-length transcript sequencing accelerates the transcriptome research of Gymnocypris namensis, an iconic fish of the Tibetan Plateau

Gymnocypris namensis, the only commercial fish in Namtso Lake of Tibet in China, is rated as nearly threatened species in the Red List of China's Vertebrates. As one of the highest-altitude schizothorax fish in China, G. namensis has strong adaptability to the plateau harsh environment. Although being an indigenous economic fish with high value in research, the biological characterization, genetic diversity, and plateau adaptability of G. namensis are still unclear. Here, we used Pacific Biosciences single molecular real time long read sequencing technology to generate full-length transcripts of G. namensis. Sequences clustering analysis and error correction with Illumina-produced short reads to obtain 319,044 polished isoforms. After removing redundant reads, 125,396 non-redundant isoforms were obtained. Among all transcripts, 103,286 were annotated to public databases. Natural selection has acted on 42 genes for G. namensis, which were enriched on the functions of mismatch repair and Glutathione metabolism. Total 89,736 open reading frames, 95,947 microsatellites, and 21,360 long non-coding RNAs were identified across all transcripts. This is the first study of transcriptome in G. namensis by using PacBio Iso-seq. The acquisition of full-length transcript isoforms might accelerate the transcriptome research of G. namensis and provide basis for further research.

Genomics-Enabled Analysis of Puroindoline b2 Genes Identifies New Alleles in Wheat and Related Triticeae Species

Kernel hardness is a key trait of wheat seeds, largely controlled by two tightly linked genes Puroindoline a and b (Pina and Pinb). Genes homologous to Pinb, namely Pinb2, have been studied. Whether these genes contribute to kernel hardness and other important seed traits remains inconclusive. Using the high-quality bread wheat reference genome, we show that PINB2 are encoded by three homoeologous loci Pinb2 not syntenic to the Hardness locus, with Pinb2-7A locus containing three tandem copies. PINB2 proteins have several features conserved for the Pin/Pinb2 phylogenetic cluster but lack a structural basis of significant impact on kernel hardness. Pinb2 are seed-specifically expressed with varied expression levels between the homoeologous copies and among wheat varieties. Using the high-quality genome information, we developed new Pinb2 allele specific markers and demonstrated their usefulness by 1) identifying new Pinb2 alleles in Triticeae species; and 2) performing an association analysis of Pinb2 with kernel hardness. The association result suggests that Pinb2 genes may have no substantial contribution to kernel hardness. Our results provide new insights into Pinb2 evolution and expression and the new allele-specific markers are useful to further explore Pinb2's contribution to seed traits in wheat.

An update on the maize zein-gene family in the post-genomics era

Maize (Zea mays) is a cereal crop of global food importance. However, the deficiency of essential amino acids, more importantly lysine, methionine and tryptophan, in the major seed storage zein proteins makes corn nutritionally of low value for human consumption. The idea of improving maize nutritional value prompted the search for maize natural mutants harboring low zein contents and higher amount of lysine. These studies resulted in the identification of more than dozens of maize opaque mutants in the previous few decades,o2mutant being the most extensively studied one. However, the high lysine contents but soft kernel texture and chalky endosperm halted the widespread application and commercial success of maize opaque mutants, which ultimately paved the way for the development of Quality Protein Maize (QPM) by modifying the soft endosperm ofo2 mutant into lysine-rich hard endosperm. The previous few decades have witnessed a marked progress in maize zein research. It includes elucidation of molecular mechanism underlying the role of different zein genes in seed endosperm development by cloning different components of zein family, exploring the general organization, function and evolution of zein family members within maize species and among other cereals, and elucidating the cis- and trans-regulatory elements modulating the regulation of different molecular players of maize seed endosperm development. The current advances in high quality reference genomes of maize lines B73 and Mo17 plus the completion of ongoing pan genome sequencing projects of more maize lines with NGS technologies are expected to revolutionize maize zein gene research in near future. This review highlights the recent advances in QPM development and its practical application in the post genomic era, genomic and physical composition and evolution of zein family, and expression, regulation and downstream role of zein genes in endosperm development. Moreover, recent genomic tools and methods developed for functional validation of maize zein genes are also discussed.

Graphical abstract

High frequency DNA rearrangement at qγ27 creates a novel allele for Quality Protein Maize breeding

Copy number variation (CNV) is a major source of genetic variation and often contributes to phenotypic variation in maize. The duplication at the 27-kDa γ-zein locus (qγ27) is essential to convert soft endosperm into hard endosperm in quality protein maize (QPM). This duplication is unstable and generally produces CNV at this locus. We conducted genetic experiments designed to directly measure DNA rearrangement frequencies occurring in males and females of different genetic backgrounds. The average frequency with which the duplication rearranges to single copies is 1.27 × 10-3 and varies among different lines. A triplication of γ27 gene was screened and showed a better potential than the duplication for the future QPM breeding. Our results highlight a novel approach to directly determine the frequency of DNA rearrangements, in this case resulting in CNV at the qγ27 locus. Furthermore, this provides a highly effective way to test suitable parents in QPM breeding.

Modelling BioNano optical data and simulation study of genome map assembly

Motivation

The launch of the BioNano next-generation mapping system has greatly enhanced the performance of physical map construction, thus rapidly expanding the application of optical mapping in genome research. Data biases have profound implications for downstream applications. However, very little is known about the properties and biases of BioNano data, and the very factors that contribute to whole-genome optical map assembly.

Results

We generated BioNano molecule data from eight organisms with diverse base compositions. We first characterized the properties/biases of BioNano molecule data, i.e. molecule length distribution, false labelling signal, variation of optical resolution and coverage distribution bias, and their inducing factors such as chimeric molecules, fragile sites and DNA molecule stretching. Second, we developed the BioNano Molecule SIMulator (BMSIM), a novel computer simulation program for optical data. BMSIM, is of great use for future genome mapping projects. Third, we evaluated the experimental variables that impact whole-genome optical map assembly. Specifically, the effects of coverage depth, molecule length, false-positive and false-negative labelling signals, chimeric molecules and nicking enzyme and nick site density were investigated. Our simulation study provides the empirical findings on how to control experimental variables and gauge analytical parameters to maximize benefit and minimize cost on whole-genome optical map assembly.

Availability and implementation

BMSIM is freely available on: https://github.com/pingchen09990102/BMSIM.

Supplementary information

Supplementary data are available at Bioinformatics online.

Plant evolution and environmental adaptation unveiled by long-read whole-genome sequencing of Spirodela

Aquatic plants have to adapt to the environments distinct from where land plants grow. A critical aspect of adaptation is the dynamics of sequence repeats, not resolved in older sequencing platforms due to incomplete and fragmented genome assemblies from short reads. Therefore, we used PacBio long-read sequencing of the Spirodela polyrhiza genome, reaching a 44-fold increase of contiguity with an N50 (a median of contig lengths) of 831 kb and filling 95.4% of gaps left from the previous version. Reconstruction of repeat regions indicates that sequentially nested long terminal repeat (LTR) retrotranspositions occur early in monocot evolution, featured with both prokaryote-like gene-rich regions and eukaryotic repeat islands. Protein-coding genes are reduced to 18,708 gene models supported by 492,435 high-quality full-length PacBio complementary DNA (cDNA) sequences. Different from land plants, the primitive architecture of Spirodela's adventitious roots and lack of lateral roots and root hairs are consistent with dispensable functions of nutrient absorption. Disease-resistant genes encoding antimicrobial peptides and dirigent proteins are expanded by tandem duplications. Remarkably, disease-resistant genes are not only amplified, but also highly expressed, consistent with low levels of 24-nucleotide (nt) small interfering RNA (siRNA) that silence the immune system of land plants, thereby protecting Spirodela against a wide spectrum of pathogens and pests. The long-read sequence information not only sheds light on plant evolution and adaptation to the environment, but also facilitates applications in bioenergy and phytoremediation.

Nonallelic homologous recombination events responsible for copy number variation within an RNA silencing locus

The structure of chalcone synthase (CHS) gene repeats in different alleles of the I (inhibitor) locus in soybean spawns endogenous RNA interference (RNAi) that leads to phenotypic change in seed coat color of this major agronomic crop. Here, we examined CHS gene copy number by digital PCR and single nucleotide polymorphisms (SNPs) through whole genome resequencing of 15 cultivars that varied in alleles of the I locus (I, ii , ik , and i) that control the pattern distribution of pigments in the seed coats. Lines homozygous for the ii allele had the highest copy number followed by the I and ik cultivars which were more related to each other than to the lines with ii alleles. Some of the recessive i alleles were spontaneous mutations, and each revealed a loss of copy number by digital PCR relative to the parent varieties. Amplicon sequencing and whole genome resequencing determined that the breakpoints of several ii to i mutations resulted from nonallelic homologous recombination (NAHR) events between CHS genes located in segmental duplications leading to large 138-kilobase deletions that erase the structure generating the CHS siRNAs along with eight other non-CHS genes. Functional hybrid CHS genes (designated CHS5:1) were formed in the process and represent rare examples of NAHR in higher plants that have been captured by examining spontaneous mutational events in isogenic mutant lines.

Mass spectrometric characterization of the zein protein composition in maize flour extracts upon protein separation by SDS-PAGE and 2D gel electrophoresis

Highly homogenous α zein protein was isolated from maize kernels in an environment-friendly process using 95% ethanol as solvent. Due to the polyploidy and genetic polymorphism of the plant source, the application of high resolution separation methods in conjunction with precise analytical methods, such as MALDI-TOF-MS, is required to accurately estimate homogeneity of products that contain natural zein protein. The α zein protein product revealed two main bands in SDS-PAGE analysis, one at 25 kDa and other at 20 kDa apparent molecular mass. Yet, high resolution 2DE revealed approximately five protein spot groups in each row, the first at ca. 25 kDa and the second at ca. 20 kDa. Peptide mass fingerprinting data of the proteins in the two dominant SDS-PAGE bands matched to 30 amino acid sequence entries out of 102 non-redundant data base entries. MALDI-TOF-MS peptide mapping of the proteins from all spots indicated the presence of only α zein proteins. The most prominent ion signals in the MALDI mass spectra of the protein mixture of the 25 kDa SDS gel band after in-gel digestion were found at m/z 1272.6 and m/z 2009.1, and the most prominent ion signals of the protein mixture of the 20 kDa band after in-gel digestion were recorded at m/z 1083.5 and m/z 1691.8. These ion signals have been found typical for α zein proteins and may serve as marker ion signals which upon chymotryptic digestion reliably indicate the presence of α zein protein in two hybrid corn products.

NAC-type transcription factors regulate accumulation of starch and protein in maize seeds

Grain starch and protein are synthesized during endosperm development, prompting the question of what regulatory mechanism underlies the synchronization of the accumulation of secondary and primary gene products. We found that two endosperm-specific NAC transcription factors, ZmNAC128 and ZmNAC130, have such a regulatory function. Knockdown of expression of ZmNAC128 and ZmNAC130 with RNA interference (RNAi) caused a shrunken kernel phenotype with significant reduction of starch and protein. We could show that ZmNAC128 and ZmNAC130 regulate the transcription of Bt2 and then reduce its protein level, a rate-limiting step in starch synthesis of maize endosperm. Lack of ZmNAC128 and ZmNAC130 also reduced accumulation of zeins and nonzeins by 18% and 24% compared with nontransgenic siblings, respectively. Although ZmNAC128 and ZmNAC130 affected expression of zein genes in general, they specifically activated transcription of the 16-kDa γ-zein gene. The two transcription factors did not dimerize with each other but exemplified redundancy, whereas individual discovery of their function was not amenable to conventional genetics but illustrated the power of RNAi. Given that both the Bt2 and the 16-kDa γ-zein genes were activated by ZmNAC128 or ZmNAC130, we could identify a core binding site ACGCAA contained within their target promoter regions by combining Dual-Luciferase Reporter and Electrophoretic Mobility Shift assays. Consistent with these properties, transcriptomic profiling uncovered that lack of ZmNAC128 and ZmNAC130 had a pleiotropic effect on the utilization of carbohydrates and amino acids.

Copy Number Variation in Domestication

Domesticated plants have long served as excellent models for studying evolution. Many genes and mutations underlying important domestication traits have been identified, and most causal mutations appear to be SNPs. Copy number variation (CNV) is an important source of genetic variation that has been largely neglected in studies of domestication. Ongoing work demonstrates the importance of CNVs as a source of genetic variation during domestication, and during the diversification of domesticated taxa. Here, we review how CNVs contribute to evolutionary processes underlying domestication, and review examples of domestication traits caused by CNVs. We draw from examples in plant species, but also highlight cases in animal systems that could illuminate the roles of CNVs in the domestication process.

Candidate gene identification of existing or induced mutations with pipelines applicable to large genomes

Bulked segregant analysis (BSA) is used to identify existing or induced variants that are linked to phenotypes. Although it is widely used in Arabidopsis and rice, it remains challenging for crops with large genomes, such as maize. Moreover, analysis of huge data sets can present a bottleneck linking phenotypes to their molecular basis, especially for geneticists without programming experience. Here, we identified two genes of maize defective kernel mutants with newly developed analysis pipelines that require no programing skills and should be applicable to any large genome. In the 1970s, Neuffer and Sheridan generated a chemically induced defective kernel (dek) mutant collection with the potential to uncover critical genes for seed development. To locate such mutations, the dek phenotypes were introgressed into two inbred lines to take advantage of maize haplotype variations and their sequenced genomes. We generated two pipelines that take fastq files derived from next-generation (nextGen) paired-end DNA and cDNA sequencing as input, call on several well established and freely available genomic analysis tools to call SNPs and INDELs, and generate lists of the most likely causal mutations together with variant index plots to locate the mutation to a specific sequence position on a chromosome. The pipelines were validated with a known strawberry mutation before cloning the dek mutants, thereby enabling phenotypic analysis of large genomes by next-generation sequencing.

Uniparental and transgressive expression of α-zeins in maize endosperm of o2 hybrid lines

The α-zein gene family encodes the most abundant storage proteins of maize (Zea mays) endosperm. Members of this family are expressed in a parent-of-origin manner. To characterize this phenomenon further, we investigated the expression of a subset of α-zein polypeptides in reciprocal crosses between o2 lines that were characterized by a simplified α-zein pattern. Maize lines that suppressed the expression of α-zeins when used as female parents were identified. The suppression was cross-specific, occurring only when specific genetic backgrounds were combined. Four α-zein sequences that were sensitive to uniparental expression were isolated. Molecular characterization of these α-zeins confirmed that their expression or suppression depended on the genetic proprieties of the endosperm tissue instead of their parental origin. DNA methylation analysis of both maternally and paternally expressed α-zeins revealed no clear correlation between this epigenetic marker and parent-of-origin allelic expression, suggesting that an additional factor(s) is involved in this process. Genetic analyses revealed that the ability of certain lines to suppress α-zein expression was unstable after one round of heterozygosity with non-suppressing lines. Interestingly, α-zeins also showed a transgressive expression pattern because unexpressed isoforms were reactivated in both F2 and backcross plants. Collectively, our results suggest that parent-of-origin expression of specific α-zein alleles depends on a complex interaction between genotypes in a manner that is reminiscent of paramutation-like phenomena.

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.

Improved Brassica rapa reference genome by single-molecule sequencing and chromosome conformation capture technologies

Brassica rapa comprises several important cultivated vegetables and oil crops. Current reference genome assemblies of Brassica rapa are quite fragmented and not highly contiguous, thereby limiting extensive genetic and genomic analyses. Here, we report an improved assembly of the B. rapa genome (v3.0) using single-molecule sequencing, optical mapping, and chromosome conformation capture technologies (Hi-C). Relative to the previous reference genomes, our assembly features a contig N50 size of 1.45 Mb, representing a ~30-fold improvement. We also identified a new event that occurred in the B. rapa genome ~1.2 million years ago, when a long terminal repeat retrotransposon (LTR-RT) expanded. Further analysis refined the relationship of genome blocks and accurately located the centromeres in the B. rapa genome. The B. rapa genome v3.0 will serve as an important community resource for future genetic and genomic studies in B. rapa. This resource will facilitate breeding efforts in B. rapa, as well as comparative genomic analysis with other Brassica species.

Demography and its effects on genomic variation in crop domestication

Over two thousand plant species have been modified morphologically through cultivation and human use. Here, we review three aspects of crop domestication that are currently undergoing marked revisions, due to analytical advancements and their application to whole genome resequencing (WGS) data. We begin by discussing the duration and demographic history of domestication. There has been debate as to whether domestication occurred quickly or slowly. The latter is tentatively supported both by fossil data and application of WGS data to sequentially Markovian coalescent methods that infer the history of effective population size. This history suggests the possibility of extended human impacts on domesticated lineages prior to their purposeful cultivation. We also make the point that demographic history matters, because it shapes patterns and levels of extant genetic diversity. We illustrate this point by discussing the evolutionary processes that contribute to the empirical observation that most crops examined to date have more putatively deleterious alleles than their wild relatives. These deleterious alleles may contribute to genetic load within crops and may be fitting targets for crop improvement. Finally, the same demographic factors are likely to shape the spectrum of structural variants (SVs) within crops. SVs are known to underlie many of the phenotypic changes associated with domestication and crop improvement, but we currently lack sufficient knowledge about the mechanisms that create SVs, their rates of origin, their population frequencies and their phenotypic effects.

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.

Overexpression of serine acetyltransferase in maize leaves increases seed-specific methionine-rich zeins

Maize kernels do not contain enough of the essential sulphur-amino acid methionine (Met) to serve as a complete diet for animals, even though maize has the genetic capacity to store Met in kernels. Prior studies indicated that the availability of the sulphur (S)-amino acids may limit their incorporation into seed storage proteins. Serine acetyltransferase (SAT) is a key control point for S-assimilation leading to Cys and Met biosynthesis, and SAT overexpression is known to enhance S-assimilation without negative impact on plant growth. Therefore, we overexpressed Arabidopsis thaliana AtSAT1 in maize under control of the leaf bundle sheath cell-specific rbcS1 promoter to determine the impact on seed storage protein expression. The transgenic events exhibited up to 12-fold higher SAT activity without negative impact on growth. S-assimilation was increased in the leaves of SAT overexpressing plants, followed by higher levels of storage protein mRNA and storage proteins, particularly the 10-kDa δ-zein, during endosperm development. This zein is known to impact the level of Met stored in kernels. The elite event with the highest expression of AtSAT1 showed 1.40-fold increase in kernel Met. When fed to chickens, transgenic AtSAT1 kernels significantly increased growth rate compared with the parent maize line. The result demonstrates the efficacy of increasing maize nutritional value by SAT overexpression without apparent yield loss. Maternal overexpression of SAT in vegetative tissues was necessary for high-Met zein accumulation. Moreover, SAT overcomes the shortage of S-amino acids that limits the expression and accumulation of high-Met zeins during kernel development.

Dynamic Evolution of α-Gliadin Prolamin Gene Family in Homeologous Genomes of Hexaploid Wheat

Wheat Gli-2 loci encode complex groups of α-gliadin prolamins that are important for breadmaking, but also major triggers of celiac disease (CD). Elucidation of α-gliadin evolution provides knowledge to produce wheat with better end-use properties and reduced immunogenic potential. The Gli-2 loci contain a large number of tandemly duplicated genes and highly repetitive DNA, making sequence assembly of their genomic regions challenging. Here, we constructed high-quality sequences spanning the three wheat homeologous α-gliadin loci by aligning PacBio-based sequence contigs with BioNano genome maps. A total of 47 α-gliadin genes were identified with only 26 encoding intact full-length protein products. Analyses of α-gliadin loci and phylogenetic tree reconstruction indicate significant duplications of α-gliadin genes in the last ~2.5 million years after the divergence of the A, B and D genomes, supporting its rapid lineage-independent expansion in different Triticeae genomes. We showed that dramatic divergence in expression of α-gliadin genes could not be attributed to sequence variations in the promoter regions. The study also provided insights into the evolution of CD epitopes and identified a single indel event in the hexaploid wheat D genome that likely resulted in the generation of the highly toxic 33-mer CD epitope.

Analysis of Extreme Phenotype Bulk Copy Number Variation (XP-CNV) Identified the Association of rp1 with Resistance to Goss's Wilt of Maize

Goss's wilt (GW) of maize is caused by the Gram-positive bacterium Clavibacter michiganensis subsp. nebraskensis (Cmn) and has spread in recent years throughout the Great Plains, posing a threat to production. The genetic basis of plant resistance is unknown. Here, a simple method for quantifying disease symptoms was developed and used to select cohorts of highly resistant and highly susceptible lines known as extreme phenotypes (XP). Copy number variation (CNV) analyses using whole genome sequences of bulked XP revealed 141 genes containing CNV between the two XP groups. The CNV genes include the previously identified common rust resistant locus rp1. Multiple Rp1 accessions with distinct rp1 haplotypes in an otherwise susceptible accession exhibited hypersensitive responses upon inoculation. GW provides an excellent system for the genetic dissection of diseases caused by closely related subspecies of C. michiganesis. Further work will facilitate breeding strategies to control GW and provide needed insight into the resistance mechanism of important related diseases such as bacterial canker of tomato and bacterial ring rot of potato.

Pan-genome analysis highlights the extent of genomic variation in cultivated and wild rice

The rich genetic diversity in Oryza sativa and Oryza rufipogon serves as the main sources in rice breeding. Large-scale resequencing has been undertaken to discover allelic variants in rice, but much of the information for genetic variation is often lost by direct mapping of short sequence reads onto the O. sativa japonica Nipponbare reference genome. Here we constructed a pan-genome dataset of the O. sativa-O. rufipogon species complex through deep sequencing and de novo assembly of 66 divergent accessions. Intergenomic comparisons identified 23 million sequence variants in the rice genome. This catalog of sequence variations includes many known quantitative trait nucleotides and will be helpful in pinpointing new causal variants that underlie complex traits. In particular, we systemically investigated the whole set of coding genes using this pan-genome data, which revealed extensive presence and absence of variation among rice accessions. This pan-genome resource will further promote evolutionary and functional studies in rice.

Gene Duplication and Evolution Dynamics in the Homeologous Regions Harboring Multiple Prolamin and Resistance Gene Families in Hexaploid Wheat

Improving end-use quality and disease resistance are important goals in wheat breeding. The genetic loci controlling these traits are highly complex, consisting of large families of prolamin and resistance genes with members present in all three homeologous A, B, and D genomes in hexaploid bread wheat. Here, orthologous regions harboring both prolamin and resistance gene loci were reconstructed and compared to understand gene duplication and evolution in different wheat genomes. Comparison of the two orthologous D regions from the hexaploid wheat Chinese Spring and the diploid progenitor Aegilops tauschii revealed their considerable difference due to the presence of five large structural variations with sizes ranging from 100 kb to 2 Mb. As a result, 44% of the Ae. tauschii and 71% of the Chinese Spring sequences in the analyzed regions, including 79 genes, are not shared. Gene rearrangement events, including differential gene duplication and deletion in the A, B, and D regions, have resulted in considerable erosion of gene collinearity in the analyzed regions, suggesting rapid evolution of prolamin and resistance gene families after the separation of the three wheat genomes. We hypothesize that this fast evolution is attributed to the co-evolution of the two gene families dispersed within a high recombination region. The identification of a full set of prolamin genes facilitated transcriptome profiling and revealed that the A genome contributes the least to prolamin expression because of its smaller number of expressed intact genes and their low expression levels, while the B and D genomes contribute similarly.

Genome Sequencing and Assembly by Long Reads in Plants

Plant genomes generated by Sanger and Next Generation Sequencing (NGS) have provided insight into species diversity and evolution. However, Sanger sequencing is limited in its applications due to high cost, labor intensity, and low throughput, while NGS reads are too short to resolve abundant repeats and polyploidy, leading to incomplete or ambiguous assemblies. The advent and improvement of long-read sequencing by Third Generation Sequencing (TGS) methods such as PacBio and Nanopore have shown promise in producing high-quality assemblies for complex genomes. Here, we review the development of sequencing, introducing the application as well as considerations of experimental design in TGS of plant genomes. We also introduce recent revolutionary scaffolding technologies including BioNano, Hi-C, and 10× Genomics. We expect that the informative guidance for genome sequencing and assembly by long reads will benefit the initiation of scientists' projects.

New insights into structural organization and gene duplication in a 1.75-Mb genomic region harboring the α-gliadin gene family in Aegilops tauschii, the source of wheat D genome

Among the wheat prolamins important for its end-use traits, α-gliadins are the most abundant, and are also a major cause of food-related allergies and intolerances. Previous studies of various wheat species estimated that between 25 and 150 α-gliadin genes reside in the Gli-2 locus regions. To better understand the evolution of this complex gene family, the DNA sequence of a 1.75-Mb genomic region spanning the Gli-2 locus was analyzed in the diploid grass, Aegilops tauschii, the ancestral source of D genome in hexaploid bread wheat. Comparison with orthologous regions from rice, sorghum, and Brachypodium revealed rapid and dynamic changes only occurring to the Ae. tauschii Gli-2 region, including insertions of high numbers of non-syntenic genes and a high rate of tandem gene duplications, the latter of which have given rise to 12 copies of α-gliadin genes clustered within a 550-kb region. Among them, five copies have undergone pseudogenization by various mutation events. Insights into the evolutionary relationship of the duplicated α-gliadin genes were obtained from their genomic organization, transcription patterns, transposable element insertions and phylogenetic analyses. An ancestral glutamate-like receptor (GLR) gene encoding putative amino acid sensor in all four grass species has duplicated only in Ae. tauschii and generated three more copies that are interspersed with the α-gliadin genes. Phylogenetic inference and different gene expression patterns support functional divergence of the Ae. tauschii GLR copies after duplication. Our results suggest that the duplicates of α-gliadin and GLR genes have likely taken different evolutionary paths; conservation for the former and neofunctionalization for the latter.

Altered nucleosome positions in maize haplotypes and mutants of a subset of SWI/SNF-like proteins

Chromatin remodelers alter DNA-histone interactions in eukaryotic organisms and have been well characterized in yeast and Arabidopsis. While there are maize proteins with similar domains as known remodelers, the ability of the maize proteins to alter nucleosome position has not been reported. Mutant alleles of several maize proteins (RMR1, CHR101, CHR106, CHR127, and CHR156) with similar functional domains to known chromatin remodelers were identified. Altered gene expression of Chr101, Chr106, Chr127, and Chr156 was demonstrated in plants homozygous for the mutant alleles. These mutant genotypes were subjected to nucleosome position analysis to determine whether misregulation of putative maize chromatin proteins would lead to altered DNA-histone interactions. Nucleosome position changes were observed in plants homozygous for chr101, chr106, chr127, and chr156 mutant alleles, suggesting that CHR101, CHR106, CHR127, and CHR156 may affect chromatin structure. The role of RNA polymerases in altering DNA-histone interactions was also tested. Changes in nucleosome position were demonstrated in homozygous mop2-1 individuals. These changes were demonstrated at the b1 tandem repeats and at newly identified loci. Additionally, differential DNA-histone interactions and altered gene expression of putative chromatin remodelers were demonstrated between different maize haplotypes.

Hybrid assembly of the large and highly repetitive genome of Aegilops tauschii, a progenitor of bread wheat, with the MaSuRCA mega-reads algorithm

Long sequencing reads generated by single-molecule sequencing technology offer the possibility of dramatically improving the contiguity of genome assemblies. The biggest challenge today is that long reads have relatively high error rates, currently around 15%. The high error rates make it difficult to use this data alone, particularly with highly repetitive plant genomes. Errors in the raw data can lead to insertion or deletion errors (indels) in the consensus genome sequence, which in turn create significant problems for downstream analysis; for example, a single indel may shift the reading frame and incorrectly truncate a protein sequence. Here, we describe an algorithm that solves the high error rate problem by combining long, high-error reads with shorter but much more accurate Illumina sequencing reads, whose error rates average <1%. Our hybrid assembly algorithm combines these two types of reads to construct mega-reads, which are both long and accurate, and then assembles the mega-reads using the CABOG assembler, which was designed for long reads. We apply this technique to a large data set of Illumina and PacBio sequences from the species Aegilops tauschii, a large and extremely repetitive plant genome that has resisted previous attempts at assembly. We show that the resulting assembled contigs are far larger than in any previous assembly, with an N50 contig size of 486,807 nucleotides. We compare the contigs to independently produced optical maps to evaluate their large-scale accuracy, and to a set of high-quality bacterial artificial chromosome (BAC)-based assemblies to evaluate base-level accuracy.

Maize defective kernel mutant generated by insertion of a Ds element in a gene encoding a highly conserved TTI2 cochaperone

We have used the newly engineered transposable element Dsg to tag a gene that gives rise to a defective kernel (dek) phenotype. Dsg requires the autonomous element Ac for transposition. Upon excision, it leaves a short DNA footprint that can create in-frame and frameshift insertions in coding sequences. Therefore, we could create alleles of the tagged gene that confirmed causation of the dek phenotype by the Dsg insertion. The mutation, designated dek38-Dsg, is embryonic lethal, has a defective basal endosperm transfer (BETL) layer, and results in a smaller seed with highly underdeveloped endosperm. The maize dek38 gene encodes a TTI2 (Tel2-interacting protein 2) molecular cochaperone. In yeast and mammals, TTI2 associates with two other cochaperones, TEL2 (Telomere maintenance 2) and TTI1 (Tel2-interacting protein 1), to form the triple T complex that regulates DNA damage response. Therefore, we cloned the maize Tel2 and Tti1 homologs and showed that TEL2 can interact with both TTI1 and TTI2 in yeast two-hybrid assays. The three proteins regulate the cellular levels of phosphatidylinositol 3-kinase-related kinases (PIKKs) and localize to the cytoplasm and the nucleus, consistent with known subcellular locations of PIKKs. dek38-Dsg displays reduced pollen transmission, indicating TTI2's importance in male reproductive cell development.

Hybrid assembly of the large and highly repetitive genome of Aegilops tauschii, a progenitor of bread wheat, with the MaSuRCA mega-reads algorithm

Long sequencing reads generated by single-molecule sequencing technology offer the possibility of dramatically improving the contiguity of genome assemblies. The biggest challenge today is that long reads have relatively high error rates, currently around 15%. The high error rates make it difficult to use this data alone, particularly with highly repetitive plant genomes. Errors in the raw data can lead to insertion or deletion errors (indels) in the consensus genome sequence, which in turn create significant problems for downstream analysis; for example, a single indel may shift the reading frame and incorrectly truncate a protein sequence. Here, we describe an algorithm that solves the high error rate problem by combining long, high-error reads with shorter but much more accurate Illumina sequencing reads, whose error rates average <1%. Our hybrid assembly algorithm combines these two types of reads to construct mega-reads, which are both long and accurate, and then assembles the mega-reads using the CABOG assembler, which was designed for long reads. We apply this technique to a large data set of Illumina and PacBio sequences from the species Aegilops tauschii, a large and extremely repetitive plant genome that has resisted previous attempts at assembly. We show that the resulting assembled contigs are far larger than in any previous assembly, with an N50 contig size of 486,807 nucleotides. We compare the contigs to independently produced optical maps to evaluate their large-scale accuracy, and to a set of high-quality bacterial artificial chromosome (BAC)-based assemblies to evaluate base-level accuracy.

CNVs into the wild: screening the genomes of conifer trees (Picea spp.) reveals fewer gene copy number variations in hybrids and links to adaptation

Background

Copy number variations (CNVs) have been linked to different phenotypes in human, including many diseases. A genome-scale understanding of CNVs is available in a few plants but none are wild species, leaving a knowledge gap regarding their genome biology and evolutionary role. We developed a reliable CNV detection method for species lacking contiguous reference genome. We selected multiple probes within 14,078 gene sequences and developed comparative genome hybridization on arrays. Gene CNVs were assessed in three full-sib families from species with 20 Gb genomes, i.e., white and black spruce, and interior spruce - a natural hybrid.

Results

We discovered hundreds of gene CNVs in each species, 3612 in total, which were enriched in functions related to stress and defense responses and narrow expression profiles, indicating a potential role in adaptation. The number of shared CNVs was in accordance with the degree of relatedness between individuals and species. The genetically mapped subset of these genes showed a wide distribution across the genome, implying numerous structural variations. The hybrid family presented significantly fewer CNVs, suggesting that the admixture of two species within one genome reduces the occurrence of CNVs.

Conclusions

The approach we developed is of particular interest in non-model species lacking a reference genome. Our findings point to a role for CNVs in adaptation. Their reduced abundance in the hybrid may limit genetic variability and evolvability of hybrids. We propose that CNVs make a qualitatively distinct contribution to adaptation which could be important for short term change.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-016-3458-8) contains supplementary material, which is available to authorized users.