Common sequence polymorphisms shaping genetic diversity in Arabidopsis thaliana

Comparative Chloroplast Genomics Reveals Intrageneric Divergence in Salix

As the most diverse genus of Salicaceae, Salix is primarily distributed in the temperate zone of the Northern Hemisphere, encompassing 350-500 species worldwide. The genus's evolutionary history is complex due to significant genetic differentiation. Chloroplast genes, being highly conserved, serve as effective tools for studying uniparental inheritance and evolution. In this study, we sequenced and assembled the chloroplast genomes of five representative Salix species. Phylogenetic relationships were constructed using chloroplast genome data, and structural differences among lineages were compared. These Salix chloroplast genomes exhibited a typical quadripartite structure, with lengths ranging from 154,444 to 155,725 bp. We successfully annotated 131 genes, including 88 protein-coding genes, 35 tRNA genes, and 8 rRNA genes. Clade I showed higher variability in the SSC region, identifying five highly variable regions: petA-psbJ, rps16-rps3, ndhD, ccsA-ndhD, and ndhG-ndhI. Two rapidly evolving genes, ndhI and ycf4, were also identified. The total length of insertions and deletions (InDels) in Clade I was 1046 bp. Clade II exhibited greater variability in the LSC region, with four highly variable regions being identified: trnK-trnQ, ndhC-trnV, trnV, and psdE-petL. Four rapidly evolving genes-infA, rpoC1, rps18, and ycf1-were identified. The total length of InDels in Clade II was 1282 bp. Therefore, this study elucidated the chloroplast genome evolution across different Salix lineages, thereby providing deeper insights into intrageneric phylogenetic relationships.

Evaluating the Roles of Drift and Selection in Trait Loss along an Elevational Gradient

Traits that have lost function sometimes persist through evolutionary time. These traits may persist if there is not enough standing genetic variation for the trait to allow a response to selection, if selection against the trait is weak relative to drift, or if the trait has a residual function. To determine the evolutionary processes shaping whether nonfunctional traits are retained or lost, we investigated short stamens in 16 populations of Arabidopsis thaliana along an elevational cline in northeast Spain. We found a cline in short stamen number from retention of short stamens in high elevation populations to incomplete loss in low elevation populations. We did not find evidence that limited genetic variation constrains the loss of short stamens at high elevations, nor evidence for divergent selection on short stamens between high and low elevations. Finally, we identified loci associated with short stamens in northeast Spain that are different from loci associated with variation in short stamen number across latitudes from a previous study. Overall, we did not identify the evolutionary mechanisms contributing to an elevational cline in short stamen number but did identify different genetic loci underlying variation in short stamen along similar phenotypic clines.

Exploring the complexity of genome size reduction in angiosperms

The genome sizes of angiosperms decreased significantly more than the genome sizes of their ancestors (pteridophytes and gymnosperms). Decreases in genome size involve a highly complex process, with remnants of the genome size reduction scattered across the genome and not directly linked to specific genomic structures. This is because the associated mechanisms operate on a much smaller scale than the mechanisms mediating increases in genome size. This review thoroughly summarizes the available literature regarding the molecular mechanisms underlying genome size reductions and introduces Utricularia gibba and Arabidopsis thaliana as model species for the examination of the effects of these molecular mechanisms. Additionally, we propose that phosphorus deficiency and drought stress are the major external factors contributing to decreases in genome size. Considering these factors affect almost all land plants, angiosperms likely gained the mechanisms for genome size reductions. These environmental factors may affect the retention rates of deletions, while also influencing the mutation rates of deletions via the functional diversification of the proteins facilitating double-strand break repair. The biased retention and mutation rates of deletions may have synergistic effects that enhance deletions in intergenic regions, introns, transposable elements, duplicates, and repeats, leading to a rapid decrease in genome size. We suggest that these selection pressures and associated molecular mechanisms may drive key changes in angiosperms during recurrent cycles of genome size decreases and increases.

Comparative phenotypic and transcriptomic analysis reveals genotypic differences in nitrogen use efficiency in sorghum

Sorghum (Sorghumbicolor L.), a model for C4 grass and an emerging biofuel crop, is known for its robust tolerance to low input field. However, the focus on enhancing nitrogen use efficiency (NUE) in sorghum under low nitrogen (N) conditions has been limited. This study conducted hydroponic experiments and field trials with two sorghum inbred lines, contrasting in their N efficiency: the N-efficient (398B) and the N-inefficient (CS3541) inbred lines. The aim was to analyze the key factors influencing NUE by integrating phenotypic, physiological, and multi-omics approaches under N deficiency conditions. The field experiments revealed that 398B displayed superior NUE and yield performance compared to CS3541. In hydroponic experiments, the growth of 398B outperformed CS3541 following N deficiency, attributing to its higher photosynthetic and sustaining activity of N metabolism-related enzymes. Genomic and transcriptomic integration highlighted fewer genomic diversities and alterations in global gene expression in 398B, which were likely contributor to its high NUE. Additionally, co-expression network analysis suggested the involvement of key genes which impact N uptake efficiency (NUpE) and N utilization efficiency (NUtE) in both lines, such as an N transporter, Sobic.003G371000.v3.2leaf(NPF5.10) and a transcription factor, Sobic.002G202800.v3.2leaf(WRKY) in bolstering NUE under low-N stress. The findings collectively suggested that 398B achieved higher NUpE and NUtE, effectively coordinating photosynthesis and N metabolism to enhance NUE. The candidate genes regulating N uptake and utilization efficiencies could provide valuable insights for developing sorghum breeds with improved NUE, contributing to sustainable agricultural practices and bioenergy crop development.

Recurrent selection shapes the genomic landscape of differentiation between a pair of host-specialized haplodiploids that diverged with gene flow

Understanding the genetics of adaptation and speciation is critical for a complete picture of how biodiversity is generated and maintained. Heterogeneous genomic differentiation between diverging taxa is commonly documented, with genomic regions of high differentiation interpreted as resulting from differential gene flow, linked selection and reduced recombination rates. Disentangling the roles of each of these non-exclusive processes in shaping genome-wide patterns of divergence is challenging but will enhance our knowledge of the repeatability of genomic landscapes across taxa. Here, we combine whole-genome resequencing and genome feature data to investigate the processes shaping the genomic landscape of differentiation for a sister-species pair of haplodiploid pine sawflies, Neodiprion lecontei and Neodiprion pinetum. We find genome-wide correlations between genome features and summary statistics are consistent with pervasive linked selection, with patterns of diversity and divergence more consistently predicted by exon density and recombination rate than the neutral mutation rate (approximated by dS). We also find that both global and local patterns of FST, dXY and π provide strong support for recurrent selection as the primary selective process shaping variation across pine sawfly genomes, with some contribution from balancing selection and lineage-specific linked selection. Because inheritance patterns for haplodiploid genomes are analogous to those of sex chromosomes, we hypothesize that haplodiploids may be especially prone to recurrent selection, even if gene flow occurred throughout divergence. Overall, our study helps fill an important taxonomic gap in the genomic landscape literature and contributes to our understanding of the processes that shape genome-wide patterns of genetic variation.

Comprehensive analysis of pepper (Capsicum annuum) RAV genes family and functional identification of CaRAV1 under chilling stress

BACKGROUND

Despite its known significance in plant abiotic stress responses, the role of the RAV gene family in the response of Capsicum annuum to chilling stress remains largely unexplored.

RESULTS

