Emerging knowledge from genome sequencing of crop species

Entailing the Next-Generation Sequencing and Metabolome for Sustainable Agriculture by Improving Plant Tolerance

Crop production is a serious challenge to provide food for the 10 billion individuals forecasted to live across the globe in 2050. The scientists' emphasize establishing an equilibrium among diversity and quality of crops by enhancing yield to fulfill the increasing demand for food supply sustainably. The exploitation of genetic resources using genomics and metabolomics strategies can help generate resilient plants against stressors in the future. The innovation of the next-generation sequencing (NGS) strategies laid the foundation to unveil various plants' genetic potential and help us to understand the domestication process to unmask the genetic potential among wild-type plants to utilize for crop improvement. Nowadays, NGS is generating massive genomic resources using wild-type and domesticated plants grown under normal and harsh environments to explore the stress regulatory factors and determine the key metabolites. Improved food nutritional value is also the key to eradicating malnutrition problems around the globe, which could be attained by employing the knowledge gained through NGS and metabolomics to achieve suitability in crop yield. Advanced technologies can further enhance our understanding in defining the strategy to obtain a specific phenotype of a crop. Integration among bioinformatic tools and molecular techniques, such as marker-assisted, QTLs mapping, creation of reference genome, de novo genome assembly, pan- and/or super-pan-genomes, etc., will boost breeding programs. The current article provides sequential progress in NGS technologies, a broad application of NGS, enhancement of genetic manipulation resources, and understanding the crop response to stress by producing plant metabolites. The NGS and metabolomics utilization in generating stress-tolerant plants/crops without deteriorating a natural ecosystem is considered a sustainable way to improve agriculture production. This highlighted knowledge also provides useful research that explores the suitable resources for agriculture sustainability.

Genome-Wide Association Study of Soybean Germplasm Derived From Canadian × Chinese Crosses to Mine for Novel Alleles to Improve Seed Yield and Seed Quality Traits

Genome-wide association study (GWAS) has emerged in the past decade as a viable tool for identifying beneficial alleles from a genomic diversity panel. In an ongoing effort to improve soybean [Glycine max (L.) Merr.], which is the third largest field crop in Canada, a GWAS was conducted to identify novel alleles underlying seed yield and seed quality and agronomic traits. The genomic panel consisted of 200 genotypes including lines derived from several generations of bi-parental crosses between modern Canadian × Chinese cultivars (CD-CH). The genomic diversity panel was field evaluated at two field locations in Ontario in 2019 and 2020. Genotyping-by-sequencing (GBS) was conducted and yielded almost 32 K high-quality SNPs. GWAS was conducted using Fixed and random model Circulating Probability Unification (FarmCPU) model on the following traits: seed yield, seed protein concentration, seed oil concentration, plant height, 100 seed weight, days to maturity, and lodging score that allowed to identify five QTL regions controlling seed yield and seed oil and protein content. A candidate gene search identified a putative gene for each of the three traits. The results of this GWAS study provide insight into potentially valuable genetic resources residing in Chinese modern cultivars that breeders may use to further improve soybean seed yield and seed quality traits.

Biotechnological Perspectives of Omics and Genetic Engineering Methods in Alfalfa

For several decades, researchers are working to develop improved major crops with better adaptability and tolerance to environmental stresses. Forage legumes have been widely spread in the world due to their great ecological and economic values. Abiotic and biotic stresses are main factors limiting legume production, however, alfalfa (Medicago sativa L.) shows relatively high level of tolerance to drought and salt stress. Efforts focused on alfalfa improvements have led to the release of cultivars with new traits of agronomic importance such as high yield, better stress tolerance or forage quality. Alfalfa has very high nutritional value due to its efficient symbiotic association with nitrogen-fixing bacteria, while deep root system can help to prevent soil water loss in dry lands. The use of modern biotechnology tools is challenging in alfalfa since full genome, unlike to its close relative barrel medic (Medicago truncatula Gaertn.), was not released yet. Identification, isolation, and improvement of genes involved in abiotic or biotic stress response significantly contributed to the progress of our understanding how crop plants cope with these environmental challenges. In this review, we provide an overview of the progress that has been made in high-throughput sequencing, characterization of genes for abiotic or biotic stress tolerance, gene editing, as well as proteomic and metabolomics techniques bearing biotechnological potential for alfalfa improvement.

Efficient QTL detection of flowering date in a soybean RIL population using the novel restricted two-stage multi-locus GWAS procedure

KEY MESSAGE

