Progress, challenges and the future of crop genomes

A Brief Overview of the Epigenetic Regulatory Mechanisms in Plants

Plants continuously adapt to their environments by responding to various intrinsic and extrinsic signals. They face numerous biotic and abiotic stresses such as extreme temperatures, drought, or pathogens, requiring complex regulatory mechanisms to control gene activity and adapt their proteome for survival. Epigenetic regulation plays a crucial role in these adaptations, potentially leading to both heritable and non-heritable changes across generations. This process enables plants to adjust their gene expression profiles and acclimate effectively. It is also vital for plant development and productivity, affecting growth, yield, and seed quality, and enabling plants to "remember" environmental stimuli and adapt accordingly. Key epigenetic mechanisms that play significant roles include DNA methylation, histone modification, and ubiquitin ligase complex activity. These processes, which have been extensively studied in the last two decades, have led to a better understanding of the underlying mechanisms and expanded the potential for improving agriculturally and economically important plant traits. DNA methylation is a fundamental process that regulates gene expression by altering chromatin structure. The addition of methyl groups to cytosines by DNA methylases leads to gene suppression, whereas DNA demethylases reverse this effect. Histone modifications, on the other hand, collectively referred to as the "histone code", influence chromatin structure and gene activity by promoting either gene transcription or gene silencing. These modifications are either recognized, added, or removed by a variety of enzymes that act practically as an environmental memory, having a significant impact on plant development and the responses of plants to environmental stimuli. Finally, ubiquitin ligase complexes, which tag specific histones or regulatory proteins with ubiquitin, are also crucial in plant epigenetic regulation. These complexes are involved in protein degradation and play important roles in regulating various cellular activities. The intricate interplay between DNA methylation, histone modifications, and ubiquitin ligases adds complexity to our understanding of epigenetic regulation. These mechanisms collectively control gene expression, generating a complex and branching network of interdependent regulatory pathways. A deeper understanding of this complex network that helps plants adapt to environmental changes and stressful conditions will provide valuable insights into the regulatory mechanisms involved. This knowledge could pave the way for new biotechnological approaches and plant breeding strategies aimed at enhancing crop resilience, productivity, and sustainable agriculture.

Transcriptomic Characterization of Genes Harboring Markers Linked to Maize Yield

BACKGROUND

It is currently believed that breeding priorities, including maize breeding, should focus on introducing varieties with greater utility value, specifically higher yields, into production. Global modern maize breeding relies on various molecular genetics techniques. Using the above mentioned technologies, we can identify regions of the genome that are associated with various phenotypic traits, including yield, which is of fundamental importance for understanding and manipulating these regions.

OBJECTIVES

The aim of the study was to analyze the expression of candidate genes associated with maize yield. To better understand the function of the analyzed genes in increasing maize yield, their expression in different organs and tissues was also assessed using publicly available transcriptome data.

METHODS

RT-qPCR analyses were performed using iTaq Universal SYBR Green Supermix (Bio-Rad, Hercules, CA, USA) and CFX96 Touch Real-Time PCR Detection System (Bio-Rad, Hercules, CA, USA). Each of the performed RT-qPCR experiments consisted of three biological replicates and three technical replicates, the results of which were averaged.

RESULTS

The research results allowed us to select three out of six candidate genes (cinnamoyl-CoA reductase 1-CCR1, aspartate aminotransferase-AAT and sucrose transporter 1-SUT1), which can significantly affect grain yield in maize. Not only our studies but also literature reports clearly indicate the participation of CCR1, AAT and SUT1 in the formation of yield. Identified molecular markers located within these genes can be used in breeding programs to select high yielding maize genotypes.

Dynamics of ER stress-induced gene regulation in plants

Endoplasmic reticulum (ER) stress is a potentially lethal condition that is induced by the abnormal accumulation of unfolded or misfolded secretory proteins in the ER. In eukaryotes, ER stress is managed by the unfolded protein response (UPR) through a tightly regulated, yet highly dynamic, reprogramming of gene transcription. Although the core principles of the UPR are similar across eukaryotes, unique features of the plant UPR reflect the adaptability of plants to their ever-changing environments and the need to balance the demands of growth and development with the response to environmental stressors. The past decades have seen notable progress in understanding the mechanisms underlying ER stress sensing and signalling transduction pathways, implicating the UPR in the effects of physiological and induced ER stress on plant growth and crop yield. Facilitated by sequencing technologies and advances in genetic and genomic resources, recent efforts have driven the discovery of transcriptional regulators and elucidated the mechanisms that mediate the dynamic and precise gene regulation in response to ER stress at the systems level.

Duckweed: Beyond an Efficient Plant Model System

Duckweed (Lemnaceae) rises as a crucial model system due to its unique characteristics and wide-ranging utility. The significance of physiological research and phytoremediation highlights the intricate potential of duckweed in the current era of plant biology. Special attention to duckweed has been brought due to its distinctive features of nutrient uptake, ion transport dynamics, detoxification, intricate signaling, and stress tolerance. In addition, duckweed can alleviate environmental pollutants and enhance sustainability by participating in bioremediation processes and wastewater treatment. Furthermore, insights into the genomic complexity of Lemnaceae species and the flourishing field of transgenic development highlight the opportunities for genetic manipulation and biotechnological innovations. Novel methods for the germplasm conservation of duckweed can be adopted to preserve genetic diversity for future research endeavors and breeding programs. This review centers around prospects in duckweed research promoting interdisciplinary collaborations and technological advancements to drive its full potential as a model organism.

Contrasting patterns of 5S rDNA repeats in European and Asian ecotypes of greater duckweed, Spirodela polyrhiza (Lemnaceae)

Ribosomal DNA (rDNA) contains highly conserved, specifically organized sequences encoding ribosomal RNAs (rRNAs) separated by variable non-transcribed intergenic spacers (NTSs) and is abundant in eukaryotic genomes. These characteristics make the rDNA an informative molecular target to study genome organization, molecular evolution, and phylogenetics. In this study, we characterized the 5S rDNA repeats in the greater duckweed Spiroldela polyrhiza, a species known for its small size, rapid growth, highly conserved genome organization, and low mutation rate. Sequence analysis of at least 12 individually cloned PCR fragments containing the 5S rDNA units for each of six ecotypes that originated from Europe (Ukraine) and Asia (China) revealed two distinct types of 5S rDNA repeats containing NTSs of different lengths and nucleotide compositions. The shorter 5S rDNA repeat units had a highly homogeneous 400-bp NTS, with few ecotype- or region-specific single-nucleotide polymorphisms (SNPs). The longer 5S rDNA units had NTSs of 1056-1084 bp with characteristic intra- and inter-genomic variants due to specific SNPs and insertions/deletions of 4-15-bp DNA elements. We also detected significant variability in the ratio of short/long 5S rDNA variants between ecotypes of S. polyrhiza. The contrasting dynamics of the two types of 5S rDNA units, combined with the unusually low repeat copy number (for plants) in S. polyrhiza (46-220 copies per genome), shows that this species could serve as an excellent model for examining the mechanisms of concerted evolution and functional significance of rDNA variability.

Current challenges and future of agricultural genomes to phenomes in the USA

Dramatic improvements in measuring genetic variation across agriculturally relevant populations (genomics) must be matched by improvements in identifying and measuring relevant trait variation in such populations across many environments (phenomics). Identifying the most critical opportunities and challenges in genome to phenome (G2P) research is the focus of this paper. Previously (Genome Biol, 23(1):1-11, 2022), we laid out how Agricultural Genome to Phenome Initiative (AG2PI) will coordinate activities with USA federal government agencies expand public-private partnerships, and engage with external stakeholders to achieve a shared vision of future the AG2PI. Acting on this latter step, AG2PI organized the "Thinking Big: Visualizing the Future of AG2PI" two-day workshop held September 9-10, 2022, in Ames, Iowa, co-hosted with the United State Department of Agriculture's National Institute of Food and Agriculture (USDA NIFA). During the meeting, attendees were asked to use their experience and curiosity to review the current status of agricultural genome to phenome (AG2P) work and envision the future of the AG2P field. The topic summaries composing this paper are distilled from two 1.5-h small group discussions. Challenges and solutions identified across multiple topics at the workshop were explored. We end our discussion with a vision for the future of agricultural progress, identifying two areas of innovation needed: (1) innovate in genetic improvement methods development and evaluation and (2) innovate in agricultural research processes to solve societal problems. To address these needs, we then provide six specific goals that we recommend be implemented immediately in support of advancing AG2P research.