In this study, we identified and characterized six members of the CaRAV gene subfamily in pepper plants through genome-wide analysis. Subsequently, the CaRAV subfamily was classified into four branches based on homology with Arabidopsis thaliana, each exhibiting relatively conserved domains within the branch. We discovered that light response elements accounted for the majority of CaRAVs, whereas low-temperature response elements were specific to the NGA gene subfamily. After pepper plants were subjected to chilling stress, qRT‒PCR analysis revealed that CaRAV1, CaRAV2 and CaNGA1 were significantly induced in response to chilling stress, indicating that CaRAVs play a role in the response to chilling stress. Using virus-induced gene silencing (VIGS) vectors, we targeted key members of the CaRAV gene family. Under normal growth conditions, the MDA content and SOD enzyme activity of the silenced plants were slightly greater than those of the control plants, and the REC activity was significantly greater than that of the control plants. The levels of MDA and electrolyte leakage were greater in the silenced plants after they were exposed to chilling stress, and the POD and CAT enzyme activities were significantly lower than those in the control, which was particularly evident under repeated chilling stress. In addition, the relative expression of CaPOD and CaCAT was greater in V2 plants upon repeated chilling stress, especially CaCAT was significantly greater in V2 plants than in the other two silenced plants, with 3.29 and 1.10 increases within 12 and 24 h. These findings suggest that CaRAV1 and CaNGA1 positively regulate the response to chilling stress.

CONCLUSIONS

Silencing of key members of the CaRAV gene family results in increased susceptibility to chilling damage and reduced antioxidant enzyme activity in plants, particularly under repeated chilling stress. This study provides valuable information for understanding the classification and putative functions of RAV transcription factors in pepper plants.

Genome-wide association identifies a BAHD acyltransferase activity that assembles an ester of glucuronosylglycerol and phenylacetic acid

Genome-wide association studies (GWAS) are an effective approach to identify new specialized metabolites and the genes involved in their biosynthesis and regulation. In this study, GWAS of Arabidopsis thaliana soluble leaf and stem metabolites identified alleles of an uncharacterized BAHD-family acyltransferase (AT5G57840) associated with natural variation in three structurally related metabolites. These metabolites were esters of glucuronosylglycerol, with one metabolite containing phenylacetic acid as the acyl component of the ester. Knockout and overexpression of AT5G57840 in Arabidopsis and heterologous overexpression in Nicotiana benthamiana and Escherichia coli demonstrated that it is capable of utilizing phenylacetyl-CoA as an acyl donor and glucuronosylglycerol as an acyl acceptor. We, thus, named the protein Glucuronosylglycerol Ester Synthase (GGES). Additionally, phenylacetyl glucuronosylglycerol increased in Arabidopsis CYP79A2 mutants that overproduce phenylacetic acid and was lost in knockout mutants of UDP-sulfoquinovosyl: diacylglycerol sulfoquinovosyl transferase, an enzyme required for glucuronosylglycerol biosynthesis and associated with glycerolipid metabolism under phosphate-starvation stress. GGES is a member of a well-supported clade of BAHD family acyltransferases that arose by duplication and neofunctionalized during the evolution of the Brassicales within a larger clade that includes HCT as well as enzymes that synthesize other plant-specialized metabolites. Together, this work extends our understanding of the catalytic diversity of BAHD acyltransferases and uncovers a pathway that involves contributions from both phenylalanine and lipid metabolism.

An NLR paralog Pit2 generated from tandem duplication of Pit1 fine-tunes Pit1 localization and function

NLR family proteins act as intracellular receptors. Gene duplication amplifies the number of NLR genes, and subsequent mutations occasionally provide modifications to the second gene that benefits immunity. However, evolutionary processes after gene duplication and functional relationships between duplicated NLRs remain largely unclear. Here, we report that the rice NLR protein Pit1 is associated with its paralogue Pit2. The two are required for the resistance to rice blast fungus but have different functions: Pit1 induces cell death, while Pit2 competitively suppresses Pit1-mediated cell death. During evolution, the suppression of Pit1 by Pit2 was probably generated through positive selection on two fate-determining residues in the NB-ARC domain of Pit2, which account for functional differences between Pit1 and Pit2. Consequently, Pit2 lost its plasma membrane localization but acquired a new function to interfere with Pit1 in the cytosol. These findings illuminate the evolutionary trajectory of tandemly duplicated NLR genes after gene duplication.

Mitochondrial genomes of four slug moths (Lepidoptera, Limacodidae): Genome description and phylogenetic implications

The family Limacodidae belongs to the superfamily Zygaenoidea, which includes 1672 species commonly referred to as slug moths. Limacodidae larvae are major pests for many economically important plant species and can cause human dermatitis. At present, the structure of the mitochondrial genome (mitogenome), phylogenetic position, and adaptive evolution of slug moths are poorly understood. Herein, the mitogenomes of Parasa lepida, Phlossa conjuncta, Thosea sinensis, and Setora sinensis were sequenced and compared with other available mitogenome sequences to better characterize the mitogenomic diversity and evolution of this moth family. The mitogenomes of P. lepida, P. conjuncta, T. sinensis, and S. sinensis were confirmed to be circular in structure with lengths of 15,575 bp, 15,553 bp, 15,535 bp, and 15,529 bp, respectively. The Limacodidae mitogenomes exhibited similar nucleotide composition, codon usage, RNA structure, and control region patterns, indicating the conservation of the mitogenome in the family Limacodidae. A sliding window, Ka/Ks, and genetic distance analyses revealed that the atp8 and nad6 genes exhibited the highest levels of variability and the most rapid evolutionary rates among the 13 protein-coding genes (PCGs) encoded in these Limacodidae mitogenomes, suggesting that they may offer value as candidate DNA markers. The phylogenetic analysis recovered the overall relationship as Tortricoidea + (Sesiidae + (Zygaenoidea + (Cossoidea/+Choreutoidea + (others)))). Within Zygaenoidea, Limacodidae was recovered as monophyletic, and the phylogenetic relationships were recovered as (Phaudidae + Zyganidae) + Limacodidae in all six phylogenetic trees. The analysis indicated that P. lepida, P. conjuncta, T. sinensis, and S. sinensis are members of the Limacodidae.

High allelic diversity in Arabidopsis NLRs is associated with distinct genomic features

Plants rely on Nucleotide-binding, Leucine-rich repeat Receptors (NLRs) for pathogen recognition. Highly variable NLRs (hvNLRs) show remarkable intraspecies diversity, while their low-variability paralogs (non-hvNLRs) are conserved between ecotypes. At a population level, hvNLRs provide new pathogen-recognition specificities, but the association between allelic diversity and genomic and epigenomic features has not been established. Our investigation of NLRs in Arabidopsis Col-0 has revealed that hvNLRs show higher expression, less gene body cytosine methylation, and closer proximity to transposable elements than non-hvNLRs. hvNLRs show elevated synonymous and nonsynonymous nucleotide diversity and are in chromatin states associated with an increased probability of mutation. Diversifying selection maintains variability at a subset of codons of hvNLRs, while purifying selection maintains conservation at non-hvNLRs. How these features are established and maintained, and whether they contribute to the observed diversity of hvNLRs is key to understanding the evolution of plant innate immune receptors.

Pangenome-Wide Association Study and Transcriptome Analysis Reveal a Novel QTL and Candidate Genes Controlling both Panicle and Leaf Blast Resistance in Rice

Cultivating rice varieties with robust blast resistance is the most effective and economical way to manage the rice blast disease. However, rice blast disease comprises leaf and panicle blast, which are different in terms of resistance mechanisms. While many blast resistant rice cultivars were bred using genes conferring resistance to only leaf or panicle blast, mining durable and effective quantitative trait loci (QTLs) for both panicle and leaf blast resistance is of paramount importance. In this study, we conducted a pangenome-wide association study (panGWAS) on 9 blast resistance related phenotypes using 414 international diverse rice accessions from an international rice panel. This approach led to the identification of 74 QTLs associated with rice blast resistance. One notable locus, qPBR1, validated in a F4:5 population and fine-mapped in a Heterogeneous Inbred Family (HIF), exhibited broad-spectrum, major and durable blast resistance throughout the growth period. Furthermore, we performed transcriptomic analysis of 3 resistant and 3 sensitive accessions at different time points after infection, revealing 3,311 differentially expressed genes (DEGs) potentially involved in blast resistance. Integration of the above results identified 6 candidate genes within the qPBR1 locus, with no significant negative effect on yield. The results of this study provide valuable germplasm resources, QTLs, blast response genes and candidate functional genes for developing rice varieties with enduring and broad-spectrum blast resistance. The qPBR1, in particular, holds significant potential for breeding new rice varieties with comprehensive and durable resistance throughout their growth period.

