Aluminum tolerance in maize is associated with higher MATE1 gene copy number

The Role of Low-Molecular-Weight Organic Acids in Metal Homeostasis in Plants

Low-molecular-weight organic acids (LMWOAs) are essential O-containing metal-binding ligands involved in maintaining metal homeostasis, various metabolic processes, and plant responses to biotic and abiotic stress. Malate, citrate, and oxalate play a crucial role in metal detoxification and transport throughout the plant. This review provides a comparative analysis of the accumulation of LMWOAs in excluders, which store metals mainly in roots, and hyperaccumulators, which accumulate metals mainly in shoots. Modern concepts of the mechanisms of LMWOA secretion by the roots of excluders and hyperaccumulators are summarized, and the formation of various metal complexes with LMWOAs in the vacuole and conducting tissues, playing an important role in the mechanisms of metal detoxification and transport, is discussed. Molecular mechanisms of transport of LMWOAs and their complexes with metals across cell membranes are reviewed. It is discussed whether different endogenous levels of LMWOAs in plants determine their metal tolerance. While playing an important role in maintaining metal homeostasis, LMWOAs apparently make a minor contribution to the mechanisms of metal hyperaccumulation, which is associated mainly with root exudates increasing metal bioavailability and enhanced xylem loading of LMWOAs. The studies of metal-binding compounds may also contribute to the development of approaches used in biofortification, phytoremediation, and phytomining.

Systematic Investigation of Aluminum Stress-Related Genes and Their Critical Roles in Plants

Aluminum (Al) stress is a dominant obstacle for plant growth in acidic soil, which accounts for approximately 40-50% of the world's potential arable land. The identification and characterization of Al stress response (Al-SR) genes in Arabidopsis, rice, and other plants have deepened our understanding of Al's molecular mechanisms. However, as a crop sensitive to acidic soil, only eight Al-SR genes have been identified and functionally characterized in maize. In this review, we summarize the Al-SR genes in plants, including their classifications, subcellular localizations, expression organs, functions, and primarily molecular regulatory networks. Moreover, we predict 166 putative Al-SR genes in maize based on orthologue analyses, facilitating a comprehensive understanding of the impact of Al stress on maize growth and development. Finally, we highlight the potential applications of alleviating Al toxicity in crop production. This review deepens our understanding of the Al response in plants and provides a blueprint for alleviating Al toxicity in crop production.

Utilizing a comparative approach to assess genome evolution during diploidization in Artemisia tridentata, a keystone species of western North America

PREMISE

Polyploidization is often followed by diploidization. Diploidization is generally studied using synthetic polyploid lines and/or crop plants, but rarely using extant diploids or nonmodel plants such as Artemisia tridentata. This threatened western North American keystone species has a large genome compared to congeneric Artemisia species; dominated by diploid and tetraploid cytotypes, with multiple origins of tetraploids with genome size reduction.

METHODS

The genome of an A. tridentata sample was resequenced to study genome evolution and compared to that of A. annua, a diploid congener. Three diploid genomes of A. tridentata were compared to test for multiple diploidization events.

RESULTS

The A. tridentata genome had many chromosomal rearrangements relative to that of A. annua, while large-scale synteny of A. tridentata chromosome 3 and A. annua chromosome 4 was conserved. The three A. tridentata genomes had similar sizes (4.19-4.2 Gbp), heterozygosity (2.24-2.25%), and sequence (98.73-99.15% similarity) across scaffolds, and in k-mer analyses, similar patterns of diploid heterozygous k-mers (AB = 41%, 47%, and 47%), triploid heterozygous k-mers (AAB = 18-21%), and tetraploid k-mers (AABB = 13-17%). Biallelic SNPs were evenly distributed across scaffolds for all individuals. Comparisons of transposable element (TE) content revealed differential enrichment of TE clades.

CONCLUSIONS

Our findings suggest population-level TE differentiation after a shared polyploidization-to-diploidization event(s) and exemplify the complex processes of genome evolution. This research approached provides new resources for exploration of abiotic stress response, especially the roles of TEs in response pathways.

The prevalence of copy number increase at multiallelic copy number variants associated with cave colonization

Copy number variation is a common contributor to phenotypic diversity, yet its involvement in ecological adaptation is not easily discerned. Instances of parallelly evolving populations of the same species in a similar environment marked by strong selective pressures present opportunities to study the role of copy number variants (CNVs) in adaptation. By identifying CNVs that repeatedly occur in multiple populations of the derived ecotype and are not (or are rarely) present in the populations of the ancestral ecotype, the association of such CNVs with adaptation to the novel environment can be inferred. We used this paradigm to identify CNVs associated with recurrent adaptation of the Mexican tetra (Astyanax mexicanus) to cave environment. Using a read-depth approach, we detected CNVs from previously re-sequenced genomes of 44 individuals belonging to two ancestral surfaces and three derived cave populations. We identified 102 genes and 292 genomic regions that repeatedly diverge in copy number between the two ecotypes and occupy 0.8% of the reference genome. Functional analysis revealed their association with processes previously recognized to be relevant for adaptation, such as vision, immunity, oxygen consumption, metabolism, and neural function and we propose that these variants have been selected for in the cave or surface waters. The majority of the ecotype-divergent CNVs are multiallelic and display copy number increases in cavefish compared to surface fish. Our findings suggest that multiallelic CNVs - including gene duplications - and divergence in copy number provide a fast route to produce novel phenotypes associated with adaptation to subterranean life.

Duplicated Copy Number Variant of the Maize 9-Lipoxygenase ZmLOX5 Improves 9,10-KODA-Mediated Resistance to Fall Armyworms

Extensive genome structure variations, such as copy number variations (CNVs) and presence/absence variations, are the basis for the remarkable genetic diversity of maize; however, the effect of CNVs on maize herbivory defense remains largely underexplored. Here, we report that the naturally occurring duplication of the maize 9-lipoxygenase gene ZmLOX5 leads to increased resistance of maize to herbivory by fall armyworms (FAWs). Previously, we showed that ZmLOX5-derived oxylipins are required for defense against chewing insect herbivores and identified several inbred lines, including Yu796, that contained duplicated CNVs of ZmLOX5, referred to as Yu796-2×LOX5. To test whether introgression of the Yu796-2×LOX5 locus into a herbivore-susceptible B73 background that contains a single ZmLOX5 gene is a feasible approach to increase resistance, we generated a series of near-isogenic lines that contained either two, one, or zero copies of the Yu796-2×LOX5 locus in the B73 background via six backcrosses (BC6). Droplet digital PCR (ddPCR) confirmed the successful introgression of the Yu796-2×LOX5 locus in B73. The resulting B73-2×LOX5 inbred line displayed increased resistance against FAW, associated with increased expression of ZmLOX5, increased wound-induced production of its primary oxylipin product, the α-ketol, 9-hydroxy-10-oxo-12(Z),15(Z)-octadecadienoic acid (9,10-KODA), and the downstream defense hormones regulated by this molecule, 12-oxo-phytodienoic acid (12-OPDA) and abscisic acid (ABA). Surprisingly, wound-induced JA-Ile production was not increased in B73-2×LOX5, resulting from the increased JA catabolism. Furthermore, B73-2×LOX5 displayed reduced water loss in response to drought stress, likely due to increased ABA and 12-OPDA content. Taken together, this study revealed that the duplicated CNV of ZmLOX5 quantitively contributes to maize antiherbivore defense and presents proof-of-concept evidence that the introgression of naturally occurring duplicated CNVs of a defensive gene into productive but susceptible crop varieties is a feasible breeding approach for enhancing plant resistance to herbivory and tolerance to abiotic stress.

Climate-resilient crops: Lessons from xerophytes

Developing climate-resilient crops is critical for future food security and sustainable agriculture under current climate scenarios. Of specific importance are drought and soil salinity. Tolerance traits to these stresses are highly complex, and the progress in improving crop tolerance is too slow to cope with the growing demand in food production unless a major paradigm shift in crop breeding occurs. In this work, we combined bioinformatics and physiological approaches to compare some of the key traits that may differentiate between xerophytes (naturally drought-tolerant plants) and mesophytes (to which the majority of the crops belong). We show that both xerophytes and salt-tolerant mesophytes have a much larger number of copies in key gene families conferring some of the key traits related to plant osmotic adjustment, abscisic acid (ABA) sensing and signalling, and stomata development. We show that drought and salt-tolerant species have (i) higher reliance on Na for osmotic adjustment via more diversified and efficient operation of Na+ /H+ tonoplast exchangers (NHXs) and vacuolar H+ - pyrophosphatase (VPPases); (ii) fewer and faster stomata; (iii) intrinsically lower ABA content; (iv) altered structure of pyrabactin resistance/pyrabactin resistance-like (PYR/PYL) ABA receptors; and (v) higher number of gene copies for protein phosphatase 2C (PP2C) and sucrose non-fermenting 1 (SNF1)-related protein kinase 2/open stomata 1 (SnRK2/OST1) ABA signalling components. We also show that the past trends in crop breeding for Na+ exclusion to improve salinity stress tolerance are counterproductive and compromise their drought tolerance. Incorporating these genetic insights into breeding practices could pave the way for more drought-tolerant and salt-resistant crops, securing agricultural yields in an era of climate unpredictability.