Characteristics of Amorphophallus konjac as indicated by its genome

Amorphophallus konjac, belonging to the genus Amorphophallus of the Araceae family, is an economically important crop widely used in health products and biomaterials. In the present work, we performed the whole-genome assembly of A. konjac based on the NovaSeq platform sequence data. The final genome assembly was 4.58 Gb with a scaffold N50 of 3212 bp. The genome includes 39,421 protein-coding genes, and 71.75% of the assemblies were repetitive sequences. Comparative genomic analysis showed 1647 gene families have expanded and 2685 contracted in the A. konjac genome. Likewise, genome evolution analysis indicated that A. konjac underwent whole-genome duplication, possibly contributing to the expansion of certain gene families. Furthermore, we identified many candidate genes involved in the tuber formation and development, cellulose and lignification synthesis. The genome of A. konjac obtained in this work provides a valuable resource for the further study of the genetics, genomics, and breeding of this economically important crop, as well as for evolutionary studies of Araceae family.

Genome resources for the elite bread wheat cultivar Aikang 58 and mining of elite homeologous haplotypes for accelerating wheat improvement

Despite recent progress in crop genomics studies, the genomic changes brought about by modern breeding selection are still poorly understood, thus hampering genomics-assisted breeding, especially in polyploid crops with compound genomes such as common wheat (Triticum aestivum). In this work, we constructed genome resources for the modern elite common wheat variety Aikang 58 (AK58). Comparative genomics between AK58 and the landrace cultivar Chinese Spring (CS) shed light on genomic changes that occurred through recent varietal improvement. We also explored subgenome diploidization and divergence in common wheat and developed a homoeologous locus-based genome-wide association study (HGWAS) approach, which was more effective than single homoeolog-based GWAS in unraveling agronomic trait-associated loci. A total of 123 major HGWAS loci were detected using a genetic population derived from AK58 and CS. Elite homoeologous haplotypes (HHs), formed by combinations of subgenomic homoeologs of the associated loci, were found in both parents and progeny, and many could substantially improve wheat yield and related traits. We built a website where users can download genome assembly sequence and annotation data for AK58, perform blast analysis, and run JBrowse. Our work enriches genome resources for wheat, provides new insights into genomic changes during modern wheat improvement, and suggests that efficient mining of elite HHs can make a substantial contribution to genomics-assisted breeding in common wheat and other polyploid crops.

A potential candidate gene associated with the angles of the ear leaf and the second leaf above the ear leaf in maize

BACKGROUND

Leaf angle is a key trait for maize plant architecture that plays a significant role in its morphological development, and ultimately impacting maize grain yield. Although many studies have been conducted on the association and localization of genes regulating leaf angle in maize, most of the candidate genes identified are associated with the regulation of ligule-ear development and phytohormone pathways, and only a few candidate genes have been reported to enhance the mechanical strength of leaf midrib and vascular tissues.

RESULTS

To address this gap, we conducted a genome-wide association study (GWAS) using the leaf angle phenotype and genotyping-by-sequencing data generated from three recombinant inbred line (RIL) populations of maize. Through GWAS analysis, we identified 156 SNPs significantly associated with the leaf angle trait and detected a total of 68 candidate genes located within 10 kb upstream and downstream of these individual SNPs. Among these candidate genes, Zm00001d045408, located on chromosome 9 emerged as a key gene controlling the angles of both the ear leaf and the second leaf above the ear leaf. Notably, this new gene's homolog in Arabidopsis promotes cell division and vascular tissue development. Further analysis revealed that a SNP transversion (G/T) at 7.536 kb downstream of the candidate gene Zm00001d045408 may have caused a reduction in leaf angles of the ear and the second leaf above the ear leaf. Our analysis of the 10 kb region downstream of this candidate gene revealed a 4.337 kb solo long-terminal reverse transcription transposon (solo LTR), located 3.112 kb downstream of Zm00001d045408, with the SNP located 87 bp upstream of the solo LTR.

CONCLUSIONS

In summary, we have identified a novel candidate gene, Zm00001d045408 and a solo LTR that are associated with the angles of both the ear leaf and the second leaf above the ear leaf. The future research holds great potential in exploring the precise role of newly identified candidate gene in leaf angle regulation. Functional characterization of this gene can help in gaining deeper insights into the complex genetic pathways underlying maize plant architecture.

The East Asian wild apples, Malus baccata (L.) Borkh and Malus hupehensis (Pamp.) Rehder., are additional contributors to the genomes of cultivated European and Chinese varieties

The domestication process in long-lived plant perennials differs dramatically from that of annuals, with a huge amount of genetic exchange between crop and wild populations. Though apple is a major fruit crop grown worldwide, the contribution of wild apple species to the genetic makeup of the cultivated apple genome remains a topic of intense study. We used population genomics approaches to investigate the contributions of several wild apple species to European and Chinese rootstock and dessert genomes, with a focus on the extent of wild-crop gene flow. Population genetic structure inferences revealed that the East Asian wild apples, Malus baccata (L.) Borkh and M. hupehensis (Pamp.), form a single panmictic group, and that the European dessert and rootstock apples form a specific gene pool whereas the Chinese dessert and rootstock apples were a mixture of three wild gene pools, suggesting different evolutionary histories of European and Chinese apple varieties. Coalescent-based inferences and gene flow estimates indicated that M. baccata - M. hupehensis contributed to the genome of both European and Chinese cultivated apples through wild-to-crop introgressions, and not as an initial contributor as previously supposed. We also confirmed the contribution through wild-to-crop introgressions of Malus sylvestris Mill. to the cultivated apple genome. Apple tree domestication is therefore one example in woody perennials that involved gene flow from several wild species from multiple geographical areas. This study provides an example of a complex protracted process of domestication in long-lived plant perennials, and is a starting point for apple breeding programmes.

Nanopore sequencing technology and its applications

Since the development of Sanger sequencing in 1977, sequencing technology has played a pivotal role in molecular biology research by enabling the interpretation of biological genetic codes. Today, nanopore sequencing is one of the leading third-generation sequencing technologies. With its long reads, portability, and low cost, nanopore sequencing is widely used in various scientific fields including epidemic prevention and control, disease diagnosis, and animal and plant breeding. Despite initial concerns about high error rates, continuous innovation in sequencing platforms and algorithm analysis technology has effectively addressed its accuracy. During the coronavirus disease (COVID-19) pandemic, nanopore sequencing played a critical role in detecting the severe acute respiratory syndrome coronavirus-2 virus genome and containing the pandemic. However, a lack of understanding of this technology may limit its popularization and application. Nanopore sequencing is poised to become the mainstream choice for preventing and controlling COVID-19 and future epidemics while creating value in other fields such as oncology and botany. This work introduces the contributions of nanopore sequencing during the COVID-19 pandemic to promote public understanding and its use in emerging outbreaks worldwide. We discuss its application in microbial detection, cancer genomes, and plant genomes and summarize strategies to improve its accuracy.

Full-length transcriptome sequencing and comparative transcriptome analysis of Eriocheir sinensis in response to infection by the microsporidian Hepatospora eriocheir

As a new generation of high-throughput sequencing technology, PacBio Iso-Seq technology (Iso-Seq) provides a better alternative sequencing method for the acquisition of full-length unigenes. In this study, a total of 22.27 gigabyte (Gb) subread bases and 128,614 non-redundant unigenes (mean length: 2,324 bp) were obtained from six main tissues of Eriocheir sinensis including the heart, nerve, intestine, muscle, gills and hepatopancreas. In addition, 74,732 unigenes were mapped to at least one of the following databases: Non-Redundant Protein Sequence Database (NR), Gene Ontology (GO), Kyoto Encyclopaedia of Genes and Genomes (KEGG), KEGG Orthology (KO) and Protein family (Pfam). In addition, 6696 transcription factors (TFs), 28,458 long non-coding RNAs (lncRNAs) and 94,230 mRNA-miRNA pairs were identified. Hepatospora eriocheir is the primary pathogen of E. sinensis and can cause hepatopancreatic necrosis disease (HPND); the intestine is the main target tissue. Here, we attempted to identify the key genes related to H. eriocheir infection in the intestines of E. sinensis. By combining Iso-Seq and Illumina RNA-seq analysis, we identified a total of 12,708 differentially expressed unigenes (DEUs; 6,696 upregulated and 6,012 downregulated) in the crab intestine following infection with H. eriocheir. Based on the biological analysis of these DEUs, several key processes were identified, including energy metabolism-related pathways, cell apoptosis and innate immune-related pathways. Twelve selected genes from these DEUs were subsequently verified by quantitative real-time PCR (qRT-PCR) analysis. Our findings enhance our understanding of the E. sinensis transcriptome and the specific association between E. sinensis and H. eriocheir infection.

