A membrane associated tandem kinase from wild emmer wheat confers broad-spectrum resistance to powdery mildew

Crop wild relatives offer natural variations of disease resistance for crop improvement. Here, we report the isolation of broad-spectrum powdery mildew resistance gene Pm36, originated from wild emmer wheat, that encodes a tandem kinase with a transmembrane domain (WTK7-TM) through the combination of map-based cloning, PacBio SMRT long-read genome sequencing, mutagenesis, and transformation. Mutagenesis assay reveals that the two kinase domains and the transmembrane domain of WTK7-TM are critical for the powdery mildew resistance function. Consistently, in vitro phosphorylation assay shows that two kinase domains are indispensable for the kinase activity of WTK7-TM. Haplotype analysis uncovers that Pm36 is an orphan gene only present in a few wild emmer wheat, indicating its single ancient origin and potential contribution to the current wheat gene pool. Overall, our findings not only provide a powdery mildew resistance gene with great potential in wheat breeding but also sheds light into the mechanism underlying broad-spectrum resistance.

Natural variation of STKc_GSK3 kinase TaSG-D1 contributes to heat stress tolerance in Indian dwarf wheat

Heat stress threatens global wheat (Triticum aestivum) production, causing dramatic yield losses worldwide. Identifying heat tolerance genes and comprehending molecular mechanisms are essential. Here, we identify a heat tolerance gene, TaSG-D1E286K, in Indian dwarf wheat (Triticum sphaerococcum), which encodes an STKc_GSK3 kinase. TaSG-D1E286K improves heat tolerance compared to TaSG-D1 by enhancing phosphorylation and stability of downstream target TaPIF4 under heat stress condition. Additionally, we reveal evolutionary footprints of TaPIF4 during wheat selective breeding in China, that is, InDels predominantly occur in the TaPIF4 promoter of Chinese modern wheat cultivars and result in decreased expression level of TaPIF4 in response to heat stress. These sequence variations with negative effect on heat tolerance are mainly introduced from European germplasm. Our study provides insight into heat stress response mechanisms and proposes a potential strategy to improve wheat heat tolerance in future.

A genomic toolkit for winged bean Psophocarpus tetragonolobus

A sustainable supply of plant protein is critical for future generations and needs to be achieved while reducing green house gas emissions from agriculture and increasing agricultural resilience in the face of climate volatility. Agricultural diversification with more nutrient-rich and stress tolerant crops could provide the solution. However, this is often hampered by the limited availability of genomic resources and the lack of understanding of the genetic structure of breeding germplasm and the inheritance of important traits. One such crop with potential is winged bean (Psophocarpus tetragonolobus), a high seed protein tropical legume which has been termed 'the soybean for the tropics'. Here, we present a chromosome level winged bean genome assembly, an investigation of the genetic diversity of 130 worldwide accessions, together with two linked genetic maps and a trait QTL analysis (and expression studies) for regions of the genome with desirable ideotype traits for breeding, namely architecture, protein content and phytonutrients.

Epiallelic variation of non-coding RNA genes and their phenotypic consequences

Epigenetic variations contribute greatly to the phenotypic plasticity and diversity. Current functional studies on epialleles have predominantly focused on protein-coding genes, leaving the epialleles of non-coding RNA (ncRNA) genes largely understudied. Here, we uncover abundant DNA methylation variations of ncRNA genes and their significant correlations with plant adaptation among 1001 natural Arabidopsis accessions. Through genome-wide association study (GWAS), we identify large numbers of methylation QTL (methylQTL) that are independent of known DNA methyltransferases and enriched in specific chromatin states. Proximal methylQTL closely located to ncRNA genes have a larger effect on DNA methylation than distal methylQTL. We ectopically tether a DNA methyltransferase MQ1v to miR157a by CRISPR-dCas9 and show de novo establishment of DNA methylation accompanied with decreased miR157a abundance and early flowering. These findings provide important insights into the genetic basis of epigenetic variations and highlight the contribution of epigenetic variations of ncRNA genes to plant phenotypes and diversity.

A mediator of OsbZIP46 deactivation and degradation negatively regulates seed dormancy in rice

Preharvest sprouting (PHS) is a deleterious phenotype that occurs frequently in rice-growing regions where the temperature and precipitation are high. It negatively affects yield, quality, and downstream grain processing. Seed dormancy is a trait related to PHS. Longer seed dormancy is preferred for rice production as it can prevent PHS. Here, we map QTLs associated with rice seed dormancy and clone Seed Dormancy 3.1 (SDR3.1) underlying one major QTL. SDR3.1 encodes a mediator of OsbZIP46 deactivation and degradation (MODD). We show that SDR3.1 negatively regulates seed dormancy by inhibiting the transcriptional activity of ABIs. In addition, we reveal two critical amino acids of SDR3.1 that are critical for the differences in seed dormancy between the Xian/indica and Geng/japonica cultivars. Further, SDR3.1 has been artificially selected during rice domestication. We propose a two-line model for the process of rice seed dormancy domestication from wild rice to modern cultivars. We believe the candidate gene and germplasm studied in this study would be beneficial for the genetic improvement of rice seed dormancy.

Comparative transcriptomics of two Salvia subg. Perovskia species contribute towards molecular background of abietane-type diterpenoid biosynthesis

Tanshinones, are a group of diterpenoid red pigments present in Danshen - an important herbal drug of Traditional Chinese Medicine which is a dried root of Salvia miltiorrhiza Bunge. Some of the tanshinones are sought after as pharmacologically active natural products. To date, the biosynthetic pathway of tanshinones has been only partially elucidated. These compounds are also present in some of the other Salvia species, i.a. from subgenus Perovskia, such as S. abrotanoides (Kar.) Sytsma and S. yangii B.T. Drew. Despite of the close genetic relationship between these species, significant qualitative differences in their diterpenoid profile have been discovered. In this work, we have used the Liquid Chromatography-Mass Spectrometry analysis to follow the content of diterpenoids during the vegetation season, which confirmed our previous observations of a diverse diterpenoid profile. As metabolic differences are reflected in different transcript profile of a species or tissues, we used metabolomics-guided transcriptomic approach to select candidate genes, which expression possibly led to observed chemical differences. Using an RNA-sequencing technology we have sequenced and de novo assembled transcriptomes of leaves and roots of S. abrotanoides and S. yangii. As a result, 134,443 transcripts were annotated by UniProt and 56,693 of them were assigned as Viridiplantae. In order to seek for differences, the differential expression analysis was performed, which revealed that 463, 362, 922 and 835 genes indicated changes in expression in four comparisons. GO enrichment analysis and KEGG functional analysis of selected DEGs were performed. The homology and expression of two gene families, associated with downstream steps of tanshinone and carnosic acid biosynthesis were studied, namely: cytochromes P-450 and 2-oxoglutarate-dependend dioxygenases. Additionally, BLAST analysis revealed existence of 39 different transcripts related to abietane diterpenoid biosynthesis in transcriptomes of S. abrotanoides and S. yangii. We have used quantitative real-time RT-PCR analysis of selected candidate genes, to follow their expression levels over the vegetative season. A hypothesis of an existence of a multifunctional CYP76AH89 in transcriptomes of S. abrotanoides and S. yangii is discussed and potential roles of other CYP450 homologs are speculated. By using the comparative transcriptomic approach, we have generated a dataset of candidate genes which provides a valuable resource for further elucidation of tanshinone biosynthesis. In a long run, our investigation may lead to optimization of diterpenoid profile in S. abrotanoides and S. yangii, which may become an alternative source of tanshinones for further research on their bioactivity and pharmacological therapy.