Rapid mini-chromosome divergence among fungal isolates causing wheat blast outbreaks in Bangladesh and Zambia

The fungal pathogen, Magnaporthe oryzae Triticum pathotype, causing wheat blast disease was first identified in South America and recently spread across continents to South Asia and Africa. Here, we studied the genetic relationship among isolates found on the three continents. Magnaporthe oryzae strains closely related to a South American field isolate B71 were found to have caused the wheat blast outbreaks in South Asia and Africa. Genomic variation among isolates from the three continents was examined using an improved B71 reference genome and whole-genome sequences. We found strong evidence to support that the outbreaks in Bangladesh and Zambia were caused by the introductions of genetically separated isolates, although they were all close to B71 and, therefore, collectively referred to as the B71 branch. In addition, B71 branch strains carried at least one supernumerary mini-chromosome. Genome assembly of a Zambian strain revealed that its mini-chromosome was similar to the B71 mini-chromosome but with a high level of structural variation. Our findings show that while core genomes of the multiple introductions are highly similar, the mini-chromosomes have undergone marked diversification. The maintenance of the mini-chromosome and rapid genomic changes suggest the mini-chromosomes may serve important virulence or niche adaptation roles under diverse environmental conditions.

CcNFYB3-CcMATE35 and LncRNA CcLTCS-CcCS modules jointly regulate the efflux and synthesis of citrate to enhance aluminium tolerance in pigeon pea

Aluminium (Al) toxicity decreases crop production in acid soils in general, but many crops have evolved complex mechanisms to resist it. However, our current understanding of how plants cope with Al stress and perform Al resistance is still at the initial stage. In this study, the citrate transporter CcMATE35 was identified to be involved in Al stress response. The release of citrate was increased substantially in CcMATE35 over-expression (OE) lines under Al stress, indicating enhanced Al resistance. It was demonstrated that transcription factor CcNFYB3 regulated the expression of CcMATE35, promoting the release of citrate from roots to increase Al resistance in pigeon pea. We also found that a Long noncoding RNA Targeting Citrate Synthase (CcLTCS) is involved in Al resistance in pigeon pea. Compared with controls, overexpression of CcLTCS elevated the expression level of the Citrate Synthase gene (CcCS), leading to increases in root citrate level and citrate release, which forms another module to regulate Al resistance in pigeon pea. Simultaneous overexpression of CcNFYB3 and CcLTCS further increased Al resistance. Taken together, these findings suggest that the two modules, CcNFYB3-CcMATE35 and CcLTCS-CcCS, jointly regulate the efflux and synthesis of citrate and may play an important role in enhancing the resistance of pigeon pea under Al stress.

k-mer-based GWAS enhances the discovery of causal variants and candidate genes in soybean

Genome-wide association studies (GWAS) are powerful statistical methods that detect associations between genotype and phenotype at genome scale. Despite their power, GWAS frequently fail to pinpoint the causal variant or the gene controlling a given trait in crop species. Assessing genetic variants other than single-nucleotide polymorphisms (SNPs) could alleviate this problem. In this study, we tested the potential of structural variant (SV)- and k-mer-based GWAS in soybean by applying these methods as well as conventional SNP/indel-based GWAS to 13 traits. We assessed the performance of each GWAS approach based on loci for which the causal genes or variants were known from previous genetic studies. We found that k-mer-based GWAS was the most versatile approach and the best at pinpointing causal variants or candidate genes. Moreover, k-mer-based analyses identified promising candidate genes for loci related to pod color, pubescence form, and resistance to Phytophthora sojae. In our dataset, SV-based GWAS did not add value compared to k-mer-based GWAS and may not be worth the time and computational resources invested. Despite promising results, significant challenges remain regarding the downstream analysis of k-mer-based GWAS. Notably, better methods are needed to associate significant k-mers with sequence variation. Our results suggest that coupling k-mer- and SNP/indel-based GWAS is a powerful approach for discovering candidate genes in crop species.

The Triticeae CBF Gene Cluster-To Frost Resistance and Beyond

The pivotal role of CBF/DREB1 transcriptional factors in Triticeae crops involved in the abiotic stress response has been highlighted. The CBFs represent an important hub in the ICE-CBF-COR pathway, which is one of the most relevant mechanisms capable of activating the adaptive response to cold and drought in wheat, barley, and rye. Understanding the intricate mechanisms and regulation of the cluster of CBF genes harbored by the homoeologous chromosome group 5 entails significant potential for the genetic improvement of small grain cereals. Triticeae crops seem to share common mechanisms characterized, however, by some peculiar aspects of the response to stress, highlighting a combined landscape of single-nucleotide variants and copy number variation involving CBF members of subgroup IV. Moreover, while chromosome 5 ploidy appears to confer species-specific levels of resistance, an important involvement of the ICE factor might explain the greater tolerance of rye. By unraveling the genetic basis of abiotic stress tolerance, researchers can develop resilient varieties better equipped to withstand extreme environmental conditions. Hence, advancing our knowledge of CBFs and their interactions represents a promising avenue for improving crop resilience and food security.

Nitrate improves aluminium resistance through SLAH-mediated citrate exudation from roots

Aluminium (Al) toxicity is one of the major constraint for crop production in acidic soil, and the inappropriate utilization of nitrogen fertilizer can accelerate soil acidification. Despite previous studies investigating the regulation of nitrogen forms in Al toxicity of plants, the underlying mechanism, particularly at the molecular level, remains unclear. This study aims to uncover the potentially regulatory mechanism of nitrate (NO3 - ) in the Al resistance of maize and Arabidopsis. NO3 - conservatively improves Al resistance in maize and Arabidopsis, with nitrate-elevated citrate synthesis and exudation potentially playing critical roles in excluding Al from the root symplast. ZmSLAH2 in maize and AtSLAH1 in Arabidopsis are essential for the regulation of citrate exudation and NO3 - -promoted Al resistance, with ZmMYB81 directly targeting the ZmSLAH2 promoter to activate its activity. Additionally, NO3 - transport is necessary for NO3 - -promoted Al resistance, with ZmNRT1.1A and AtNRT1.1 potentially playing vital roles. The suppression of NO3 - transport in roots by ammonium (NH4 + ) may inhibit NO3 - -promoted Al resistance. This study provides novel insights into the understanding of the crucial role of NO3 - -mediated signalling in the Al resistance of plants and offers guidance for nitrogen fertilization on acid soils.

Accurate prediction of quantitative traits with failed SNP calls in canola and maize

In modern plant breeding, genomic selection is becoming the gold standard to select superior genotypes in large breeding populations that are only partially phenotyped. Many breeding programs commonly rely on single-nucleotide polymorphism (SNP) markers to capture genome-wide data for selection candidates. For this purpose, SNP arrays with moderate to high marker density represent a robust and cost-effective tool to generate reproducible, easy-to-handle, high-throughput genotype data from large-scale breeding populations. However, SNP arrays are prone to technical errors that lead to failed allele calls. To overcome this problem, failed calls are often imputed, based on the assumption that failed SNP calls are purely technical. However, this ignores the biological causes for failed calls-for example: deletions-and there is increasing evidence that gene presence-absence and other kinds of genome structural variants can play a role in phenotypic expression. Because deletions are frequently not in linkage disequilibrium with their flanking SNPs, permutation of missing SNP calls can potentially obscure valuable marker-trait associations. In this study, we analyze published datasets for canola and maize using four parametric and two machine learning models and demonstrate that failed allele calls in genomic prediction are highly predictive for important agronomic traits. We present two statistical pipelines, based on population structure and linkage disequilibrium, that enable the filtering of failed SNP calls that are likely caused by biological reasons. For the population and trait examined, prediction accuracy based on these filtered failed allele calls was competitive to standard SNP-based prediction, underlying the potential value of missing data in genomic prediction approaches. The combination of SNPs with all failed allele calls or the filtered allele calls did not outperform predictions with only SNP-based prediction due to redundancy in genomic relationship estimates.