Eighty-six R1 QTLs accounting for 89.92% phenotypic variance in a soybean RIL population were identified using RTM-GWAS with SNPLDB marker which performed superior over CIM and MLM-GWAS with BIN/SNPLDB marker. A population (NJRIKY) composed of 427 recombinant inbred lines (RILs) derived from Kefeng-1 × NN1138-2 (MGII × MGV, MG maturity group) was applied for detecting flowering date (R1) quantitative trait locus (QTL) system in soybean. From a low-depth re-sequencing (~ 0.75 ×), 576,874 SNPs were detected and organized into 4737 BINs (recombination breakpoint determinations) and 3683 SNP linkage disequilibrium blocks (SNPLDBs), respectively. Using the association mapping procedures "Restricted Two-stage Multi-locus Genome-wide Association Study" (RTM-GWAS), "Mixed Linear Model Genome-wide Association Study" (MLM-GWAS) and the linkage mapping procedure "Composite Interval Mapping" (CIM), 67, 36 and 10 BIN-QTLs and 86, 14 and 23 SNPLDB-QTLs were detected with their phenotypic variance explained (PVE) 88.70-89.92% (within heritability 98.2%), 146.41-353.62% (overflowing) and 88.29-172.34% (overflowing), respectively. The RTM-GWAS with SNPLDBs which showed to be more efficient and reasonable than the others was used to identify the R1 QTL system in NJRIKY. The detected 86 SNPLDB-QTLs with their PVE from 0.02 to 30.66% in a total of 89.92% covered 51 out of 104 R1 QTLs in 18 crosses in SoyBase and 26 out of 139 QTLs in a nested association mapping population, while the rest 29 QTLs were novel ones. From the QTL system, 52 candidate genes were annotated, including the verified gene E1, E2, E9 and J, and grouped into 3 categories of biological processes, among which 24 genes were enriched into three protein-protein interaction networks, suggesting gene networks working together. Since NJRIKY involves only MGII and MGV, the QTL/gene system among MG000-MGX should be explored further.

Computational pan-genomics: status, promises and challenges

Many disciplines, from human genetics and oncology to plant breeding, microbiology and virology, commonly face the challenge of analyzing rapidly increasing numbers of genomes. In case of Homo sapiens, the number of sequenced genomes will approach hundreds of thousands in the next few years. Simply scaling up established bioinformatics pipelines will not be sufficient for leveraging the full potential of such rich genomic data sets. Instead, novel, qualitatively different computational methods and paradigms are needed. We will witness the rapid extension of computational pan-genomics, a new sub-area of research in computational biology. In this article, we generalize existing definitions and understand a pan-genome as any collection of genomic sequences to be analyzed jointly or to be used as a reference. We examine already available approaches to construct and use pan-genomes, discuss the potential benefits of future technologies and methodologies and review open challenges from the vantage point of the above-mentioned biological disciplines. As a prominent example for a computational paradigm shift, we particularly highlight the transition from the representation of reference genomes as strings to representations as graphs. We outline how this and other challenges from different application domains translate into common computational problems, point out relevant bioinformatics techniques and identify open problems in computer science. With this review, we aim to increase awareness that a joint approach to computational pan-genomics can help address many of the problems currently faced in various domains.

Advances in genome editing for improved animal breeding: A review

Since centuries, the traits for production and disease resistance are being targeted while improving the genetic merit of domestic animals, using conventional breeding programs such as inbreeding, outbreeding, or introduction of marker-assisted selection. The arrival of new scientific concepts, such as cloning and genome engineering, has added a new and promising research dimension to the existing animal breeding programs. Development of genome editing technologies such as transcription activator-like effector nuclease, zinc finger nuclease, and clustered regularly interspaced short palindromic repeats systems begun a fresh era of genome editing, through which any change in the genome, including specific DNA sequence or indels, can be made with unprecedented precision and specificity. Furthermore, it offers an opportunity of intensification in the frequency of desirable alleles in an animal population through gene-edited individuals more rapidly than conventional breeding. The specific research is evolving swiftly with a focus on improvement of economically important animal species or their traits all of which form an important subject of this review. It also discusses the hurdles to commercialization of these techniques despite several patent applications owing to the ambiguous legal status of genome-editing methods on account of their disputed classification. Nonetheless, barring ethical concerns gene-editing entailing economically important genes offers a tremendous potential for breeding animals with desirable traits.

Toward understanding of the high number of plant aquaporin isoforms and multiple regulation mechanisms