Comparative epigenetic and genetic spatial structure in Mediterranean mountain plants: a multispecies study

Changes in epigenetic states can allow individuals to cope with environmental changes. If such changes are heritable, this may lead to epigenetic adaptation. Thus, it is likely that in sessile organisms such as plants, part of the spatial epigenetic variation found across individuals will reflect the environmental heterogeneity within populations. The departure of the spatial epigenetic structure from the baseline genetic variation can help in understanding the value of epigenetic regulation in species with different breadth of optimal environmental requirements. Here, we hypothesise that in plants with narrow environmental requirements, epigenetic variability should be less structured in space given the lower variability in suitable environmental conditions. We performed a multispecies study that considered seven pairs of congeneric plant species, each encompassing a narrow endemic with habitat specialisation and a widespread species. In three populations per species we used AFLP and methylation-sensitive AFLP markers to characterise the spatial genetic and epigenetic structures. Narrow endemics showed a significantly lower epigenetic than genetic differentiation between populations. Within populations, epigenetic variation was less spatially structured than genetic variation, mainly in narrow endemics. In these species, structural equation models revealed that such pattern was associated to a lack of correlation between epigenetic and genetic information. Altogether, these results show a greater decoupling of the spatial epigenetic variation from the baseline spatial genetic pattern in endemic species. These findings highlight the value of studying genetic and epigenetic spatial variation to better understand habitat specialisation in plants.

Root system architecture associated zinc variability in wheat (Triticum aestivum L.)

Root system architecture (RSA) plays a fundamental role in nutrient uptake, including zinc (Zn). Wheat grains are inheritably low in Zn. As Zn is an essential nutrient for plants, improving its uptake will not only improve their growth and yield but also the nutritional quality of staple grains. A rhizobox study followed by a pot study was conducted to evaluate Zn variability with respect to RSA and its impact on grain Zn concentration. The grain Zn content of one hundred wheat varieties was determined and grown in rhizoboxes with differential Zn (no Zn and 0.05 mg L-1 ZnSO4). Seedlings were harvested 12 days after sowing, and root images were taken and analyzed by SmartRoot software. Using principal component analysis, twelve varieties were screened out based on vigorous and weaker RSA with high and low grain Zn content. The screened varieties were grown in pots with (11 mg ZnSO4 kg-1 soil) and without Zn application to the soil. Zinc translocation, localization, and agronomic parameters were recorded after harvesting at maturity. In the rhizobox experiment, 4% and 8% varieties showed higher grain Zn content with vigorous and weaker RSA, respectively, while 45% and 43% varieties had lower grain Zn content with vigorous and weaker RSA. However, the pot experiment revealed that varieties with vigorous root system led to higher grain yield, though the grain Zn concentration were variable, while all varieties with weaker root system had lower yield as well as grain Zn concentration. Zincol-16 revealed the highest Zn concentration (28.07 mg kg-1) and grain weight (47.9 g). Comparatively higher level of Zn was localized in the aleurone layer than in the embryonic region and endosperm. It is concluded that genetic variability exists among wheat varieties for RSA and grain Zn content, with a significant correlation. Therefore, RSA attributes are promising targets for the Zn biofortification breeding program. However, Zn localization in endosperm needs to be further investigated to achieve the goal of reducing Zn malnutrition.

Changes in microbial community structure and yield responses with the use of nano-fertilizers of nitrogen and zinc in wheat-maize system

The growing popularity of nano-fertilization around the world for enhancing yield and nutrient use efficiency has been realized, however its influence on soil microbial structure is not fully understood. The purpose of carrying out this study was to assess the combined effect of nano and conventional fertilizers on the soil biological indicators and crop yield in a wheat-maize system. The results indicate that the at par grain yield of wheat and maize was obtained with application of 75% of recommended nitrogen (N) with full dose of phosphorus (P) and potassium (K) through conventional fertilizers along with nano-N (nano-urea) or nano-N plus nano-Zn sprays and N100PK i.e. business as usual (recommended dose of fertilizer). Important soil microbial property like microbial biomass carbon was found statistically similar with nano fertilizer-based management (N75PK + nano-N, and N75PK + nano-N + nano-Zn) and conventional management (N100PK), during both wheat and maize seasons. The experimental data indicated that the application of foliar spray of nano-fertilizers along with 75% N as basal is a sustainable nutrient management approach with respect to growth, yield and rhizosphere biological activity. Furthermore, two foliar sprays of nano-N or nano-N + nano-Zn curtailed N requirement by 25%, furthermore enhanced soil microbial diversity and the microbial community structure. The specific microbial groups, including Actinobacteria, Bacteroidia, and Proteobacteria, were present in abundance and were positively correlated with wheat and maize yield and soil microbial biomass carbon. Thus, one of the best nutrient management approaches for sustaining productivity and maintaining sound microbial diversity in wheat-maize rotation is the combined use of nano-fertilizers and conventional fertilizers.

Proline-mediated redox regulation in wheat for mitigating nickel-induced stress and soil decontamination

Nickel (Ni) is known as a plant micronutrient and serves as a component of many significant enzymes, however, it can be extremely toxic to plants when present in excess concentration. Scientists are looking for natural compounds that can influence the development processes of plants. Therefore, it was decided to use proline as a protective agent against Ni toxicity. Proline (Pro) is a popularly known osmoprotectant to regulate the biomass and developmental processes of plants under a variety of environmental stresses, but its role in the modulation of Ni-induced toxicity in wheat is very little explored. This investigation indicated the role of exogenously applied proline (10 mM) on two wheat varieties (V1 = Punjab-11, V2 = Ghazi-11) exposed to Ni (100 mg/kg) stress. Proline mediated a positive rejoinder on morphological, photosynthetic indices, antioxidant enzymes, oxidative stress markers, ion uptake were analyzed with and without Ni stress. Proline alone and in combination with Ni improved the growth, photosynthetic performance, and antioxidant capacity of wheat plants. However, Ni application alone exhibited strong oxidative damage through increased H2O2 (V1 = 28.96, V2 = 55.20) accumulation, lipid peroxidation (V1 = 26.09, V2 = 38.26%), and reduced translocation of macronutrients from root to shoot. Application of Pro to Ni-stressed wheat plants enhanced actions of catalase (CAT), peroxidase (POD), superoxide dismutase (SOD), and total soluble protein (TSP) contents by 45.70, 44.06, 43.40, and 25.11% in V1, and 39.32, 46.46, 42.22, 55.29% in V2, compared to control plants. The upregulation of antioxidant enzymes, proline accumulation, and uptake of essential mineral ions has maintained the equilibrium of Ni in both wheat cultivars, indicating Ni detoxification. This trial insight into an awareness that foliar application of proline can be utilized as a potent biochemical method in mitigating Ni-induced stress and might serve as a strong remedial technique for the decontamination of polluted soil particularly with metals.

