Population Genomic Analysis and De Novo Assembly Reveal the Origin of Weedy Rice as an Evolutionary Game

Weed biology and management in the multi-omics era: Progress and perspectives

Weeds pose a significant threat to crop production, resulting in substantial yield reduction. In addition, they possess robust weedy traits that enable them to survive in extreme environments and evade human control. In recent years, the application of multi-omics biotechnologies has helped to reveal the molecular mechanisms underlying these weedy traits. In this review, we systematically describe diverse applications of multi-omics platforms for characterizing key aspects of weed biology, including the origins of weed species, weed classification, and the underlying genetic and molecular bases of important weedy traits such as crop-weed interactions, adaptability to different environments, photoperiodic flowering responses, and herbicide resistance. In addition, we discuss limitations to the application of multi-omics techniques in weed science, particularly compared with their extensive use in model plants and crops. In this regard, we provide a forward-looking perspective on the future application of multi-omics technologies to weed science research. These powerful tools hold great promise for comprehensively and efficiently unraveling the intricate molecular genetic mechanisms that underlie weedy traits. The resulting advances will facilitate the development of sustainable and highly effective weed management strategies, promoting greener practices in agriculture.

Harnessing weedy rice as functional food and source of novel traits for crop improvement

A relative of cultivated rice (Oryza sativa L.), weedy or red rice (Oryza spp.) is currently recognized as the dominant weed, leading to a drastic loss of yield of cultivated rice due to its highly competitive abilities like producing more tillers, panicles, and biomass with better nutrient uptake. Due to its high nutritional value, antioxidant properties (anthocyanin and proanthocyanin), and nutrient absorption ability, weedy rice is gaining immense research attentions to understand its genetic constitution to augment future breeding strategies and to develop nutrition-rich functional foods. Consequently, this review focuses on the unique gene source of weedy rice to enhance the cultivated rice for its crucial features like water use efficiency, abiotic and biotic stress tolerance, early flowering, and the red pericarp of the seed. It explores the debating issues on the origin and evolution of weedy rice, including its high diversity, signalling aspects, quantitative trait loci (QTL) mapping under stress conditions, the intricacy of the mechanism in the expression of the gene flow, and ecological challenges of nutrient removal by weedy rice. This review may create a foundation for future researchers to understand the gene flow between cultivated crops and weedy traits and support an improved approach for the applicability of several models in predicting multiomics variables.

A derived weedy rice × ancestral cultivar cross identifies evolutionarily relevant weediness QTLs

Weedy rice (Oryza spp.) is a weedy relative of the cultivated rice that competes with the crop and causes significant production loss. The BHA (blackhull awned) US weedy rice group has evolved from aus cultivated rice and differs from its ancestors in several important weediness traits, including flowering time, plant height and seed shattering. Prior attempts to determine the genetic basis of weediness traits in plants using linkage mapping approaches have not often considered weed origins. However, the timing of divergence between crossed parents can affect the detection of quantitative trait loci (QTL) relevant to the evolution of weediness. Here, we used a QTL-seq approach that combines bulked segregant analysis and high-throughput whole genome resequencing to map the three important weediness traits in an F2 population derived from a cross between BHA weedy rice with an ancestral aus cultivar. We compared these QTLs with those previously detected in a cross of BHA with a more distantly related crop, indica. We identified multiple QTLs that overlapped with regions under selection during the evolution of weedy BHA rice and some candidate genes possibly underlying the evolution weediness traits in BHA. We showed that QTLs detected with ancestor-descendant crosses are more likely to be involved in the evolution of weediness traits than those detected from crosses of more diverged taxa.

Integrated genome-wide differentiation and association analyses identify causal genes underlying breeding-selected grain quality traits in japonica rice

Improving grain quality is a primary objective in contemporary rice breeding. Japanese modern rice breeding has developed two different types of rice, eating and sake-brewing rice, with different grain characteristics, indicating the selection of variant gene alleles during the breeding process. Given the critical importance of promptly and efficiently identifying genes selected in past breeding for future molecular breeding, we conducted genome scans for divergence, genome-wide association studies, and map-based cloning. Consequently, we successfully identified two genes, OsMnS and OsWOX9D, both contributing to rice grain traits. OsMnS encodes a mannan synthase that increases the white core frequency in the endosperm, a desirable trait for sake brewing but decreases the grain appearance quality. OsWOX9D encodes a grass-specific homeobox-containing transcription factor, which enhances grain width for better sake brewing. Furthermore, haplotype analysis revealed that their defective alleles were selected in East Asia, but not Europe, during modern improvement. In addition, our analyses indicate that a reduction in grain mannan content during African rice domestication may also be caused a defective OsMnS allele due to breeding selection. This study not only reveals the delicate balance between grain appearance quality and nutrition in rice but also provides a new strategy for isolating causal genes underlying complex traits, based on the concept of "breeding-assisted genomics" in plants.

A syntelog-based pan-genome provides insights into rice domestication and de-domestication

BACKGROUND

Asian rice is one of the world's most widely cultivated crops. Large-scale resequencing analyses have been undertaken to explore the domestication and de-domestication genomic history of Asian rice, but the evolution of rice is still under debate.

RESULTS

Here, we construct a syntelog-based rice pan-genome by integrating and merging 74 high-accuracy genomes based on long-read sequencing, encompassing all ecotypes and taxa of Oryza sativa and Oryza rufipogon. Analyses of syntelog groups illustrate subspecies divergence in gene presence-and-absence and haplotype composition and identify massive genomic regions putatively introgressed from ancient Geng/japonica to ancient Xian/indica or its wild ancestor, including almost all well-known domestication genes and a 4.5-Mbp centromere-spanning block, supporting a single domestication event in main rice subspecies. Genomic comparisons between weedy and cultivated rice highlight the contribution from wild introgression to the emergence of de-domestication syndromes in weedy rice.

CONCLUSIONS

This work highlights the significance of inter-taxa introgression in shaping diversification and divergence in rice evolution and provides an exploratory attempt by utilizing the advantages of pan-genomes in evolutionary studies.

Advances in Rice Seed Shattering