Since the discovery of the first plant aquaporin (AQP) in 1993, our conception of the way plants control cell water homeostasis as well as their global water balance has been revisited. Plant AQPs constitute a large family of evolutionarily related channels that, in addition to water, can also facilitate the membrane diffusion of a number of small solutes, such as urea, CO₂, H₂O₂, ammonia, metalloids, and even ions, indicating a wide range of cellular functions. At the cellular level, AQPs are subject to various regulation mechanisms leading to active/inactive channels in their target membranes. In this review, we discuss several specific questions that need to be addressed in future research. Why are so many different AQPs simultaneously expressed in specific cellular types? How is their selectivity to different solutes controlled (in particular in the case of multiple permeation properties)? What does the molecular interaction between AQPs and other molecules tell us about their regulation and their involvement in specific cellular and physiological processes? Resolving these questions will definitely help us better understand the physiological advantages that plants have to express and regulate so many AQP isoforms.

Leaf Transcriptome Sequencing for Identifying Genic-SSR Markers and SNP Heterozygosity in Crossbred Mango Variety 'Amrapali' (Mangifera indica L.)

Mango (Mangifera indica L.) is called “king of fruits” due to its sweetness, richness of taste, diversity, large production volume and a variety of end usage. Despite its huge economic importance genomic resources in mango are scarce and genetics of useful horticultural traits are poorly understood. Here we generated deep coverage leaf RNA sequence data for mango parental varieties ‘Neelam’, ‘Dashehari’ and their hybrid ‘Amrapali’ using next generation sequencing technologies. De-novo sequence assembly generated 27,528, 20,771 and 35,182 transcripts for the three genotypes, respectively. The transcripts were further assembled into a non-redundant set of 70,057 unigenes that were used for SSR and SNP identification and annotation. Total 5,465 SSR loci were identified in 4,912 unigenes with 288 type I SSR (n ≥ 20 bp). One hundred type I SSR markers were randomly selected of which 43 yielded PCR amplicons of expected size in the first round of validation and were designated as validated genic-SSR markers. Further, 22,306 SNPs were identified by aligning high quality sequence reads of the three mango varieties to the reference unigene set, revealing significantly enhanced SNP heterozygosity in the hybrid Amrapali. The present study on leaf RNA sequencing of mango varieties and their hybrid provides useful genomic resource for genetic improvement of mango.

Next generation breeding

The genomic revolution of the past decade has greatly improved our understanding of the genetic make-up of living organisms. The sequencing of crop genomes has completely changed our vision and interpretation of genome organization and evolution. Re-sequencing allows the identification of an unlimited number of markers as well as the analysis of germplasm allelic diversity based on allele mining approaches. High throughput marker technologies coupled with advanced phenotyping platforms provide new opportunities for discovering marker-trait associations which can sustain genomic-assisted breeding. The availability of genome sequencing information is enabling genome editing (site-specific mutagenesis), to obtain gene sequences desired by breeders. This review illustrates how next generation sequencing-derived information can be used to tailor genomic tools for different breeders' needs to revolutionize crop improvement.

Genetic structure of Miscanthus sinensis and Miscanthus sacchariflorus in Japan indicates a gradient of bidirectional but asymmetric introgression

Unilateral introgression from diploids to tetraploids has been hypothesized to be an important evolutionary mechanism in plants. However, few examples have been definitively identified, perhaps because data of sufficient depth and breadth were difficult to obtain before the advent of affordable high-density genotyping. Throughout Japan, tetraploid Miscanthus sacchariflorus and diploid Miscanthus sinensis are common, and occasionally hybridize. In this study, 667 M. sinensis and 78 M. sacchariflorus genotypes from Japan were characterized using 20 704 SNPs and ten plastid microsatellites. Similarity of SNP genotypes between diploid and tetraploid M. sacchariflorus indicated that the tetraploids originated through autopolyploidy. Structure analysis indicated a gradient of introgression from diploid M. sinensis into tetraploid M. sacchariflorus throughout Japan; most tetraploids had some M. sinensis DNA. Among phenotypically M. sacchariflorus tetraploids, M. sinensis ancestry averaged 7% and ranged from 1-39%, with introgression greatest in southern Japan. Unexpectedly, rare (~1%) diploid M. sinensis individuals from northern Japan were found with 6-27% M. sacchariflorus ancestry. Population structure of M. sinensis in Japan included three groups, and was driven primarily by distance, and secondarily by geographic barriers such as mountains and straits. Miscanthus speciation is a complex and dynamic process. In contrast to limited introgression between diploid M. sacchariflorus and M. sinensis in northern China, selection for adaptation to a moderate maritime climate probably favoured cross-ploidy introgressants in southern Japan. These results will help guide the selection of Miscanthus accessions for the breeding of biomass cultivars.