Development of tobacco rattle virus-based platform for dual heterologous gene expression and CRISPR/Cas reagent delivery

A large number of viral delivery systems have been developed for characterizing functional genes and producing heterologous recombinant proteins in plants, and but most of them are unable to co-express two fusion-free foreign proteins in the whole plant for extended periods of time. In this study, we modified tobacco rattle virus (TRV) as a TRVe dual delivery vector, using the strategy of gene substitution. The reconstructed TRVe had the capability to simultaneously produce two fusion-free foreign proteins at the whole level of Nicotiana benthamiana, and maintained the genetic stability for the insert of double foreign genes. Moreover, TRVe allowed systemic expression of two foreign proteins with the total lengths up to ∼900 aa residues. In addition, Cas12a protein and crRNA were delivered by the TRVe expression system for site-directed editing of genomic DNA in N. benthamiana 16c line constitutively expressing green fluorescent protein (GFP). Taker together, the TRV-based delivery system will be a simple and powerful means to rapidly co-express two non-fused foreign proteins at the whole level and facilitate functional genomics studies in plants.

Insights from the genomes of 4 diploid Camelina spp

Plant evolution has been a complex process involving hybridization and polyploidization making understanding the origin and evolution of a plant's genome challenging even once a published genome is available. The oilseed crop, Camelina sativa (Brassicaceae), has a fully sequenced allohexaploid genome with 3 unknown ancestors. To better understand which extant species best represent the ancestral genomes that contributed to C. sativa's formation, we sequenced and assembled chromosome level draft genomes for 4 diploid members of Camelina: C. neglecta C. hispida var. hispida, C. hispida var. grandiflora, and C. laxa using long and short read data scaffolded with proximity data. We then conducted phylogenetic analyses on regions of synteny and on genes described for Arabidopsis thaliana, from across each nuclear genome and the chloroplasts to examine evolutionary relationships within Camelina and Camelineae. We conclude that C. neglecta is closely related to C. sativa's sub-genome 1 and that C. hispida var. hispida and C. hispida var. grandiflora are most closely related to C. sativa's sub-genome 3. Further, the abundance and density of transposable elements, specifically Helitrons, suggest that the progenitor genome that contributed C. sativa's sub-genome 3 maybe more similar to the genome of C. hispida var. hispida than that of C. hispida var. grandiflora. These diploid genomes show few structural differences when compared to C. sativa's genome indicating little change to chromosome structure following allopolyploidization. This work also indicates that C. neglecta and C. hispida are important resources for understanding the genetics of C. sativa and potential resources for crop improvement.

Sequencing and de novo assembly of the Koshihikari genome and identification of the genomic region related to the eating quality of cooked rice

The japonica rice (Oryza sativa L.) cultivar Koshihikari is considered an important breeding material with good eating quality (EQ). To effectively utilize Koshihikari in molecular breeding programs, determining its whole genome sequence including cultivar-specific segment is crucial. Here, the Koshihikari genome was sequenced using Nanopore and Illumina platforms, and de novo assembly was performed. A highly contiguous Koshihikari genome sequence was compared with Nipponbare, the reference genome of japonica. Genome-wide synteny was observed, as expected, without large structural variations. However, several gaps in alignment were detected on chromosomes 3, 4, 9, and 11. It was notable that previously identified EQ-related QTLs were found in these gaps. Moreover, sequence variations were identified in chromosome 11 at a region flanking the P5 marker, one of the significant markers of good EQ. The Koshihikari-specific P5 region was found to be transmitted through the lineage. High EQ cultivars derived from Koshihikari possessed P5 sequences; on the other hand, Koshihikari-derived low EQ cultivars didn't contain the P5 region, which implies that the P5 genomic region affects the EQ of Koshihikari progenies. The EQ of near-isogenic lines (NILs) of Samnam (a low EQ cultivar) genetic background harboring the P5 segment was improved compared to that of Samnam in Toyo taste value. The structure of the Koshihikari-specific P5 genomic region associated with good EQ was analyzed, which is expected to facilitate the molecular breeding of rice cultivars with superior EQ.

Supplementary Information

The online version contains supplementary material available at 10.1007/s11032-022-01335-3.

A novel full-length transcriptome resource from multiple immune-related tissues in turbot (Scophthalmus maximus) using Pacbio SMART sequencing

Turbot (Scophthalmus maximus) is an important cold-water economic fish. However, the production and development of turbot industry has been constantly hindered by the frequent occurrence of some diseases. Lacking full-length transcriptome for turbot limits immune gene discoveries and gene structures analysis. Therefore, we generated a full-length transcriptome using mixed immune-related tissues of turbot with PacBio Sequel platform. In this study, a total of 31.7 Gb high quality data were generated with the average subreads length of 2618 bp. According to the presence of 5' and 3' primers as well as poly (A) tails, FL (Full-length) and NFL (Non-full-length) isoforms were obtained. Meanwhile, we identified 32,003 non-redundant transcripts, 76.02% of which was novel isoforms of known genes. In addition, 12,176 alternative splicing (AS) events, 6614 polyadenylation (APA) events, 1905 transcription factors, and 2703 lncRNAs were identified. This work is a comprehensive report on the full-length transcriptome of immune-related tissues of turbot, and it also provides valuable molecular resources for future research on the adaptation mechanisms and functional genomics of turbot.

Descriptive Genomic Analysis and Sequence Genotyping of the Two Papaya Species (Vasconcellea pubescens and Vasconcellea chilensis) Using GBS Tools

A genotyping by sequencing (GBS) approach was used to analyze the organization of genetic diversity in V. pubescens and V. chilensis. GBS identified 4675 and 4451 SNPs/INDELs in two papaya species. The cultivated orchards of V. pubescens exhibited scarce genetic diversity and low but significant genetic differentiation. The neutrality test yielded a negative and significant result, suggesting that V. pubescens suffered a selective sweep or a rapid expansion after a bottleneck during domestication. In contrast, V. chilensis exhibited a high level of genetic diversity. The genetic differentiation among the populations was slight, but it was possible to distinguish the two genetic groups. The neutrality test indicated no evidence that natural selection and genetic drift affect the natural population of V. chilensis. Using the Carica papaya genome as a reference, we identified critical SNPs/INDELs associated with putative genes. Most of the identified genes are related to stress responses (salt and nematode) and vegetative and reproductive development. These results will be helpful for future breeding and conservation programs of the Caricaceae family.

Breeding of Vegetable Cowpea for Nutrition and Climate Resilience in Sub-Saharan Africa: Progress, Opportunities, and Challenges

Currently, the world population is increasing, and humanity is facing food and nutritional scarcity. Climate change and variability are a major threat to global food and nutritional security, reducing crop productivity in the tropical and subtropical regions of the globe. Cowpea has the potential to make a significant contribution to global food and nutritional security. In addition, it can be part of a sustainable food system, being a genetic resource for future crop improvement, contributing to resilience and improving agricultural sustainability under climate change conditions. In malnutrition prone regions of sub-Saharan Africa (SSA) countries, cowpea has become a strategic dryland legume crop for addressing food insecurity and malnutrition. Therefore, this review aims to assess the contribution of cowpea to SSA countries as a climate-resilient crop and the existing production challenges and perspectives. Cowpea leaves and immature pods are rich in diverse nutrients, with high levels of protein, vitamins, macro and micronutrients, minerals, fiber, and carbohydrates compared to its grain. In addition, cowpea is truly a multifunctional crop for maintaining good health and for reducing non-communicable human diseases. However, as a leafy vegetable, cowpea has not been researched and promoted sufficiently because it has not been promoted as a food security crop due to its low yield potential, susceptibility to biotic and abiotic stresses, quality assurance issues, policy regulation, and cultural beliefs (it is considered a livestock feed). The development of superior cowpea as a leafy vegetable can be approached in different ways, such as conventional breeding and gene stacking, speed breeding, mutation breeding, space breeding, demand-led breeding, a pan-omics approach, and local government policies. The successful breeding of cowpea genotypes that are high-yielding with a good nutritional value as well as having resistance to biotics and tolerant to abiotic stress could also be used to address food security and malnutrition-related challenges in sub-Saharan Africa.

Plant Proteome Dynamics