Seed shattering is an important trait that wild rice uses to adapt to the natural environment and maintain population reproduction, and weedy rice also uses it to compete with the rice crop. The loss of shattering is a key event in rice domestication. The degree of shattering is not only one of the main reasons for rice yield reduction but also affects its adaptability to modern mechanical harvesting methods. Therefore, it is important to cultivate rice varieties with a moderate shattering degree. In this paper, the research progress on rice seed shattering in recent years is reviewed, including the physiological basis, morphological and anatomical characteristics of rice seed shattering, inheritance and QTL/gene mapping of rice seed shattering, the molecular mechanism regulating rice seed shattering, the application of seed-shattering genes, and the relationship between seed-shattering genes and domestication.

Weed genomics: yielding insights into the genetics of weedy traits for crop improvement

Weeds cause tremendous economic and ecological damage worldwide. The number of genomes established for weed species has sharply increased during the recent decade, with some 26 weed species having been sequenced and de novo genomes assembled. These genomes range from 270 Mb (Barbarea vulgaris) to almost 4.4 Gb (Aegilops tauschii). Importantly, chromosome-level assemblies are now available for 17 of these 26 species, and genomic investigations on weed populations have been conducted in at least 12 species. The resulting genomic data have greatly facilitated studies of weed management and biology, especially origin and evolution. Available weed genomes have indeed revealed valuable weed-derived genetic materials for crop improvement. In this review, we summarize the recent progress made in weed genomics and provide a perspective for further exploitation in this emerging field.

Transcriptomic and QTL Analysis of Seed Germination Vigor under Low Temperature in Weedy Rice WR04-6

Low temperature is one of the major factors affecting rice germination, and low temperature germination (LTG) is an important agronomic trait. Although significant progress has been made in the study of rice LTG, the molecular mechanism of LTG remains poorly understood. To explore more rice LTG gene resources, we first demonstrated that weedy rice WR04-6 (Oryza sativa f. spontanea) had significantly higher LTG ability at 10 °C than the cultivated rice Qishanzhan (QSZ Oryza sativa L. ssp. indica). RNA-seq was used to investigate the gene expression of WR04-6 and QSZ at 10 °C for 10, 12 and 14 days after imbibition (DAI) of seed germination. The results of Gene Ontology (GO) enrichment and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment revealed that the differentially expressed genes (DEGs) between WR04-6 and QSZ were mainly concentrated on the response to starch catabolic processes and the response to abscisic acid (ABA). This is consistent with the results of α-amylase activity, ABA and gibberellins (GA) treatment. A recombinant inbred line (RIL) population derived from a cross between WR04-6 and QSZ and its high-density SNP genetic map were used to detect quantitative trait loci (QTL) for LTG rates. The results showed that two new QTLs were located on chromosome 3 and chromosome 12. Combined with the mapped QTLs and RNA-seq DEGs, sixteen candidate genes potentially associated with LTG were identified. Validation of the expression of the candidates by qRT-PCR were consistent with the RNA-seq data. These results will enable us to understand the genetic basis of LTG in weedy rice and provide new genetic resources for the generation of rice germplasm with improved LTG.

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.

Global whole-genome comparison and analysis to classify subpopulations and identify resistance genes in weedy rice relevant for improving crops

Common weedy rice plants are important genetic resources for modern breeding programs because they are the closest relatives to rice cultivars and their genomes contain elite genes. Determining the utility and copy numbers of WRKY and nucleotide-binding site (NBS) resistance-related genes may help to clarify their variation patterns and lead to crop improvements. In this study, the weedy rice line LM8 was examined at the whole-genome level. To identify the Oryza sativa japonica subpopulation that LM8 belongs to, the single nucleotide polymorphisms (SNPs) of 180 cultivated and 23 weedy rice varieties were used to construct a phylogenetic tree and a principal component analysis and STRUCTURE analysis were performed. The results indicated that LM8 with admixture components from japonica (GJ) and indica (XI) belonged to GJ-admixture (GJ-adm), with more than 60% of its genetic background derived from XI-2 (22.98%), GJ-tropical (22.86%), and GJ-subtropical (17.76%). Less than 9% of its genetic background was introgressed from weedy rice. Our results also suggested LM8 may have originated in a subtropical or tropical geographic region. Moreover, the comparisons with Nipponbare (NIP) and Shuhui498 (R498) revealed many specific structure variations (SVs) in the LM8 genome and fewer SVs between LM8 and NIP than between LM8 and R498. Next, 96 WRKY and 464 NBS genes were identified and mapped on LM8 chromosomes to eliminate redundancies. Three WRKY genes (ORUFILM02g002693, ORUFILM05g002725, and ORUFILM05g001757) in group III and one RNL [including the resistance to powdery mildew 8 (RPW8) domain, NBS, and leucine rich repeats (LRRs)] type NBS gene (ORUFILM12g000772) were detected in LM8. Among the NBS genes, the RPW8 domain was detected only in ORUFILM12g000772. This gene may improve plant resistance to pathogens as previously reported. Its classification and potential utility imply LM8 should be considered as a germplasm resource relevant for rice breeding programs.

Genetic basis of the early heading of high-latitude weedy rice

Japonica rice (Oryza sativa L.) is an important staple food in high-latitude regions and is widely distributed in northern China, Japan, Korea, and Europe. However, the genetic diversity of japonica rice is relatively narrow and poorly adapted. Weedy rice (Oryza sativa f. spontanea) is a semi-domesticated rice. Its headings are earlier than the accompanied japonica rice, making it a potential new genetic resource, which can make up for the defects of wild rice that are difficult to be directly applied to japonica rice improvement caused by reproductive isolation. In this study, we applied a natural population consisting of weedy rice, japonica landrace, and japonica cultivar to conduct a genome-wide association study (GWAS) of the heading date and found four loci that could explain the natural variation of the heading date in this population. At the same time, we developed recombinant inbred lines (RILs) crossed by the early-heading weedy rice WR04-6 and its accompanied japonica cultivar ShenNong 265 (SN265) to carry out a QTL mapping analysis of the heading date and mapped four quantitative trait locus (QTLs) and three epistatic effect gene pairs. The major locus on chromosome 6 overlapped with the GWAS result. Further analysis found that two genes, Hd1 and OsCCT22, on chromosome 6 (Locus 2 and Locus 3) may be the key points of the early-heading character of weedy rice. As minor effect genes, Dth7 and Hd16 also have genetic contributions to the early heading of weedy rice. In the process of developing the RIL population, we introduced fragments of Locus 2 and Locus 3 from the weedy rice into super-high-yielding japonica rice, which successfully promoted its heading date by at least 10 days and expanded the rice suitable cultivation area northward by about 400 km. This study successfully revealed the genetic basis of the early heading of weedy rice and provided a new idea for the genetic improvement of cultivated rice by weedy rice.