The NLR family of innate immune and cell death sensors

Nucleotide-binding oligomerization domain (NOD)-like receptors, also known as nucleotide-binding leucine-rich repeat receptors (NLRs), are a family of cytosolic pattern recognition receptors that detect a wide variety of pathogenic and sterile triggers. Activation of specific NLRs initiates pro- or anti-inflammatory signaling cascades and the formation of inflammasomes-multi-protein complexes that induce caspase-1 activation to drive inflammatory cytokine maturation and lytic cell death, pyroptosis. Certain NLRs and inflammasomes act as integral components of larger cell death complexes-PANoptosomes-driving another form of lytic cell death, PANoptosis. Here, we review the current understanding of the evolution, structure, and function of NLRs in health and disease. We discuss the concept of NLR networks and their roles in driving cell death and immunity. An improved mechanistic understanding of NLRs may provide therapeutic strategies applicable across infectious and inflammatory diseases and in cancer.

Population structure and genetic diversity in Eucalyptus pellita based on SNP markers

Eucalyptus pellita has the characteristics of rapid growth and high resistance. However, there is little research on molecular breeding of E. pellita, which is essential to shortening breeding life and selecting quality varieties. Therefore, a crucial step before selective breeding can be carried out to increase the wood quality of E. pellita is identifying genetic diversity and population structure using single nucleotide polymorphism (SNP) markers. In this study, the genetic diversity of 1st generation 196 E. pellita families from 23 geographically defined was assessed using 1,677,732 SNP markers identified by whole genome resequencing. SNP annotation showed that the ratio of non-synonymous to synonymous coding mutations was 0.83. Principal component analysis (PCA), phylogenetic tree, and population structure analysis permitted the families to be categorized into three groups, one of which (G2) contains most of the Indonesian (IDN) and Papua New Guinea (PNG) families. Genetic relationship analysis showed that IDN was closely related to PNG. Genetic diversity analysis showed that He, PIC, I, and H mean values were 0.2502, 0.2027, 0.3815, and 0.2680, respectively. PCA analysis classified various provenances in QLD into two categories (G1 and G3). The genetic diversity of G3 was higher than that of G2. The results of genetic differentiation (Fst) showed that PNG region was divided into two groups (PNG1 and PNG2), the Fst (0.172) between QLD and PNG2 region was higher than QLD and PNG1, and the Fst (0.024) between IDN and PNG1 is smaller than IDN and PNG2. A Mantel test revealed a positive correlation between the genetic and geographic distance of E. pellita. This study has a certain reference value for genetic identification, germplasm preservation, and breeding of E. pellita. Also, it provides a basis for subsequent association analysis to explore excellent alleles and introduction.

NLR receptors in plant immunity: making sense of the alphabet soup

Plants coordinately use cell-surface and intracellular immune receptors to perceive pathogens and mount an immune response. Intracellular events of pathogen recognition are largely mediated by immune receptors of the nucleotide binding and leucine rich-repeat (NLR) classes. Upon pathogen perception, NLRs trigger a potent broad-spectrum immune reaction, usually accompanied by a form of programmed cell death termed the hypersensitive response. Some plant NLRs act as multifunctional singleton receptors which combine pathogen detection and immune signaling. However, NLRs can also function in higher order pairs and networks of functionally specialized interconnected receptors. In this article, we cover the basic aspects of plant NLR biology with an emphasis on NLR networks. We highlight some of the recent advances in NLR structure, function, and activation and discuss emerging topics such as modulator NLRs, pathogen suppression of NLRs, and NLR bioengineering. Multi-disciplinary approaches are required to disentangle how these NLR immune receptor pairs and networks function and evolve. Answering these questions holds the potential to deepen our understanding of the plant immune system and unlock a new era of disease resistance breeding.

Pan-European study of genotypes and phenotypes in the Arabidopsis relative Cardamine hirsuta reveals how adaptation, demography, and development shape diversity patterns

We study natural DNA polymorphisms and associated phenotypes in the Arabidopsis relative Cardamine hirsuta. We observed strong genetic differentiation among several ancestry groups and broader distribution of Iberian relict strains in European C. hirsuta compared to Arabidopsis. We found synchronization between vegetative and reproductive development and a pervasive role for heterochronic pathways in shaping C. hirsuta natural variation. A single, fast-cycling ChFRIGIDA allele evolved adaptively allowing range expansion from glacial refugia, unlike Arabidopsis where multiple FRIGIDA haplotypes were involved. The Azores islands, where Arabidopsis is scarce, are a hotspot for C. hirsuta diversity. We identified a quantitative trait locus (QTL) in the heterochronic SPL9 transcription factor as a determinant of an Azorean morphotype. This QTL shows evidence for positive selection, and its distribution mirrors a climate gradient that broadly shaped the Azorean flora. Overall, we establish a framework to explore how the interplay of adaptation, demography, and development shaped diversity patterns of 2 related plant species.

Differences in the intraspecies copy number variation of Arabidopsis thaliana conserved and nonconserved miRNA genes

MicroRNAs (miRNAs) regulate gene expression by RNA interference mechanism. In plants, miRNA genes (MIRs) which are grouped into conserved families, i.e. they are present among the different plant taxa, are involved in the regulation of many developmental and physiological processes. The roles of the nonconserved MIRs-which are MIRs restricted to one plant family, genus, or even species-are less recognized; however, many of them participate in the responses to biotic and abiotic stresses. Both over- and underproduction of miRNAs may influence various biological processes. Consequently, maintaining intracellular miRNA homeostasis seems to be crucial for the organism. Deletions and duplications in the genomic sequence may alter gene dosage and/or activity. We evaluated the extent of copy number variations (CNVs) among Arabidopsis thaliana (Arabidopsis) MIRs in over 1000 natural accessions, using population-based analysis of the short-read sequencing data. We showed that the conserved MIRs were unlikely to display CNVs and their deletions were extremely rare, whereas nonconserved MIRs presented moderate variation. Transposon-derived MIRs displayed exceptionally high diversity. Conversely, MIRs involved in the epigenetic control of transposons reactivated during development were mostly invariable. MIR overlap with the protein-coding genes also limited their variability. At the expression level, a higher rate of nonvariable, nonconserved miRNAs was detectable in Col-0 leaves, inflorescence, and siliques compared to nonconserved variable miRNAs, although the expression of both groups was much lower than that of the conserved MIRs. Our data indicate that CNV rate of Arabidopsis MIRs is related with their age, function, and genomic localization.

Exploring the sorghum race level diversity utilizing 272 sorghum accessions genomic resources

Due to evolutionary divergence, sorghum race populations exhibit significant genetic and morphological variation. A k-mer-based sorghum race sequence comparison identified the conserved k-mers of all 272 accessions from sorghum and the race-specific genetic signatures identified the gene variability in 10,321 genes (PAVs). To understand sorghum race structure, diversity and domestication, a deep learning-based variant calling approach was employed in a set of genotypic data derived from a diverse panel of 272 sorghum accessions. The data resulted in 1.7 million high-quality genome-wide SNPs and identified selective signature (both positive and negative) regions through a genome-wide scan with different (iHS and XP-EHH) statistical methods. We discovered 2,370 genes associated with selection signatures including 179 selective sweep regions distributed over 10 chromosomes. Co-localization of these regions undergoing selective pressure with previously reported QTLs and genes revealed that the signatures of selection could be related to the domestication of important agronomic traits such as biomass and plant height. The developed k-mer signatures will be useful in the future to identify the sorghum race and for trait and SNP markers for assisting in plant breeding programs.

Molecular mechanisms of adaptive evolution in wild animals and plants