Proteins are intimately involved in executing and controlling virtually all cellular processes. To understand the molecular mechanisms that underlie plant phenotypes, it is essential to investigate protein expression, interactions, and modifications, to name a few. The proteome is highly dynamic in time and space, and a plethora of protein modifications, protein interactions, and network constellations are at play under specific conditions and developmental stages. Analysis of proteomes aims to characterize the entire protein complement of a particular cell type, tissue, or organism-a challenging task, given the dynamic nature of the proteome. Modern mass spectrometry-based proteomics technology can be used to address this complexity at a system-wide scale by the global identification and quantification of thousands of proteins. In this review, we present current methods and technologies employed in mass spectrometry-based proteomics and provide examples of dynamic changes in the plant proteome elucidated by proteomic approaches.

The Use of DArTseq Technology to Identify New SNP and SilicoDArT Markers Related to the Yield-Related Traits Components in Maize

In the last decade, many scientists have used molecular biology methods in their research to locate the grain-yield-determining loci and yield structure characteristics in maize. Large-scale molecular analyses in maize do not only focus on the identification of new markers and quantitative trait locus (QTL) regions. DNA analysis in the selection of parental components for heterotic crosses is a very important tool for breeders. The aim of this research was to identify and select new markers for maize (SNP and SilicoDArT) linked to genes influencing the size of the yield components in maize. The plant material used for the research was 186 inbred maize lines. The field experiment was established in twolocations. The yield and six yield components were analyzed. For identification of SNP and SilicoDArT markers related to the yield and yield components, next-generation sequencing was used. As a result of the biometric measurements analysis, differentiation in the average elevation of the analyzed traits for the lines in both locations was found. The above-mentioned results indicate the existence of genotype–environment interactions. The analysis of variance for the observed quality between genotypes indicated a statistically significant differentiation between genotypes and a statistically significant differentiation for all the observed properties betweenlocations. A canonical variable analysis was applied to present a multi-trait assessment of the similarity of the tested maize genotypes in a lower number of dimensions with the lowest possible loss of information. No grouping of lines due to the analyzed was observed. As a result of next-generation sequencing, the molecular markers SilicoDArT (53,031) and SNP (28,571) were obtained. The genetic distance between the analyzed lines was estimated on the basis of these markers. Out of 81,602 identified SilicoDArT and SNP markers, 15,409 (1559 SilicoDArT and 13,850 SNPs) significantly related to the analyzed yield components were selected as a result of association mapping. The greatest numbers of molecular markers were associated with cob length (1203), cob diameter (1759), core length (1201) and core diameter (2326). From 15,409 markers significantly related to the analyzed traits of the yield components, 18 DArT markers were selected, which were significant for the same four traits (cob length, cob diameter, core length, core diameter) in both Kobierzyce and Smolice. These markers were used for physical mapping. As a result of the analyses, it was found that 6 out of 18 (1818; 14,506; 2317; 3233; 11,657; 12,812) identified markers are located inside genes. These markers are located on chromosomes 8, 9, 7, 3, 5, and 1, respectively.

Dissection and ultramicroscopic observation of an apical pollen tube of Pyrus

The pollen tube is ideal for studying cell polar growth, and observing the ultrastructure of the pollen tube tip using transmission electron microscopy (TEM) is the primary method for studying pollen tube growth. The preparation of ultrathin sections of the pollen tube tip sample is important for its successful microscopic observation. The direction of pollen tube growth in vitro is irregular, and it is difficult to dissect the tip of the pollen tube during ultrathin sectioning. Here, we used two methods to efficiently obtain an ultrathin section of the pollen tube tip of Pyrus. In the first method, laser micro-cutting was used to obtain the pollen tube tip, followed by ultrathin sectioning. In the other method, the pollen tubes were cultured in the same growth direction, followed by ultrathin sectioning. Ultrathin sections, which were observed via TEM, showed typical characteristics of the pollen tube tip, such as dense vesicles, numerous mitochondria, and secretory vesicles of the Golgi. We concluded that these two methods are effective in pollen tube tip sample preparation for ultrathin sectioning and provide the foundation for observing the ultrastructure of pollen tube tips.

Diversification, spread, and admixture of octoploid strawberry in the Western Hemisphere

PREMISE

Polyploid species often have complex evolutionary histories that have, until recently, been intractable due to limitations of genomic resources. While recent work has further uncovered the evolutionary history of the octoploid strawberry (Fragaria L.), there are still open questions. Much is unknown about the evolutionary relationship of the wild octoploid species, Fragaria virginiana and Fragaria chiloensis, and gene flow within and among species after the formation of the octoploid genome.

METHODS

We leveraged a collection of wild octoploid ecotypes of strawberry representing the recognized subspecies and ranging from Alaska to southern Chile, and a high-density SNP array to investigate wild octoploid strawberry evolution. Evolutionary relationships were interrogated with phylogenetic analysis and genetic clustering algorithms. Additionally, admixture among and within species is assessed with model-based and tree-based approaches.

RESULTS

Phylogenetic analysis revealed that the two octoploid strawberry species are monophyletic sister lineages. The genetic clustering results show substructure between North and South American F. chiloensis populations. Additionally, model-based and tree-based methods support gene flow within and among the two octoploid species, including newly identified admixture in the Hawaiian F. chiloensis subsp. sandwicensis population.

CONCLUSIONS

F. virginiana and F. chiloensis are supported as monophyletic and sister lineages. All but one of the subspecies show extensive paraphyly. Furthermore, phylogenetic relationships among F. chiloensis populations supports a single population range expansion southward from North America. The inter- and intraspecific relationships of octoploid strawberry are complex and suggest substantial gene flow between sympatric populations among and within species.

A genome sequence resource for the genus Passiflora, the genome of the wild diploid species Passiflora organensis

The genus Passiflora comprises a large group of plants popularly known as passionfruit, much appreciated for their exotic flowers and edible fruits. The species (∼500) are morphologically variable (e.g., growth habit, size, and color of flowers) and are adapted to distinct tropical ecosystems. In this study, we generated the genome of the wild diploid species Passiflora organensis Gardner by adopting a hybrid assembly approach. Passiflora organensis has a small genome of 259 Mbp and a heterozygosity rate of 81%, consistent with its reproductive system. Most of the genome sequences could be integrated into its chromosomes with cytogenomic markers (satellite DNA) as references. The repeated sequences accounted for 58.55% of the total DNA analyzed, and the Tekay lineage was the prevalent retrotransposon. In total, 25,327 coding genes were predicted. Passiflora organensis retains 5,609 singletons and 15,671 gene families. We focused on the genes potentially involved in the locus determining self-incompatibility and the MADS-box gene family, allowing us to infer expansions and contractions within specific subfamilies. Finally, we recovered the organellar DNA. Structural rearrangements and two mitoviruses, besides relics of other mobile elements, were found in the chloroplast and mt-DNA molecules, respectively. This study presents the first draft genome assembly of a wild Passiflora species, providing a valuable sequence resource for genomic and evolutionary studies on the genus, and support for breeding cropped passionfruit species.

Maize Transformation: From Plant Material to the Release of Genetically Modified and Edited Varieties

Over the past decades, advances in plant biotechnology have allowed the development of genetically modified maize varieties that have significantly impacted agricultural management and improved the grain yield worldwide. To date, genetically modified varieties represent 30% of the world's maize cultivated area and incorporate traits such as herbicide, insect and disease resistance, abiotic stress tolerance, high yield, and improved nutritional quality. Maize transformation, which is a prerequisite for genetically modified maize development, is no longer a major bottleneck. Protocols using morphogenic regulators have evolved significantly towards increasing transformation frequency and genotype independence. Emerging technologies using either stable or transient expression and tissue culture-independent methods, such as direct genome editing using RNA-guided endonuclease system as an in vivo desired-target mutator, simultaneous double haploid production and editing/haploid-inducer-mediated genome editing, and pollen transformation, are expected to lead significant progress in maize biotechnology. This review summarises the significant advances in maize transformation protocols, technologies, and applications and discusses the current status, including a pipeline for trait development and regulatory issues related to current and future genetically modified and genetically edited maize varieties.

Gapless indica rice genome reveals synergistic contributions of active transposable elements and segmental duplications to rice genome evolution

