Going broad and deep: sequencing-driven insights into plant physiology, evolution, and crop domestication.

Deep sequencing is a term that has become embedded in the plant genomic literature in recent years and with good reason. A torrent of (largely) high-quality genomic and transcriptomic data has been collected and most of this has been publicly released. Indeed, almost 1000 plant genomes have been reported (www.plabipd.de) and the 2000 Plant Transcriptomes Project has long been completed. The EarthBioGenome project will dwarf even these milestones. That said, massive progress in understanding plant physiology, evolution, and crop domestication has been made by sequencing broadly (across a species) as well as deeply (within a single individual). We will outline the current state of the art in genome and transcriptome sequencing before we briefly review the most visible of these broad approaches, namely genome-wide association and transcriptome-wide association studies, as well as the compilation of pangenomes. This will include both (i) the most commonly used methods reliant on single nucleotide polymorphisms and short InDels and (ii) more recent examples which consider structural variants. We will subsequently present case studies exemplifying how their application has brought insight into either plant physiology or evolution and crop domestication. Finally, we will provide conclusions and an outlook as to the perspective for the extension of such approaches to different species, tissues, and biological processes.

ZmWAK02 encoding an RD-WAK protein confers maize resistance against gray leaf spot.

Gray leaf spot (GLS) caused by Cercospora zeina or C. zeae-maydis is a major maize disease throughout the world. Although more than 100 QTLs resistant against GLS have been identified, very few of them have been cloned. Here, we identified a major resistance QTL against GLS, qRglsSB, explaining 58.42% phenotypic variation in SB12×SA101 BC1 F1 population. By fine-mapping, it was narrowed down into a 928 kb region. By using transgenic lines, mutants and complementation lines, it was confirmed that the ZmWAK02 gene, encoding an RD wall-associated kinase, is the responsible gene in qRglsSB resistant against GLS. The introgression of the ZmWAK02 gene into hybrid lines significantly improves their grain yield in the presence of GLS pressure and does not reduce their grain yield in the absence of GLS. In summary, we cloned a gene, ZmWAK02, conferring large effect of GLS resistance and confirmed its great value in maize breeding.

The transcription factor GATA10 regulates fertility conversion of a two-line hybrid tms5 mutant rice via the modulation of UbL40 expression.

The thermosensitive genic male sterile 5 (tms5) mutation causes thermosensitive genic male sterility in rice (Oryza sativa) through loss of RNase ZS1 function, which influences ubiquitin fusion ribosomal protein L40 (UbL40 ) messenger RNA levels during male development. Here, we used ATAC-seq, combined with analysis of H3K9ac and H3K4me2, to identify changes in accessible chromatin during fertility conversion of the two-line hybrid rice Wuxiang S (WXS) derived from a mutant tms5 allele. Furthermore, RNA-seq and bioinformatic analyses identified specific transcription factors (TFs) in differentially accessible chromatin regions. Among these TFs, only GATA10 targeted UbL40 . Osgata10 knockout mutations, which resulted in low expression of UbL40 and a tendency toward male fertility, confirmed that GATA10 regulated fertility conversion via the modulation of UbL40 . Meanwhile, GATA10 acted as a mediator for interactions with ERF65, which revealed that transcriptional regulation is a complex process involving multiple complexes of TFs, namely TF modules. It appears that the ERF141/MADS7/MADS50/MYB modules affect metabolic processes that control anther and pollen development, especially cell wall formation. Our analysis revealed that these modules directly or indirectly affect metabolic pathway-related genes to coordinate plant growth with proper anther development, and furthermore, that GATA10 regulates fertility conversion via the modulation of UbL40 expression.

Improving Nelumbo nucifera genome assemblies using high-resolution genetic maps and BioNano genome mapping reveals ancient chromosome rearrangements.