Importance of genetic architecture in marker selection decisions for genomic prediction

KEY MESSAGE

We demonstrate potential for improved multi-environment genomic prediction accuracy using structural variant markers. However, the degree of observed improvement is highly dependent on the genetic architecture of the trait. Breeders commonly use genetic markers to predict the performance of untested individuals as a way to improve the efficiency of breeding programs. These genomic prediction models have almost exclusively used single nucleotide polymorphisms (SNPs) as their source of genetic information, even though other types of markers exist, such as structural variants (SVs). Given that SVs are associated with environmental adaptation and not all of them are in linkage disequilibrium to SNPs, SVs have the potential to bring additional information to multi-environment prediction models that are not captured by SNPs alone. Here, we evaluated different marker types (SNPs and/or SVs) on prediction accuracy across a range of genetic architectures for simulated traits across multiple environments. Our results show that SVs can improve prediction accuracy, but it is highly dependent on the genetic architecture of the trait and the relative gain in accuracy is minimal. When SVs are the only causative variant type, 70% of the time SV predictors outperform SNP predictors. However, the improvement in accuracy in these instances is only 1.5% on average. Further simulations with predictors in varying degrees of LD with causative variants of different types (e.g., SNPs, SVs, SNPs and SVs) showed that prediction accuracy increased as linkage disequilibrium between causative variants and predictors increased regardless of the marker type. This study demonstrates that knowing the genetic architecture of a trait in deciding what markers to use in large-scale genomic prediction modeling in a breeding program is more important than what types of markers to use.

Harnessing the power of exogenous factors to enhance plant resistance to aluminum toxicity; a critical review

Aluminum (Al) is the most prevalent element in the earth crust and is toxic to plants in acidic soils. However, plants can address Al toxicity through external exclusion (which prevents Al from entering roots) and internal detoxification (which counterbalances the toxic-Al absorbed by roots). Nowadays, certain categories of exogenously added regulatory factors (EARF), such as nutritional elements, organic acids, amino acids, phytohormones, or biochar, etc. play a critical role in reducing the bioavailability/toxicity of Al in plants. Numerous studies suggest that regulating factors against Al toxicity mediate the expression of Al-responsive genes and transcription factors, thereby regulating the secretion of organic acids, alkalizing rhizosphere pH, modulating cell wall (CW) modifications, improving antioxidant defense systems, and promoting the compartmentalization of non-toxic Al within intracellular. This review primarily discusses recent and older published papers to demonstrate the basic concepts of Al phytotoxicity. Furthermore, we provide a comprehensive explanation of the crucial roles of EARF-induced responses against Al toxicity in plants. This information may serve as a foundation for improving plant resistance to Al and enhancing the growth of susceptible species in acidic soils. And this review holds significant theoretical significance for EARF to improve the quality of acidic soils cultivated land, increase crop yield and quality, and ensure food security.

Exploring Pan-Genomes: An Overview of Resources and Tools for Unraveling Structure, Function, and Evolution of Crop Genes and Genomes

The availability of multiple sequenced genomes from a single species made it possible to explore intra- and inter-specific genomic comparisons at higher resolution and build clade-specific pan-genomes of several crops. The pan-genomes of crops constructed from various cultivars, accessions, landraces, and wild ancestral species represent a compendium of genes and structural variations and allow researchers to search for the novel genes and alleles that were inadvertently lost in domesticated crops during the historical process of crop domestication or in the process of extensive plant breeding. Fortunately, many valuable genes and alleles associated with desirable traits like disease resistance, abiotic stress tolerance, plant architecture, and nutrition qualities exist in landraces, ancestral species, and crop wild relatives. The novel genes from the wild ancestors and landraces can be introduced back to high-yielding varieties of modern crops by implementing classical plant breeding, genomic selection, and transgenic/gene editing approaches. Thus, pan-genomic represents a great leap in plant research and offers new avenues for targeted breeding to mitigate the impact of global climate change. Here, we summarize the tools used for pan-genome assembly and annotations, web-portals hosting plant pan-genomes, etc. Furthermore, we highlight a few discoveries made in crops using the pan-genomic approach and future potential of this emerging field of study.

The Role of TCA Cycle Enzymes in Plants

As one of the iconic pathways in plant metabolism, the tricarboxylic acid (TCA) cycle is commonly thought to not only be responsible for the oxidization of respiratory substrate to drive ATP synthesis but also provide carbon skeletons to anabolic processes and contribute to carbon-nitrogen interaction and biotic stress responses. The functions of the TCA cycle enzymes are characterized by a saturation transgenesis approach, whereby the constituent expression of proteins is knocked out or reduced in order to investigate their function in vivo. The alteration of TCA cycle enzyme expression results in changed plant growth and photosynthesis under controlled conditions. Moreover, improvements in plant performance and postharvest properties are reported by overexpression of either endogenous forms or heterologous genes of a number of the enzymes. Given the importance of the TCA cycle in plant metabolism regulation, here, the function of each enzyme and its roles in different tissues are discussed. This article additionally highlights the recent finding that the plant TCA cycle, like that of mammals and microbes, dynamically assembles functional substrate channels or metabolons and discusses the implications of this finding to the current understanding of the metabolic regulation of the plant TCA cycle.

k-mer-Based Genome-Wide Association Studies in Plants: Advances, Challenges, and Perspectives

Genome-wide association studies (GWAS) have allowed the discovery of marker-trait associations in crops over recent decades. However, their power is hampered by a number of limitations, with the key one among them being an overreliance on single-nucleotide polymorphisms (SNPs) as molecular markers. Indeed, SNPs represent only one type of genetic variation and are usually derived from alignment to a single genome assembly that may be poorly representative of the population under study. To overcome this, k-mer-based GWAS approaches have recently been developed. k-mer-based GWAS provide a universal way to assess variation due to SNPs, insertions/deletions, and structural variations without having to specifically detect and genotype these variants. In addition, k-mer-based analyses can be used in species that lack a reference genome. However, the use of k-mers for GWAS presents challenges such as data size and complexity, lack of standard tools, and potential detection of false associations. Nevertheless, efforts are being made to overcome these challenges and a general analysis workflow has started to emerge. We identify the priorities for k-mer-based GWAS in years to come, notably in the development of user-friendly programs for their analysis and approaches for linking significant k-mers to sequence variation.

A comprehensive analysis of copy number variations in diverse apple populations

BACKGROUND

As an important source of genetic variation, copy number variation (CNV) can alter the dosage of DNA segments, which in turn may affect gene expression level and phenotype. However, our knowledge of CNV in apple is still limited. Here, we obtained high-confidence CNVs and investigated their functional impact based on genome resequencing data of two apple populations, cultivars and wild relatives.

RESULTS

In this study, we identified 914,610 CNVs comprising 14,839 CNV regions (CNVRs) from 346 apple accessions, including 289 cultivars and 57 wild relatives. CNVRs summed to 71.19 Mb, accounting for 10.03% of the apple genome. Under the low linkage disequilibrium (LD) with nearby SNPs, they could also accurately reflect the population structure of apple independent of SNPs. Furthermore, A total of 3,621 genes were covered by CNVRs and functionally involved in biological processes such as defense response, reproduction and metabolic processes. In addition, the population differentiation index ([Formula: see text]) analysis between cultivars and wild relatives revealed 127 CN-differentiated genes, which may contribute to trait differences in these two populations.

CONCLUSIONS

This study was based on identification of CNVs from 346 diverse apple accessions, which to our knowledge was the largest dataset for CNV analysis in apple. Our work presented the first comprehensive CNV map and provided valuable resources for understanding genomic variations in apple.

Linking plant functional genes to rhizosphere microbes: a review

The importance of rhizomicrobiome in plant development, nutrition acquisition and stress tolerance is unquestionable. Relevant plant genes corresponding to the above functions also regulate rhizomicrobiome construction. Deciphering the molecular regulatory network of plant-microbe interactions could substantially contribute to improving crop yield and quality. Here, the plant gene-related nutrient uptake, biotic and abiotic stress resistance, which may influence the composition and function of microbial communities, are discussed in this review. In turn, the influence of microbes on the expression of functional plant genes, and thereby plant growth and immunity, is also reviewed. Moreover, we have specifically paid attention to techniques and methods used to link plant functional genes and rhizomicrobiome. Finally, we propose to further explore the molecular mechanisms and signalling pathways of microbe-host gene interactions, which could potentially be used for managing plant health in agricultural systems.

Structural variation of GL1 gene determines the trichome formation in Brassica juncea