Identification of candidate genes for drought tolerance by whole-genome resequencing in maize

BACKGROUND

Drought stress is one of the major limiting factors for maize production. With the availability of maize B73 reference genome and whole-genome resequencing of 15 maize inbreds, common variants (CV) and clustering analyses were applied to identify non-synonymous SNPs (nsSNPs) and corresponding candidate genes for drought tolerance.

RESULTS

A total of 524 nsSNPs that were associated with 271 candidate genes involved in plant hormone regulation, carbohydrate and sugar metabolism, signaling molecules regulation, redox reaction and acclimation of photosynthesis to environment were detected by CV and cluster analyses. Most of the nsSNPs identified were clustered in bin 1.07 region that harbored six previously reported QTL with relatively high phenotypic variation explained for drought tolerance. Genes Ontology (GO) analysis of candidate genes revealed that there were 35 GO terms related to biotic stimulus and membrane-bounded organelle, showing significant differences between the candidate genes and the reference B73 background. Changes of expression level in these candidate genes for drought tolerance were detected using RNA sequencing for fertilized ovary, basal leaf meristem tissue and roots collected under drought stressed and well-watered conditions. The results indicated that 70% of candidate genes showed significantly expression changes under two water treatments and our strategies for mining candidate genes are feasible and relatively efficient.

CONCLUSIONS

Our results successfully revealed candidate nsSNPs and associated genes for drought tolerance by comparative sequence analysis of 16 maize inbred lines. Both methods we applied were proved to be efficient for identifying candidate genes for complex traits through the next-generation sequencing technologies (NGS). These selected genes will not only facilitate understanding of genetic basis of drought stress response, but also accelerate genetic improvement through marker-assisted selection in maize.

Genomic resources in fruit plants: an assessment of current status

The availability of many genomic resources such as genome sequences, functional genomics resources including microarrays and RNA-seq, sufficient numbers of molecular markers, express sequence tags (ESTs) and high-density genetic maps is causing a rapid acceleration of genetics and genomic research of many fruit plants. This is leading to an increase in our knowledge of the genes that are linked to many horticultural and agronomically important traits. Recently, some progress has also been made on the identification and functional analysis of miRNAs in some fruit plants. This is one of the most active research fields in plant sciences. The last decade has witnessed development of genomic resources in many fruit plants such as apple, banana, citrus, grapes, papaya, pears, strawberry etc.; however, many of them are still not being exploited. Furthermore, owing to lack of resources, infrastructure and research facilities in many lesser-developed countries, development of genomic resources in many underutilized or less-studied fruit crops, which grow in these countries, is limited. Thus, research emphasis should be given to those fruit crops for which genomic resources are relatively scarce. The development of genomic databases of these less-studied fruit crops will enable biotechnologists to identify target genes that underlie key horticultural and agronomical traits. This review presents an overview of the current status of the development of genomic resources in fruit plants with the main emphasis being on genome sequencing, EST resources, functional genomics resources including microarray and RNA-seq, identification of quantitative trait loci and construction of genetic maps as well as efforts made on the identification and functional analysis of miRNAs in fruit plants.

DOG1-like genes in cereals: investigation of their function by means of ectopic expression in Arabidopsis

The Arabidopsis gene DOG1 (AtDOG1) functions in seed dormancy and in sugar signaling. Little is known about the structural and functional features of plant genes homologous to AtDOG1, except for one type (clade 1) of Triticeae AtDOG1-like genes, which was previously demonstrated to be functionally orthologous to AtDOG1. Here, through phylogenetic, structural, and functional analyses of cereal AtDOG1-like genes, we characterized their features: these genes exist as a gene family that can be classified into five distinct clades (1-5). Of these, AtDOG1-like genes in clades 1-4 have a similar architecture to AtDOG1: they encode proteins with three conserved regions. In contrast, the clade 5 genes are distinct; their encoded proteins lack these conserved regions, but harbor domains that interact with DNA. Ectopic expression of the cereal AtDOG1-like genes of clades 2-4 in Arabidopsis demonstrated that like the clade 1 genes, they performed the same function as AtDOG1. The correlation between the depth of seed dormancy and the efficiency of sugar signaling in transgenic Arabidopsis conferred by genes in clades 1-4 suggests a close link in the underlying mechanisms between the seed dormancy and sugar signaling functions of AtDOG1.