Wild animals and plants have developed a variety of adaptive traits driven by adaptive evolution, an important strategy for species survival and persistence. Uncovering the molecular mechanisms of adaptive evolution is the key to understanding species diversification, phenotypic convergence, and inter-species interaction. As the genome sequences of more and more non-model organisms are becoming available, the focus of studies on molecular mechanisms of adaptive evolution has shifted from the candidate gene method to genetic mapping based on genome-wide scanning. In this study, we reviewed the latest research advances in wild animals and plants, focusing on adaptive traits, convergent evolution, and coevolution. Firstly, we focused on the adaptive evolution of morphological, behavioral, and physiological traits. Secondly, we reviewed the phenotypic convergences of life history traits and responding to environmental pressures, and the underlying molecular convergence mechanisms. Thirdly, we summarized the advances of coevolution, including the four main types: mutualism, parasitism, predation and competition. Overall, these latest advances greatly increase our understanding of the underlying molecular mechanisms for diverse adaptive traits and species interaction, demonstrating that the development of evolutionary biology has been greatly accelerated by multi-omics technologies. Finally, we highlighted the emerging trends and future prospects around the above three aspects of adaptive evolution.

Allelic variation in the Arabidopsis TNL CHS3/CSA1 immune receptor pair reveals two functional cell-death regulatory modes

Some plant NLR immune receptors are encoded in head-to-head "sensor-executor" pairs that function together. Alleles of the NLR pair CHS3/CSA1 form three clades. The clade 1 sensor CHS3 contains an integrated domain (ID) with homology to regulatory domains, which is lacking in clades 2 and 3. In this study, we defined two cell-death regulatory modes for CHS3/CSA1 pairs. One is mediated by ID domain on clade 1 CHS3, and the other relies on CHS3/CSA1 pairs from all clades detecting perturbation of an associated pattern-recognition receptor (PRR) co-receptor. Our data support the hypothesis that an ancestral Arabidopsis CHS3/CSA1 pair gained a second recognition specificity and regulatory mechanism through ID acquisition while retaining its original specificity as a "guard" against PRR co-receptor perturbation. This likely comes with a cost, since both ID and non-ID alleles of the pair persist in diverse Arabidopsis populations through balancing selection.

Adaptation in Outbred Sexual Yeast is Repeatable, Polygenic and Favors Rare Haplotypes

We carried out a 200 generation Evolve and Resequence (E&R) experiment initiated from an outbred diploid recombined 18-way synthetic base population. Replicate populations were evolved at large effective population sizes (>105 individuals), exposed to several different chemical challenges over 12 weeks of evolution, and whole-genome resequenced. Weekly forced outcrossing resulted in an average between adjacent-gene per cell division recombination rate of ∼0.0008. Despite attempts to force weekly sex, roughly half of our populations evolved cheaters and appear to be evolving asexually. Focusing on seven chemical stressors and 55 total evolved populations that remained sexual we observed large fitness gains and highly repeatable patterns of genome-wide haplotype change within chemical challenges, with limited levels of repeatability across chemical treatments. Adaptation appears highly polygenic with almost the entire genome showing significant and consistent patterns of haplotype change with little evidence for long-range linkage disequilibrium in a subset of populations for which we sequenced haploid clones. That is, almost the entire genome is under selection or drafting with selected sites. At any given locus adaptation was almost always dominated by one of the 18 founder's alleles, with that allele varying spatially and between treatments, suggesting that selection acts primarily on rare variants private to a founder or haplotype blocks harboring multiple mutations.

Strategies for Intraspecific Genotyping of Duckweed: Comparison of Five Orthogonal Methods Applied to the Giant Duckweed Spirodela polyrhiza

The predominantly vegetative propagating duckweeds are of growing commercial interest. Since clonal accessions within a respective species can vary considerably with respect to their physiological as well as biochemical traits, it is critical to be able to track the clones of species of interest after their characterization. Here, we compared the efficacy of five different genotyping methods for Spirodela polyrhiza, a species with very low intraspecific sequence variations, including polymorphic NB-ARC-related loci, tubulin-gene-based polymorphism (TBP), simple sequence repeat variations (SSR), multiplexed ISSR genotyping by sequencing (MIG-seq), and low-coverage, reduced-representation genome sequencing (GBS). Four of the five approaches could distinguish 20 to 22 genotypes out of the 23 investigated clones, while TBP resolved just seven genotypes. The choice for a particular method for intraspecific genotyping can depend on the research question and the project budget, while the combination of orthogonal methods may increase the confidence and resolution for the results obtained.

The Boechera model system for evolutionary ecology

Model systems in biology expand the research capacity of individuals and the community. Closely related to Arabidopsis, the genus Boechera has emerged as an important ecological model owing to the ability to integrate across molecular, functional, and eco-evolutionary approaches. Boechera species are broadly distributed in relatively undisturbed habitats predominantly in western North America and provide one of the few experimental systems for identification of ecologically important genes through genome-wide association studies and investigations of selection with plants in their native habitats. The ecologically, evolutionarily, and agriculturally important trait of apomixis (asexual reproduction via seeds) is common in the genus, and field experiments suggest that abiotic and biotic environments shape the evolution of sex. To date, population genetic studies have focused on the widespread species B. stricta, detailing population divergence and demographic history. Molecular and ecological studies show that balancing selection maintains genetic variation in ~10% of the genome, and ecological trade-offs contribute to complex trait variation for herbivore resistance, flowering phenology, and drought tolerance. Microbiome analyses have shown that host genotypes influence leaf and root microbiome composition, and the soil microbiome influences flowering phenology and natural selection. Furthermore, Boechera offers numerous opportunities for investigating biological responses to global change. In B. stricta, climate change has induced a shift of >2 weeks in the timing of first flowering since the 1970s, altered patterns of natural selection, generated maladaptation in previously locally-adapted populations, and disrupted life history trade-offs. Here we review resources and results for this eco-evolutionary model system and discuss future research directions.

A genetically linked pair of NLR immune receptors shows contrasting patterns of evolution

Throughout their evolution, plant nucleotide-binding leucine-rich-repeat receptors (NLRs) have acquired widely divergent unconventional integrated domains that enhance their ability to detect pathogen effectors. However, the functional dynamics that drive the evolution of NLRs with integrated domains (NLR-IDs) remain poorly understood. Here, we reconstructed the evolutionary history of an NLR locus prone to unconventional domain integration and experimentally tested hypotheses about the evolution of NLR-IDs. We show that the rice (Oryza sativa) NLR Pias recognizes the effector AVR-Pias of the blast fungal pathogen Magnaporthe oryzae. Pias consists of a functionally specialized NLR pair, the helper Pias-1 and the sensor Pias-2, that is allelic to the previously characterized Pia pair of NLRs: the helper RGA4 and the sensor RGA5. Remarkably, Pias-2 carries a C-terminal DUF761 domain at a similar position to the heavy metal-associated (HMA) domain of RGA5. Phylogenomic analysis showed that Pias-2/RGA5 sensor NLRs have undergone recurrent genomic recombination within the genus Oryza, resulting in up to six sequence-divergent domain integrations. Allelic NLRs with divergent functions have been maintained transspecies in different Oryza lineages to detect sequence-divergent pathogen effectors. By contrast, Pias-1 has retained its NLR helper activity throughout evolution and is capable of functioning together with the divergent sensor-NLR RGA5 to respond to AVR-Pia. These results suggest that opposite selective forces have driven the evolution of paired NLRs: highly dynamic domain integration events maintained by balancing selection for sensor NLRs, in sharp contrast to purifying selection and functional conservation of immune signaling for helper NLRs.

BLSSpeller to discover novel regulatory motifs in maize

With the decreasing cost of sequencing and availability of larger numbers of sequenced genomes, comparative genomics is becoming increasingly attractive to complement experimental techniques for the task of transcription factor (TF) binding site identification. In this study, we redesigned BLSSpeller, a motif discovery algorithm, to cope with larger sequence datasets. BLSSpeller was used to identify novel motifs in Zea mays in a comparative genomics setting with 16 monocot lineages. We discovered 61 motifs of which 20 matched previously described motif models in Arabidopsis. In addition, novel, yet uncharacterized motifs were detected, several of which are supported by available sequence-based and/or functional data. Instances of the predicted motifs were enriched around transcription start sites and contained signatures of selection. Moreover, the enrichment of the predicted motif instances in open chromatin and TF binding sites indicates their functionality, supported by the fact that genes carrying instances of these motifs were often found to be co-expressed and/or enriched in similar GO functions. Overall, our study unveiled several novel candidate motifs that might help our understanding of the genotype to phenotype association in crops.