The ultimate goal of genome assembly is a high-accuracy gapless genome. Here, we report a new assembly pipeline that is used to produce a gapless genome for the indica rice cultivar Minghui 63. The resulting 397.71-Mb final assembly is composed of 12 contigs with a contig N50 size of 31.93 Mb. Each chromosome is represented by a single contig and the genomic sequences of all chromosomes are gapless. Quality evaluation of this gapless genome assembly showed that gene regions in our assembly have the highest completeness compared with the other 15 reported high-quality rice genomes. Further comparison with the japonica rice genome revealed that the gapless indica genome assembly contains more transposable elements (TEs) and segmental duplications (SDs), the latter of which produce many duplicated genes that can affect agronomic traits through dose effect or sub-/neo-functionalization. The insertion of TEs can also affect the expression of duplicated genes, which may drive the evolution of these genes. Furthermore, we found the expansion of nucleotide-binding site with leucine-rich repeat disease-resistance genes and cis-zeatin-O-glucosyltransferase growth-related genes in SDs in the gapless indica genome assembly, suggesting that SDs contribute to the adaptive evolution of rice disease resistance and developmental processes. Collectively, our findings suggest that active TEs and SDs synergistically contribute to rice genome evolution.

Estimation of Genome Size in the Endemic Species Reseda pentagyna and the Locally Rare Species Reseda lutea Using comparative Analyses of Flow Cytometry and K-Mer Approaches

Genome size is one of the fundamental cytogenetic features of a species, which is critical for the design and initiation of any genome sequencing projects and can provide essential insights in studying taxonomy, cytogenetics, phylogenesis, and evolutionary studies. However, this key cytogenetic information is almost lacking in the endemic species Reseda pentagyna and the locally rare species Reseda lutea in Saudi Arabia. Therefore, genome size was analyzed by propidium iodide PI flow cytometry and compared to k-mer analysis methods. The standard method for genome size measures (flow cytometry) estimated the genome size of R. lutea and R. pentagyna with nuclei isolation MB01 buffer were found to be 1.91 ± 0.02 and 2.09 ± 0.03 pg/2 °C, respectively, which corresponded approximately to a haploid genome size of 934 and 1.022 Mbp, respectively. For validation, K-mer analysis was performed on both species' Illumina paired-end sequencing data from both species. Five k-mer analysis approaches were examined for biocomputational estimation of genome size: A general formula and four well-known programs (CovEST, Kmergenie, FindGSE, and GenomeScope). The parameter preferences had a significant impact on GenomeScope and Kmergenie estimates. While the general formula estimations did not differ considerably, with an average genome size of 867.7 and 896. Mbp. The differences across flow cytometry and biocomputational predictions may be due to the high repeat content, particularly long repetitive regions in both genomes, 71% and 57%, which interfered with k-mer analysis. GenomeScope allowed quantification of high heterozygosity levels (1.04 and 1.37%) of R. lutea and R. pentagyna genomes, respectively. Based on our observations, R. lutea may have a tetraploid genome or higher. Our results revealed fundamental cytogenetic information for R. lutea and R. pentagyna, which should be used in future taxonomic studies and whole-genome sequencing.

Advances in Cereal Crop Genomics for Resilience under Climate Change

Adapting to climate change, providing sufficient human food and nutritional needs, and securing sufficient energy supplies will call for a radical transformation from the current conventional adaptation approaches to more broad-based and transformative alternatives. This entails diversifying the agricultural system and boosting productivity of major cereal crops through development of climate-resilient cultivars that can sustainably maintain higher yields under climate change conditions, expanding our focus to crop wild relatives, and better exploitation of underutilized crop species. This is facilitated by the recent developments in plant genomics, such as advances in genome sequencing, assembly, and annotation, as well as gene editing technologies, which have increased the availability of high-quality reference genomes for various model and non-model plant species. This has necessitated genomics-assisted breeding of crops, including underutilized species, consequently broadening genetic variation of the available germplasm; improving the discovery of novel alleles controlling important agronomic traits; and enhancing creation of new crop cultivars with improved tolerance to biotic and abiotic stresses and superior nutritive quality. Here, therefore, we summarize these recent developments in plant genomics and their application, with particular reference to cereal crops (including underutilized species). Particularly, we discuss genome sequencing approaches, quantitative trait loci (QTL) mapping and genome-wide association (GWAS) studies, directed mutagenesis, plant non-coding RNAs, precise gene editing technologies such as CRISPR-Cas9, and complementation of crop genotyping by crop phenotyping. We then conclude by providing an outlook that, as we step into the future, high-throughput phenotyping, pan-genomics, transposable elements analysis, and machine learning hold much promise for crop improvements related to climate resilience and nutritional superiority.

Full-length transcriptome sequencing from multiple immune-related tissues of Paralichthys olivaceus

Olive flounder (Paralichthys olivaceus) is an important economical flatfish in Japan, Korea and China, but its production has been greatly threatened by various of diseases. Although RNA-seq has provided valuable insights into the host-pathogen interaction, there are still some disadvantages, such as a short sequencing length, the incomplete or inaccurate splicing. Therefore, we generated a full-length transcriptome using mixed immune-related tissues of P. olivaceus with PacBio Sequel platform. In this study, 379,671 full-length non-chimeric (flnc) reads were generated with average length of 2482 bp, which is longer than any previously reported in P. olivaceus. A total of 66,420 isoforms of transcript were identified, 46,850 of which were novel isoforms of known genes accounting for 70.54%. In addition, 7720 novel genes, 12,540 alternative splicing (AS) events, 9296 alternative polyadenylation (APA) events, 2298 transcription factors (TFs), 10,270 lncRNAs and 5400 fusion transcripts were identified. Furthermore, functional annotation showed that most of the full-length transcripts were enriched in immune-related signaling pathways. Otherwise, the mRNA-miRNA interacting networks confirmed that 28.5% of mRNAs were predicted to be targeted by more than one miRNA. These results facilitate the understanding of gene structure, post-transcriptional regulatory networks, and subsequently proteomic diversity. In conclusion, our study provides the full-length transcriptome from multiple immune-related tissues of P. olivaceus, which is valuable for exploring its immune responses.

The high-quality genome of diploid strawberry (Fragaria nilgerrensis) provides new insights into anthocyanin accumulation

Fragaria nilgerrensis is a wild diploid strawberry species endemic to east and southeast region in Asia and provides a rich source of genetic variations for strawberry improvement. Here, we present a chromosome-scale assembly of F. nilgerrensis using single-molecule real-time (SMRT) Pacific Biosciences sequencing and chromosome conformation capture (Hi-C) genome scaffolding. The genome assembly size was 270.3 Mb, with a contig N50 of ∼8.5 Mb. A total of 28 780 genes and 117.2 Mb of transposable elements were annotated for this genome. Next, detailed comparative genomics with the high-quality F. vesca reference genome was conducted to obtain the difference among transposable elements, SNPs, Indels, and so on. The genome size of F. nilgerrensis was enhanced by around 50 Mb relatively to F. vesca, which is mainly due to expansion of transposable elements. In comparison with the F. vesca genome, we identified 4 561 825 SNPs, 846 301 Indels, 4243 inversions, 35 498 translocations and 10 099 relocations. We also found a marked expansion of genes involved in phenylpropanoid biosynthesis, starch and sucrose metabolism, cyanoamino acid metabolism, plant-pathogen interaction, brassinosteroid biosynthesis and plant hormone signal transduction in F. nilgerrensis, which may account for its specific phenotypes and considerable environmental adaptability. Interestingly, we found sequence variations in the upstream regulatory region of FnMYB10, a core transcriptional activator of anthocyanin biosynthesis, resulted in the low expression level of the FnMYB10 gene, which is likely responsible for white fruit phenotype of F. nilgerrensis. The high-quality F. nilgerrensis genome will be a valuable resource for biological research and comparative genomics research.

The Role of Genetic Resources in Breeding for Climate Change: The Case of Public Breeding Programmes in Eighteen Developing Countries

The role of plant breeding in adapting crops to climate changes that affect food production in developing countries is recognized as extremely important and urgent, alongside other agronomic, socio-economic and policy adaptation pathways. To enhance plant breeders' capacity to respond to climate challenges, it is acknowledged that they need to be able to access and use as much genetic diversity as they can get. Through an analysis of data from a global survey, we explore if and how public breeders in selected developing countries are responding to climate challenges through a renewed or innovative use of plant genetic resources, particularly in terms of types of material incorporated into their breeding work as well as sources of such germplasm. It also looks at the possible limitations breeders encounter in their efforts towards exploring diversity for adaptation. Breeders are clearly considering climate challenges. In general, their efforts are aimed at intensifying their breeding work on traits that they were already working on before climate change was so widely discussed. Similarly, the kinds of germplasm they use, and the sources from which they obtain it, do not appear to have changed significantly over the course of recent years. The main challenges breeders faced in accessing germplasm were linked to administrative/legal factors, particularly related to obtaining genetic resources across national borders. They also underscore technical challenges such as a lack of appropriate technologies to exploit germplasm sets such as crop wild relatives and landraces. Addressing these limitations will be crucial to fully enhance the role of public sector breeders in helping to adapt vulnerable agricultural systems to the challenges of climate change.