Microbiome homeostasis on rice leaves is regulated by a precursor molecule of lignin biosynthesis

In terrestrial ecosystems, plant leaves provide the largest biological habitat for highly diverse microbial communities, known as the phyllosphere microbiota. However, the underlying mechanisms of host-driven assembly of these ubiquitous communities remain largely elusive. Here, we conduct a large-scale and in-depth assessment of the rice phyllosphere microbiome aimed at identifying specific host-microbe links. A genome-wide association study reveals a strong association between the plant genotype and members of four bacterial orders, Pseudomonadales, Burkholderiales, Enterobacterales and Xanthomonadales. Some of the associations are specific to a distinct host genomic locus, pathway or even gene. The compound 4-hydroxycinnamic acid (4-HCA) is identified as the main driver for enrichment of bacteria belonging to Pseudomonadales. 4-HCA can be synthesized by the host plant's OsPAL02 from the phenylpropanoid biosynthesis pathway. A knockout mutant of OsPAL02 results in reduced Pseudomonadales abundance, dysbiosis of the phyllosphere microbiota and consequently higher susceptibility of rice plants to disease. Our study provides a direct link between a specific plant metabolite and rice phyllosphere homeostasis opening possibilities for new breeding strategies.

A pilot study on non-invasive in situ detection of phytochemicals and plant endogenous status using fiber optic infrared spectroscopy

Traditional methods for assessing plant health often lack the necessary attributes for continuous and non-destructive monitoring. In this pilot study, we present a novel technique utilizing a customized fiber optic probe based on attenuated total reflection Fourier transform infrared spectroscopy (ATR-FTIR) with a contact force control unit for non-invasive and continuous plant health monitoring. We also developed a normalized difference mid-infrared reflectance index through statistical analysis of spectral features, enabling differentiation of drought and age conditions in plants. Our research aims to characterize phytochemicals and plant endogenous status optically, addressing the need for improved analytical measurement methods for in situ plant health assessment. The probe configuration was optimized with a triple-loop tip and a 3 N contact force, allowing sensitive measurements while minimizing leaf damage. By combining polycrystalline and chalcogenide fiber probes, a comprehensive wavenumber range analysis (4000-900 cm-1) was achieved. Results revealed significant variations in phytochemical composition among plant species, for example, red spinach with the highest polyphenolic content and green kale with the highest lignin content. Petioles displayed higher lignin and cellulose absorbance values compared to veins. The technique effectively monitored drought stress on potted green bok choy plants in situ, facilitating the quantification of changes in water content, antioxidant activity, lignin, and cellulose levels. This research represents the first demonstration of the potential of fiber optic ATR-FTIR probes for non-invasive and rapid plant health measurements, providing insights into plant health and advancements in quantitative monitoring for indoor farming practices, bioanalytical chemistry, and environmental sciences.

Control of compound leaf patterning by MULTI-PINNATE LEAF1 (MPL1) in chickpea

Plant lateral organs are often elaborated through repetitive formation of developmental units, which progress robustly in predetermined patterns along their axes. Leaflets in compound leaves provide an example of such units that are generated sequentially along the longitudinal axis, in species-specific patterns. In this context, we explored the molecular mechanisms underlying an acropetal mode of leaflet initiation in chickpea pinnate compound leaf patterning. By analyzing naturally occurring mutants multi-pinnate leaf1 (mpl1) that develop higher-ordered pinnate leaves with more than forty leaflets, we show that MPL1 encoding a C2H2-zinc finger protein sculpts a morphogenetic gradient along the proximodistal axis of the early leaf primordium, thereby conferring the acropetal leaflet formation. This is achieved by defining the spatiotemporal expression pattern of CaLEAFY, a key regulator of leaflet initiation, and also perhaps by modulating the auxin signaling pathway. Our work provides novel molecular insights into the sequential progression of leaflet formation.

Copy number variation of the restorer Rf4 underlies human selection of three-line hybrid rice breeding

Cytoplasmic male sterility (CMS) lines are important for breeding hybrid crops, and utilization of CMS lines requires strong fertility restorer (Rf) genes. Rf4, a major Rf for Wild-Abortive CMS (CMS-WA), has been cloned in rice. However, the Rf4 evolution and formation of CMS-WA/Rf system remain elusive. Here, we show that the Rf4 locus emerges earlier than the CMS-WA gene WA352 in wild rice, and 69 haplotypes of the Rf4 locus are generated in the Oryza genus through the copy number and sequence variations. Eight of these haplotypes of the Rf4 locus are enriched in modern rice cultivars during natural and human selections, whereas non-functional rf4i is preferentially selected for breeding current CMS-WA lines. We further verify that varieties carrying two-copy Rf4 haplotype have stronger fertility restoration ability and are widely used in three-line hybrid rice breeding. Our findings increase our understanding of CMS/Rf systems and will likely benefit crop breeding.

In vitro and in silico prediction of antibacterial interaction between essential oils via graph embedding approach

Essential oils contain a variety of volatile metabolites, and are expected to be utilized in wide fields such as antimicrobials, insect repellents and herbicides. However, it is difficult to foresee the effect of oil combinations because hundreds of compounds can be involved in synergistic and antagonistic interactions. In this research, it was developed and evaluated a machine learning method to classify types of (synergistic/antagonistic/no) antibacterial interaction between essential oils. Graph embedding was employed to capture structural features of the interaction network from literature data, and was found to improve in silico predicting performances to classify synergistic interactions. Furthermore, in vitro antibacterial assay against a standard strain of Staphylococcus aureus revealed that four essential oil pairs (Origanum compactum-Trachyspermum ammi, Cymbopogon citratus-Thujopsis dolabrata, Cinnamomum verum-Cymbopogon citratus and Trachyspermum ammi-Zingiber officinale) exhibited synergistic interaction as predicted. These results indicate that graph embedding approach can efficiently find synergistic interactions between antibacterial essential oils.

Spatial genetic structure and seed quality of a southernmost Abies nephrolepis population

Abies nephrolepis (Trautv. ex Maxim.) Maxim. has its southernmost populations in South Korea and they are expected to decline under climate change. To establish a strategic conservation plan, this study aimed to investigate the spatial genetic structure and seed characteristics of A. nephrolepis. We used nine microsatellite markers on 165 individuals of A. nephrolepis and sampled seeds in a southernmost population at Mt. Hambaeksan, South Korea. We observed a high level of heterozygosity, and a simulation study found that sampling 20 individuals was enough to secure sufficient genetic diversity on average. Spatial autocorrelation analysis revealed that individuals had a positive genetic relationship until 30 m. Bayesian clustering models, STRUCTURE and GENELAND, failed to achieve a consensus in the optimal number of population (K), estimating K = 1 and K = 2, respectively. Principal coordinate analysis supported the absence of genetic substructure within the study population. There was a large variance in seed production among mother trees. On average, seeds of A. nephrolepis from Mt. Hambaeksan had a purity of 70.4% and a germination percentage of 32.2%. We found that seed weight was the most effective indicator of seed quality. Mother trees at higher altitudes had poorer purity which is threatening to A. nephrolepis considering the upslope retreat of subalpine species under climate change. Our results provide insights into the interactions among spatial processes, genetic structure, and seed quality within a population of A. nephrolepis.