Genomic selection and genetic architecture of agronomic traits during modern rapeseed breeding

Rapeseed (Brassica napus L.) is an important oil-producing crop for the world. Its adaptation, yield and quality have been considerably improved in recent decades, but the genomic basis underlying successful breeding selection remains unclear. Hence, we conducted a comprehensive genomic assessment of rapeseed in the breeding process based on the whole-genome resequencing of 418 diverse rapeseed accessions. We unraveled the genomic basis for the selection of adaptation and agronomic traits. Genome-wide association studies identified 628 associated loci-related causative candidate genes for 56 agronomically important traits, including plant architecture and yield traits. Furthermore, we uncovered nonsynonymous mutations in plausible candidate genes for agronomic traits with significant differences in allele frequency distributions across the improvement process, including the ribosome recycling factor (BnRRF) gene for seed weight. This study provides insights into the genomic basis for improving rapeseed varieties and a valuable genomic resource for genome-assisted rapeseed breeding.

Rapid intraspecies evolution of fitness effects of yeast genes

Organisms within species have numerous genetic and phenotypic variations. Growing evidences show intraspecies variation of mutant phenotypes may be more complicated than expected. Current studies on intraspecies variations of mutant phenotypes are limited to just a few strains. This study investigated the intraspecies variation of fitness effects of 5,630 gene mutants in ten Saccharomyces cerevisiae strains using CRISPR–Cas9 screening. We found that the variability of fitness effects induced by gene disruptions is very large across different strains. Over 75% of genes affected cell fitness in a strain-specific manner to varying degrees. The strain specificity of the fitness effect of a gene is related to its evolutionary and functional properties. Subsequent analysis revealed that younger genes, especially those newly acquired in S. cerevisiae species, are more likely to be strongly strain-specific. Intriguingly, there seems to exist a ceiling of fitness effect size for strong strain-specific genes, and among them, the newly acquired genes are still evolving and have yet to reach this ceiling. Additionally, for a large proportion of protein complexes, the strain specificity profile is inconsistent among genes encoding the same complex. Taken together, these results offer a genome-wide map of intraspecies variation for fitness effect as a mutant phenotype and provide an updated insight on intraspecies phenotypic evolution.

Global Identification and Characterization of C2 Domain-Containing Proteins Associated with Abiotic Stress Response in Rice (Oryza sativa L.)

C2 domain-containing proteins (C2DPs) have been identified in different genomes that contain single or multiple C2 domains in their C- or N-terminal. It possesses higher functional activity in the transmembrane regions. The identification of C2 domains were reported in a previous study, such as multiple C2 domains and transmembrane-region proteins (MCTPs) and N-terminal-TM-C2 domain proteins (NTMC2s) of rice, Arabidopsis thaliana, and cotton, whereas the C2DP gene family in rice has not been comprehensively studied, and the role of the C2DP gene in rice in response to abiotic stress is not yet fully understood. In this study, we identified 82 C2DPs in the rice genome and divided them into seven groups through phylogenetic analysis. The synteny analysis revealed that duplication events were either exhibited within the genome of rice or between the genomes of rice and other species. Through the analysis of cis-acting elements in promoters, expression profiles, and qRT-PCR results, the functions of OsC2DPs were found to be widely distributed in diverse tissues and were extensively involved in phytohormones-related and abiotic stresses response in rice. The prediction of the microRNA (miRNA) targets of OsC2DPs revealed the possibility of regulation by consistent miRNAs. Notably, OsC2DP50/51/52 as a co-tandem duplication exhibited similar expression variations and involved the coincident miRNA-regulation pathway. Moreover, the results of the genotypic variation and haplotype analysis revealed that OsC2DP17, OsC2DP29, and OsC2DP49 were associated with cold stress responses. These findings provided comprehensive insights for characterizations of OsC2DPs in rice as well as for their roles for abiotic stress.

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.

Analysis of 427 genomes reveals moso bamboo population structure and genetic basis of property traits

Moso bamboo (Phyllostachys edulis) is an economically and ecologically important nontimber forestry species. Further development of this species as a sustainable bamboo resource has been hindered by a lack of population genome information. Here, we report a moso bamboo genomic variation atlas of 5.45 million single-nucleotide polymorphisms (SNPs) from whole-genome resequencing of 427 individuals covering 15 representative geographic areas. We uncover low genetic diversity, high genotype heterozygosity, and genes under balancing selection underlying moso bamboo population adaptation. We infer its demographic history with one bottleneck and its recently small population without a rebound. We define five phylogenetic groups and infer that one group probably originated by a single-origin event from East China. Finally, we conduct genome-wide association analysis of nine important property-related traits to identify candidate genes, many of which are involved in cell wall, carbohydrate metabolism, and environmental adaptation. These results provide a foundation and resources for understanding moso bamboo evolution and the genetic mechanisms of agriculturally important traits.

Zea mays RNA-seq estimated transcript abundances are strongly affected by read mapping bias

Background

Genetic variation for gene expression is a source of phenotypic variation for natural and agricultural species. The common approach to map and to quantify gene expression from genetically distinct individuals is to assign their RNA-seq reads to a single reference genome. However, RNA-seq reads from alleles dissimilar to this reference genome may fail to map correctly, causing transcript levels to be underestimated. Presently, the extent of this mapping problem is not clear, particularly in highly diverse species. We investigated if mapping bias occurred and if chromosomal features associated with mapping bias. Zea mays presents a model species to assess these questions, given it has genotypically distinct and well-studied genetic lines.

Results

In Zea mays, the inbred B73 genome is the standard reference genome and template for RNA-seq read assignments. In the absence of mapping bias, B73 and a second inbred line, Mo17, would each have an approximately equal number of regulatory alleles that increase gene expression. Remarkably, Mo17 had 2–4 times fewer such positively acting alleles than did B73 when RNA-seq reads were aligned to the B73 reference genome. Reciprocally, over one-half of the B73 alleles that increased gene expression were not detected when reads were aligned to the Mo17 genome template. Genes at dissimilar chromosomal ends were strongly affected by mapping bias, and genes at more similar pericentromeric regions were less affected. Biased transcript estimates were higher in untranslated regions and lower in splice junctions. Bias occurred across software and alignment parameters.

Conclusions

Mapping bias very strongly affects gene transcript abundance estimates in maize, and bias varies across chromosomal features. Individual genome or transcriptome templates are likely necessary for accurate transcript estimation across genetically variable individuals in maize and other species.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07577-3.

Genomic analyses of diverse wild and cultivated accessions provide insights into the evolutionary history of jujube

The Chinese jujube (Ziziphus jujuba Mill.), a member of the Rhamnaceae family, is an important perennial fruit tree crop of substantial economic, ecological and nutritional value, and is also used as a traditional herbal medicine. Here, we report the resequencing of 493 jujube accessions, including 202 wild and 291 cultivated accessions at >16× depth. Our population genomic analyses revealed that the Shanxi-Shaanxi area of China was jujube's primary domestication centre and that jujube was then disseminated into East China before finally extending into South China. Divergence events analysis indicated that Ziziphus acidojujuba and Ziziphus jujuba diverged around 2.7 Mya, suggesting the interesting possibility that a long pre-domestication period may have occurred prior to human intervention. Using the large genetic polymorphism data set, we identified a 15-bp tandem insertion in the promoter of the jujube ortholog of the POLLEN DEFECTIVE IN GUIDANCE 1 (POD1) gene, which was strongly associated with seed-setting rate. Integrating genome-wide association study (GWAS), transcriptome data, expression analysis and transgenic validation in tomato, we identified a DA3/UBIQUITIN-SPECIFIC PROTEASE 14 (UBP14) ortholog, which negatively regulate fruit weight in jujube. We also identified candidate genes, which have likely influenced the selection of fruit sweetness and crispness texture traits among fresh and dry jujubes. Our study not only illuminates the genetic basis of jujube evolution and domestication and provides a deep and rich genomic resource to facilitate both crop improvement and hypothesis-driven basic research, but also identifies multiple agriculturally important genes for this unique perennial tree fruit species.

