Maca (Lepidium meyenii Walp, 2n = 8x = 64), belonging to the Brassicaceae family, is an economic plant cultivated in the central Andes sierra in Peru (4000-4500 m). Considering that the rapid uplift of the central Andes occurred 5-10 million years ago (Ma), an evolutionary question arises regarding how plants such as maca acquire high-altitude adaptation within a short geological period. Here, we report the high-quality genome assembly of maca, in which two closely spaced maca-specific whole-genome duplications (WGDs; ∼6.7 Ma) were identified. Comparative genomic analysis between maca and closely related Brassicaceae species revealed expansions of maca genes and gene families involved in abiotic stress response, hormone signaling pathway, and secondary metabolite biosynthesis via WGDs. The retention and subsequent functional divergence of many duplicated genes may account for the morphological and physiological changes (i.e., small leaf shape and self-fertility) in maca in a high-altitude environment. In addition, some duplicated maca genes were identified with functions in morphological adaptation (i.e., LEAF CURLING RESPONSIVENESS) and abiotic stress response (i.e., GLYCINE-RICH RNA-BINDING PROTEINS and DNA-DAMAGE-REPAIR/TOLERATION 2) under positive selection. Collectively, the maca genome provides useful information to understand the important roles of WGDs in the high-altitude adaptation of plants in the Andes.

Many plants are facultatively asexual, balancing short-term benefits with long-term costs of asexuality. During range expansion, natural selection likely influences the genetic controls of asexuality in these organisms. However, evidence of natural selection driving asexuality is limited, and the evolutionary consequences of asexuality on the genomic and epigenomic diversity remain controversial. We analyzed population genomes and epigenomes of Spirodela polyrhiza, (L.) Schleid., a facultatively asexual plant that flowers rarely, revealing remarkably low genomic diversity and DNA methylation levels. Within species, demographic history and the frequency of asexual reproduction jointly determined intra-specific variations of genomic diversity and DNA methylation levels. Genome-wide scans revealed that genes associated with stress adaptations, flowering and embryogenesis were under positive selection. These data are consistent with the hypothesize that natural selection can shape the evolution of asexuality during habitat expansions, which alters genomic and epigenomic diversity levels.

Many plants emit diverse floral scents that mediate plant–environment interactions and attain reproductive success. However, how plants evolve novel and adaptive biosynthetic pathways for floral volatiles remains unclear. Here, we show that in the wild tobacco, Nicotiana attenuata, a dominant species-specific floral volatile (benzyl acetone, BA) that attracts pollinators and deters florivore is synthesized by phenylalanine ammonia-lyase 4 (NaPAL4), isoflavone reductase 3 (NaIFR3), and chalcone synthase 3 (NaCHAL3). Transient expression of NaFIR3 alone in N. attenuata leaves is sufficient and necessary for ectopic foliar BA emissions, and coexpressing NaIFR3 with NaPAL4 and NaCHAL3 increased the BA emission levels. Independent changes in transcription of NaPAL4 and NaCHAL3 contributed to intraspecific variations of floral BA emission. However, among species, the gain of expression of NaIFR3 resulted in the biosynthesis of BA, which was only found in N. attenuata. This study suggests that novel metabolic pathways associated with adaptation can arise via reconfigurations of gene expression.

The Colorado potato beetle (CPB) is a major pest of potato crops, known for its remarkable ability to develop resistance to more than 50 pesticides. For decades, CPB has served as a model species for studying insecticide resistance, insect physiology, diapause, reproduction, and evolution. However, research progress on CPB has been hindered by the lack of comprehensive genomic and transcriptomic resources. Here, leveraging a recently established chromosome-level genome assembly, we constructed a gene expression atlas of CPB using transcriptomic data from 61 samples representing major organs and developmental stages. By integrating short- and long-read sequencing technologies, we enhanced the genome annotation and identified 6,623 additional genes that were previously undetected. Furthermore, we developed a web portal to facilitate the search and visualization of the gene expression atlas, providing an accessible resource for the research community. The CPB gene expression atlas offers valuable tools and comprehensive data that will accelerate future research in pest control and insect biology.

The immune systems are fundamentally vital for evolution and survival of species; as such, selection patterns in innate immune loci are of special interest in molecular evolutionary research. The interferon regulatory factor (IRF) gene family control many different aspects of the innate and adaptive immune responses in vertebrates. Among these, IRF3 is known to take active part in very many biological processes. We assembled and evaluated 1356 base pairs of the IRF3 gene coding region in domesticated goats from Africa (Nigeria, Ethiopia and South Africa) and Asia (Iran and China) and the wild goat (Capra aegagrus). Five segregating sites with θ value of 0.0009 for this gene demonstrated a low diversity across the goats’ populations. Fu and Li tests were significantly positive but Tajima’s D test was significantly negative, suggesting its deviation from neutrality. Neighbor joining tree of IRF3 gene in domesticated goats, wild goat and sheep showed that all domesticated goats have a closer relationship than with the wild goat and sheep. Maximum likelihood tree of the gene showed that different domesticated goats share a common ancestor and suggest single origin. Four unique haplotypes were observed across all the sequences, of which, one was particularly common to African goats (MOCH-K14-0425, Poitou and WAD). In assessing the evolution mode of the gene, we found that the codon model d_N/d_S ratio for all goats was greater than one. Phylogenetic Analysis by Maximum Likelihood (PAML) gave a ω₀ (d_N/d_S) value of 0.067 with LnL value of -6900.3 for the first Model (M1) while ω2 = 1.667 in model M2 with LnL value of -6900.3 with positive selection inferred in 3 codon sites. Mechanistic empirical combination (MEC) model for evaluating adaptive selection pressure on particular codons also confirmed adaptive selection pressure in three codons (207, 358 and 408) in IRF3 gene. Positive diversifying selection inferred with recent evolutionary changes in domesticated goat IRF3 led us to conclude that the gene evolution may have been influenced by domestication processes in goats.

Understanding non-target-site resistance (NTSR) to herbicides represents a pressing challenge as NTSR is widespread in many weeds. Using giant duckweed (Spirodela polyrhiza) as a model, we systematically investigated genetic and molecular mechanisms of diquat resistance, which can only be achieved via NTSR. Quantifying the diquat resistance of 138 genotypes, we revealed an 8.5-fold difference in resistance levels between the most resistant and most susceptible genotypes. Further experiments suggested that diquat uptake and antioxidant-related processes jointly contributed to diquat resistance in S. polyrhiza. Using a genome-wide association approach, we identified several candidate genes, including a homolog of dienelactone hydrolase, that are associated with diquat resistance in S. polyrhiza. Together, these results provide new insights into the mechanisms and evolution of NTSR in plants.