Hybrid-derived weedy rice maintains adaptive combinations of alleles associated with seed dormancy

Plant hybridization is a pathway for the evolution of adaptive traits. However, hybridization between adapted and nonadapted populations may affect the persistence of combinations of adaptive alleles evolved through natural selection. Seed dormancy is an adaptive trait for weedy rice because it regulates the timing of seed germination and the persistence of the soil seed bank. Hybridization between weedy and cultivated rice has been confirmed with an adaptive introgression of deep seed dormancy alleles from cultivated rice. Here, we explored the influence of hybridization on the conservation of adaptive allele combinations by evaluating natural variation and genetic structure in seed dormancy-associated genomic regions. Based on sequence variation in the genomic regions associated with seed dormancy, hybrid-derived weedy rice strains maintained most of the adaptive combinations for this trait observed in the parental weedy rice, despite equal representation of the parental weedy and cultivated rice in the whole genome sequence. Moreover, hybrid-derived weedy rice strains were more dormant than their parental weedy rice strains, and this trait was strongly influenced by the environment. This study suggests that hybridization between weedy rice (adaptive allelic combinations for seed dormancy) and cultivated rice (nonadaptive combinations) generates weedy rice strains expressing deep seed dormancy caused by genome stabilization through the removal of alleles derived from cultivated rice, in addition to the adaptive introgression of deep seed dormancy alleles derived from cultivated rice. Thus, hybridization between adapted and nonadapted populations appears to be reinforcing the trajectory towards the evolution of adaptive traits.

Prospects of Feral Crop De Novo Redomestication

Modern agriculture depends on a narrow variety of crop species, leaving global food and nutritional security highly vulnerable to the adverse effects of climate change and population expansion. Crop improvement using conventional and molecular breeding approaches leveraging plant genetic diversity using crop wild relatives (CWRs) has been one approach to address these issues. However, the rapid pace of the global change requires additional innovative solutions to adapt agriculture to meet global needs. Neodomestication-the rapid and targeted introduction of domestication traits using introgression or genome editing of CWRs-is being explored as a supplementary approach. These methods show promise; however, they have so far been limited in efficiency and applicability. We propose expanding the scope of neodomestication beyond truly wild CWRs to include feral crops as a source of genetic diversity for novel crop development, in this case 'redomestication'. Feral crops are plants that have escaped cultivation and evolved independently, typically adapting to their local environments. Thus, feral crops potentially contain valuable adaptive features while retaining some domestication traits. Due to their genetic proximity to crop species, feral crops may be easier targets for de novo domestication (i.e. neodomestication via genome editing techniques). In this review, we explore the potential of de novo redomestication as an application for novel crop development by genome editing of feral crops. This approach to efficiently exploit plant genetic diversity would access an underutilized reservoir of genetic diversity that could prove important in support of global food insecurity in the face of the climate change.

Compared analysis with a high-quality genome of weedy rice reveals the evolutionary game of de-domestication

The weedy rice (Oryza sativa f. spontanea) harbors large numbers of excellent traits and genetic diversities, which serves as a valuable germplasm resource and has been considered as a typical material for research about de-domestication. However, there are relatively few reference genomes on weedy rice that severely limit exploiting these genetic resources and revealing more details about de-domestication events. In this study, a high-quality genome (~376.4 Mb) of weedy rice A02 was assembled based on Nanopore ultra-long platform with a coverage depth of about 79.3× and 35,423 genes were predicted. Compared to Nipponbare genome, 5,574 structural variations (SVs) were found in A02. Based on super pan-genome graph, population SVs of 238 weedy rice and cultivated rice accessions were identified using public resequencing data. Furthermore, the de-domestication sites of weedy rice and domestication sites of wild rice were analyzed and compared based on SVs and single-nucleotide polymorphisms (SNPs). Interestingly, an average of 2,198 genes about de-domestication could only be found by F _ST analysis based on SVs (SV-F _ST) while not by F _ST analysis based on SNPs (SNP-F _ST) in divergent region. Additionally, there was a low overlap between domestication and de-domestication intervals, which demonstrated that two different mechanisms existed in these events. Our finding could facilitate pinpointing of the evolutionary events that had shaped the genomic architecture of wild, cultivated, and weedy rice, and provide a good foundation for cloning of the superior alleles for breeding.

The rhizosphere bacterial community contributes to the nutritional competitive advantage of weedy rice over cultivated rice in paddy soil

Background

Weedy rice competes for nutrients and living space with cultivated rice, which results in serious reductions in rice production. The rhizosphere bacterial community plays an important role in nutrient competition between species. It is therefore important to clarify the differences in the diversities of the inter rhizosphere bacterial community between cultivated rice and weedy rice. The differences in compositions and co-occurrence networks of the rhizosphere bacterial community of cultivated rice and weedy rice are largely unknown and thus the aim of our study.

Results

In our study, the different rhizosphere bacterial community structures in weedy rice (AW), cultivated rice (AY) and cultivated rice surrounded by weedy rice (WY) were determined based on 16S rRNA gene sequencing. The majority of the WY rhizosphere was enriched with unique types of microorganisms belonging to Burkholderia. The rhizosphere bacterial community showed differences in relative abundance among the three groups. Network analysis revealed a more complex co-occurrence network structure in the rhizosphere bacterial community of AW than in those of AY and WY due to a higher degree of Microbacteriaceae and Micrococcaceae in the network. Both network analysis and functional predictions reveal that weedy rice contamination dramatically impacts the iron respiration of the rhizosphere bacterial community of cultivated rice.

Conclusions

Our study shows that there are many differences in the rhizosphere bacterial community of weedy rice and cultivated rice. When cultivated rice was disturbed by weedy rice, the rhizosphere bacterial community and co-occurrence network also changed. The above differences tend to lead to a nutritional competitive advantage for weedy rice in paddy soils.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12866-022-02648-1.

Regain flood adaptation in rice through a 14-3-3 protein OsGF14h