PH13 improves soybean shade traits and enhances yield for high-density planting at high latitudes

Shading in combination with extended photoperiods can cause exaggerated stem elongation (ESE) in soybean, leading to lodging and reduced yields when planted at high-density in high-latitude regions. However, the genetic basis of plant height in adaptation to these regions remains unclear. Here, through a genome-wide association study, we identify a plant height regulating gene on chromosome 13 (PH13) encoding a WD40 protein with three main haplotypes in natural populations. We find that an insertion of a Ty1/Copia-like retrotransposon in the haplotype 3 leads to a truncated PH13H3 with reduced interaction with GmCOP1s, resulting in accumulation of STF1/2, and reduced plant height. In addition, PH13H3 allele has been strongly selected for genetic improvement at high latitudes. Deletion of both PH13 and its paralogue PHP can prevent shade-induced ESE and allow high-density planting. This study provides insights into the mechanism of shade-resistance and offers potential solutions for breeding high-yielding soybean cultivar for high-latitude regions.

Diversity of gene expression responses to light quality in barley

Light quality influence on barley development is poorly understood. We exposed three barley genotypes with either sensitive or insensitive response to two light sources producing different light spectra, fluorescent bulbs, and metal halide lamps, keeping constant light intensity, duration, and temperature. Through RNA-seq, we identified the main genes and pathways involved in the genotypic responses. A first analysis identified genotypic differences in gene expression of development-related genes, including photoreceptors and flowering time genes. Genes from the vernalization pathway of light quality-sensitive genotypes were affected by fluorescent light. In particular, vernalization-related repressors reacted differently: HvVRN2 did not experience relevant changes, whereas HvOS2 expression increased under fluorescent light. To identify the genes primarily related to light quality responses, and avoid the confounding effect of plant developmental stage, genes influenced by development were masked in a second analysis. Quantitative expression levels of PPD-H1, which influenced HvVRN1 and HvFT1, explained genotypic differences in development. Upstream mechanisms (light signaling and circadian clock) were also altered, but no specific genes linking photoreceptors and the photoperiod pathway were identified. The variety of light-quality sensitivities reveals the presence of possible mechanisms of adaptation of winter and facultative barley to latitudinal variation in light quality, which deserves further research.

The pan-genome and local adaptation of Arabidopsis thaliana

Arabidopsis thaliana serves as a model species for investigating various aspects of plant biology. However, the contribution of genomic structural variations (SVs) and their associate genes to the local adaptation of this widely distribute species remains unclear. Here, we de novo assemble chromosome-level genomes of 32 A. thaliana ecotypes and determine that variable genes expand the gene pool in different ecotypes and thus assist local adaptation. We develop a graph-based pan-genome and identify 61,332 SVs that overlap with 18,883 genes, some of which are highly involved in ecological adaptation of this species. For instance, we observe a specific 332 bp insertion in the promoter region of the HPCA1 gene in the Tibet-0 ecotype that enhances gene expression, thereby promotes adaptation to alpine environments. These findings augment our understanding of the molecular mechanisms underlying the local adaptation of A. thaliana across diverse habitats.

Evolution of phenotypic disparity in the plant kingdom

The plant kingdom exhibits diverse bodyplans, from single-celled algae to complex multicellular land plants, but it is unclear how this phenotypic disparity was achieved. Here we show that the living divisions comprise discrete clusters within morphospace, separated largely by reproductive innovations, the extinction of evolutionary intermediates and lineage-specific evolution. Phenotypic complexity correlates not with disparity but with ploidy history, reflecting the role of genome duplication in plant macroevolution. Overall, the plant kingdom exhibits a pattern of episodically increasing disparity throughout its evolutionary history that mirrors the evolutionary floras and reflects ecological expansion facilitated by reproductive innovations. This pattern also parallels that seen in the animal and fungal kingdoms, suggesting a general pattern for the evolution of multicellular bodyplans.

A de novo evolved gene contributes to rice grain shape difference between indica and japonica

The role of de novo evolved genes from non-coding sequences in regulating morphological differentiation between species/subspecies remains largely unknown. Here, we show that a rice de novo gene GSE9 contributes to grain shape difference between indica/xian and japonica/geng varieties. GSE9 evolves from a previous non-coding region of wild rice Oryza rufipogon through the acquisition of start codon. This gene is inherited by most japonica varieties, while the original sequence (absence of start codon, gse9) is present in majority of indica varieties. Knockout of GSE9 in japonica varieties leads to slender grains, whereas introgression to indica background results in round grains. Population evolutionary analyses reveal that gse9 and GSE9 are derived from wild rice Or-I and Or-III groups, respectively. Our findings uncover that the de novo GSE9 gene contributes to the genetic and morphological divergence between indica and japonica subspecies, and provide a target for precise manipulation of rice grain shape.

The phased telomere-to-telomere reference genome of Musa acuminata, a main contributor to banana cultivars

Musa acuminata is a main wild contributor to banana cultivars. Here, we reported a haplotype-resolved and telomere-to-telomere reference genome of M. acuminata by incorporating PacBio HiFi reads, Nanopore ultra-long reads, and Hi-C data. The genome size of the two haploid assemblies was estimated to be 469.83 Mb and 470.21 Mb, respectively. Multiple assessments confirmed the contiguity (contig N50: 16.53 Mb and 18.58 Mb; LAI: 20.18 and 19.48), completeness (BUSCOs: 98.57% and 98.57%), and correctness (QV: 45.97 and 46.12) of the genome. The repetitive sequences accounted for about half of the genome size. In total, 40,889 and 38,269 protein-coding genes were annotated in the two haploid assemblies, respectively, of which 9.56% and 3.37% were newly predicted. Genome comparison identified a large reciprocal translocation involving 3 Mb and 10 Mb from chromosomes 01 and 04 within M. acuminata. This reference genome of M. acuminata provides a valuable resource for further understanding of subgenome evolution of Musa species, and precise genetic improvement of banana.

A pathogen-induced putative NAC transcription factor mediates leaf rust resistance in barley

Leaf rust, caused by Puccinia hordei, is one of the most widespread and damaging foliar diseases affecting barley. The barley leaf rust resistance locus Rph7 has been shown to have unusually high sequence and haplotype divergence. In this study, we isolate the Rph7 gene using a fine mapping and RNA-Seq approach that is confirmed by mutational analysis and transgenic complementation. Rph7 is a pathogen-induced, non-canonical resistance gene encoding a protein that is distinct from other known plant disease resistance proteins in the Triticeae. Structural analysis using an AlphaFold2 protein model suggests that Rph7 encodes a putative NAC transcription factor with a zinc-finger BED domain with structural similarity to the N-terminal DNA-binding domain of the NAC transcription factor (ANAC019) from Arabidopsis. A global gene expression analysis suggests Rph7 mediates the activation and strength of the basal defence response. The isolation of Rph7 highlights the diversification of resistance mechanisms available for engineering disease control in crops.