Detection of 30 bp DNA fragments with a sensitive modified Southern blot analysis

To evaluate crops generated by new breeding techniques, it is important to confirm the removal of recombinant DNAs (rDNAs) derived from foreign genes including unintentionally introduced short rDNA(s). We attempted to develop a sensitive detection method for such short rDNAs using Southern blot analysis and performed a model study targeting single-copy endogenous genes in plants. To increase the detection sensitivity, the general protocol for Southern blot analysis was modified. In the model study, we used endogenous-gene-targeting probes in which complementary sequences were serially replaced by dummy sequences, and detected complementary sequences as well as 30 bp. We further evaluated the sensitivity using short rDNAs derived from GM sequences as pseudoinsertions, and the results demonstrated that rDNA-insertions as small as 30 bp could be detected. The results suggested that unintentionally introduced rDNA-insertions were 30 bp or more in length could be detected by the Southern blot analysis.

Evaluation and management of fungal-infected carrot seeds

Carrot (Daucus carota L.), which is one of the 10 most important vegetable crops worldwide, is an edible root vegetable desired for its taste as well as its medicinal uses. However, a fungus isolated from carrot seeds was observed to substantially decrease the germination rate. The isolate was identified as Alternaria alternata based on morphological and molecular characteristics as well as a phylogenetic tree. The maximum seed infection rate of selected carrot cultivars was approximately 60%, with the main infection site just underneath the seed shell. Additionally, the germination rate of infected seeds decreased by 28.7%. However, the seed infection rate varied among the examined carrot cultivars. Regarding the effects of chemical fungicides, the optimal treatment involved immersing seeds in amistar top suspension concentrate (SC) (effective concentration of 0.65 g/L) for 6 h, which effectively killed the fungi inside the carrot seeds. The results of this study provide a theoretical basis for the development of efficient methods for preventing the infection of carrot seeds by specific fungi and increasing the germination rate and vigour index.

A transcriptome atlas of silkworm silk glands revealed by PacBio single-molecule long-read sequencing

The silk gland of the silkworm Bombyx mori is a specialized organ where silk proteins are efficiently synthesized under precise regulation that largely determines the properties of silk fibers. To understand the genes involved in the regulation of silk protein synthesis, considerable research has focused on the transcripts expressed in silk glands; however, the complete transcriptome profile of this organ has yet to be elucidated. Here, we report a full-length silk gland transcriptome obtained by PacBio single-molecule long-read sequencing technology. In total, 11,697 non-redundant transcripts were identified in mixed samples of silk glands dissected from larvae at five developmental stages. When compared with the published reference, the full-length transcripts optimized the structures of 3002 known genes, and a total of 9061 novel transcripts with an average length of 2171 bp were detected. Among these, 1403 (15.5%) novel transcripts were computationally revealed to be lncRNAs, 8135 (89.8%) novel transcripts were annotated to different protein and nucleotide databases, and 5655 (62.4%) novel transcripts were predicted to have complete ORFs. Furthermore, we found 1867 alternative splicing events, 2529 alternative polyadenylation events, 784 fusion events and 6596 SSRs. This study provides a comprehensive set of reference transcripts and greatly revises and expands the available silkworm transcript data. In addition, these data will be very useful for studying the regulatory mechanisms of silk protein synthesis.

Kmasker plants - a tool for assessing complex sequence space in plant species

Many plant genomes display high levels of repetitive sequences. The assembly of these complex genomes using short high-throughput sequence reads is still a challenging task. Underestimation or disregard of repeat complexity in these datasets can easily misguide downstream analysis. Detection of repetitive regions by k-mer counting methods has proved to be reliable. Easy-to-use applications utilizing k-mer counting are in high demand, especially in the domain of plants. We present Kmasker plants, a tool that uses k-mer count information as an assistant throughout the analytical workflow of genome data that is provided as a command-line and web-based solution. Beside its core competence to screen and mask repetitive sequences, we have integrated features that enable comparative studies between different cultivars or closely related species and methods that estimate target specificity of guide RNAs for application of site-directed mutagenesis using Cas9 endonuclease. In addition, we have set up a web service for Kmasker plants that maintains pre-computed indices for 10 of the economically most important cultivated plants. Source code for Kmasker plants has been made publically available at https://github.com/tschmutzer/kmasker. The web service is accessible at https://kmasker.ipk-gatersleben.de.

Modeling Plant Metabolism: From Network Reconstruction to Mechanistic Models

Mathematical modeling of plant metabolism enables the plant science community to understand the organization of plant metabolism, obtain quantitative insights into metabolic functions, and derive engineering strategies for manipulation of metabolism. Among the various modeling approaches, metabolic pathway analysis can dissect the basic functional modes of subsections of core metabolism, such as photorespiration, and reveal how classical definitions of metabolic pathways have overlapping functionality. In the many studies using constraint-based modeling in plants, numerous computational tools are currently available to analyze large-scale and genome-scale metabolic networks. For ¹³C-metabolic flux analysis, principles of isotopic steady state have been used to study heterotrophic plant tissues, while nonstationary isotope labeling approaches are amenable to the study of photoautotrophic and secondary metabolism. Enzyme kinetic models explore pathways in mechanistic detail, and we discuss different approaches to determine or estimate kinetic parameters. In this review, we describe recent advances and challenges in modeling plant metabolism.

UNBRANCHED3 Expression and Inflorescence Development is Mediated by UNBRANCHED2 and the Distal Enhancer, KRN4, in Maize

Enhancers are cis-acting DNA segments with the ability to increase target gene expression. They show high sensitivity to DNase and contain specific DNA elements in an open chromatin state that allows the binding of transcription factors (TFs). While numerous enhancers are annotated in the maize genome, few have been characterized genetically. KERNEL ROW NUMBER4 (KRN4), an intergenic quantitative trait locus for kernel row number, is assumed to be a cis-regulatory element of UNBRANCHED3 (UB3), a key inflorescence gene. However, the mechanism by which KRN4 controls UB3 expression remains unclear. Here, we found that KRN4 exhibits an open chromatin state, harboring sequences that showed high enhancer activity toward the 35S and UB3 promoters. KRN4 is bound by UB2-centered transcription complexes and interacts with the UB3 promoter by three duplex interactions to affect UB3 expression. Sequence variation at KRN4 enhances ub2 and ub3 mutant ear fasciation. Therefore, we suggest that KRN4 functions as a distal enhancer of the UB3 promoter via chromatin interactions and recruitment of UB2-centered transcription complexes for the fine-tuning of UB3 expression in meristems of ear inflorescences. These results provide evidence that an intergenic region helps to finely tune gene expression, providing a new perspective on the genetic control of quantitative traits.

Genome-wide DNA polymorphisms in four Actinidia arguta genotypes based on whole-genome re-sequencing

Among the genus Actinidia, Actinidia arguta possesses the strongest cold resistance and produces fresh fruit with an intense flavor. To investigate genomic variation that may contribute to variation in phenotypic traits, we performed whole-genome re-sequencing of four A. arguta genotypes originating from different regions in China and identified the polymorphisms using InDel markers. In total, 4,710,650, 4,787,750, 4,646,026, and 4,590,616 SNPs and 1,481,002, 1,534,198, 1,471,304, and 1,425,393 InDels were detected in the ‘Ruby-3’, ‘Yongfeng male’, ‘Kuilv male’, and ‘Hongbei male’ genomes, respectively, compared with the reference genome sequence of cv ‘Hongyang’. A subset of 120 InDels were selected for re-sequencing validation. Additionally, genes related to non-synonymous SNPs and InDels in coding domain sequences were screened for functional analysis. The analysis of GO and KEGG showed that genes involved in cellular responses to water deprivation, sucrose transport, decreased oxygen levels and plant hormone signal transduction were significantly enriched in A. arguta. The results of this study provide insight into the genomic variation of kiwifruit and can inform future research on molecular breeding to improve cold resistance in kiwifruit.

Building near-complete plant genomes