Contemporary climatic stress seriously affects rice production. Unfortunately, long-term domestication and improvement modified the phytohormones network to achieve the production needs of cultivated rice, thus leading to a decrease in adaptation. Here, we identify a 14-3-3 protein-coding gene OsGF14h in weedy rice that confers anaerobic germination and anaerobic seedling development tolerance. OsGF14h acts as a signal switch to balance ABA signaling and GA biosynthesis by interacting with the transcription factors OsHOX3 and OsVP1, thereby boosting the seeding rate from 13.5% to 60.5% for anaerobic sensitive variety under flooded direct-seeded conditions. Meanwhile, OsGF14h co-inheritance with the Rc (red pericarp gene) promotes divergence between temperate japonica cultivated rice and temperate japonica weedy rice through artificial and natural selection. Our study retrieves a superior allele that has been lost during modern japonica rice improvement and provides a fine-tuning tool to improve flood adaptation for elite rice varieties.

Waterlogging tolerance is important in direct seeding rice cultivation practice. Here, the authors identify a 14-3-3 protein-coding gene OsGF14h in weedy rice that confers anaerobic germination and anaerobic seedling development tolerance by balancing ABA signaling and GA biosynthesis.

The complete chloroplast genome of japonica type weedy rice (Oryza sativa f. spontanea)

As a noxious weed, weedy rice (Oryza sativa f. spontanea Roshev. 1931) has threatened global food security and sustainable crop production. On the other hand, weedy rice has a strong tolerance for abiotic stresses and the potential to provide rich resources for rice genetic improvement. Thus, for a more comprehensive understanding of its speciation, we sequenced and assembled the first complete chloroplast genome of Oryza sativa f. spontanea (japonica type). The complete chloroplast genome was 134,555 bp in length and encoded 133 genes, including 83 protein-coding genes, 42 tRNA genes and 8 rRNA genes. Phylogenetic analysis revealed that the indica-japonica differentiation of weedy rice was closely related to cultivated rice, and Oryza sativa f. spontanea (japonica type) was genetically more closely clustered with cultivated rice O. sativa (japonica type) than O. nivara or other wild rice.

Identification of candidate genes and clarification of the maintenance of the green pericarp of weedy rice grains

The weedy rice (Oryza sativa f. spontanea) pericarp has diverse colors (e.g., purple, red, light-red, and white). However, research on pericarp colors has focused on red and purple, but not green. Unlike many other common weedy rice resources, LM8 has a green pericarp at maturity. In this study, the coloration of the LM8 pericarp was evaluated at the cellular and genetic levels. First, an examination of their ultrastructure indicated that LM8 chloroplasts were normal regarding plastid development and they contained many plastoglobules from the early immature stage to maturity. Analyses of transcriptome profiles and differentially expressed genes revealed that most chlorophyll (Chl) degradation-related genes in LM8 were expressed at lower levels than Chl a/b cycle-related genes in mature pericarps, suggesting that the green LM8 pericarp was associated with inhibited Chl degradation in intact chloroplasts. Second, the F₂ generation derived from a cross between LM8 (green pericarp) and SLG (white pericarp) had a pericarp color segregation ratio of 9:3:4 (green:brown:white). The bulked segregant analysis of the F₂ populations resulted in the identification of 12 known genes in the chromosome 3 and 4 hotspot regions as candidate genes related to Chl metabolism in the rice pericarp. The RNA-seq and sqRT-PCR assays indicated that the expression of the Chl a/b cycle-related structural gene DVR (encoding divinyl reductase) was sharply up-regulated. Moreover, genes encoding magnesium-chelatase subunit D and the light-harvesting Chl a/b-binding protein were transcriptionally active in the fully ripened dry pericarp. Regarding the ethylene signal transduction pathway, the CTR (encoding an ethylene-responsive protein kinase) and ERF (encoding an ethylene-responsive factor) genes expression profiles were determined. The findings of this study highlight the regulatory roles of Chl biosynthesis- and degradation-related genes influencing Chl accumulation during the maturation of the LM8 pericarp.

Beyond a reference genome: pangenomes and population genomics of underutilized and orphan crops for future food and nutrition security

Underutilized crops are, by definition, under-researched compared to staple crops yet come with traits that may be especially important given climate change and the need to feed a globally increasing population. These crops are often stress-tolerant, and this combined with unique and beneficial nutritional profiles. Whilst progress is being made by generating reference genome sequences, in this Tansley Review, we show how this is only the very first step. We advocate that going 'beyond a reference genome' should be a priority, as it is only at this stage one can identify the specific genes and the adaptive alleles that underpin the valuable traits. We sum up how population genomic and pangenomic approaches have led to the identification of stress- and disease-tolerant alleles in staple crops and compare this to the small number of examples from underutilized crops. We also demonstrate how previously underutilized crops have benefitted from genomic advances and that many breeding targets in underutilized crops are often well studied in staple crops. This cross-crop population-level resequencing could lead to an understanding of the genetic basis of adaptive traits in underutilized crops. This level of investment may be crucial for fully understanding the value of these crops before they are lost.

Evidence for evolution and selection of drought-resistant genes based on high-throughput resequencing in weedy rice

Weedy rice (Oryza sativa f. spontanea) is a relative of cultivated rice that propagates in paddy fields and has strong drought resistance. In this study, we used 501 rice accessions to reveal the selection mechanism of drought resistance in weedy rice through a combination of selection analysis, genome-wide association studies, gene knockout and overexpression analysis, and Ca2+ and K+ ion flux assays. The results showed that the weedy rice species investigated have gene introgression with cultivated rice, which is consistent with the hypothesis that weedy rice originated from de-domestication of cultivated rice. Regions related to tolerance have particularly diversified during de-domestication and three drought-tolerance genes were identified. Of these, Os01g0800500 was also identified using an assay of the degree of leaf withering under drought, and it was named as PAPH1, encoding a PAP family protein. The drought-resistance capacity of PAPH1-knockout lines was much lower than that of the wild type, while that of overexpression lines was much higher. Concentrations of Ca2+ and K+ were lower in the knockout lines and higher in the overexpression lines compared with those of the wild type, suggesting that PAPH1 plays important roles in coping with drought stress. Our study therefore provides new insights into the genetic mechanisms underlying adaptive tolerance to drought in wild rice and highlights potential new resistance genes for future breeding programs in cultivated rice.