QTG-Miner aids rapid dissection of the genetic base of tassel branch number in maize

Genetic dissection of agronomic traits is important for crop improvement and global food security. Phenotypic variation of tassel branch number (TBN), a major breeding target, is controlled by many quantitative trait loci (QTLs). The lack of large-scale QTL cloning methodology constrains the systematic dissection of TBN, which hinders modern maize breeding. Here, we devise QTG-Miner, a multi-omics data-based technique for large-scale and rapid cloning of quantitative trait genes (QTGs) in maize. Using QTG-Miner, we clone and verify seven genes underlying seven TBN QTLs. Compared to conventional methods, QTG-Miner performs well for both major- and minor-effect TBN QTLs. Selection analysis indicates that a substantial number of genes and network modules have been subjected to selection during maize improvement. Selection signatures are significantly enriched in multiple biological pathways between female heterotic groups and male heterotic groups. In summary, QTG-Miner provides a large-scale approach for rapid cloning of QTGs in crops and dissects the genetic base of TBN for further maize breeding.

Divergent responses of native and invasive macroalgae to submarine groundwater discharge

Marine macroalgae are important indicators of healthy nearshore groundwater dependent ecosystems (GDEs), which are emergent global conservation priorities. Submarine groundwater discharge (SGD) supports abundant native algal communities in GDEs via elevated but naturally derived nutrients. GDEs are threatened by anthropogenic nutrient inputs that pollute SGD above ambient levels, favoring invasive algae. Accordingly, this case study draws on the GDE conditions of Kona, Hawai'i where we evaluated daily photosynthetic production and growth for two macroalgae; a culturally valued native (Ulva lactuca) and an invasive (Hypnea musciformis). Manipulative experiments-devised to address future land-use, climate change, and water-use scenarios for Kona-tested algal responses under a natural range of SGD nutrient and salinity levels. Our analyses demonstrate that photosynthesis and growth in U. lactuca are optimal in low-salinity, high-nutrient waters, whereas productivity for H. musciformis appears limited to higher salinities despite elevated nutrient subsidies. These findings suggest that reductions in SGD via climate change decreases in rainfall or increased water-use from the aquifer may relax physiological constraints on H. musciformis. Collectively, this study reveals divergent physiologies of a native and an invasive macroalga to SGD and highlights the importance of maintaining SGD quantity and quality to protect nearshore GDEs.

Effect of cutting management on the forage production and quality of tepary bean (Phaseolus acutifolius A. Gray)

Tepary bean (Phaseolus acutifolius A. Gray) is an underutilized drought tolerant annual legume, originating from the Sonoran Desert, that may be a beneficial forage/hay for beef cattle in the Southern Great Plains of the US (SGP). The SGP has erratic rainfall and periods of intermittent drought exacerbated by high summer temperatures. In 2020 and 2021, a split-plot design was used to evaluate 13 genotypes of tepary bean and a forage soybean (control) at El Reno, OK, USA to compare production of plant biomass and forage nutritive value parameters under seven harvest regimes. Genotypes were used as the main plot and cutting management as the sub-plot. Biomass production of all tepary bean genotypes equaled that of soybean (p > 0.05), while several genotypes had superior forage nutritive value traits (p ≤ 0.05). Overall, a 15-cm cutting height and 30-day harvest interval produced the best overall product (average dry biomass of 5.8 Mg ha-1 with average relative feed values (RFV) of 165). Although all harvest regimes reduced total seasonal biomass, forage nutritive value increased. However, the tradeoff between forage production and nutritive value may be unacceptable to most producers. Further agronomic and breeding research is needed to encourage producers to grow tepary bean as a forage/hay in the SGP.

Allelic variation of TaWD40-4B.1 contributes to drought tolerance by modulating catalase activity in wheat

Drought drastically restricts wheat production, so to dissect allelic variations of drought tolerant genes without imposing trade-offs between tolerance and yield is essential to cope with the circumstance. Here, we identify a drought tolerant WD40 protein encoding gene TaWD40-4B.1 of wheat via the genome-wide association study. The full-length allele TaWD40-4B.1^C but not the truncated allele TaWD40-4B.1^T possessing a nonsense nucleotide variation enhances drought tolerance and grain yield of wheat under drought. TaWD40-4B.1^C interacts with canonical catalases, promotes their oligomerization and activities, and reduces H₂O₂ levels under drought. The knock-down of catalase genes erases the role of TaWD40-4B.1^C in drought tolerance. TaWD40-4B.1^C proportion in wheat accessions is negatively correlative with the annual rainfall, suggesting this allele may be selected during wheat breeding. The introgression of TaWD40-4B.1^C enhances drought tolerance of the cultivar harboring TaWD40-4B.1^T. Therefore, TaWD40-4B.1^C could be useful for molecular breeding of drought tolerant wheat.

The giant diploid faba genome unlocks variation in a global protein crop

Increasing the proportion of locally produced plant protein in currently meat-rich diets could substantially reduce greenhouse gas emissions and loss of biodiversity1. However, plant protein production is hampered by the lack of a cool-season legume equivalent to soybean in agronomic value2. Faba bean (Vicia faba L.) has a high yield potential and is well suited for cultivation in temperate regions, but genomic resources are scarce. Here, we report a high-quality chromosome-scale assembly of the faba bean genome and show that it has expanded to a massive 13 Gb in size through an imbalance between the rates of amplification and elimination of retrotransposons and satellite repeats. Genes and recombination events are evenly dispersed across chromosomes and the gene space is remarkably compact considering the genome size, although with substantial copy number variation driven by tandem duplication. Demonstrating practical application of the genome sequence, we develop a targeted genotyping assay and use high-resolution genome-wide association analysis to dissect the genetic basis of seed size and hilum colour. The resources presented constitute a genomics-based breeding platform for faba bean, enabling breeders and geneticists to accelerate the improvement of sustainable protein production across the Mediterranean, subtropical and northern temperate agroecological zones.

Common evolutionary trajectory of short life-cycle in Brassicaceae ruderal weeds

Weed species are detrimental to crop yield. An understanding of how weeds originate and adapt to field environments is needed for successful crop management and reduction of herbicide use. Although early flowering is one of the weed trait syndromes that enable ruderal weeds to overcome frequent disturbances, the underlying genetic basis is poorly understood. Here, we establish Cardamine occulta as a model to study weed ruderality. By genome assembly and QTL mapping, we identify impairment of the vernalization response regulator gene FLC and a subsequent dominant mutation in the blue-light receptor gene CRY2 as genetic drivers for the establishment of short life cycle in ruderal weeds. Population genomics study further suggests that the mutations in these two genes enable individuals to overcome human disturbances through early deposition of seeds into the soil seed bank and quickly dominate local populations, thereby facilitating their spread in East China. Notably, functionally equivalent dominant mutations in CRY2 are shared by another weed species, Rorippa palustris, suggesting a common evolutionary trajectory of early flowering in ruderal weeds in Brassicaceae.

Graph-based pangenomics maximizes genotyping density and reveals structural impacts on fungal resistance in melon