Plant genomes span several orders of magnitude in size, vary in levels of ploidy and heterozygosity, and contain old and recent bursts of transposable elements, which render them challenging but interesting to assemble. Recent advances in single molecule sequencing and physical mapping technologies have enabled high-quality, chromosome scale assemblies of plant species with increasing complexity and size. Single molecule reads can now exceed megabases in length, providing unprecedented opportunities to untangle genomic regions missed by short read technologies. However, polyploid and heterozygous plant genomes are still difficult to assemble but provide opportunities for new tools and approaches. Haplotype phasing, structural variant analysis and de novo pan-genomics are the emerging frontiers in plant genome assembly.

Phased genome sequence of an interspecific hybrid flowering cherry, 'Somei-Yoshino' (Cerasus × yedoensis)

We report the phased genome sequence of an interspecific hybrid, the flowering cherry ‘Somei-Yoshino’ (Cerasus × yedoensis). The sequence data were obtained by single-molecule real-time sequencing technology, split into two subsets based on genome information of the two probable ancestors, and assembled to obtain two haplotype phased genome sequences of the interspecific hybrid. The resultant genome assembly consisting of the two haplotype sequences spanned 690.1 Mb with 4,552 contigs and an N50 length of 1.0 Mb. We predicted 95,076 high-confidence genes, including 94.9% of the core eukaryotic genes. Based on a high-density genetic map, we established a pair of eight pseudomolecule sequences, with highly conserved structures between the two haplotype sequences with 2.4 million sequence variants. A whole genome resequencing analysis of flowering cherries suggested that ‘Somei-Yoshino’ might be derived from a cross between C. spachiana and either C. speciosa or its relatives. A time-course transcriptome analysis of floral buds and flowers suggested comprehensive changes in gene expression in floral bud development towards flowering. These genome and transcriptome data are expected to provide insights into the evolution and cultivation of flowering cherry and the molecular mechanism underlying flowering.

Improved hybrid de novo genome assembly and annotation of African wild rice, Oryza longistaminata, from Illumina and PacBio sequencing reads

African wild rice Oryza longistaminata, one of the eight AA- genome species in the genus Oryza, possesses highly valued traits, such as the rhizomatousness for perennial rice breeding, strong tolerance to biotic and abiotic stresses, and high biomass production on poor soils. To obtain the high-quality reference genome for O. longistaminata we employed a hybrid assembly approach through incorporating Illumina and PacBio sequencing datasets. The final genome assembly comprised only 107 scaffolds and was approximately ∼363.5 Mb, representing ∼92.7% of the estimated African wild rice genome (∼392 Mb). The N50 lengths of the assembled contigs and scaffolds were ∼46.49 Kb and ∼6.83 Mb, indicating ∼3.72-fold and ∼18.8-fold improvement in length compared to the earlier released assembly (∼12.5 Kb and 364 Kb, respectively). Aided with Hi-C data and syntenic relationship with O. sativa, these assembled scaffolds were anchored into 12 pseudo-chromosomes. Genome annotation and comparative genomic analysis reveal that lineage-specific expansion of gene families that respond to biotic- and abiotic stresses are of great potential for mining novel alleles to overcome major diseases and abiotic adaptation in rice breeding programs. This reference genome of African wild rice will greatly enlarge the existing database of rice genome resources and unquestionably form a solid base to understand genomic basis underlying highly valued phenotypic traits and search for novel gene sources in O. longistaminata for the future rice breeding programs.

The genome of jojoba (Simmondsia chinensis): A taxonomically isolated species that directs wax ester accumulation in its seeds

Seeds of the desert shrub, jojoba (Simmondsia chinensis), are an abundant, renewable source of liquid wax esters, which are valued additives in cosmetic products and industrial lubricants. Jojoba is relegated to its own taxonomic family, and there is little genetic information available to elucidate its phylogeny. Here, we report the high-quality, 887-Mb genome of jojoba assembled into 26 chromosomes with 23,490 protein-coding genes. The jojoba genome has only the whole-genome triplication (γ) shared among eudicots and no recent duplications. These genomic resources coupled with extensive transcriptome, proteome, and lipidome data helped to define heterogeneous pathways and machinery for lipid synthesis and storage, provided missing evolutionary history information for this taxonomically segregated dioecious plant species, and will support efforts to improve the agronomic properties of jojoba.

Nanopore sequencing-based genome assembly and evolutionary genomics of circum-basmati rice

Background

The circum-basmati group of cultivated Asian rice (Oryza sativa) contains many iconic varieties and is widespread in the Indian subcontinent. Despite its economic and cultural importance, a high-quality reference genome is currently lacking, and the group’s evolutionary history is not fully resolved. To address these gaps, we use long-read nanopore sequencing and assemble the genomes of two circum-basmati rice varieties.

Results

We generate two high-quality, chromosome-level reference genomes that represent the 12 chromosomes of Oryza. The assemblies show a contig N50 of 6.32 Mb and 10.53 Mb for Basmati 334 and Dom Sufid, respectively. Using our highly contiguous assemblies, we characterize structural variations segregating across circum-basmati genomes. We discover repeat expansions not observed in japonica—the rice group most closely related to circum-basmati—as well as the presence and absence variants of over 20 Mb, one of which is a circum-basmati-specific deletion of a gene regulating awn length. We further detect strong evidence of admixture between the circum-basmati and circum-aus groups. This gene flow has its greatest effect on chromosome 10, causing both structural variation and single-nucleotide polymorphism to deviate from genome-wide history. Lastly, population genomic analysis of 78 circum-basmati varieties shows three major geographically structured genetic groups: Bhutan/Nepal, India/Bangladesh/Myanmar, and Iran/Pakistan.

Conclusion

The availability of high-quality reference genomes allows functional and evolutionary genomic analyses providing genome-wide evidence for gene flow between circum-aus and circum-basmati, describes the nature of circum-basmati structural variation, and reveals the presence/absence variation in this important and iconic rice variety group.

Genome size, chromosome number determination, and analysis of the repetitive elements in Cissus quadrangularis

Cissus quadrangularis (Vitaceae) is a perennial climber endemic to Africa and is characterized by succulent angular stems. The plant grows in arid and semi-arid regions of Africa especially in the African savanna. The stem of C. quadrangularis has a wide range of applications in both human and animal medicine, but there is limited cytogenetic information available for this species. In this study, the chromosome number, genome size, and genome composition for C. quadrangularis were determined. Flow cytometry results indicated that the genome size of C. quadrangularis is approximately 2C = 1.410 pg. Fluorescence microscopy combined with DAPI stain showed the chromosome numbers to be 2n = 48. It is likely that C. quadrangularis has a tetraploid genome after considering the basic chromosome numbers in Cissus genus (n = 10, 11, or 12). A combination of low-throughput genome sequencing and bioinformatics analysis allowed identification and quantification of repetitive elements that make up about 52% of the C. quadrangularis genome, which was dominated by LTR-retrotransposons. Two LTR superfamilies were identified as Copia and Gypsy, with 24% and 15% of the annotated clusters, respectively. The comparison of repeat elements for C. quadrangularis, Vitis vinifera, and four other selected members in the Cissus genus revealed a high diversity in the repetitive element components, which could suggest recent amplification events in the Cissus genus. Our data provides a platform for further studies on the phylogeny and karyotype evolution in this genus and in the family Vitaceae.

Draft genome of Tanacetum cinerariifolium, the natural source of mosquito coil

Pyrethrum (Tanacetum cinerariifolium), which is a perennial Asteraceae plant with white daisy-like flowers, is the original source of mosquito coils and is known for the biosynthesis of the pyrethrin class of natural insecticides. However, the molecular basis of the production of pyrethrins by T. cinerariifolium has yet to be fully elucidated. Here, we present the 7.1-Gb draft genome of T. cinerariifolium, consisting of 2,016,451 scaffolds and 60,080 genes predicted with high confidence. Notably, analyses of transposable elements (TEs) indicated that TEs occupy 33.84% of the genome sequence. Furthermore, TEs of the sire and oryco clades were found to be enriched in the T. cinerariifolium-specific evolutionary lineage, occupying a total of 13% of the genome sequence, a proportion approximately 8-fold higher than that in other plants. InterProScan analysis demonstrated that biodefense-related toxic proteins (e.g., ribosome inactivating proteins), signal transduction-related proteins (e.g., histidine kinases), and metabolic enzymes (e.g., lipoxygenases, acyl-CoA dehydrogenases/oxygenases, and P450s) are also highly enriched in the T. cinerariifolium genome. Molecular phylogenetic analysis detected a variety of enzymes with genus-specific multiplication, including both common enzymes and others that appear to be specific to pyrethrin biosynthesis. Together, these data identify possible novel components of the pyrethrin biosynthesis pathway and provide new insights into the unique genomic features of T. cinerariifolium.