Population differentiation and epidemic tracking of Bursaphelenchus xylophilus in China based on chromosome-level assembly and whole-genome sequencing data

BACKGROUND

Bursaphelenchus xylophilus, the pinewood nematode, kills millions of pine trees worldwide every year, and causes enormous economic and ecological losses. Despite extensive research on population variation, there is little understanding of the population-wide variation spectrum in China.

RESULTS

We sequenced an inbred B. xylophilus strain using Pacbio+Illumina+Bionano+Hi-C and generated a chromosome-level assembly (AH1) with six chromosomes of 77.1 Mb (chromosome N50: 12 Mb). The AH1 assembly shows very high continuity and completeness, and contains novel genes with potentially important functions compared with previous assemblies. Subsequently, we sequenced 181 strains from China and the USA and found ~7.8 million single nucleotide polymorphisms (SNPs). Analysis shows that the B. xylophilus population in China can be divided into geographically bounded subpopulations with severe cross-infection and potential migrations. In addition, distribution of B. xylophilus is dominated by temperature zones while geographically associated SNPs are mainly located on adaptation related GPCR gene families, suggesting the nematode has been evolving to adapt to different temperatures. A machine-learning based epidemic tracking method has been established to predict their geographical origins, which can be applied to any other species.

CONCLUSION

Our study provides the community with the first high-quality chromosome-level assembly which includes a comprehensive catalogue of genetic variations. It provides insights into population structure and effective tracking method for this invasive species, which facilitates future studies to address a variety of applied, genomic and evolutionary questions in B. xylophilus as well as related species.

Whole-Genome Sequencing and RNA-Seq Reveal Differences in Genetic Mechanism for Flowering Response between Weedy Rice and Cultivated Rice

Flowering is a key agronomic trait that influences adaptation and productivity. Previous studies have indicated the genetic complexity associated with the flowering response in a photoinsensitive weedy rice accession PSRR-1 despite the presence of a photosensitive allele of a key flowering gene Hd1. In this study, we used whole-genome and RNA sequencing data from both cultivated and weedy rice to add further insights. The de novo assembly of unaligned sequences predicted 225 genes, in which 45 were specific to PSRR-1, including two genes associated with flowering. Comparison of the variants in PSRR-1 with the 3K rice genome (RG) dataset identified unique variants within the heading date QTLs. Analyses of the RNA-Seq result under both short-day (SD) and long-day (LD) conditions revealed that many differentially expressed genes (DEGs) colocalized with the flowering QTLs, and some DEGs such as Hd1, OsMADS56, Hd3a, and RFT1 had unique variants in PSRR-1. Ehd1, Hd1, OsMADS15, and OsMADS56 showed different alternate splicing (AS) events between genotypes and day length conditions. OsMADS56 was expressed in PSRR-1 but not in Cypress under both LD and SD conditions. Based on variations in both sequence and expression, the unique flowering response in PSRR-1 may be due to the high-impact variants of flowering genes, and OsMADS56 is proposed as a key regulator for its day-neutral flowering response.

Artificial selection in the expansion of rice cultivation

Gene distributions and population genomics suggest artificial selection of ghd7 osprr37, for extremely early heading date of rice, in the Tohoku region of Japan. The ranges of cultivated crops expanded into various environmental conditions around the world after their domestication. Hokkaido, Japan, lies at the northern limit of cultivation of rice, which originated in the tropics. Novel genotypes for extremely early heading date in Hokkaido are controlled by loss-of-function of both Grain number, plant height and heading date 7 (Ghd7) and Oryza sativa Pseudo-Response Regulator 37 (OsPRR37). We traced genotypes for extremely early heading date and analyzed the phylogeny of rice varieties grown historically in Japan. The mutations in Ghd7 and OsPRR37 had distinct local distributions. Population genomics revealed that varieties collected from the Tohoku region of northern Japan formed three clusters. Mutant alleles of Ghd7 and OsPRR37 appear to have allowed rice cultivation to spread into Hokkaido. Our results show that the mutations of two genes might be occurred in the process of artificial selection during early rice cultivation in the Tohoku region.

Rice functional genomics: decades' efforts and roads ahead

Rice (Oryza sativa L.) is one of the most important crops in the world. Since the completion of rice reference genome sequences, tremendous progress has been achieved in understanding the molecular mechanisms on various rice traits and dissecting the underlying regulatory networks. In this review, we summarize the research progress of rice biology over past decades, including omics, genome-wide association study, phytohormone action, nutrient use, biotic and abiotic responses, photoperiodic flowering, and reproductive development (fertility and sterility). For the roads ahead, cutting-edge technologies such as new genomics methods, high-throughput phenotyping platforms, precise genome-editing tools, environmental microbiome optimization, and synthetic methods will further extend our understanding of unsolved molecular biology questions in rice, and facilitate integrations of the knowledge for agricultural applications.

High-Quality Genomes and High-Density Genetic Map Facilitate the Identification of Genes From a Weedy Rice

Genes have been lost or weakened from cultivated rice during rice domestication and breeding. Weedy rice (Oryza sativa f. spontanea) is usually recognized as the progeny between cultivated rice and wild rice and is also known to harbor an gene pool for rice breeding. Therefore, identifying genes from weedy rice germplasms is an important way to break the bottleneck of rice breeding. To discover genes from weedy rice germplasms, we constructed a genetic map based on w-hole-genome sequencing of a F₂ population derived from the cross between LM8 and a cultivated rice variety. We further identified 31 QTLs associated with 12 important agronomic traits and revealed that ORUFILM03g000095 gene may play an important role in grain length regulation and participate in grain formation. To clarify the genomic characteristics from weedy rice germplasms of LM8, we generated a high-quality genome assembly using single-molecule sequencing, Bionano optical mapping, and Hi-C technologies. The genome harbored a total size of 375.8 Mb, a scaffold N50 of 24.1 Mb, and originated approximately 0.32 million years ago (Mya) and was more closely related to Oryza sativa ssp. japonica. and contained 672 unique genes. It is related to the formation of grain shape, heading date and tillering. This study generated a high-quality reference genome of weedy rice and high-density genetic map that would benefit the analysis of genome evolution for related species and suggested an effective way to identify genes related to important agronomic traits for further rice breeding.

Genomic divergence during feralization reveals both conserved and distinct mechanisms of parallel weediness evolution