The genomic sequences segregating in experimental populations are often highly divergent from the community reference and from one another. Such divergence is problematic under various short-read-based genotyping strategies. In addition, large structural differences are often invisible despite being strong candidates for causal variation. These issues are exacerbated in specialty crop breeding programs with fewer, lower-quality sequence resources. Here, we examine the benefits of complete genomic information, based on long-read assemblies, in a biparental mapping experiment segregating at numerous disease resistance loci in the non-model crop, melon (Cucumis melo). We find that a graph-based approach, which uses both parental genomes, results in 19% more variants callable across the population and raw allele calls with a 2 to 3-fold error-rate reduction, even relative to single reference approaches using a parent genome. We show that structural variation has played a substantial role in shaping two Fusarium wilt resistance loci with known causal genes. We also report on the genetics of powdery mildew resistance, where copy number variation and local recombination suppression are directly interpretable via parental genome alignments. Benefits observed, even in this low-resolution biparental experiment, will inevitably be amplified in more complex populations.

Genomic basis of the giga-chromosomes and giga-genome of tree peony Paeonia ostii

Tree peony (Paeonia ostii) is an economically important ornamental plant native to China. It is also notable for its seed oil, which is abundant in unsaturated fatty acids such as α-linolenic acid (ALA). Here, we report chromosome-level genome assembly (12.28 Gb) of P. ostii. In contrast to monocots with giant genomes, tree peony does not appear to have undergone lineage-specific whole-genome duplication. Instead, explosive LTR expansion in the intergenic regions within a short period (~ two million years) may have contributed to the formation of its giga-genome. In addition, expansion of five types of histone encoding genes may have helped maintain the giga-chromosomes. Further, we conduct genome-wide association studies (GWAS) on 448 accessions and show expansion and high expression of several genes in the key nodes of fatty acid biosynthetic pathway, including SAD, FAD2 and FAD3, may function in high level of ALAs synthesis in tree peony seeds. Moreover, by comparing with cultivated tree peony (P. suffruticosa), we show that ectopic expression of class A gene AP1 and reduced expression of class C gene AG may contribute to the formation of petaloid stamens. Genomic resources reported in this study will be valuable for studying chromosome/genome evolution and tree peony breeding.

Chemical composition of heartwood and sapwood of Tectona grandis characterized by CG/MS-PY

Teak wood has chemical compounds that can be used for pharmaceutical and textile industries, in addition, this compounds are related to resistance to biodeterioration, color and modification processes. Heartwood and sapwood of T. grandis (teak), 15 years-old, were characterized by Py-CG/MS analysis and syringyl (S)/guaiacyl (G) ratio was evaluated. Heartwood and sapwood were pyrolyzed at 550 °C and 62 and 51 compounds were identified from them, respectively. The acetic acid (10%) and levoglucosan (26.65%) were the most abundant compound in the sapwood and heartwood, respectively. The high acetic acid content enhances the use of teak wood to production of artificial essences for perfumery, paints, dyes. While levoglucosan can be used in the manufacture of epoxy resins, antiparasitic and insecticides. The organic compounds identified include 2-methylanthraquinone as one of the main component responsible for the resistance of the teak wood to biological factors (fungi and termites). The syringyl (S)/guaiacyl (G) ratio of heartwood and sapwood was 0.51 and 0.50, respectively.

Natural variation of Dt2 determines branching in soybean

Shoot branching is fundamentally important in determining soybean yield. Here, through genome-wide association study, we identify one predominant association locus on chromosome 18 that confers soybean branch number in the natural population. Further analyses determine that Dt2 is the corresponding gene and the natural variations in Dt2 result in significant differential transcriptional levels between the two major haplotypes. Functional characterization reveals that Dt2 interacts with GmAgl22 and GmSoc1a to physically bind to the promoters of GmAp1a and GmAp1d and to activate their transcription. Population genetic investigation show that the genetic differentiation of Dt2 display significant geographic structure. Our study provides a predominant gene for soybean branch number and may facilitate the breeding of high-yield soybean varieties.

Novel metallomic profiling and non-carcinogenic risk assessment of botanical ingredients for use in herbal, phytopharmaceutical and dietary products using HR-ICP-SFMS

Knowledge of element concentrations in botanical extracts is relevant to assure consumer protection given the increased interest in plant-based ingredients. This study demonstrates successful multi-element investigations in order to address the lack of comprehensive profiling data for botanical extracts, while reporting for the first time the metallomic profile(s) of arnica, bush vetch, sweet cicely, yellow rattle, bogbean, rock-tea and tufted catchfly. Key element compositions were quantified using a validated HR-ICP-SFMS method (µg kg-1) and were found highly variable between the different plants: Lithium (18-3964); Beryllium (3-121); Molybdenum (75-4505); Cadmium (5-325); Tin (6-165); Barium (747-4646); Platinum (2-33); Mercury (5-30); Thallium (3-91); Lead (12-4248); Bismuth (2-30); Titanium (131-5827); Vanadium (15-1758); Chromium (100-4534); Cobalt (21-652); Nickel (230-6060) and Copper (1910-6340). Compendial permissible limits were not exceeded. Overall, no evidence of a health risk to consumers could be determined from consumption of the investigated plants at reasonable intake rates. Mathematical risk modelling (EDI, CDI, HQ, HI) estimated levels above safe oral thresholds only for Cd (16%) and Pb (8%) from higher intakes of the respective plant-derived material. Following high consumption of certain plants, 42% of the samples were categorised as potentially unsafe due to cumulative exposure to Cu, Cd, Hg and Pb. PCA suggested a potential influence of post-harvest processing on Cr, Ti and V levels in commercially-acquired plant material compared to wild-collected and farm-grown plants. Moreover, a strong correlation was observed between Pb-Bi, Be-V, Bi-Sn, and Tl-Mo occurrence. This study may support future research by providing both robust methodology and accompanying reference profile(s) suitable for the quality evaluation of essential elements and/or metal contaminants in botanical ingredients.

Concerted expansion and contraction of immune receptor gene repertoires in plant genomes

Recent reports suggest that cell-surface and intracellular immune receptors function synergistically to activate robust defence against pathogens, but whether they co-evolve is unclear. Here we determined the numbers of cell-surface and intracellular immune receptors in 350 species. Surprisingly, the number of receptor genes that are predicted to encode cell-surface and intracellular immune receptors is strongly correlated. We suggest this is consistent with mutual potentiation of immunity initiated by cell-surface and intracellular receptors being reflected in the concerted co-evolution of the size of their repertoires across plant species.

A natural uORF variant confers phosphorus acquisition diversity in soybean

Phosphorus (P) is an essential element for all organisms. Because P fertilizers are a non-renewable resource and high fixation in soils, sustainable agriculture requires researchers to improve crop P acquisition efficiency. Here, we report a strong association signal at a locus of CPU1 (component of phosphorus uptake 1), from a genome-wide association study of P acquisition efficiency in a soybean core collection grown in the field. A SEC12-like gene, GmPHF1, is identified as the causal gene for CPU1. GmPHF1 facilitates the ER (endoplasmic reticulum) exit of the phosphate transporter, GmPT4, to the plasma membrane of root epidermal cells. A common SNP in an upstream open reading frame (uORF) of GmPHF1, which alters the abundance of GmPHF1 in a tissue-specific manner, contributes to P acquisition diversity in soybean. A natural genetic variation conditions diversity in soybean P acquisition, which can be used to develop P-efficient soybean genotypes.