A 448 kb region on chromosome B02 was delimited to be associated with trichome trait in Brassica juncea, in which the BjuVB02G54610 gene with a structural variation of 3 kb structure variation (SV) encoding a MYB transcription factor was predicted as the possible candidate gene. Mustards (Brassica juncea) are allopolyploid crops in the worldwide, and trichomes are essential quality attributes that significantly influence its taste and palpability in vegetable-use cultivars. As important accessory tissues from specialized epidermal cells, trichomes also play an important role in mitigating biotic and abiotic stresses. In this study, we constructed a F₂ segregating population using YJ27 with intensive trichome leaves and 03B0307 with glabrous leaves as parents. By bulked segregant analysis (BSA-seq), we obtained a 2.1 Mb candidate region on B02 chromosome associated with the trichome or glabrous trait formation. Then, we used 13 Kompetitive Allele Specific PCR (KASP) markers for fine mapping and finally narrowed down the candidate region to about 448 kb in length. Interestingly, among the region, there was a 3 kb sequence deletion that located on the BjuVB02G54610 gene in the F₂ individuals with trichome leaves. Genotyping results of F₂ populations confirmed this deletion (R² = 81.44%) as a major QTL. Natural population re-sequencing analysis and genotyping results further validated the key role of the 3 kb structure variation (SV) of insertion/deletion type in trichome development in B. juncea. Our findings provide important information on the formation of trichomes and potential target gene for breeding vegetable mustards.

Identification of structural variations related to drought tolerance in wheat (Triticum aestivum L.)

Structural variations are common in plant genomes, affecting meiotic recombination and distorted segregation in wheat. And presence/absence variations can significantly affect drought tolerance in wheat. Drought is a major abiotic stress limiting wheat production. Common wheat has a complex genome with three sub-genomes, which host large numbers of structural variations (SVs). SVs play critical roles in understanding the genetic contributions of plant domestication and phenotypic plasticity, but little is known about their genomic characteristics and their effects on drought tolerance. In the present study, high-resolution karyotypes of 180 doubled haploids (DHs) were developed. Signal polymorphisms between the parents involved with 8 presence-absence variations (PAVs) of tandem repeats (TR) distributed on the 7 (2A, 4A, 5A, 7A, 3B, 7B, and 2D) of 21 chromosomes. Among them, PAV on chromosome 2D showed distorted segregation, others transmit normal conforming to a 1:1 segregation ration in the population; and a PAVs recombination occurred on chromosome 2A. Association analysis of PAV and phenotypic traits under different water regimes, we found PAVs on chromosomes 4A, 5A, and 7B showed negative effect on grain length (GL) and grain width (GW); PAV.7A had opposite effect on grain thickness (GT) and spike length (SL), with the effect on traits differing under different water regimes. PAVs on linkage group 2A, 4A, 7A, 2D, and 7B associated with the drought tolerance coefficients (DTCs), and significant negative effect on drought resistance values (D values) were detected in PAV.7B. Additionally, quantitative trait loci (QTL) associated with phenotypic traits using the 90 K SNP array showed QTL for DTCs and grain-related traits in chromosomes 4A, and 5A, 3B were co-localized in differential regions of PAVs. These PAVs can cause the differentiation of the target region of SNP and could be used for genetic improvement of agronomic traits under drought stress via marker-assisted selection (MAS) breeding.

The triterpenoid saponin content difference is associated with the two type oxidosqualene cyclase gene copy numbers of Pulsatilla chinensis and Pulsatilla cernua

Pulsatilla chinensis is an important medicinal herb, its dried radix is used to treat the inflammation since ancient China. Triterpenoid saponins are proved to be the main active compounds of Pulsatilla genus. The triterpenoid saponin contents vary widely in different Pulsatilla species. But no enzyme involved in the triterpenoid saponin biosynthetic pathway was identified in Pulsitilla genus. This seriously limits the explanation of the triterpene content difference of Pulsatilla species. In this article, we obtained two oxidosqualene cyclase (OSC) genes from P. chinensis and P. cernua by touchdown PCR and anchored PCR. These two OSCs converted 2,3-oxidosqualene into different triterpenoids. The OSC from P. cernua is a monofunctional enzyme for β-amyrin synthesis, while the OSC from P. chinensis is a multifunctional enzyme for lupeol and β-amyrin synthesis, and the lupeol is the main product. Then we identified the 260th amino acid residue was the key site for the product difference by gene fusion and site-directed mutant technology. When the 260th amino acid residue was tryptophan (W260) and phenylalanine (F260), the main catalysate was β-amyrin and lupeol, respectively. Then we found that the expression of these two genes was strongly correlated with the lupeol-type and β-amyrin-type triterpenoid contents in P. cernua and P. chinensis. Finally, we found the gene copy number difference of these two genotypes leaded to the triterpenoid diversity in P. cernua and P. chinensis. This study provides useful information for the molecular breeding and quality improvement of P. chinensis and a molecular marker to identify the P. chinensis decoction pieces.

Era of gapless plant genomes: innovations in sequencing and mapping technologies revolutionize genomics and breeding

Whole-genome sequencing and assembly have revolutionized plant genetics and molecular biology over the last two decades. However, significant shortcomings in first- and second-generation technology resulted in imperfect reference genomes: numerous and large gaps of low quality or undeterminable sequence in areas of highly repetitive DNA along with limited chromosomal phasing restricted the ability of researchers to characterize regulatory noncoding elements and genic regions that underwent recent duplication events. Recently, advances in long-read sequencing have resulted in the first gapless, telomere-to-telomere (T2T) assemblies of plant genomes. This leap forward has the potential to increase the speed and confidence of genomics and molecular experimentation while reducing costs for the research community.

High-quality Fagopyrum esculentum genome provides insights into the flavonoid accumulation among different tissues and self-incompatibility

Common buckwheat (Fagopyrum esculentum) and Tartary buckwheat (Fagopyrum tataricum), the two most widely cultivated buckwheat species, differ greatly in flavonoid content and reproductive mode. Here, we report the first high-quality and chromosome-level genome assembly of common buckwheat with 1.2 Gb. Comparative genomic analysis revealed that common buckwheat underwent a burst of long terminal repeat retrotransposons insertion accompanied by numerous large chromosome rearrangements after divergence from Tartary buckwheat. Moreover, multiple gene families involved in stress tolerance and flavonoid biosynthesis such as multidrug and toxic compound extrusion (MATE) and chalcone synthase (CHS) underwent significant expansion in buckwheat, especially in common buckwheat. Integrated multi-omics analysis identified high expression of catechin biosynthesis-related genes in flower and seed in common buckwheat and high expression of rutin biosynthesis-related genes in seed in Tartary buckwheat as being important for the differences in flavonoid type and content between these buckwheat species. We also identified a candidate key rutin-degrading enzyme gene (Ft8.2377) that was highly expressed in Tartary buckwheat seed. In addition, we identified a haplotype-resolved candidate locus containing many genes reportedly associated with the development of flower and pollen, which was potentially related to self-incompatibility in common buckwheat. Our study provides important resources facilitating future functional genomics-related research of flavonoid biosynthesis and self-incompatibility in buckwheat.

Aluminum in plant: Benefits, toxicity and tolerance mechanisms

Aluminum (Al) is the third most ubiquitous metal in the earth's crust. A decrease in soil pH below 5 increases its solubility and availability. However, its impact on plants depends largely on concentration, exposure time, plant species, developmental age, and growing conditions. Although Al can be beneficial to plants by stimulating growth and mitigating biotic and abiotic stresses, it remains unknown how Al mediates these effects since its biological significance in cellular systems is still unidentified. Al is considered a major limiting factor restricting plant growth and productivity in acidic soils. It instigates a series of phytotoxic symptoms in several Al-sensitive crops with inhibition of root growth and restriction of water and nutrient uptake as the obvious symptoms. This review explores advances in Al benefits, toxicity and tolerance mechanisms employed by plants on acidic soils. These insights will provide directions and future prospects for potential crop improvement.

Rice Pangenome Genotyping Array: an efficient genotyping solution for pangenome-based accelerated genetic improvement in rice