Agricultural weeds are the most important biotic constraints to global crop production, and chief among these is weedy rice. Despite increasing yield losses from weedy rice in recent years worldwide, the genetic basis of weediness evolution remains unclear. Using whole-genome sequence analyses, we examined the origins and adaptation of Japanese weedy rice. We find evidence for a weed origin from tropical japonica crop ancestry, which has not previously been documented in surveys of weedy rice worldwide. We further show that adaptation occurs largely through different genetic mechanisms between independently-evolved temperate japonica- and tropical japonica-derived strains; most genomic signatures of positive selection are unique within weed types. In addition, some weedy rice strains have evolved through hybridization between weedy and cultivated rice with adaptive introgression from the crop. Surprisingly, introgression from cultivated rice confers not only crop-like adaptive traits (such as shorter plant height, facilitating crop mimicry) but also weedy-like traits (such as seed dormancy). These findings reveal how hybridization with cultivated rice can promote persistence and proliferation of weedy rice.

Comparative multi-omics analysis of hypoxic germination tolerance in weedy rice embryos and coleoptiles

Hypoxic germination tolerance is an important trait for seedling establishment of direct-seeded rice. Our comparative metabolomics analysis revealed that weedy rice accumulated more sugar and amino acids than cultivated rice accumulated in the embryo and coleoptile tissues under hypoxic stress. At the transcriptional level, oxidative phosphorylation activity in weedy rice was higher than in cultivated rice that likely led to more efficient energy metabolism during hypoxic stress. Based on our comparative proteomics analysis, enriched proteins related to cell wall implied that the advantages in energy metabolism of weedy rice were ultimately reflected in the formation of tissue structures. In this study, we found that most of key hypoxic germination tolerance (HGT) genes shared the same genetic backgrounds with Oryza japonica, however, several of them genetically similar to other Oryza plant also play important roles. Our findings suggest weedy rice can serve as genetic resources for the improvement of direct-seeding rice.

Uncovering phylogenetic relationships and genetic diversity of water dropwort using phenotypic traits and SNP markers

The water dropworts Oenanthe linearis Wall. ex DC. and O. javanica (Blume) DC. are aquatic perennial herbs that have been used in China as vegetables and traditional medicines. However, their phylogenetic relationships and genetic diversity are poorly understood. Here, we presented the phenotypic traits and genome-wide DNA marker-based analysis of 158 water dropwort accessions representing both species. The analysis revealed that Oenanthe linearis was readily segregated into linear-leaf and deep-cleft leaf water dropworts according to their leaf shapes at flowering. Oenanthe javanica was classified by clustering analysis into two clusters based mainly on the morphological characteristics of their ultimate segments (leaflets). A set of 11 493 high-quality single-nucleotide polymorphisms was identified and used to construct a phylogenetic tree. There was strong discrimination between O. linearis and O. javanica, which was consistent with their phenotype diversification. The population structure and phylogenetic tree analyses suggested that the O. linearis accessions formed two major groups, corresponding to the linear-leaf and deep-cleft leaf types. The most obvious phenotypic differences between them were fully expressed at the reproductive growth stage. A single-nucleotide polymorphism-based analysis revealed that the O. javanica accessions could be categorized into groups I andII. However, this finding did not entirely align with the clusters revealed by morphological classification. Landraces were clustered into one group along with the remaining wild accessions. Hence, water dropwort domestication was short in duration. The level of genetic diversity for O. linearis (π = 0.1902) was slightly lower than that which was estimated for O. javanica (π = 0.2174). There was a low level of genetic differentiation between O. linearis and O. javanica (Fst = 0.0471). The mean genetic diversity among accessions ranged from 0.1818 for the linear-leaf types to 0.2318 for the groupII accessions. The phenotypic traits and the single-nucleotide polymorphism markers identified here lay empirical foundation for future genomic studies on water dropwort.

De-Domestication: An Extension of Crop Evolution

De-domestication or feralization is an interesting phenomenon in crops and livestock. Previously, evidence for crop de-domestication was based mainly on studies using phenotypic and genotypic data from limited molecular markers or gene segments. Recent genomic studies in rice, barley, and wheat provide comprehensive landscapes of de-domestication on a whole-genome scale. Here, we summarize crop de-domestication processes, ecological roles of de-domesticates, mechanisms underlying crop de-domestication syndromes, and conditions potentially favoring de-domestication events. We further explain how recent de-domestication studies have expanded our understanding of the complexity of crop evolution, and highlight the genetic novelties of de-domesticates beneficial for modern crop breeding.

Explore the genetics of weedy traits using rice 3K database

BACKGROUND

Weedy rice, a conspecific weedy counterpart of the cultivated rice (Oryza sativa L.), has been problematic in rice-production area worldwide. Although we started to know about the origin of some weedy traits for some rice-growing regions, an overall assessment of weedy trait-related loci was not yet available. On the other hand, the advances in sequencing technologies, together with community efforts, have made publicly available a large amount of genomic data. Given the availability of public data and the need of "weedy" allele mining for a better management of weedy rice, the objective of the present study was to explore the genetic architecture of weedy traits based on publicly available data, mainly from the 3000 Rice Genome Project (3K-RGP).

RESULTS

Based on the results of population structure analysis, we have selected 1378 individuals from four sub-populations (aus, indica, temperate japonica, tropical japonica) without admixed genomic composition for genome-wide association analysis (GWAS). Five traits were investigated: awn color, seed shattering, seed threshability, seed coat color, and seedling height. GWAS was conducted for each sub-population × trait combination and we have identified 66 population-specific trait-associated SNPs. Eleven significant SNPs fell into an annotated gene and four other SNPs were close to a putative candidate gene (± 25 kb). SNPs located in or close to Rc were particularly predictive of the occurrence of seed coat color and our results showed that different sub-populations required different SNPs for a better seed coat color prediction. We compared the data of 3K-RGP to a publicly available weedy rice dataset. The profile of allele frequency, phenotype-genotype segregation of target SNP, as well as GWAS results for the presence and absence of awns diverged between the two sets of data.

CONCLUSIONS

The genotype of trait-associated SNPs identified in this study, especially those located in or close to Rc, can be developed to diagnostic SNPs to trace the origin of weedy trait occurred in the field. The difference of results from the two publicly available datasets used in this study emphasized the importance of laboratory experiments to confirm the allele mining results based on publicly available data.