Soybean (Glycine max) Haplotype Map (GmHapMap): a universal resource for soybean translational and functional genomics

Here, we describe a worldwide haplotype map for soybean (GmHapMap) constructed using whole-genome sequence data for 1007 Glycine max accessions and yielding 14.9 million variants as well as 4.3 M tag single-nucleotide polymorphisms (SNPs). When sampling random subsets of these accessions, the number of variants and tag SNPs plateaued beyond approximately 800 and 600 accessions, respectively. This suggests extensive coverage of diversity within the cultivated soybean. GmHapMap variants were imputed onto 21 618 previously genotyped accessions with up to 96% success for common alleles. A local association analysis was performed with the imputed data using markers located in a 1-Mb region known to contribute to seed oil content and enabled us to identify a candidate causal SNP residing in the NPC1 gene. We determined gene-centric haplotypes (407 867 GCHs) for the 55 589 genes and showed that such haplotypes can help to identify alleles that differ in the resulting phenotype. Finally, we predicted 18 031 putative loss-of-function (LOF) mutations in 10 662 genes and illustrated how such a resource can be used to explore gene function. The GmHapMap provides a unique worldwide resource for applied soybean genomics and breeding.

Genetics of autoimmunity in plants: an evolutionary genetics perspective

Autoimmunity in plants has been found in numerous hybrids as a form of hybrid necrosis and mutant panels. Uncontrolled cell death is a main cellular outcome of autoimmunity, which negatively impacts growth. Its occurrence highlights the vulnerable nature of the plant immune system. Genetic investigation of autoimmunity in hybrid plants revealed that extreme variation in the immune receptor repertoire is a major contributor, reflecting an evolutionary conundrum that plants face in nature. In this review, we discuss natural variation in the plant immune system and its contribution to fitness. The value of autoimmunity genetics lies in its ability to identify combinations of a natural immune receptor and its partner that are predisposed to triggering autoimmunity. The network of immune components for autoimmunity becomes instrumental in revealing mechanistic details of how immune receptors recognize cellular invasion and activate signaling. The list of autoimmunity-risk variants also allows us to infer evolutionary processes contributing to their maintenance in the natural population. Our approach to autoimmunity, which integrates mechanistic understanding and evolutionary genetics, has the potential to serve as a prognosis tool to optimize immunity in crops.

Gramene 2021: harnessing the power of comparative genomics and pathways for plant research