Asymmetric functional divergence of young, dispersed gene duplicates in Arabidopsis thaliana

One prediction of the classic Ohno model of gene duplication predicts that new genes form from the asymmetric functional divergence of a newly arisen, redundant duplicate locus. In order to understand the mechanisms which give rise to functional divergence of newly formed dispersed duplicates, we assessed the expression and molecular evolutionary divergence of a suite of 19 highly similar dispersed duplicates in Arabidopsis thaliana. These duplicates have a K sil equal to or less than 5 % and are specific to the A. thaliana lineage; thus, they predictably represent some of the youngest duplicates in the A. thaliana genome. We found that the majority of young duplicate loci exhibit asymmetric expression patterns, with the daughter locus exhibiting reduced expression across all tissues analyzed relative to the progenitor locus or simply not expressed. Furthermore, daughter loci, on the whole, have significantly more nonsynonymous substitutions than the progenitor loci. We also identified four pairs of loci which exhibit significant (P < 0.05) evolutionary rate asymmetry, three of which exhibit elevated dN/dS in the duplicate copy. We suggest, based on these data, that functional diversification initially takes the form of asymmetric regulatory divergence that can be a direct consequence of the mode of duplication. The reduced and/or absence of expression in the daughter copy relaxes functional constraint on its protein coding sequence leading to the asymmetric accumulation of nonsynonymous mutations. Thus, our data both affirm Ohno's prediction while explaining the mechanism by which functional divergence initially occurs following duplication for dispersed gene duplicates.

Multiple origins of the determinate growth habit in domesticated common bean (Phaseolus vulgaris)

BACKGROUND AND AIMS

The actual number of domestications of a crop is one of the key questions in domestication studies. Answers to this question have generally been based on relationships between wild progenitors and domesticated descendants determined with anonymous molecular markers. In this study, this question was investigated by determining the number of instances a domestication phenotype had been selected in a crop species. One of the traits that appeared during domestication of common bean (Phaseolus vulgaris) is determinacy, in which stems end with a terminal inflorescence. It has been shown earlier that a homologue of the arabidopsis TFL1 gene - PvTFL1y - controls determinacy in a naturally occurring variation of common bean.

METHODS

Sequence variation was analysed for PvTFL1y in a sample of 46 wild and domesticated accessions that included determinate and indeterminate accessions.

KEY RESULTS

Indeterminate types - wild and domesticated - showed only synonymous nucleotide substitutions. Determinate types - observed only among domesticated accessions - showed, in addition to synonymous substitutions, non-synonymous substitutions, indels, a putative intron-splicing failure, a retrotransposon insertion and a deletion of the entire locus. The retrotransposon insertion was observed in 70 % of determinate cultivars, in the Americas and elsewhere. Other determinate mutants had a more restricted distribution in the Americas only, either in the Andean or in the Mesoamerican gene pool of common bean.

CONCLUSIONS

Although each of the determinacy haplotypes probably does not represent distinct domestication events, they are consistent with the multiple (seven) domestication pattern in the genus Phaseolus. The predominance of determinacy in the Andean gene pool may reflect domestication of common bean prior to maize introduction in the Andes.

Characterization of transposable elements in the ectomycorrhizal fungus Laccaria bicolor

Background

The publicly available Laccaria bicolor genome sequence has provided a considerable genomic resource allowing systematic identification of transposable elements (TEs) in this symbiotic ectomycorrhizal fungus. Using a TE-specific annotation pipeline we have characterized and analyzed TEs in the L. bicolor S238N-H82 genome.

Methodology/Principal Findings

TEs occupy 24% of the 60 Mb L. bicolor genome and represent 25,787 full-length and partial copy elements distributed within 171 families. The most abundant elements were the Copia-like. TEs are not randomly distributed across the genome, but are tightly nested or clustered. The majority of TEs exhibits signs of ancient transposition except some intact copies of terminal inverted repeats (TIRS), long terminal repeats (LTRs) and a large retrotransposon derivative (LARD) element. There were three main periods of TE expansion in L. bicolor: the first from 57 to 10 Mya, the second from 5 to 1 Mya and the most recent from 0.5 Mya ago until now. LTR retrotransposons are closely related to retrotransposons found in another basidiomycete, Coprinopsis cinerea.

Conclusions

This analysis 1) represents an initial characterization of TEs in the L. bicolor genome, 2) contributes to improve genome annotation and a greater understanding of the role TEs played in genome organization and evolution and 3) provides a valuable resource for future research on the genome evolution within the Laccaria genus.