Genetic and physical maps are powerful tools to anchor fragmented draft genome assemblies generated from next-generation sequencing. Currently, two draft assemblies of Nelumbo nucifera, the genomes of 'China Antique' and 'Chinese Tai-zi', have been released. However, there is presently no information on how the sequences are assembled into chromosomes in N. nucifera. The lack of physical maps and inadequate resolution of available genetic maps hindered the assembly of N. nucifera chromosomes. Here, a linkage map of N. nucifera containing 2371 bin markers [217 577 single nucleotide polymorphisms (SNPs)] was constructed using restriction-site associated DNA sequencing data of 181 F2 individuals and validated by adding 197 simple sequence repeat (SSR) markers. Additionally, a BioNano optical map covering 86.20% of the 'Chinese Tai-zi' genome was constructed. The draft assembly of 'Chinese Tai-zi' was improved based on the BioNano optical map, showing an increase of the scaffold N50 from 0.989 to 1.48 Mb. Using a combination of multiple maps, 97.9% of the scaffolds in the 'Chinese Tai-zi' draft assembly and 97.6% of the scaffolds in the 'China Antique' draft assembly were anchored into pseudo-chromosomes, and the centromere regions along the pseudo-chromosomes were identified. An evolutionary scenario was proposed to reach the modern N. nucifera karyotype from the seven ancestral eudicot chromosomes. The present study provides the highest-resolution linkage map, the optical map and chromosome level genome assemblies for N. nucifera, which are valuable for the breeding and cultivation of N. nucifera and future studies of comparative and evolutionary genomics in angiosperms.

The NnCenH3 protein and centromeric DNA sequence profiles of Nelumbo nucifera Gaertn. (sacred lotus) reveal the DNA structures and dynamics of centromeres in basal eudicots.

Centromeres on eukaryotic chromosomes consist of large arrays of DNA repeats that undergo very rapid evolution. Nelumbo nucifera Gaertn. (sacred lotus) is a phylogenetic relict and an aquatic perennial basal eudicot. Studies concerning the centromeres of this basal eudicot species could provide ancient evolutionary perspectives. In this study, we characterized the centromeric marker protein NnCenH3 (sacred lotus centromere-specific histone H3 variant), and used a chromatin immunoprecipitation (ChIP)-based technique to recover the NnCenH3 nucleosome-associated sequences of sacred lotus. The properties of the centromere-binding protein and DNA sequences revealed notable divergence between sacred lotus and other flowering plants, including the following factors: (i) an NnCenH3 alternative splicing variant comprising only a partial centromere-targeting domain, (ii) active genes with low transcription levels in the NnCenH3 nucleosomal regions, and (iii) the prevalence of the Ty1/copia class of long terminal repeat (LTR) retrotransposons in the centromeres of sacred lotus chromosomes. In addition, the dynamic natures of the centromeric region showed that some of the centromeric repeat DNA sequences originated from telomeric repeats, and a pair of centromeres on the dicentric chromosome 1 was inactive in the metaphase cells of sacred lotus. Our characterization of the properties of centromeric DNA structure within the sacred lotus genome describes a centromeric profile in ancient basal eudicots and might provide evidence of the origins and evolution of centromeres. Furthermore, the identification of centromeric DNA sequences is of great significance for the assembly of the sacred lotus genome.

Molecular evolution and functional characterisation of an ancient phenylalanine ammonia-lyase gene (NnPAL1) from Nelumbo nucifera: novel insight into the evolution of the PAL family in angiosperms.

BACKGROUND: Phenylalanine ammonia-lyase (PAL; E.C.4.3.1.5) is a key enzyme of the phenylpropanoid pathway in plant development, and it catalyses the deamination of phenylalanine to trans-cinnamic acid, leading to the production of secondary metabolites. This enzyme has been identified in many organisms, ranging from prokaryotes to higher plants. Because Nelumbo nucifera is a basal dicot rich in many secondary metabolites, it is a suitable candidate for research on the phenylpropanoid pathway. RESULTS: Three PAL members, NnPAL1, NnPAL2 and NnPAL3, have been identified in N. nucifera using genome-wide analysis. NnPAL1 contains two introns; however, both NnPAL2 and NnPAL3 have only one intron. Molecular and evolutionary analysis of NnPAL1 confirms that it is an ancient PAL member of the angiosperms and may have a different origin. However, PAL clusters, except NnPAL1, are monophyletic after the split between dicots and monocots. These observations suggest that duplication events remain an important occurrence in the evolution of the PAL gene family. Molecular assays demonstrate that the mRNA of the NnPAL1 gene is 2343 bp in size and encodes a 717 amino acid polypeptide. The optimal pH and temperature of the recombinant NnPAL1 protein are 9.0 and 55°C, respectively. The NnPAL1 protein retains both PAL and weak TAL catalytic activities with Km values of 1.07 mM for L-phenylalanine and 3.43 mM for L-tyrosine, respectively. Cis-elements response to environmental stress are identified and confirmed using real-time PCR for treatments with abscisic acid (ABA), indoleacetic acid (IAA), ultraviolet light, Neurospora crassa (fungi) and drought. CONCLUSIONS: We conclude that the angiosperm PAL genes are not derived from a single gene in an ancestral angiosperm genome; therefore, there may be another ancestral duplication and vertical inheritance from the gymnosperms. The different evolutionary histories for PAL genes in angiosperms suggest different mechanisms of functional regulation. The expression patterns of NnPAL1 in response to stress may be necessary for the survival of N. nucifera since the Cretaceous Period. The discovery and characterisation of the ancient NnPAL1 help to elucidate PAL evolution in angiosperms.