Gramene (http://www.gramene.org), a knowledgebase founded on comparative functional analyses of genomic and pathway data for model plants and major crops, supports agricultural researchers worldwide. The resource is committed to open access and reproducible science based on the FAIR data principles. Since the last NAR update, we made nine releases; doubled the genome portal's content; expanded curated genes, pathways and expression sets; and implemented the Domain Informational Vocabulary Extraction (DIVE) algorithm for extracting gene function information from publications. The current release, #63 (October 2020), hosts 93 reference genomes-over 3.9 million genes in 122 947 families with orthologous and paralogous classifications. Plant Reactome portrays pathway networks using a combination of manual biocuration in rice (320 reference pathways) and orthology-based projections to 106 species. The Reactome platform facilitates comparison between reference and projected pathways, gene expression analyses and overlays of gene-gene interactions. Gramene integrates ontology-based protein structure-function annotation; information on genetic, epigenetic, expression, and phenotypic diversity; and gene functional annotations extracted from plant-focused journals using DIVE. We train plant researchers in biocuration of genes and pathways; host curated maize gene structures as tracks in the maize genome browser; and integrate curated rice genes and pathways in the Plant Reactome.

The patterns of deleterious mutations during the domestication of soybean

Globally, soybean is a major protein and oil crop. Enhancing our understanding of the soybean domestication and improvement process helps boost genomics-assisted breeding efforts. Here we present a genome-wide variation map of 10.6 million single-nucleotide polymorphisms and 1.4 million indels for 781 soybean individuals which includes 418 domesticated (Glycine max), 345 wild (Glycine soja), and 18 natural hybrid (G. max/G. soja) accessions. We describe the enhanced detection of 183 domestication-selective sweeps and the patterns of putative deleterious mutations during domestication and improvement. This predominantly selfing species shows 7.1% reduction of overall deleterious mutations in domesticated soybean relative to wild soybean and a further 1.4% reduction from landrace to improved accessions. The detected domestication-selective sweeps also show reduced levels of deleterious alleles. Importantly, genotype imputation with this resource increases the mapping resolution of genome-wide association studies for seed protein and oil traits in a soybean diversity panel.

RNA-Seq transcriptome analysis to identify candidate genes involved in non-target site-based mesosulfuron-methyl resistance in Beckmannia syzigachne

American sloughgrass (Beckmannia syzigachne Steud.) has become a dominant weed in fields with rice-wheat rotation. Moreover, herbicide resistance has rendered weed control difficult. We identified a biotype showing resistance to ALS inhibitor mesosulfuron-methyl with a resistant index 3.3, but without any ALS mutation. This study aims to identify and confirm the factors associated with non-target site resistance of this biotype to mesosulfuron-methyl using RNA-Seq. 118,111 unigenes were assembled, and 50.9% of these were annotated across seven databases. Eleven contigs related to metabolic resistance were identified based on differential expression via RNA-Seq which include a novel resistance-related transcription factor (MYC3) and two disease resistance proteins were also identified (At1g58602 and At1g15890). Fold changes in expression of these genes in comparison M-R vs. M-S ranged from 3.9 to 11.6, as confirmed by qPCR. The expression of a contig annotated as cytochrome P450 (CYP86B1) in resistant individuals was over 3 times higher than that in sensitive individuals at 0-72 h after mesosulfuron-methyl treatment. A similar trend was noted for three other genes annotated as glutathione S-transferase (GST), namely GST-T3, GST-U6, and GST-U14; the expression of GST-U6 in resistant individuals was up to 142.3 times higher than that in sensitive individuals at 24 h after mesosulfuron-methyl treatment. In addition, GST activity in resistant individuals was 2.1 to 5.3 times higher than that in sensitive individuals. The GR₅₀ of resistant biotype decreased from 24.4 to 11.3 g a.i. ha^-1 after P450 inhibitor malathion treatment. This study identified a cytochrome P450 gene CYP86B1 and three GST genes GST-T3, GST-U6, and GST-U14 that have higher expression in mesosulfuron-methyl resistant B. syzigachne, suggesting that both P450- and GST-based activities could be involved in resistance.

Advancement of research on plant NLRs evolution, biochemical activity, structural association, and engineering

In this review, we have included evolution of plant intracellular immune receptors, oligomeric complex formation, enzymatic action, engineering, and mechanisms of immune inspection for appropriate defense outcomes. NLR (Nucleotide binding oligomerization domain containing leucine-rich repeat) proteins are the intracellular immune receptors that recognize pathogen-derived virulence factors to confer effector-triggered immunity (ETI). Activation of plant defense by the NLRs are often conveyed through N-terminal Toll-like/ IL-1 receptor (TIR) or non-TIR (coiled-coils or CC) domains. Homodimerization or self-association property of CC/ TIR domains of plant NLRs contribute to their auto-activity and induction of in planta ectopic cell death. High resolution crystal structures of Arabidopsis thaliana RPS4^TIR, L6^TIR, SNC1^TIR, RPP1^TIR and Muscadinia rotundifolia RPV1^TIR showed that interaction is mediated through one or two distinct interfaces i.e., αA and αE helices comprise AE interface and αD and αE helices were found to form DE interface. By contrast, conserved helical regions were determined for CC domains of plant NLRs. Evolutionary history of NLRs diversification has shown that paired forms were originated from NLR singletons. Plant TIRs executed NAD⁺ hydrolysis activity for cell death promotion. Plant NLRs were found to form large oligomeric complexes as observed in animal inflammasomes. We have also discussed different protein engineering methods includes domain shuffling, and decoy modification that increase effector recognition spectrum of plant NLRs. In summary, our review highlights structural basis of perception of the virulence factors by NLRs or NLR pairs to design novel classes of plant immune receptors.

High quality genome of Erigeron breviscapus provides a reference for herbal plants in Asteraceae

Erigeron breviscapus is an important medicinal plant in Compositae and the first species to realize the whole process from the decoding of the draft genome sequence to scutellarin biosynthesis in yeast. However, the previous low‐quality genome assembly has hindered the optimization of candidate genes involved in scutellarin synthesis and the development of molecular‐assisted breeding based on the genome. Here, the E. breviscapus genome was updated using PacBio RSII sequencing data and Hi‐C data, and increased in size from 1.2 Gb to 1.43 Gb, with a scaffold N50 of 156.82 Mb and contig N50 of 140.95 kb, and a total of 43,514 protein‐coding genes were obtained and oriented onto nine pseudo‐chromosomes, thus becoming the third plant species assembled to chromosome level after sunflower and lettuce in Compositae. Fourteen genes with evidence for positive selection were identified and found to be related to leaf morphology, flowering and secondary metabolism. The number of genes in some gene families involved in flavonoid biosynthesis in E. breviscapus have been significantly expanded. In particular, additional candidate genes involved in scutellarin biosynthesis, such as flavonoid‐7‐O‐glucuronosyltransferase genes (F7GATs) were identified using updated genome. In addition, three candidate genes encoding indole‐3‐pyruvate monooxygenase YUCCA2 (YUC2), serine carboxypeptidase‐like 18 (SCPL18), and F‐box protein (FBP), respectively, were identified to be probably related to leaf development and flowering by resequencing 99 individuals. These results provided a substantial genetic basis for improving agronomic and quality traits of E. breviscapus, and provided a platform for improving other draft genome assemblies to chromosome‐level.

Pro-197-Ser Mutation in ALS and High-Level GST Activities: Multiple Resistance to ALS and ACCase Inhibitors in Beckmannia syzigachne

American sloughgrass (Beckmannia syzigachne Steud.) is one of the most troublesome weeds infesting wheat and canola fields in China. Some biotypes cannot be controlled, either by acetolactate synthase (ALS) or acetyl coenzyme A carboxylase (ACCase) inhibitors, which are the main herbicides for controlling this weed. However, very few studies have investigated multiple resistance mechanism in B. syzigachne. In this study, a B. syzigachne biotype with a high resistance to ALS inhibitors we have reported was also showed relatively lower resistance to ACCase inhibitors, with a resistance index around 7. RNA-seq analysis was used to investigate the factors responsible for multiple resistance, and 60,108 unigenes were assembled by de novo transcriptome assembly and then annotated across eight databases. A Pro-197-Ser mutation was identified in the ALS gene by SNPs analysis and validated by PCR, while no mutation was identified in the ACCase gene. Nineteen candidate metabolic genes were screened and their overexpression was confirmed by qPCR. The expression of GST-T3 and GST-U6 in resistant plants ranged from 7.5- to 109.4-folds than that in susceptible ones at different times after two kinds of herbicide treatment. In addition, GST activities in resistant plants were 3.0-5.0 times higher than that in susceptible plants. Other novel resistance factors also showed high correlation with multiple resistance which included four genes encoding disease resistance proteins, a transcription factor (MYC3), and one gene conferring blight resistance. In this research, a B. syzigachne biotype was confirmed to have evolved multiple resistance to ACCase and ALS inhibitors. The Pro-197-Ser mutation in ALS gene and high-level GST activities were confirmed responsible for the multiple resistance. Characterized disease-resistance proteins, transcription factor, and blight-resistance proteins may play an essential role in these multiple herbicide resistance.

Degradation of the Repetitive Genomic Landscape in a Close Relative of Caenorhabditis elegans

The abundance, diversity, and genomic distribution of repetitive elements is highly variable among species. These patterns are thought to be driven in part by reproductive mode and the interaction of selection and recombination, and recombination rates typically vary by chromosomal position. In the nematode Caenorhabditis elegans, repetitive elements are enriched at chromosome arms and depleted on centers, and this mirrors the chromosomal distributions of other genomic features such as recombination rate. How conserved is this genomic landscape of repeats, and what evolutionary forces maintain it? To address this, we compared the genomic organization of repetitive elements across five Caenorhabditis species with chromosome-level assemblies. As previously reported, repeat content is enriched on chromosome arms in most Caenorhabditis species, and no obvious patterns of repeat content associated with reproductive mode were observed. However, the fig-associated C. inopinata has experienced repetitive element expansion and reveals no association of global repeat density with chromosome position. Patterns of repeat superfamily specific distributions reveal this global pattern is driven largely by a few repeat superfamilies that in C. inopinata have expanded in number and have weak associations with chromosome position. Additionally, 15% of predicted protein-coding genes in C. inopinata align to transposon-related proteins. When these are excluded, C. inopinata has no enrichment of genes in chromosome centers, in contrast to its close relatives who all have such clusters. Forward evolutionary simulations reveal that chromosomal heterogeneity in recombination rate alone can generate structured repetitive genomic landscapes when insertions are weakly deleterious, whereas chromosomal heterogeneity in the fitness effects of transposon insertion can promote such landscapes across a variety of evolutionary scenarios. Thus, patterns of gene density along chromosomes likely contribute to global repetitive landscapes in this group, although other historical or genomic factors are needed to explain the idiosyncrasy of genomic organization of various transposable element taxa within C. inopinata. Taken together, these results highlight the power of comparative genomics and evolutionary simulations in testing hypotheses regarding the causes of genome organization.

The differential expression patterns of paralogs in response to stresses indicate expression and sequence divergences

BACKGROUND

Theoretically, paralogous genes generated through whole genome duplications should share identical expression levels due to their identical sequences and chromatin environments. However, functional divergences and expression differences have arisen due to selective pressures throughout evolution. A comprehensive investigation of the expression patterns of paralogous gene pairs in response to various stresses and a study of correlations between the expression levels and sequence divergences of the paralogs are needed.

RESULTS

In this study, we analyzed the expression patterns of paralogous genes under different types of stress and investigated the correlations between the expression levels and sequence divergences of the paralogs. We analyzed the differential expression patterns of the paralogs under four different types of stress (drought, cold, infection, and herbivory) and classified them into three main types according to their expression patterns. We then further analyzed the differential expression patterns under various degrees of stress and constructed corresponding co-expression networks of differentially expressed paralogs and transcription factors. Finally, we investigated the correlations between the expression levels and sequence divergences of the paralogs and identified positive correlations between expression level and sequence divergence. With regard to sequence divergence, we identified correlations between selective pressures and phylogenetic relationships.

CONCLUSIONS

These results shed light on differential expression patterns of paralogs in response to environmental stresses and are helpful for understanding the relationships between expression levels and sequences divergences.

AthCNV: A Map of DNA Copy Number Variations in the Arabidopsis Genome

Copy number variations (CNVs) greatly contribute to intraspecies genetic polymorphism and phenotypic diversity. Recent analyses of sequencing data for >1000 Arabidopsis (Arabidopsis thaliana) accessions focused on small variations and did not include CNVs. Here, we performed genome-wide analysis and identified large indels (50 to 499 bp) and CNVs (500 bp and larger) in these accessions. The CNVs fully overlap with 18.3% of protein-coding genes, with enrichment for evolutionarily young genes and genes involved in stress and defense. By combining analysis of both genes and transposable elements (TEs) affected by CNVs, we revealed that the variation statuses of genes and TEs are tightly linked and jointly contribute to the unequal distribution of these elements in the genome. We also determined the gene copy numbers in a set of 1060 accessions and experimentally validated the accuracy of our predictions by multiplex ligation-dependent probe amplification assays. We then successfully used the CNVs as markers to analyze population structure and migration patterns. Finally, we examined the impact of gene dosage variation triggered by a CNV spanning the SEC10 gene on SEC10 expression at both the transcript and protein levels. The catalog of CNVs, CNV-overlapping genes, and their genotypes in a top model dicot will stimulate the exploration of the genetic basis of phenotypic variation.

Genome-wide variations analysis of sorghum cultivar Hongyingzi for brewing Moutai liquor

Background

Hongyingzi is a sorghum (Sorghum bicolor L. Moench) cultivar for brewing Moutai liquor. For an overall understanding of the whole genome of Hongyingzi, we performed whole-genome resequencing technology to reveal its comprehensive variations.

Results

Compared with the BTx623 reference genome, we uncovered 1,885,774 single nucleotide polymorphisms (SNPs), 309,381 small fragments insertions and deletions (Indels), 31,966 structural variations (SVs), and 217,273 copy number variations (CNVs). These alterations conferred 29,614 gene variations. It was also predicted that 35 gene variations were related to the multidrug and toxic efflux (MATE) transporter, chalcone synthase (CHS), ATPase isoform 10 (AHA10) transporter, dihydroflavonol-4-reductase (DFR), the laccase 15 (LAC15), flavonol 3′-hydroxylase (F3′H), flavanone 3-hydroxylase (F3H), O-methyltransferase (OMT), flavonoid 3′5′ hydroxylase (F3′5′H), UDP-glucose:sterol-glucosyltransferase (SGT), flavonol synthase (FLS), and chalcone isomerase (CHI) involved in the tannin synthesis.

Conclusions

These results would provide theoretical supports for the molecular markers developments and gene function studies related to the tannin synthesis, and the genetic improvement of liquor-making sorghum based on the genome editing technology.

Natural variation in the Arabidopsis AGO2 gene is associated with susceptibility to potato virus X

RNA silencing functions as an anti-viral defence in plants through the action of DICER-like (DCL) and ARGONAUTE (AGO) proteins. Despite the importance of this mechanism, little is known about the functional consequences of variation in genes encoding RNA silencing components. The AGO2 protein has been shown to be important for defense against multiple viruses, and we investigated how naturally occurring differences in AGO2 between and within species affects its antiviral activities. We find that the AGO2 protein from Arabidopsis thaliana, but not Nicotiana benthamiana, effectively limits potato virus X (PVX). Consistent with this, we find that the A. thaliana AGO2 gene shows a high incidence of polymorphisms between accessions, with evidence of selective pressure. Using functional analyses, we identify polymorphisms that specifically affect AGO2 antiviral activity, without interfering with other AGO2-associated functions such as anti-bacterial resistance or DNA methylation. Our results suggest that viruses adapt to overcome RNA silencing in their hosts. Furthermore, they indicate that plant-virus interactions have influenced natural variation in RNA-silencing components and that the latter may be a source of genetically encoded virus resistance.

Advances in genome-wide association studies of complex traits in rice

Genome-wide association studies (GWAS), genetic surveys of the whole genome to detect variants associated with a trait in natural populations, are a powerful approach for dissecting complex traits. This genetic mapping approach has been applied in rice over the last 10 years. During the last decade, GWAS was used to identify the loci underlying tens of rice traits, and several important genes were detected in GWAS and further confirmed in follow-up functional experiments. In this review, we present an overview of the whole process in a typical GWAS, including population design, genotyping, phenotyping and analysis methods. Recent advances in rice GWAS are also provided, including several examples of the functional characterization of candidate genes. The possible breakthroughs of rice GWAS in the next decade are discussed with regard to their application in breeding, the consideration of epistatic interactions and in-depth functional annotations of DNA elements and genetic variants throughout the rice genome.

The Earth BioGenome project: opportunities and challenges for plant genomics and conservation

Sequencing them all. That is the ambitious goal of the recently launched Earth BioGenome project (Proceedings of the National Academy of Sciences of the United States of America, 115, 4325-4333), which aims to produce reference genomes for all eukaryotic species within the next decade. In this perspective, we discuss the opportunities of this project with a plant focus, but highlight also potential limitations. This includes the question of how to best capture all plant diversity, as the green taxon is one of the most complex clades in the tree of life, with over 300 000 species. For this, we highlight four key points: (i) the unique biological insights that could be gained from studying plants, (ii) their apparent underrepresentation in sequencing efforts given the number of threatened species, (iii) the necessity of phylogenomic methods that are aware of differences in genome complexity and quality, and (iv) the accounting for within-species genetic diversity and the historical aspect of conservation genetics.

Chromosome-level assemblies of multiple Arabidopsis genomes reveal hotspots of rearrangements with altered evolutionary dynamics

Despite hundreds of sequenced Arabidopsis genomes, very little is known about the degree of genomic collinearity within single species, due to the low number of chromosome-level assemblies. Here, we report chromosome-level reference-quality assemblies of seven Arabidopsis thaliana accessions selected across its global range. Each genome reveals between 13–17 Mb rearranged, and 5–6 Mb non-reference sequences introducing copy-number changes in ~5000 genes, including ~1900 non-reference genes. Quantifying the collinearity between the genomes reveals ~350 euchromatic regions, where accession-specific tandem duplications destroy the collinearity between the genomes. These hotspots of rearrangements are characterized by reduced meiotic recombination in hybrids and genes implicated in biotic stress response. This suggests that hotspots of rearrangements undergo altered evolutionary dynamics, as compared to the rest of the genome, which are mostly based on the accumulation of new mutations and not on the recombination of existing variation, and thereby enable a quick response to the biotic stress.

Despite tremendous genomic resources in the Arabidopsis community, only a few whole genome de novo assemblies are available. Here, the authors report chromosome-level reference-quality assemblies of seven A. thaliana accessions and reveal hotspots of rearrangements with altered evolutionary dynamics.

Comparative population genomics reveals genetic divergence and selection in lotus, Nelumbo nucifera

Background

Lotus (Nelumbo nucifera) is an aquatic plant with important agronomic, horticulture, art and religion values. It was the basal eudicot species occupying a critical phylogenetic position in flowering plants. After the domestication for thousands of years, lotus has differentiated into three cultivated types -flower lotus, seed lotus and rhizome lotus. Although the phenotypic and genetic differentiations based on molecular markers have been reported, the variation on whole-genome level among the different lotus types is still ambiguous.

Results

In order to reveal the evolution and domestication characteristics of lotus, a total of 69 lotus accessions were selected, including 45 cultivated accessions, 22 wild sacred lotus accessions, and 2 wild American lotus accessions. With Illumina technology, the genomes of these lotus accessions were resequenced to > 13× raw data coverage. On the basis of these genomic data, 25 million single-nucleotide polymorphisms (SNPs) were identified in lotus. Population analysis showed that the rhizome and seed lotus were monophyletic and genetically homogeneous, whereas the flower lotus was biphyletic and genetically heterogeneous. Using population SNP data, we identified 1214 selected regions in seed lotus, 95 in rhizome lotus, and 37 in flower lotus. Some of the genes in these regions contributed to the essential domestication traits of lotus. The selected genes of seed lotus mainly affected lotus seed weight, size and nutritional quality. While the selected genes were responsible for insect resistance, antibacterial immunity and freezing and heat stress resistance in flower lotus, and improved the size of rhizome in rhizome lotus, respectively.

Conclusions

The genome differentiation and a set of domestication genes were identified from three types of cultivated lotus- flower lotus, seed lotus and rhizome lotus, respectively. Among cultivated lotus, flower lotus showed the greatest variation. The domestication genes may show agronomic importance via enhancing insect resistance, improving seed weight and size, or regulating lotus rhizome size. The domestication history of lotus enhances our knowledge of perennial aquatic crop evolution, and the obtained dataset provides a basis for future genomics-enabled breeding.