Studies of rice Hd1 haplotypes worldwide reveal adaptation of flowering time to different environments

Rice domestication/adaptation is a good model for studies of the development and spread of this important crop. Mutations that caused morphological and physiological change, followed by human selection/expansion, finally led to the improvement of phenotypes suitable for different kinds of environments. We used the sequence information for Heading date 1 (Hd1) gene to reveal the association between sequence changes and flowering phenotypes of rice in different regions. Seven loss-of-function hd1 haplotypes had been reported. By data-mining the genome sequencing information in the public domain, we discovered 3 other types. These loss-of-function allele haplotypes are present in subtropical and tropical regions, which indicates human selection. Some of these haplotypes are present locally. However, types 7 and 13 are present in more than one-third of the world's rice accessions, including landraces and modern varieties. In the present study, phylogenetic, allele network and selection pressure analyses revealed that these two haplotypes might have occurred early in Southeastern Asia and then were introgressed in many local landraces in nearby regions. We also demonstrate that these haplotypes are present in weedy rice populations, which again indicates that these alleles were present in rice cultivation for long time. In comparing the wild rice sequence information, these loss-of-function haplotypes occurred in agro but were not from wild rice.

Whole-genome sequencing reveals uniqueness of black-hulled and straw-hulled weedy rice genomes

Both SH and BHA weedy rice genotypes evolved independently and have distinct genomic composition. Different genetic mechanisms may be responsible for their competitiveness and adaptation to diverse environmental conditions. Two major types of weedy rice are recognized in the USA based on morphology: straw-hull (SH) and black-hull awned (BHA) weedy rice. We performed whole-genome resequencing of a SH weedy rice 'PSRR-1', a BHA weedy rice 'BHA1115', and a japonica cultivar 'Cypress' to delineate genome-wide differences and their relevance to genetics and evolution of weedy attributes. The high-quality reads were uniformly distributed with 82-88% genome coverage. The number of genotype-specific SNPs and InDels was highest in Cypress, followed by BHA1115 and PSRR-1. However, more genes were affected in BHA1115 compared with other two genotypes which is evident from the number of high-impact SNPs and InDels. Haplotype analysis of selected genes involved in domestication, adaptation, and agronomic performance not only differentiated SH from BHA weedy rice and supported evolution of weedy rice through de-domestication, but also validated the function of several genes such as qAn-1, qAn-2, Bh4, Rc, SD1, OsLG1, and OsC1. Several candidate genes were identified for previously reported seed dormancy and seed shattering QTLs. The SH and BHA weedy rice have distinct genomic composition, and the BHA weedy rice likely diverged earlier than SH weedy rice. The accumulation of plant development, reproduction, and defense-related genes in weedy rice possibly helped them to compete, survive, and spread under a wide range of environmental conditions by employing novel and diverse mechanisms. The genomic resources will be useful for both weed management and rice improvement by exploring the molecular basis of key agronomic, adaptive, and domestication attributes.

Weedy Rice as a Novel Gene Resource: A Genome-Wide Association Study of Anthocyanin Biosynthesis and an Evaluation of Nutritional Quality

The pericarp color of rice grains is an important agronomic trait affected by domestication, and the color pigment, anthocyanin, is one of the key determinants of rice nutritional quality. Weedy rice, also called red rice because its pericarp is often red, may be a novel gene resource for the development of new rice. However, the genetic basis and nutritional quality of anthocyanin are poorly known. In this study, we used a genome-wide association study (GWAS) to find novel and specific QTLs of red pericarp in weedy rice. The known key gene site of red pericarp Rc was detected as the common genetic basis of both weedy and cultivated rice, and another 13 associated signals of pericarp color that were identified may contribute specifically to weedy rice pericarp color. We then nominated three pericarp color genes that may contribute to weedy rice divergence from cultivated rice based on selection sweep analysis. After clarifying the distribution and growth dynamics of pigment in weedy rice caryopsis, we compared its nutritional quality with cultivated rice. We found that sampled weedy rice pericarps had much greater quantities of anthocyanin, beneficial trace elements, free amino acids, and unsaturated fatty acids than the cultivated rice. In conclusion, the gene resources and novel genetic systems of rice anthocyanin biosynthesis explored in this study are of great value for the development of nutritious, high anthocyanin content rice.

Using demographic model selection to untangle allopatric divergence and diversification mechanisms in the Rheum palmatum complex in the Eastern Asiatic Region

Allopatric divergence is often initiated by geological uplift and restriction to sky-islands, climate oscillations, or river capture. However, it can be difficult to establish which mechanism was the most likely to generate the current phylogeographical structure of a species. Recently, genomic data in conjunction with a model testing framework have been applied to address this issue in animals. To test whether such an approach is also likely to be successful in plants, we used population genomic data of the Rheum palmatum complex from the Eastern Asiatic Region, in conjunction with biogeographical reconstruction and demographic model selection, to identify the potential mechanism(s) which have led to the current level of divergence. Our results indicate that the R. palmatum complex originated in the central Hengduan Mts and possibly in regions further east, and then dispersed westward and eastward resulting in genetically distinct lineages. Populations are likely to have diverged in refugia during climate oscillations followed by subsequent expansion and secondary contact. However, model simulations within the western lineage of the R. palmatum complex cannot reject a restriction to sky-islands as a possible mechanism of diversification due to the genetically ambiguous position of one population. This highlights that genetically mixed populations might introduce ambiguity regarding the best diversification model in some cases. Although it might be possible to resolve this ambiguity using other data, sometimes this could prove to be difficult in complex biogeographical areas.

Diverse genetic mechanisms underlie worldwide convergent rice feralization

Background

Worldwide feralization of crop species into agricultural weeds threatens global food security. Weedy rice is a feral form of rice that infests paddies worldwide and aggressively outcompetes cultivated varieties. Despite increasing attention in recent years, a comprehensive understanding of the origins of weedy crop relatives and how a universal feralization process acts at the genomic and molecular level to allow the rapid adaptation to weediness are still yet to be explored.

Results