Benchmarking variant identification tools for plant diversity discovery

BACKGROUND

The ability to accurately and comprehensively identify genomic variations is critical for plant studies utilizing high-throughput sequencing. Most bioinformatics tools for processing next-generation sequencing data were originally developed and tested in human studies, raising questions as to their efficacy for plant research. A detailed evaluation of the entire variant calling pipeline, including alignment, variant calling, variant filtering, and imputation was performed on different programs using both simulated and real plant genomic datasets.

RESULTS

A comparison of SOAP2, Bowtie2, and BWA-MEM found that BWA-MEM was consistently able to align the most reads with high accuracy, whereas Bowtie2 had the highest overall accuracy. Comparative results of GATK HaplotypCaller versus SAMtools mpileup indicated that the choice of variant caller affected precision and recall differentially depending on the levels of diversity, sequence coverage and genome complexity. A cross-reference experiment of S. lycopersicum and S. pennellii reference genomes revealed the inadequacy of single reference genome for variant discovery that includes distantly-related plant individuals. Machine-learning-based variant filtering strategy outperformed the traditional hard-cutoff strategy resulting in higher number of true positive variants and fewer false positive variants. A 2-step imputation method, which utilized a set of high-confidence SNPs as the reference panel, showed up to 60% higher accuracy than direct LD-based imputation.

CONCLUSIONS

Programs in the variant discovery pipeline have different performance on plant genomic dataset. Choice of the programs is subjected to the goal of the study and available resources. This study serves as an important guiding information for plant biologists utilizing next-generation sequencing data for diversity characterization and crop improvement.

Are buffalograss (Buchloë dactyloides) cytotypes spatially and ecologically differentiated?

PREMISE

Although autopolyploidy is common among dominant Great Plains grasses, the distribution of cytotypes within a given species is typically poorly understood. This study aims to establish the geographic distribution of cytotypes within buffalograss (Buchloë dactyloides) and to assess whether individual cytotypes have differing ecological tolerances.

METHODS

A range-wide set of 578 B. dactyloides individuals was obtained through field collecting and sampling from herbarium specimens. The cytotype of each sample was estimated by determining allele numbers at 13 simple sequence repeat loci, a strategy that was assessed by comparing estimated to known cytotype in 79 chromosome-counted samples. Ecological differentiation between the dominant tetraploid and hexaploid cytotypes was assessed with analyses of macroclimatic variables.

RESULTS

Simple sequence repeat variation accurately estimated cytotype in 89% of samples from which a chromosome count had been obtained. Applying this approach to samples of unknown ploidy established that diploids and pentaploids are rare, with the common tetraploid and hexaploid cytotypes generally occurring in sites to the north/west (tetraploid) or south/east (hexaploid) portions of the species range. Both MANOVA and niche modeling approaches identified significant but subtle differences in macroclimatic conditions at the set of locations occupied by these two dominant cytotypes.

CONCLUSIONS

Incorporating chromosome count vouchers and cytotype-estimated herbarium records allowed us to perform the largest study of cytotype niche differentiation to date. Buffalograss cytotypes differ greatly in frequency, the common tetraploid and hexaploid cytotypes are non-randomly distributed, and these two cytotypes are subtly ecologically differentiated.

Combining genome-wide and transcriptome-wide analyses reveal the evolutionary conservation and functional diversity of aquaporins in cotton

BACKGROUND

Aquaporins (AQPs) are integral membrane proteins from a larger family of major intrinsic proteins (MIPs) and function in a huge variety of processes such as water transport, plant growth and stress response. The availability of the whole-genome data of different cotton species allows us to study systematic evolution and function of cotton AQPs on a genome-wide level.

RESULTS

Here, a total of 53, 58, 113 and 111 AQP genes were identified in G. arboreum, G. raimondii, G. hirsutum and G. barbadense, respectively. A comprehensive analysis of cotton AQPs, involved in exon/intron structure, functional domains, phylogenetic relationships and gene duplications, divided these AQPs into five subfamilies (PIP, NIP, SIP, TIP and XIP). Comparative genome analysis among 30 species from algae to angiosperm as well as common tandem duplication events in 24 well-studied plants further revealed the evolutionary conservation of AQP family in the organism kingdom. Combining transcriptome analysis and Quantitative Real-time PCR (qRT-PCR) verification, most AQPs exhibited tissue-specific expression patterns both in G. raimondii and G. hirsutum. Meanwhile, a bias of time to peak expression of several AQPs was also detected after treating G. davidsonii and G. hirsutum with 200 mM NaCl. It is interesting that both PIP1;4 h/i/j and PIP2;2a/e showed the highly conserved tandem structure, but differentially contributed to tissue development and stress response in different cotton species.

CONCLUSIONS

These results demonstrated that cotton AQPs were structural conservation while experienced the functional differentiation during the process of evolution and domestication. This study will further broaden our insights into the evolution and functional elucidation of AQP gene family in cotton.

Advances in research on the carrot, an important root vegetable in the Apiaceae family

Carrots (Daucus carota L.), among the most important root vegetables in the Apiaceae family, are cultivated worldwide. The storage root is widely utilized due to its richness in carotenoids, anthocyanins, dietary fiber, vitamins and other nutrients. Carrot extracts, which serve as sources of antioxidants, have important functions in preventing many diseases. The biosynthesis, metabolism, and medicinal properties of carotenoids in carrots have been widely studied. Research on hormone regulation in the growth and development of carrots has also been widely performed. Recently, with the development of high-throughput sequencing technology, many efficient tools have been adopted in carrot research. A large amount of sequence data has been produced and applied to improve carrot breeding. A genome editing system based on CRISPR/Cas9 was also constructed for carrot research. In this review, we will briefly summarize the origins, genetic breeding, resistance breeding, genome editing, omics research, hormone regulation, and nutritional composition of carrots. Perspectives about future research work on carrots are also briefly provided.

Apple whole genome sequences: recent advances and new prospects

In 2010, a major scientific milestone was achieved for tree fruit crops: publication of the first draft whole genome sequence (WGS) for apple (Malus domestica). This WGS, v1.0, was valuable as the initial reference for sequence information, fine mapping, gene discovery, variant discovery, and tool development. A new, high quality apple WGS, GDDH13 v1.1, was released in 2017 and now serves as the reference genome for apple. Over the past decade, these apple WGSs have had an enormous impact on our understanding of apple biological functioning, trait physiology and inheritance, leading to practical applications for improving this highly valued crop. Causal gene identities for phenotypes of fundamental and practical interest can today be discovered much more rapidly. Genome-wide polymorphisms at high genetic resolution are screened efficiently over hundreds to thousands of individuals with new insights into genetic relationships and pedigrees. High-density genetic maps are constructed efficiently and quantitative trait loci for valuable traits are readily associated with positional candidate genes and/or converted into diagnostic tests for breeders. We understand the species, geographical, and genomic origins of domesticated apple more precisely, as well as its relationship to wild relatives. The WGS has turbo-charged application of these classical research steps to crop improvement and drives innovative methods to achieve more durable, environmentally sound, productive, and consumer-desirable apple production. This review includes examples of basic and practical breakthroughs and challenges in using the apple WGSs. Recommendations for "what's next" focus on necessary upgrades to the genome sequence data pool, as well as for use of the data, to reach new frontiers in genomics-based scientific understanding of apple.

Directions for research and training in plant omics: Big Questions and Big Data

A key remit of the NSF-funded "Arabidopsis Research and Training for the 21st Century" (ART-21) Research Coordination Network has been to convene a series of workshops with community members to explore issues concerning research and training in plant biology, including the role that research using Arabidopsis thaliana can play in addressing those issues. A first workshop focused on training needs for bioinformatic and computational approaches in plant biology was held in 2016, and recommendations from that workshop have been published (Friesner et al., Plant Physiology, 175, 2017, 1499). In this white paper, we provide a summary of the discussions and insights arising from the second ART-21 workshop. The second workshop focused on experimental aspects of omics data acquisition and analysis and involved a broad spectrum of participants from academics and industry, ranging from graduate students through post-doctorates, early career and established investigators. Our hope is that this article will inspire beginning and established scientists, corporations, and funding agencies to pursue directions in research and training identified by this workshop, capitalizing on the reference species Arabidopsis thaliana and other valuable plant systems.