A precise chloroplast genome of Nelumbo nucifera (Nelumbonaceae) evaluated with Sanger, Illumina MiSeq, and PacBio RS II sequencing platforms: insight into the plastid evolution of basal eudicots.

BACKGROUND: The chloroplast genome is important for plant development and plant evolution. Nelumbo nucifera is one member of relict plants surviving from the late Cretaceous. Recently, a new sequencing platform PacBio RS II, known as 'SMRT (Single Molecule, Real-Time) sequencing', has been developed. Using the SMRT sequencing to investigate the chloroplast genome of N. nucifera will help to elucidate the plastid evolution of basal eudicots. RESULTS: The sizes of the de novo assembled complete chloroplast genome of N. nucifera were 163,307 bp, 163,747 bp and 163,600 bp with average depths of coverage of 7×, 712× and 105× sequenced by Sanger, Illumina MiSeq and PacBio RS II, respectively. The precise chloroplast genome of N. nucifera was obtained from PacBio RS II data proofread by Illumina MiSeq reads, with a quadripartite structure containing a large single copy region (91,846 bp) and a small single copy region (19,626 bp) separated by two inverted repeat regions (26,064 bp). The genome contains 113 different genes, including four distinct rRNAs, 30 distinct tRNAs and 79 distinct peptide-coding genes. A phylogenetic analysis of 133 taxa from 56 orders indicated that Nelumbo with an age of 177 million years is a sister clade to Platanus, which belongs to the basal eudicots. Basal eudicots began to emerge during the early Jurassic with estimated divergence times at 197 million years using MCMCTree. IR expansions/contractions within the basal eudicots seem to have occurred independently. CONCLUSIONS: Because of long reads and lack of bias in coverage of AT-rich regions, PacBio RS II showed a great promise for highly accurate 'finished' genomes, especially for a de novo assembly of genomes. N. nucifera is one member of basal eudicots, however, evolutionary analyses of IR structural variations of N. nucifera and other basal eudicots suggested that IR expansions/contractions occurred independently in these basal eudicots or were caused by independent insertions and deletions. The precise chloroplast genome of N. nucifera will present new information for structural variation of chloroplast genomes and provide new insight into the evolution of basal eudicots at the primary sequence and structural level.

LAZY4 acts additively with the starch-statolith-dependent gravity-sensing pathway to regulate shoot gravitropism and tiller angle in rice.

Rice tiller angle is a key agronomic trait that has significant effects on the establishment of a high-yield rice population. However, the molecular mechanism underlying the control of rice tiller angle remains to be clarified. Here, we characterized the novel tiller-angle gene LAZY4 (LA4) in rice through map-based cloning. LA4 encodes a C3H2C3-type RING zinc-finger E3 ligase localized in the nucleus, and an in vitro ubiquitination assay revealed that the conserved RING finger domain is essential for its E3 ligase activity. We found that expression of LA4 can be induced by gravistimulation and that loss of LA4 function leads to defective shoot gravitropism caused by impaired asymmetric auxin redistribution upon gravistimulation. Genetic analysis demonstrated that LA4 acts in a distinct pathway from the starch biosynthesis regulators LA2 and LA3, which function in the starch-statolith-dependent pathway. Further genetic analysis showed that LA4 regulates shoot gravitropism and tiller angle by acting upstream of LA1 to mediate lateral auxin transport upon gravistimulation. Our studies reveal that LA4 regulates shoot gravitropism and tiller angle upstream of LA1 through a novel pathway independent of the LA2-LA3-mediated gravity-sensing mechanism, providing new insights into the rice tiller-angle regulatory network.

Two teosintes made modern maize.