The advent of the pangenome era has unraveled previously unknown genetic variation existing within diverse crop plants, including rice. This untapped genetic variation is believed to account for a major portion of phenotypic variation existing in crop plants. However, the use of conventional single reference-guided genotyping often fails to capture a large portion of this genetic variation leading to a reference bias. This makes it difficult to identify and utilize novel population/cultivar-specific genes for crop improvement. Thus, we developed a Rice Pangenome Genotyping Array (RPGA) harboring probes assaying 80K single-nucleotide polymorphisms (SNPs) and presence-absence variants spanning the entire 3K rice pangenome. This array provides a simple, user-friendly and cost-effective (60-80 USD per sample) solution for rapid pangenome-based genotyping in rice. The genome-wide association study (GWAS) conducted using RPGA-SNP genotyping data of a rice diversity panel detected a total of 42 loci, including previously known as well as novel genomic loci regulating grain size/weight traits in rice. Eight of these identified trait-associated loci (dispensable loci) could not be detected with conventional single reference genome-based GWAS. A WD repeat-containing PROTEIN 12 gene underlying one of such dispensable locus on chromosome 7 (qLWR7) along with other non-dispensable loci were subsequently detected using high-resolution quantitative trait loci mapping confirming authenticity of RPGA-led GWAS. This demonstrates the potential of RPGA-based genotyping to overcome reference bias. The application of RPGA-based genotyping for population structure analysis, hybridity testing, ultra-high-density genetic map construction and chromosome-level genome assembly, and marker-assisted selection was also demonstrated. A web application (http://www.rpgaweb.com) was further developed to provide an easy to use platform for the imputation of RPGA-based genotyping data using 3K rice reference panel and subsequent GWAS.

Genome-wide association study reveals structural chromosome variations with phenotypic effects in wheat (Triticum aestivum L.)

Structural chromosome variations (SCVs) are large-scale genomic variations that can be detected by fluorescence in situ hybridization (FISH). SCVs have played important roles in the genome evolution of wheat (Triticum aestivum L.), but little is known about their genetic effects. In this study, a total of 543 wheat accessions from the Chinese wheat mini-core collection and the Shanxi Province wheat collection were used for chromosome analysis using oligonucleotide probe multiplex FISH. A total of 139 SCVs including translocations, pericentric inversions, presence/absence variations (PAVs), and copy number variations (CNVs) in heterochromatin were identified at 230 loci. The distribution frequency of SCVs varied between ecological regions and between landraces and modern cultivars. Structural analysis using SCVs as markers clearly divided the landraces and modern cultivars into different groups. There are very clear instances illustrating alien introgression and wide application of foreign germplasms improved the chromosome diversity of Chinese modern wheat cultivars. A genome-wide association study (GWAS) identified 29 SCVs associated with 12 phenotypic traits, and five (RT4AS•4AL-1DS/1DL•1DS-4AL, Mg2A-3, Mr3B-10, Mr7B-13, and Mr4A-7) of them were further validated using a doubled haploid population and advanced sib-lines, implying the potential value of these SCVs. Importantly, the number of favored SCVs that were associated with agronomic trait improvement was significantly higher in modern cultivars compared to landraces, indicating positive selection in wheat breeding. This study demonstrates the significant effects of SCVs during wheat breeding and provides an efficient method of mining favored SCVs in wheat and other crops.

Recent advances and current challenges in population genomics of structural variation in animals and plants

The field of population genomics has seen a surge of studies on genomic structural variation over the past two decades. These studies witnessed that structural variation is taxonomically ubiquitous and represent a dominant form of genetic variation within species. Recent advances in technology, especially the development of long-read sequencing platforms, have enabled the discovery of structural variants (SVs) in previously inaccessible genomic regions which unlocked additional structural variation for population studies and revealed that more SVs contribute to evolution than previously perceived. An increasing number of studies suggest that SVs of all types and sizes may have a large effect on phenotype and consequently major impact on rapid adaptation, population divergence, and speciation. However, the functional effect of the vast majority of SVs is unknown and the field generally lacks evidence on the phenotypic consequences of most SVs that are suggested to have adaptive potential. Non-human genomes are heavily under-represented in population-scale studies of SVs. We argue that more research on other species is needed to objectively estimate the contribution of SVs to evolution. We discuss technical challenges associated with SV detection and outline the most recent advances towards more representative reference genomes, which opens a new era in population-scale studies of structural variation.

Global analysis of common bean multidrug and toxic compound extrusion transporters (PvMATEs): PvMATE8 and pinto bean seed coat darkening

In common bean (Phaseolus vulgaris L.), postharvest seed coat darkening is an undesirable trait that affects crop value. The increased accumulation of proanthocyanidins (PAs) in the seed coat results in darker seeds in many market classes of colored beans after harvest. The precursors of PAs are synthesized in the cytoplasm, and subsequently get glycosylated and then transported to the vacuoles where polymerization occurs. Thus, vacuolar transporters play an important role in the accumulation of PAs. Here, we report that common bean genome contains 59 multidrug and toxic compound extrusion genes (PvMATEs). Phylogenetic analysis of putative PvMATEs with functionally characterized MATEs from other plant species categorized them into substrate-specific clades. Our data demonstrate that a vacuolar transporter PvMATE8 is expressed at a higher level in the pinto bean cultivar CDC Pintium (regular darkening) compared to 1533-15 (slow darkening). PvMATE8 localizes in the vacuolar membrane and rescues the PA deficient (tt12) mutant phenotype in Arabidopsis thaliana. Analysis of PA monomers in transgenic seeds together with wild-type and mutants suggests a possible feedback regulation of PA biosynthesis and accumulation. Identification of PvMATE8 will help better understand the mechanism of PA accumulation in common bean.

Characterization of GmMATE13 in its contribution of citrate efflux and aluminum resistance in soybeans

Citrate exudation mediated by a citrate transporter of the MATE protein family is critical for resisting aluminum (Al) toxicity in soybeans. However, the expression patterns of citrate transporter genes differ under Al stress. Thus, exploring the responsive pattern of GmMATEs in response to Al stress is of great importance to understand the Al resistance mechanism in soybeans. In the present study, the phylogenetic analysis, transcriptionally expressed pattern, and function of GmMATE13 were investigated. The results show that soybean GmMATE13 is highly homologous to known citrate transporter proteins from other plants. Under Al exposure, the transcript abundance of GmMATE13 was increased during a 24 h Al treatment period. The expression of GmMATE13 is specifically induced by Al exposure, but not by the status of Fe, Cu, Cd, or La. Moreover, it was also highly increased when soybean seedlings were grown on acidic soil with a high Al content. Subcellular localization showed that GmMATE13 was localized on the plasma membrane when it was transiently expressed in Arabidopsis protoplasts. Investigation of tissue localization of GmMATE13 expression by investigating GUS activity staining under control of the GmMATE13 promoter showed that it was mainly expressed in the central cylinder in the root tips of the soybean under Al-free conditions, yet extended to cortical and epidermis cells under Al stress. Finally, overexpressing GmMATE13 in soybean hairy roots enhanced Al resistance by increasing citrate efflux. Collectively, we conclude that GmMATE13 is a promising candidate to improve the resistance of soybean to Al toxicity in acidic soil.

Structural variants in the barley gene pool: precision and sensitivity to detect them using short-read sequencing and their association with gene expression and phenotypic variation

Structural variants (SV) of 23 barley inbreds, detected by the best combination of SV callers based on short-read sequencing, were associated with genome-wide and gene-specific gene expression and, thus, were evaluated to predict agronomic traits. In human genetics, several studies have shown that phenotypic variation is more likely to be caused by structural variants (SV) than by single nucleotide variants. However, accurate while cost-efficient discovery of SV in complex genomes remains challenging. The objectives of our study were to (i) facilitate SV discovery studies by benchmarking SV callers and their combinations with respect to their sensitivity and precision to detect SV in the barley genome, (ii) characterize the occurrence and distribution of SV clusters in the genomes of 23 barley inbreds that are the parents of a unique resource for mapping quantitative traits, the double round robin population, (iii) quantify the association of SV clusters with transcript abundance, and (iv) evaluate the use of SV clusters for the prediction of phenotypic traits. In our computer simulations based on a sequencing coverage of 25x, a sensitivity > 70% and precision > 95% was observed for all combinations of SV types and SV length categories if the best combination of SV callers was used. We observed a significant (P < 0.05) association of gene-associated SV clusters with global gene-specific gene expression. Furthermore, about 9% of all SV clusters that were within 5 kb of a gene were significantly (P < 0.05) associated with the gene expression of the corresponding gene. The prediction ability of SV clusters was higher compared to that of single-nucleotide polymorphisms from an array across the seven studied phenotypic traits. These findings suggest the usefulness of exploiting SV information when fine mapping and cloning the causal genes underlying quantitative traits as well as the high potential of using SV clusters for the prediction of phenotypes in diverse germplasm sets.

Identification of two chickpea multidrug and toxic compound extrusion transporter genes transcriptionally upregulated upon aluminum treatment in root tips

Aluminum (Al) toxicity poses a significant challenge for the yield improvement of chickpea, which is an economically important legume crop with high nutritional value in human diets. The genetic basis of Al-tolerance in chickpea remains unclear. Here, we assessed the Al-tolerance of 8 wild Cicer and one cultivated chickpea (PBA Pistol) accessions by measuring the root elongation in solution culture under control (0 μM Al³⁺) and Al treatments (15, 30 μM Al³⁺). Compared to PBA Pistol, the wild Cicer accessions displayed both tolerant and sensitive phenotypes, supporting wild Cicer as a potential genetic pool for Al-tolerance improvement. To identify potential genes related to Al-tolerance in chickpea, genome-wide screening of multidrug and toxic compound extrusion (MATE) encoding genes was performed. Fifty-six MATE genes were identified in total, which can be divided into 4 major phylogenetic groups. Four chickpea MATE genes (CaMATE1-4) were clustered with the previously characterized citrate transporters MtMATE66 and MtMATE69 in Medicago truncatula. Transcriptome data showed that CaMATE1-4 have diverse expression profiles, with CaMATE2 being root-specific. qRT-PCR analyses confirmed that CaMATE2 and CaMATE4 were highly expressed in root tips and were up-regulated upon Al treatment in all chickpea lines. Further measurement of carboxylic acids showed that malonic acid, instead of malate or citrate, is the major extruded acid by Cicer spp. root. Protein structural modeling analyses revealed that CaMATE2 has a divergent substrate-binding cavity from Arabidopsis AtFRD3, which may explain the different acid-secretion profile for chickpea. Pangenome survey showed that CaMATE1-4 have much higher genetic diversity in wild Cicer than that in cultivated chickpea. This first identification of CaMATE2 and CaMATE4 responsive to Al³⁺ treatment in Cicer paves the way for future functional characterization of MATE genes in Cicer spp., and to facilitate future design of gene-specific markers for Al-tolerant line selection in chickpea breeding programs.

microRNAs: Key Players in Plant Response to Metal Toxicity

Environmental metal pollution is a common problem threatening sustainable and safe crop production. Heavy metals (HMs) cause toxicity by targeting key molecules and life processes in plant cells. Plants counteract excess metals in the environment by enhancing defense responses, such as metal chelation, isolation to vacuoles, regulating metal intake through transporters, and strengthening antioxidant mechanisms. In recent years, microRNAs (miRNAs), as a small non-coding RNA, have become the central regulator of a variety of abiotic stresses, including HMs. With the introduction of the latest technologies such as next-generation sequencing (NGS), more and more miRNAs have been widely recognized in several plants due to their diverse roles. Metal-regulated miRNAs and their target genes are part of a complex regulatory network. Known miRNAs coordinate plant responses to metal stress through antioxidant functions, root growth, hormone signals, transcription factors (TF), and metal transporters. This article reviews the research progress of miRNAs in the stress response of plants to the accumulation of HMs, such as Cu, Cd, Hg, Cr, and Al, and the toxicity of heavy metal ions.

The NAC transcription factor ANAC017 regulates aluminum tolerance by regulating the cell wall-modifying genes

Aluminum (Al) toxicity is one of the key factors limiting crop production in acid soils; however, little is known about its transcriptional regulation in plants. In this study, we characterized the role of a NAM, ATAF1/2, and cup-shaped cotyledon 2 (NAC) transcription factors (TFs), ANAC017, in the regulation of Al tolerance in Arabidopsis (Arabidopsis thaliana). ANAC017 was localized in the nucleus and exhibited constitutive expression in the root, stem, leaf, flower, and silique, although its expression and protein accumulation were repressed by Al stress. Loss of function of ANAC017 enhanced Al tolerance when compared with wild-type Col-0 and was accompanied by lower root and root cell wall Al content. Furthermore, both hemicellulose and xyloglucan content decreased in the anac017 mutants, indicating the possible interaction between ANAC017 and xyloglucan endotransglucosylase/hydrolase (XTH). Interestingly, the expression of XTH31, which is responsible for xyloglucan modification, was downregulated in the anac017 mutants regardless of Al supply, supporting the possible interaction between ANAC017 and XTH31. Yeast one-hybrid, dual-luciferase reporter assay, and chromatin immunoprecipitation-quantitative PCR analysis revealed that ANAC017 positively regulated the expression of XTH31 through directly binding to the XTH31 promoter region, and overexpression of XTH31 in the anac017 mutant background rescued its Al-tolerance phenotype. In conclusion, we identified that the tTF ANAC017 acts upstream of XTH31 to regulate Al tolerance in Arabidopsis.

Identification of Heat-Tolerant Genes in Non-Reference Sequences in Rice by Integrating Pan-Genome, Transcriptomics, and QTLs

The availability of large-scale genomic data resources makes it very convenient to mine and analyze genes that are related to important agricultural traits in rice. Pan-genomes have been constructed to provide insight into the genome diversity and functionality of different plants, which can be used in genome-assisted crop improvement. Thus, a pan-genome comprising all genetic elements is crucial for comprehensive variation study among the heat-resistant and -susceptible rice varieties. In this study, a rice pan-genome was firstly constructed by using 45 heat-tolerant and 15 heat-sensitive rice varieties. A total of 38,998 pan-genome genes were identified, including 37,859 genes in the reference and 1141 in the non-reference contigs. Genomic variation analysis demonstrated that a total of 76,435 SNPs were detected and identified as the heat-tolerance-related SNPs, which were specifically present in the highly heat-resistant rice cultivars and located in the genic regions or within 2 kbp upstream and downstream of the genes. Meanwhile, 3214 upregulated and 2212 downregulated genes with heat stress tolerance-related SNPs were detected in one or multiple RNA-seq datasets of rice under heat stress, among which 24 were located in the non-reference contigs of the rice pan-genome. We then mapped the DEGs with heat stress tolerance-related SNPs to the heat stress-resistant QTL regions. A total of 1677 DEGs, including 990 upregulated and 687 downregulated genes, were mapped to the 46 heat stress-resistant QTL regions, in which 2 upregulated genes with heat stress tolerance-related SNPs were identified in the non-reference sequences. This pan-genome resource is an important step towards the effective and efficient genetic improvement of heat stress resistance in rice to help meet the rapidly growing needs for improved rice productivity under different environmental stresses. These findings provide further insight into the functional validation of a number of non-reference genes and, especially, the two genes identified in the heat stress-resistant QTLs in rice.

Metabolic and DNA checkpoints for the enhancement of Al tolerance

Acidic soils are a major limiting factor for food production in many developing countries. High concentrations of soluble Al cations, particularly Al³⁺, inhibit cell division and root elongation in plants. Al³⁺ damages several biomolecules, including DNA, impairing gene expression and cell cycle progression. Notably, the loss-of-function mutants of DNA checkpoints may mediate Al tolerance. Furthermore, mitochondrial organic acids play key roles in neutralizing Al³⁺ within the cell and around the rhizosphere. Here, we provide knowledge synthesis on interactions between checkpoints related to mitochondrial organic acid homeostasis and DNA integrity mediating Al tolerance in land plants. These interactions, coupled with remarkable advances in tools related to metabolism and cell cycle, may facilitate the development of next-generation productive crops under Al toxicity.

Genome-Wide Identification and Characterisation of Wheat MATE Genes Reveals Their Roles in Aluminium Tolerance

The Multidrug and toxin efflux (MATE) gene family plays crucial roles in plant growth and development and response to adverse stresses. This work investigated the structural and evolutionary characteristics, expression profiling and potential functions involved in aluminium (Al) tolerance from a genome-wide level. In total, 211 wheat MATE genes were identified, which were classified into four subfamilies and unevenly distributed on chromosomes. Duplication analysis showed that fragments and tandem repeats played the main roles in the amplification of TaMATEs, and Type II functional disproportionation had a leading role in the differentiation of TaMATEs. TaMATEs had abundant Al resistance and environmental stress-related elements, and generally had a high expression level in roots and leaves and in response to Al stress. The 3D structure prediction by AlphaFold and molecular docking showed that six TaMATE proteins localised in the plasmalemma could combine with citrate via amino acids in the citrate exuding motif and other sites, and then transport citrate to soil to form citrate aluminium. Meanwhile, citrate aluminium formed in root cells might be transported to leaves by TaMATEs to deposit in vacuoles, thereby alleviating Al toxicity.

Genomic structural variation in tomato and its role in plant immunity

It is well known that large genomic variations can greatly impact the phenotype of an organism. Structural Variants (SVs) encompass any genomic variation larger than 30 base pairs, and include changes caused by deletions, inversions, duplications, transversions, and other genome modifications. Due to their size and complex nature, until recently, it has been difficult to truly capture these variations. Recent advances in sequencing technology and computational analyses now permit more extensive studies of SVs in plant genomes. In tomato, advances in sequencing technology have allowed researchers to sequence hundreds of genomes from tomatoes, and tomato relatives. These studies have identified SVs related to fruit size and flavor, as well as plant disease response, resistance/susceptibility, and the ability of plants to detect pathogens (immunity). In this review, we discuss the implications for genomic structural variation in plants with a focus on its role in tomato immunity. We also discuss how advances in sequencing technology have led to new discoveries of SVs in more complex genomes, the current evidence for the role of SVs in biotic and abiotic stress responses, and the outlook for genetic modification of SVs to advance plant breeding objectives.

Combined use of Oxford Nanopore and Illumina sequencing yields insights into soybean structural variation biology

BACKGROUND

Structural variants (SVs), including deletions, insertions, duplications, and inversions, are relatively long genomic variations implicated in a diverse range of processes from human disease to ecology and evolution. Given their complex signatures, tendency to occur in repeated regions, and large size, discovering SVs based on short reads is challenging compared to single-nucleotide variants. The increasing availability of long-read technologies has greatly facilitated SV discovery; however, these technologies remain too costly to apply routinely to population-level studies. Here, we combined short-read and long-read sequencing technologies to provide a comprehensive population-scale assessment of structural variation in a panel of Canadian soybean cultivars.

RESULTS

We used Oxford Nanopore long-read sequencing data (~12× mean coverage) for 17 samples to both benchmark SV calls made from Illumina short-read data and predict SVs that were subsequently genotyped in a population of 102 samples using Illumina data. Benchmarking results show that variants discovered using Oxford Nanopore can be accurately genotyped from the Illumina data. We first use the genotyped deletions and insertions for population genetics analyses and show that results are comparable to those based on single-nucleotide variants. We observe that the population frequency and distribution within the genome of deletions and insertions are constrained by the location of genes. Gene Ontology and PFAM domain enrichment analyses also confirm previous reports that genes harboring high-frequency deletions and insertions are enriched for functions in defense response. Finally, we discover polymorphic transposable elements from the deletions and insertions and report evidence of the recent activity of a Stowaway MITE.

CONCLUSIONS

We show that structural variants discovered using Oxford Nanopore data can be genotyped with high accuracy from Illumina data. Our results demonstrate that long-read and short-read sequencing technologies can be efficiently combined to enhance SV analysis in large populations, providing a reusable framework for their study in a wider range of samples and non-model species.

Melatonin alleviates aluminum-induced growth inhibition by modulating carbon and nitrogen metabolism, and reestablishing redox homeostasis in Zea mays L

Melatonin, a regulatory molecule, performs pleiotropic functions in plants, including aluminum (Al) stress mitigation. Here, we conducted transcriptomic and physiological analyses to identify metabolic processes associated with the alleviated Al-induced growth inhibition of the melatonin-treated (MT) maize (Zea mays L.) seedlings. Melatonin decreased Al concentration in maize roots and leaves under Al stress. Al stress reduced the total dry weight (DW) by 41.2% after 7 days of treatment. By contrast, the total DW was decreased by only 19.4% in MT plants. According to RNA-Seq, enzyme activity, and metabolite content data, MT plants exhibited a higher level of relatively stable carbon and nitrogen metabolism than non-treated (NT) plants. Under Al stress, MT plants showed higher photosynthetic rate and sucrose content by 29.9% and 20.5% than NT plants, respectively. Similarly, the nitrate reductase activity and protein content of MT plants were 34.0% and 15.0% higher than those of NT plants, respectively. Furthermore, exogenous supply of melatonin mitigated Al-induced oxidative stress. Overall, our results suggest that melatonin alleviates aluminum-induced growth inhibition through modulating carbon and nitrogen metabolism, and reestablishing redox homeostasis in maize. Graphical Abstarct.

Duplication and specialization of NUDX1 in Rosaceae led to geraniol production in rose petals

Nudix hydrolases are conserved enzymes ubiquitously present in all kingdoms of life. Recent research revealed that several Nudix hydrolases are involved in terpenoid metabolism in plants. In modern roses, RhNUDX1 is responsible for formation of geraniol, a major compound of rose scent. Nevertheless, this compound is produced by monoterpene synthases in many geraniol-producing plants. As a consequence, this raised the question about the origin of RhNUDX1 function and the NUDX1 gene evolution in Rosaceae, in wild roses or/and during the domestication process. Here, we showed that three distinct clades of NUDX1 emerged in the Rosoidae subfamily (Nudx1-1 to Nudx1-3 clades), and two subclades evolved in the Rosa genus (Nudx1-1a and Nudx1-1b subclades). We also showed that the Nudx1-1b subclade was more ancient than the Nudx1-1a subclade, and that the NUDX1-1a gene emerged by a trans-duplication of the more ancient NUDX1-1b gene. After the transposition, NUDX1-1a was cis-duplicated, leading to a gene dosage effect on the production of geraniol in different species. Furthermore, the NUDX1-1a appearance was accompanied by the evolution of its promoter, most likely from a Copia retrotransposon origin, leading to its petal-specific expression. Thus, our data strongly suggest that the unique function of NUDX1-1a in geraniol formation was evolved naturally in the genus Rosa before domestication.

Genome-wide analysis of deletions in maize population reveals abundant genetic diversity and functional impact

Two read depth methods were jointly used in next-generation sequencing data to identify deletions in maize population. GWAS by deletions were analyzed for gene expression pattern and classical traits, respectively. Many studies have confirmed that structural variation (SV) is pervasive throughout the maize genome. Deletion is one type of SV that may impact gene expression and cause phenotypic changes in quantitative traits. In this study, two read count approaches were used to analyze the deletions in the whole-genome sequencing data of 270 maize inbred lines. A total of 19,754 deletion windows overlapped 12,751 genes, which were unevenly distributed across the genome. The deletions explained population structure well and correlated with genomic features. The deletion proportion of genes was determined to be negatively correlated with its expression. The detection of gene expression quantitative trait loci (eQTL) indicated that local eQTL were fewer but had larger effects than distant ones. The common associated genes were related to basic metabolic processes, whereas unique associated genes with eQTL played a role in the stress or stimulus responses in multiple tissues. Compared with the eQTL detected by SNPs derived from the same sequencing data, 89.4% of the associated genes could be detected by both markers. The effect of top eQTL detected by SNPs was usually larger than that detected by deletions for the same gene. A genome-wide association study (GWAS) on flowering time and plant height illustrated that only a few loci could be consistently captured by SNPs, suggesting that combining deletion and SNP for GWAS was an excellent strategy to dissect trait architecture. Our findings will provide insights into characteristic and biological function of genome-wide deletions in maize.

Bacteriophage self-counting in the presence of viral replication

When host cells are in low abundance, temperate bacteriophages opt for dormant (lysogenic) infection. Phage lambda implements this strategy by increasing the frequency of lysogeny at higher multiplicity of infection (MOI). However, it remains unclear how the phage reliably counts infecting viral genomes even as their intracellular number increases because of replication. By combining theoretical modeling with single-cell measurements of viral copy number and gene expression, we find that instead of hindering lambda's decision, replication facilitates it. In a nonreplicating mutant, viral gene expression simply scales with MOI rather than diverging into lytic (virulent) and lysogenic trajectories. A similar pattern is followed during early infection by wild-type phage. However, later in the infection, the modulation of viral replication by the decision genes amplifies the initially modest gene expression differences into divergent trajectories. Replication thus ensures the optimal decision-lysis upon single-phage infection and lysogeny at higher MOI.

Aluminium occurrence in plant feed from Northwestern Italy

BACKGROUND

Aluminium (Al) is widespread in nature, but the main route of exposure for humans and animals is through food. Cereals and vegetables give and important contribute to dietary intake. This investigation focused on the occurrence aluminium in raw materials intended for animal consumption. In fact, Al is not included in the official monitoring programs, and little or no data are still available about its presence in feed or plant material.

METHODS

Rice, wheat, corn, barley, rye, triticale, oats, and soy samples were collected in Northwestern Italy in an area characterized by clayey soil rich in aluminium silicates. Samples were subjected to homogenization and microwave acid digestion and Al was quantified by Inductively Coupled Plasma-Mass Spectrometer (ICP-MS).

RESULTS

Different aluminium concentrations were found in the analyzed raw materials: soy (89 mg kg^-1) > rye > rice > oat > triticale > barley > wheat > corn (2.9 mg kg^-1). In this study, Al concentrations were found higher than those reported for same matrices in other UE and non-EU countries, especially in rice and soy.

CONCLUSIONS

Al concentrations were related to the soil geological feature and to the physiological characteristics of the plant species. These two factors must be carefully considered for a reliable evaluation of aluminium exposure through feed and food.

Deciphering the major metabolic pathways associated with aluminum tolerance in popcorn roots using label-free quantitative proteomics

Al responsive proteins are associated with starch, sucrose, and other carbohydrate metabolic pathways. Sucrose synthase is a candidate to Al tolerance. Al responses are regulated at transcriptional and post-transcriptional levels. Aluminum toxicity is one of the important abiotic stresses that affects worldwide crop production. The soluble form of aluminum (Al³⁺) inhibits root growth by altering water and nutrient uptake, a process that also reduces plant growth and development. Under long-term Al³⁺ exposure, plants can activate several tolerance mechanisms. To date, no reports of large-scale proteomic data concerning maize responses to this ion have been published. To investigate the post-transcriptional regulation in response to Al toxicity, we performed label-free quantitative proteomics for comparative analysis of two Al-contrasting popcorn inbred lines and an Al-tolerant commercial hybrid during 72 h under Al-stress conditions. A total of 489 differentially accumulated proteins (DAPs) were identified in the Al-sensitive inbred line, 491 in the Al-tolerant inbred line, and 277 in the commercial hybrid. Among them, 120 DAPs were co-expressed in both Al tolerant genotypes. Bioinformatics analysis indicated that starch, sucrose, and other components of carbohydrate metabolism and glycolysis/gluconeogenesis are the biochemical processes regulated in response to Al toxicity. Sucrose synthase accumulation and an increase in sucrose content and starch degradation suggest that these components may enhance popcorn tolerance to Al stress. The accumulation of citrate synthase suggests a key role for this enzyme in the detoxification process in the Al-tolerant inbred line. The integration of transcriptomic and proteomic data indicates that the Al tolerance response presents a complex regulatory network into the transcription and translation dynamics of popcorn root development.

Calling large indels in 1047 Arabidopsis with IndelEnsembler

Large indels greatly impact the observable phenotypes in different organisms including plants and human. Hence, extracting large indels with high precision and sensitivity is important. Here, we developed IndelEnsembler to detect large indels in 1047 Arabidopsis whole-genome sequencing data. IndelEnsembler identified 34 093 deletions, 12 913 tandem duplications and 9773 insertions. Our large indel dataset was more comprehensive and accurate compared with the previous dataset of AthCNV (1). We captured nearly twice of the ground truth deletions and on average 27% more ground truth duplications compared with AthCNV, though our dataset has less number of large indels compared with AthCNV. Our large indels were positively correlated with transposon elements across the Arabidopsis genome. The non-homologous recombination events were the major formation mechanism of deletions in Arabidopsis genome. The Neighbor joining (NJ) tree constructed based on IndelEnsembler's deletions clearly divided the geographic subgroups of 1047 Arabidopsis. More importantly, our large indels represent a previously unassessed source of genetic variation. Approximately 49% of the deletions have low linkage disequilibrium (LD) with surrounding single nucleotide polymorphisms. Some of them could affect trait performance. For instance, using deletion-based genome-wide association study (DEL-GWAS), the accessions containing a 182-bp deletion in AT1G11520 had delayed flowering time and all accessions in north Sweden had the 182-bp deletion. We also found the accessions with 65-bp deletion in the first exon of AT4G00650 (FRI) flowered earlier than those without it. These two deletions cannot be detected in AthCNV and, interestingly, they do not co-occur in any Arabidopsis thaliana accession. By SNP-GWAS, surrounding SNPs of these two deletions do not correlate with flowering time. This example demonstrated that existing large indel datasets miss phenotypic variations and our large indel dataset filled in the gap.

LD-CNV: rapid and simple discovery of chromosomal translocations using linkage disequilibrium between copy number variable loci

Large-scale structural variations, such as chromosomal translocations, can have profound effects on fitness and phenotype, but are difficult to identify and characterize. Here, we describe a simple and effective method aimed at identifying translocations using only the dosage of sequence reads mapped on the reference genome. We binned reads on genomic segments sized according to sequencing coverage and identified instances when copy number segregated in populations. For each dosage-polymorphic 1 Mb bin, we tested independence, effectively an apparent linkage disequilibrium (LD), with other variable bins. In nine potato (Solanum tuberosum) dihaploid families translocations affecting pericentromeric regions were common and in two cases were due to genomic misassembly. In two populations, we found evidence for translocation affecting euchromatic arms. In cv. PI 310467, a nonreciprocal translocation between chromosomes (chr.) 7 and 8 resulted in a 5-3 copy number change affecting several Mb at the respective chromosome tips. In cv. "Alca Tarma," the terminal arm of chr. 4 translocated to the tip of chr. 1. Using oligonucleotide-based fluorescent in situ hybridization painting probes (oligo-FISH), we tested and confirmed the predicted arrangement in PI 310467. In 192 natural accessions of Arabidopsis thaliana, dosage haplotypes tended to vary continuously and resulted in higher noise, while apparent LD between pericentromeric regions suggested the effect of repeats. This method, LD-CNV, should be useful in species where translocations are suspected because it tests linkage without the need for genotyping.

Uncovering the transcriptional response of popcorn (Zea mays L. var. everta) under long-term aluminum toxicity

To date, the investigation of genes involved in Al resistance has focused mainly on microarrays and short periods of Al exposure. We investigated genes involved in the global response under Al stress by tracking the expression profile of two inbred popcorn lines with different Al sensitivity during 72 h of Al stress. A total of 1003 differentially expressed genes were identified in the Al-sensitive line, and 1751 were identified in the Al-resistant line, of which 273 were shared in both lines. Genes in the category of “response to abiotic stress” were present in both lines, but there was a higher number in the Al-resistant line. Transcription factors, genes involved in fatty acid biosynthesis, and genes involved in cell wall modifications were also detected. In the Al-resistant line, GST6 was identified as one of the key hub genes by co-expression network analysis, and ABC6 may play a role in the downstream regulation of CASP-like 5. In addition, we suggest a class of SWEET transporters that might be involved in the regulation of vacuolar sugar storage and may serve as mechanisms for Al resistance. The results and conclusions expand our understanding of the complex mechanisms involved in Al toxicity and provide a platform for future functional analyses and genomic studies of Al stress in popcorn.

ZmMATE6 from maize encodes a citrate transporter that enhances aluminum tolerance in transgenic Arabidopsis thaliana

The yields of cereal crops grown on acidic soils are often reduced by aluminum (Al) toxicity because the prevalence of toxic Al³⁺ cations increases as pH falls below 5.0. The Al-dependent release of citrate from resistant lines of maize is controlled by ZmMATE1 which encodes a multidrug and toxic compound extrusion (MATE) transporter protein. ZmMATE6 is another member of this family in maize whose expression is also increased by Al treatment. We investigated the function of this gene in more detail to determine whether it also contributes to Al resistance. Quantitative RT-PCR measurements found that ZmMATE6 was expressed in the roots and leaves of Al-resistant and sensitive inbred lines. Treatment with Al induced ZmMATE6 expression in all tissues but several other divalent or trivalent cations tested had no effect on expression. This expression pattern and the induction by Al treatment was confirmed in ZmMATE6 promoter-β-glucuronidase fusion lines. Heterogeneous expression of ZmMATE6 displayed a greater Al-activated release of citrate from the roots and was significantly resistant to Al toxicity than controls. This was associated with reduced accumulation of Al in the root tissues. Our results demonstrated that ZmMATE6 expression is induced by Al and functions as a citrate transporter.

MicroRNAs: Tiny, powerful players of metal stress responses in plants

Metal contamination of the environment is a widespread problem threatening sustainable and safe crop production. Physio-biochemical and molecular mechanisms of plant responses to metal exposure have been studied to establish the best possible agronomical or biotechnological methods to tackle metal contamination. Metal stress tolerance is regulated by several molecular effectors among which microRNAs are one of the key master regulators of plant growth and stress responses in plants. MicroRNAs are known to coordinate multitude of plant responses to metal stress through antioxidant functions, root growth, hormonal signalling, transcription factors and metal transporters. The present review discusses integrative functions of microRNAs in the regulation of metal stress in plants, which will be useful for engineering stress tolerance traits for improved plant growth and productivity in metal stressed situations.