We use whole-genome sequencing to examine the origin and adaptation of 524 global weedy rice samples representing all major regions of rice cultivation. Weed populations have evolved multiple times from cultivated rice, and a strikingly high proportion of contemporary Asian weed strains can be traced to a few Green Revolution cultivars that were widely grown in the late twentieth century. Latin American weedy rice stands out in having originated through extensive hybridization. Selection scans indicate that most genomic regions underlying weedy adaptations do not overlap with domestication targets of selection, suggesting that feralization occurs largely through changes at loci unrelated to domestication.

Conclusions

This is the first investigation to provide detailed genomic characterizations of weedy rice on a global scale, and the results reveal diverse genetic mechanisms underlying worldwide convergent rice feralization.

Electronic supplementary material

The online version of this article (10.1186/s13059-020-01980-x) contains supplementary material, which is available to authorized users.

Exploitation of Hi-C sequencing for improvement of genome assembly and in-vitro validation of differentially expressing genes in Jatropha curcas L

Jatropha curcas is one of the major sources of renewable energy due to potential use of its oil as a biofuel. The genome of this crop is constituted by the high content of repetitive elements. We employed the Hi-C proximity ligation technique to re-scaffold our existing hybrid genome assembly of an elite genotype (RJC1) developed using Illumina and Pacbio technologies. We assembled 99.81% of non-truncated reads to achieve 266.80 Mbp of the genome with an N50 value of 1.58 Mb. Furthermore, we compared the efficiency of Hi-C-augmented genome assembly with the hybrid genome assembly and observed a ~ 50% reduction in scaffolds and a tenfold increase in the N50 value. The gene ontology analysis revealed the identification of terms for molecular function (45.52%), cellular component (33.47%), and biological function (20.99%). Comparative genomic analysis of 13-plant species showed the conservation of 414 lipid metabolizing genes identified in the KEGG pathway analysis. Differential gene expression (DGE) studies were conducted in the healthy and Jatropha mosaic virus-infected leaves via RNA-seq analysis and observed gene expression changes for 2185 genes. Out of these, we observed 546 genes having more than two-fold change of transcript level and among these 259 genes were down-regulated and 287 genes were up-regulated. To validate RNA-seq data, two DEGs were selected for gene expression analysis using qRT-PCR and the data was in correlation with in silico results. RNA-seq analysis further shows the identification of some of the candidate genes and may be useful to develop JMV resistant plants after functional validation. This Hi-C genome assembly provides a detailed accurate reference genome which could be utilized to improve Jatropha and other economically important Euphorbiaceae family members.

Quantitative Trait Locus (QTLs) Mapping for Quality Traits of Wheat Based on High Density Genetic Map Combined With Bulked Segregant Analysis RNA-seq (BSR-Seq) Indicates That the Basic 7S Globulin Gene Is Related to Falling Number

Numerous quantitative trait loci (QTLs) have been identified for wheat quality; however, most are confined to low-density genetic maps. In this study, based on specific-locus amplified fragment sequencing (SLAF-seq), a high-density genetic map was constructed with 193 recombinant inbred lines derived from Chuanmai 42 and Chuanmai 39. In total, 30 QTLs with phenotypic variance explained (PVE) up to 47.99% were identified for falling number (FN), grain protein content (GPC), grain hardness (GH), and starch pasting properties across three environments. Five NAM genes closely adjacent to QGPC.cib-4A probably have effects on GPC. QGH.cib-5D was the only one detected for GH with high PVE of 33.31-47.99% across the three environments and was assumed to be related to the nearest pina-D1 and pinb-D1genes. Three QTLs were identified for FN in at least two environments, of which QFN.cib-3D had relatively higher PVE of 16.58-25.74%. The positive effect of QFN.cib-3D for high FN was verified in a double-haploid population derived from Chuanmai 42 × Kechengmai 4. The combination of these QTLs has a considerable effect on increasing FN. The transcript levels of Basic 7S globulin and Basic 7S globulin 2 in QFN.cib-3D were significantly different between low FN and high FN bulks, as observed through bulk segregant RNA-seq (BSR). These QTLs and candidate genes based on the high-density genetic map would be beneficial for further understanding of the genetic mechanism of quality traits and molecular breeding of wheat.

Analysis of evolutionary relationships provides new clues to the origins of weedy rice

Weedy rice (WR) (Oryza sativa f. spontanea) is considered to be a pest in modern rice production systems because it competes for resources, has poor yield characteristics, and subsequently has a negative effect on rice grain yield. The evolutionary relationships among WR, landrace rice (LR), improved rice (IR) cultivars, and wild rice are largely unknown. In this study, we conducted a population genetic analysis based on neutral markers and gene haplotypes in 524 rice accessions and a comparative transcriptomic analysis using 15 representative samples. The results showed that WR populations have the highest level of genetic diversity (H e = 0.8386) and can be divided into two groups (japonica-type and indica-type). The japonica-type WR accessions from Heilongjiang province (HLJ), Jilin province (JL), Liaoning province (LN), and NX provinces clustered with the landraces grown in these same provinces. The indica-types from Jiangsu province (JS) also clustered with the indica-type landraces from JS province. Comparative transcriptome analysis of WR' IR and LR from HLJ, JL, and LN provinces showed that the WR still clustered with the LR, and that the IR lines comprise a single population. Thirty-two differentially expressed genes were shared by the IR and LR groups as well as between the IR and WR groups. Using Gene ontology (GO) analysis, we identified 19 shared GO terms in the IR and LR groups as well as between the IR and WR groups. Our results suggest that WR populations in China have diverse origins, and comparative transcriptome analysis of different types of rice from HLJ, JL, and LN provinces suggests that IR populations have become a end point in the evolution of WR, which provides a new perspective for the study of WR origins and lays a solid foundation for rice breeding.

Getting Back to Nature: Feralization in Animals and Plants

Formerly domesticated organisms and artificially selected genes often escape controlled cultivation, but their subsequent evolution is not well studied. In this review, we examine plant and animal feralization through an evolutionary lens, including how natural selection, artificial selection, and gene flow shape feral genomes, traits, and fitness. Available evidence shows that feralization is not a mere reversal of domestication. Instead, it is shaped by the varied and complex histories of feral populations, and by novel selection pressures. To stimulate further insight we outline several future directions. These include testing how 'domestication genes' act in wild settings, studying the brains and behaviors of feral animals, and comparative analyses of feral populations and taxa. This work offers feasible and exciting research opportunities with both theoretical and practical applications.