The origins of maize were the topic of vigorous debate for nearly a century, but neither the current genetic model nor earlier archaeological models account for the totality of available data, and recent work has highlighted the potential contribution of a wild relative, Zea mays ssp. mexicana. Our population genetic analysis reveals that the origin of modern maize can be traced to an admixture between ancient maize and Zea mays ssp. mexicana in the highlands of Mexico some 4000 years after domestication began. We show that variation in admixture is a key component of maize diversity, both at individual loci and for additive genetic variation underlying agronomic traits. Our results clarify the origin of modern maize and raise new questions about the anthropogenic mechanisms underlying dispersal throughout the Americas.

A pan-Zea genome map for enhancing maize improvement.

BACKGROUND: Maize (Zea mays L.) is at the vanguard facing the upcoming breeding challenges. However, both a super pan-genome for the Zea genus and a comprehensive genetic variation map for maize breeding are still lacking. RESULTS: Here, we construct an approximately 6.71-Gb pan-Zea genome that contains around 4.57-Gb non-B73 reference sequences from fragmented de novo assemblies of 721 pan-Zea individuals. We annotate a total of 58,944 pan-Zea genes and find around 44.34% of them are dispensable in the pan-Zea population. Moreover, 255,821 common structural variations are identified and genotyped in a maize association mapping panel. Further analyses reveal gene presence/absence variants and their potential roles during domestication of maize. Combining genetic analyses with multi-omics data, we demonstrate how structural variants are associated with complex agronomic traits. CONCLUSIONS: Our results highlight the underexplored role of the pan-Zea genome and structural variations to further understand domestication of maize and explore their potential utilization in crop improvement.

Three types of genes underlying the Gametophyte factor1 locus cause unilateral cross incompatibility in maize.

Unilateral cross incompatibility (UCI) occurs between popcorn and dent corn, and represents a critical step towards speciation. It has been reported that ZmGa1P, encoding a pectin methylesterase (PME), is a male determinant of the Ga1 locus. However, the female determinant and the genetic relationship between male and female determinants at this locus are unclear. Here, we report three different types, a total of seven linked genes underlying the Ga1 locus, which control UCI phenotype by independently affecting pollen tube growth in both antagonistic and synergistic manners. These include five pollen-expressed PME genes (ZmGa1Ps-m), a silk-expressed PME gene (ZmPME3), and another silk-expressed gene (ZmPRP3), encoding a pathogenesis-related (PR) proteins. ZmGa1Ps-m confer pollen compatibility. Presence of ZmPME3 causes silk to reject incompatible pollen. ZmPRP3 promotes incompatibility pollen tube growth and thereby breaks the blocking effect of ZmPME3. In addition, evolutionary genomics analyses suggest that the divergence of the Ga1 locus existed before maize domestication and continued during breeding improvement. The knowledge gained here deepen our understanding of the complex regulation of cross incompatibility.

Genome sequencing reveals evidence of adaptive variation in the genus Zea.

Maize is a globally valuable commodity and one of the most extensively studied genetic model organisms. However, we know surprisingly little about the extent and potential utility of the genetic variation found in wild relatives of maize. Here, we characterize a high-density genomic variation map from 744 genomes encompassing maize and all wild taxa of the genus Zea, identifying over 70 million single-nucleotide polymorphisms. The variation map reveals evidence of selection within taxa displaying novel adaptations. We focus on adaptive alleles in highland teosinte and temperate maize, highlighting the key role of flowering-time-related pathways in their adaptation. To show the utility of variants in these data, we generate mutant alleles for two flowering-time candidate genes. This work provides an extensive sampling of the genetic diversity of Zea, resolving questions on evolution and identifying adaptive variants for direct use in modern breeding.

ZEAMAP, a Comprehensive Database Adapted to the Maize Multi-Omics Era.

As one of the most extensively cultivated crops, maize (Zea mays L.) has been extensively studied by researchers and breeders for over a century. With advances in high-throughput detection of various omics data, a wealth of multi-dimensional and multi-omics information has been accumulated for maize and its wild relative, teosinte. Integration of this information has the potential to accelerate genetic research and generate improvements in maize agronomic traits. To this end, we constructed ZEAMAP, a comprehensive database incorporating multiple reference genomes, annotations, comparative genomics, transcriptomes, open chromatin regions, chromatin interactions, high-quality genetic variants, phenotypes, metabolomics, genetic maps, genetic mapping loci, population structures, and populational DNA methylation signals within maize inbred lines. ZEAMAP is user friendly, with the ability to interactively integrate, visualize, and cross-reference multiple different omics datasets.