POWR1 is a domestication gene pleiotropically regulating seed quality and yield in soybean

Seed protein, oil content and yield are highly correlated agronomically important traits that essentially account for the economic value of soybean. The underlying molecular mechanisms and selection of these correlated seed traits during soybean domestication are, however, less known. Here, we demonstrate that a CCT gene, POWR1, underlies a large-effect protein/oil QTL. A causative TE insertion truncates its CCT domain and substantially increases seed oil content, weight, and yield while decreasing protein content. POWR1 pleiotropically controls these traits likely through regulating seed nutrient transport and lipid metabolism genes. POWR1 is also a domestication gene. We hypothesize that the TE insertion allele is exclusively fixed in cultivated soybean due to selection for larger seeds during domestication, which significantly contributes to shaping soybean with increased yield/seed weight/oil but reduced protein content. This study provides insights into soybean domestication and is significant in improving seed quality and yield in soybean and other crop species.

Identification of positive and negative regulators of antiviral RNA interference in Arabidopsis thaliana

Virus-host coevolution often drives virus immune escape. However, it remains unknown whether natural variations of plant virus resistance are enriched in genes of RNA interference (RNAi) pathway known to confer essential antiviral defense in plants. Here, we report two genome-wide association study screens to interrogate natural variation among wild-collected Arabidopsis thaliana accessions in quantitative resistance to the endemic cucumber mosaic virus (CMV). We demonstrate that the highest-ranked gene significantly associated with resistance from both screens acts to regulate antiviral RNAi in ecotype Columbia-0. One gene, corresponding to Reduced Dormancy 5 (RDO5), enhances resistance by promoting amplification of the virus-derived small interfering RNAs (vsiRNAs). Interestingly, the second gene, designated Antiviral RNAi Regulator 1 (VIR1), dampens antiviral RNAi so its genetic inactivation by CRISPR/Cas9 editing enhances both vsiRNA production and CMV resistance. Our findings identify positive and negative regulators of the antiviral RNAi defense that may play important roles in virus-host coevolution.

Faster than expected Rubisco deactivation in shade reduces cowpea photosynthetic potential in variable light conditions

Cowpea is the major source of vegetable protein for rural populations in sub-Saharan Africa and average yields are not keeping pace with population growth. Each day, crop leaves experience many shade events and the speed of photosynthetic adjustment to this dynamic environment strongly affects daily carbon gain. Rubisco activity is particularly important because it depends on the speed and extent of deactivation in shade and recovers slowly on return to sun. Here, direct biochemical measurements showed a much faster rate of Rubisco deactivation in cowpea than prior estimates inferred from dynamics of leaf gas exchange in other species^1-3. Shade-induced deactivation was driven by decarbamylation, and half-times for both deactivation in shade and activation in saturating light were shorter than estimates from gas exchange (≤53% and 79%, respectively). Incorporating these half-times into a model of diurnal canopy photosynthesis predicted a 21% diurnal loss of productivity and suggests slowing Rubisco deactivation during shade is an unexploited opportunity for improving crop productivity.

Pepper Fruit Elongation Is Controlled by Capsicum annuum Ovate Family Protein 20

Fruit shape is one of the most important quality traits of pepper (Capsicum spp.) and is used as a major attribute for the classification of fruit types. Wide natural variation in fruit shape exists among the major cultivated species Capsicum annuum, allowing the identification of several QTLs controlling the trait. However, to date, no genes underlying fruit shape QTLs have been conclusively identified, nor has their function been verified in pepper. We constructed a mapping population from a cross of round- and elongated-fruited C. annuum parents and identified a single major QTL on chromosome 10, termed fs10, explaining 68 and 70% of the phenotypic variation for fruit shape index and for distal fruit end angle, respectively. The QTL was mapped in several generations and was localized to a 5 Mbp region containing the ortholog of SlOFP20 that suppresses fruit elongation in tomato. Virus-induced gene silencing of the pepper ortholog CaOFP20 resulted in increased fruit elongation on two independent backgrounds. Furthermore, CaOFP20 exhibited differential expression in fs10 near-isogenic lines, as well as in an association panel of elongated- and round-fruited accessions. A 42-bp deletion in the upstream region of CaOFP20 was most strongly associated with fruit shape variation within the locus. Histological observations in ovaries and fruit pericarps indicated that fs10 exerts its effect on fruit elongation by controlling cell expansion and replication. Our results indicate that CaOFP20 functions as a suppressor of fruit elongation in C. annuum and is the most likely candidate gene underlying fs10.

Branch water uptake and redistribution in two conifers at the alpine treeline

During winter, conifers at the alpine treeline suffer dramatic losses of hydraulic conductivity, which are successfully recovered during late winter. Previous studies indicated branch water uptake to support hydraulic recovery. We analyzed water absorption and redistribution in Picea abies and Larix decidua growing at the treeline by in situ exposure of branches to δ2H-labelled water. Both species suffered high winter embolism rates (> 40-60% loss of conductivity) and recovered in late winter (< 20%). Isotopic analysis showed water to be absorbed over branches and redistributed within the crown during late winter. Labelled water was redistributed over 425 ± 5 cm within the axes system and shifted to the trunk, lower and higher branches (tree height 330 ± 40 cm). This demonstrated relevant branch water uptake and re-distribution in treeline conifers. The extent of water absorption and re-distribution was species-specific, with L. decidua showing higher rates. In natura, melting snow might be the prime source for absorbed and redistributed water, enabling embolism repair and restoration of water reservoirs prior to the vegetation period. Pronounced water uptake in the deciduous L. decidua indicated bark to participate in the process of water absorption.

A giant NLR gene confers broad-spectrum resistance to Phytophthora sojae in soybean

Phytophthora root and stem rot caused by P. sojae is a destructive soybean soil-borne disease found worldwide. Discovery of genes conferring broad-spectrum resistance to the pathogen is a need to prevent the outbreak of the disease. Here, we show that soybean Rps11 is a 27.7-kb nucleotide-binding site-leucine-rich repeat (NBS-LRR or NLR) gene conferring broad-spectrum resistance to the pathogen. Rps11 is located in a genomic region harboring a cluster of large NLR genes of a single origin in soybean, and is derived from rounds of unequal recombination. Such events result in promoter fusion and LRR expansion that may contribute to the broad resistance spectrum. The NLR gene cluster exhibits drastic structural diversification among phylogenetically representative varieties, including gene copy number variation ranging from five to 23 copies, and absence of allelic copies of Rps11 in any of the non-Rps11-donor varieties examined, exemplifying innovative evolution of NLR genes and NLR gene clusters.

The APETALA2 homolog CaFFN regulates flowering time in pepper