Genome assembly of a tropical maize inbred line provides insights into structural variation and crop improvement.

Maize is one of the most important crops globally, and it shows remarkable genetic diversity. Knowledge of this diversity could help in crop improvement; however, gold-standard genomes have been elucidated only for modern temperate varieties. Here, we present a high-quality reference genome (contig N50 of 15.78 megabases) of the maize small-kernel inbred line, which is derived from a tropical landrace. Using haplotype maps derived from B73, Mo17 and SK, we identified 80,614 polymorphic structural variants across 521 diverse lines. Approximately 22% of these variants could not be detected by traditional single-nucleotide-polymorphism-based approaches, and some of them could affect gene expression and trait performance. To illustrate the utility of the diverse SK line, we used it to perform map-based cloning of a major effect quantitative trait locus controlling kernel weight-a key trait selected during maize improvement. The underlying candidate gene ZmBARELY ANY MERISTEM1d provides a target for increasing crop yields.

A High Density Genetic Map Derived from RAD Sequencing and Its Application in QTL Analysis of Yield-Related Traits in Vigna unguiculata.

Cowpea [Vigna unguiculata (L.) Walp.] is an annual legume of economic importance and widely grown in the semi-arid tropics. However, high-density genetic maps of cowpea are still lacking. Here, we identified 34,868 SNPs (single nucleotide polymorphisms) that were distributed in the cowpea genome based on the RAD sequencing (restriction-site associated DNA sequencing) technique using a population of 170 individuals (two cowpea parents and 168 F2:3 progenies). Of these, 17,996 reliable SNPs were allotted to 11 consensus linkage groups (LGs). The length of the genetic map was 1,194.25 cM in total with a mean distance of 0.066 cM/SNP marker locus. Using this map and the F2:3 population, combined with the CIM (composite interval mapping) method, eleven quantitative trait loci (QTL) of yield-related trait were detected on seven LGs (LG4, 5, 6, 7, 9, 10, and 11) in cowpea. These QTL explained 0.05-17.32% of the total phenotypic variation. Among these, four QTL were for pod length, four QTL for thousand-grain weight (TGW), two QTL for grain number per pod, and one QTL for carpopodium length. Our results will provide a foundation for understanding genes related to grain yield in the cowpea and genus Vigna.

Development and characterization of polymorphic microRNA-based microsatellite markers in Nelumbo nucifera (Nelumbonaceae).

PREMISE OF THE STUDY: Polymorphic microRNA (miRNA)-based microsatellite markers were developed to investigate the genetic diversity and population structure of Nelumbo nucifera (Nelumbonaceae). METHODS AND RESULTS: A total of 485 miRNA-based microsatellites were found from the genomic DNA sequences of N. nucifera. After several rounds of screening, 21 primer pairs flanking di-, tri-, or pentanucleotide repeats were identified that revealed high levels of genetic diversity in four populations with two to five alleles per locus. The observed and expected heterozygosity per locus ranged from 0.000 to 1.000 and from 0.000 to 0.803, respectively. CONCLUSIONS: The polymorphic microsatellite markers will be useful for studying the genetic diversity and population structure of N. nucifera.

The mitochondrial genome map of Nelumbo nucifera reveals ancient evolutionary features.

Nelumbo nucifera is an evolutionary relic from the Late Cretaceous period. Sequencing the N. nucifera mitochondrial genome is important for elucidating the evolutionary characteristics of basal eudicots. Here, the N. nucifera mitochondrial genome was sequenced using single molecule real-time sequencing technology (SMRT), and the mitochondrial genome map was constructed after de novo assembly and annotation. The results showed that the 524,797-bp N. nucifera mitochondrial genome has a total of 63 genes, including 40 protein-coding genes, three rRNA genes and 20 tRNA genes. Fifteen collinear gene clusters were conserved across different plant species. Approximately 700 RNA editing sites in the protein-coding genes were identified. Positively selected genes were identified with selection pressure analysis. Nineteen chloroplast-derived fragments were identified, and seven tRNAs were derived from the chloroplast. These results suggest that the N. nucifera mitochondrial genome retains evolutionarily conserved characteristics, including ancient gene content and gene clusters, high levels of RNA editing, and low levels of chloroplast-derived fragment insertions. As the first publicly available basal eudicot mitochondrial genome, the N. nucifera mitochondrial genome facilitates further analysis of the characteristics of basal eudicots and provides clues of the evolutionary trajectory from basal angiosperms to advanced eudicots.