Flowering time is an important agronomic trait that contributes to fitness in plants. However, the genetic basis of flowering time has not been extensively studied in pepper. To understand the genetics underlying flowering time, we constructed an F₂ population by crossing a spontaneous early flowering mutant and a late-flowering pepper line. Using bulked segregant RNA-seq, a major locus controlling flowering time in this population was mapped to the end of chromosome 2. An APETALA2 (AP2) homolog (CaFFN) cosegregated with flowering time in 297 individuals of the F₂ population. A comparison between the parents revealed a naturally occurring rare SNP (SNP2T > C) that resulted in the loss of a start codon in CaFFN in the early flowering mutant. Transgenic Nicotiana benthamiana plants with high CaFFN expression exhibited a delay in flowering time and floral patterning defects. On the other hand, pepper plants with CaFFN silencing flowered early. Therefore, the CaFFN gene acts as a flowering repressor in pepper. CaFFN may function as a transcriptional activator to activate the expression of CaAGL15 and miR156e and as a transcriptional repressor to repress the expression of CaAG, CaAP1, CaSEP3, CaSOC1, and miR172b based on a qRT-PCR assay. Direct activation of CaAGL15 by CaFFN was detected using yeast one-hybrid and dual-luciferase reporter assays, consistent with the hypothesis that CaFFN regulates flowering time. Moreover, the CaFFN gene association analysis revealed a significant association with flowering time in a natural pepper population, indicating that the CaFFN gene has a broad effect on flowering time in pepper. Finally, the phylogeny, evolutionary expansion and expression patterns of CaFFN/AP2 homologs were analyzed to provide valuable insight into CaFFN. This study increases our understanding of the involvement of CaFFN in controlling flowering time in pepper, thus making CaFFN a target gene for breeding early maturing pepper.

Genomic and transcriptomic analyses of Citrus sinensis varieties provide insights into Valencia orange fruit mastication trait formation

Valencia orange (Citrus sinensis Osbeck) (VO) is a type of late-ripening sweet orange whose ripening occurs 4 to 5 months later than that of the mid-ripening common sweet orange (CO). Notably, the mastication trait of VO fruit is inferior to that of CO fruit. To date, how inferior pulp mastication trait forms in VO has not been determined. In this study, 13 VO varieties and 12 CO varieties were subjected to whole-genome resequencing. A total of 2.98 million SNPs were identified from 25 varieties, and a SNP molecular marker was developed to distinguish VO and CO. Moreover, 144 and 141 genes identified by selective sweep analysis were selected during VO and CO evolution, respectively. Based on gene functional enrichment analysis, most of the selected VO genes were related to the stress response and lignin biosynthesis. Simultaneously, we comparatively analyzed the transcriptome profiles of peel and pulp tissues among three VO varieties and three CO varieties, and the results demonstrated differences in lignin biosynthesis between VO and CO fruits. Furthermore, coexpression network analysis was performed to identify hub genes of lignin-related and variety-specific networks, which included CsERF74, CsNAC25, CsHSFB3, CsSPL4/13, etc. Overall, this study provides important insights into the mastication trait formation of Valencia orange fruit.

Cannabis labelling is associated with genetic variation in terpene synthase genes

Analysis of over 100 Cannabis samples quantified for terpene and cannabinoid content and genotyped for over 100,000 single nucleotide polymorphisms indicated that Sativa- and Indica-labelled samples were genetically indistinct on a genome-wide scale. Instead, we found that Cannabis labelling was associated with variation in a small number of terpenes whose concentrations are controlled by genetic variation at tandem arrays of terpene synthase genes.

The identification of grain size genes by RapMap reveals directional selection during rice domestication

Cloning quantitative trait locus (QTL) is time consuming and laborious, which hinders the understanding of natural variation and genetic diversity. Here, we introduce RapMap, a method for rapid multi-QTL mapping by employing F2 gradient populations (F2GPs) constructed by minor-phenotypic-difference accessions. The co-segregation standard of the single-locus genetic models ensures simultaneous integration of a three-in-one framework in RapMap i.e. detecting a real QTL, confirming its effect, and obtaining its near-isogenic line-like line (NIL-LL). We demonstrate the feasibility of RapMap by cloning eight rice grain-size genes using 15 F2GPs in three years. These genes explain a total of 75% of grain shape variation. Allele frequency analysis of these genes using a large germplasm collection reveals directional selection of the slender and long grains in indica rice domestication. In addition, major grain-size genes have been strongly selected during rice domestication. We think application of RapMap in crops will accelerate gene discovery and genomic breeding.

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.

Multi-tissue transcriptome analysis of two Begonia species reveals dynamic patterns of evolution in the chalcone synthase gene family

Begonia is an important horticultural plant group, as well as one of the most speciose Angiosperm genera, with over 2000 described species. Genus wide studies of genome size have shown that Begonia has a highly variable genome size, and analysis of paralog pairs has previously suggested that Begonia underwent a whole genome duplication. We address the contribution of gene duplication to the generation of diversity in Begonia using a multi-tissue RNA-seq approach. We chose to focus on chalcone synthase (CHS), a gene family having been shown to be involved in biotic and abiotic stress responses in other plant species, in particular its importance in maximising the use of variable light levels in tropical plants. We used RNA-seq to sample six tissues across two closely related but ecologically and morphologically divergent species, Begonia conchifolia and B. plebeja, yielding 17,012 and 19,969 annotated unigenes respectively. We identified the chalcone synthase gene family members in our Begonia study species, as well as in Hillebrandia sandwicensis, the monotypic sister genus to Begonia, Cucumis sativus, Arabidopsis thaliana, and Zea mays. Phylogenetic analysis suggested the CHS gene family has high duplicate turnover, all members of CHS identified in Begonia arising recently, after the divergence of Begonia and Cucumis. Expression profiles were similar within orthologous pairs, but we saw high inter-ortholog expression variation. Sequence analysis showed relaxed selective constraints on some ortholog pairs, with substitutions at conserved sites. Evidence of pseudogenisation and species specific duplication indicate that lineage specific differences are already beginning to accumulate since the divergence of our study species. We conclude that there is evidence for a role of gene duplication in generating diversity through sequence and expression divergence in Begonia.

Region-level epimutation rates in Arabidopsis thaliana

Failure to maintain DNA methylation patterns during plant development can occasionally give rise to so-called "spontaneous epimutations". These stochastic methylation changes are sometimes heritable across generations and thus accumulate in plant genomes over time. Recent evidence indicates that spontaneous epimutations have a major role in shaping patterns of methylation diversity in plant populations. Using single CG dinucleotides as units of analysis, previous work has shown that the epimutation rate is several orders of magnitude higher than the genetic mutation rate. While these large rate differences have obvious implications for understanding genome-methylome co-evolution, the functional relevance of single CG methylation changes remains questionable. In contrast to single CG, solid experimental evidence has linked methylation gains and losses in larger genomic regions with transcriptional variation and heritable phenotypic effects. Here we show that such region-level changes arise stochastically at about the same rate as those at individual CG sites, are only marginal dependent on region size and cytosine density, but strongly dependent on chromosomal location. We also find consistent evidence that region-level epimutations are not restricted to CG contexts but also frequently occur in non-CG regions at the genome-wide scale. Taken together, our results support the view that many differentially methylated regions (DMRs) in natural populations originate from epimutation events and may not be effectively tagged by proximal SNPs. This possibility reinforces the need for epigenome-wide association studies (EWAS) in plants as a way to identify the epigenetic basis of complex traits.