A Principal Components Analysis and Functional Annotation of Differentially Expressed Genes in Brain Regions of Gray Rats Selected for Tame or Aggressive Behavior.

The process of domestication, despite its short duration as it compared with the time scale of the natural evolutionary process, has caused rapid and substantial changes in the phenotype of domestic animal species. Nonetheless, the genetic mechanisms underlying these changes remain poorly understood. The present study deals with an analysis of the transcriptomes from four brain regions of gray rats (Rattus norvegicus), serving as an experimental model object of domestication. We compared gene expression profiles in the hypothalamus, hippocampus, periaqueductal gray matter, and the midbrain tegmental region between tame domesticated and aggressive gray rats and revealed subdivisions of differentially expressed genes by principal components analysis that explain the main part of differentially gene expression variance. Functional analysis (in the DAVID (Database for Annotation, Visualization and Integrated Discovery) Bioinformatics Resources database) of the differentially expressed genes allowed us to identify and describe the key biological processes that can participate in the formation of the different behavioral patterns seen in the two groups of gray rats. Using the STRING- DB (search tool for recurring instances of neighboring genes) web service, we built a gene association network. The genes engaged in broad network interactions have been identified. Our study offers data on the genes whose expression levels change in response to artificial selection for behavior during animal domestication.

Behavioural synchronisation between different groups of dogs and wolves and their owners/handlers: Exploring the effect of breed and human interaction.

Dogs have previously been shown to synchronise their behaviour with their owner and the aim of this study was to test the effect of immediate interactions, breed, and the effects of domestication. The behavioural synchronisation test was conducted in outdoor enclosures and consisted of 30 s where the owner/handler was walking and 30 s of standing still. Three studies were conducted to explore the effect of immediate interaction (study A), the effect of breed group (study B), and the effect of domestication (study C). In study A, a group of twenty companion dogs of various breeds were tested after three different human interaction treatments: Ignore, Pet, and Play. The results showed that dogs adjusted their movement pattern to align with their owner's actions regardless of treatment. Furthermore, exploration, eye contact, and movement were all influenced by the owners moving pattern, and exploration also decreased after the Play treatment. In study B, the synchronisation test was performed after the Ignore treatment on three groups: 24 dogs of ancient dog breeds, 17 solitary hunting dogs, and 20 companion dogs (data from study A). Irrespective of the group, all dogs synchronised their moving behaviour with their owner. In addition, human walking positively influenced eye contact behaviour while simultaneously decreasing exploration behaviour. In study C, a group of six socialised pack-living wolves and six similarly socialised pack-living dogs were tested after the Ignore treatment. Interestingly, these animals did not alter their moving behaviour in response to their handler. In conclusion, dogs living together with humans synchronise with their owner's moving behaviour, while wolves and dogs living in packs do not. Hence, the degree of interspecies behavioural synchronisation may be influenced by the extent to which the dogs are immersed in everyday life with humans.

Domestication constrains the ability of dogs to convey emotions via facial expressions in comparison to their wolf ancestors.

Dogs (Canis lupus familiaris) are the domestically bred descendant of wolves (Canis lupus). However, selective breeding has profoundly altered facial morphologies of dogs compared to their wolf ancestors. We demonstrate that these morphological differences limit the abilities of dogs to successfully produce the same affective facial expressions as wolves. We decoded facial movements of captive wolves during social interactions involving nine separate affective states. We used linear discriminant analyses to predict affective states based on combinations of facial movements. The resulting confusion matrix demonstrates that specific combinations of facial movements predict nine distinct affective states in wolves; the first assessment of this many affective facial expressions in wolves. However, comparative analyses with kennelled rescue dogs revealed reduced ability to predict affective states. Critically, there was a very low predictive power for specific affective states, with confusion occurring between negative and positive states, such as Friendly and Fear. We show that the varying facial morphologies of dogs (specifically non-wolf-like morphologies) limit their ability to produce the same range of affective facial expressions as wolves. Confusion among positive and negative states could be detrimental to human-dog interactions, although our analyses also suggest dogs likely use vocalisations to compensate for limitations in facial communication.

Low frequency of the wild-type freezing-tolerance LsCBF7 allele among lettuce population suggests a negative selection during domestication and breeding.

KEY MESSAGE: Sustainable winter production in lettuce requires freezing tolerant varieties. This study identified a wild-type allele of LsCBF7 that could contribute to freezing tolerance improvement in lettuce. Lettuce is one of the most consumed vegetables globally. While ideally grown in 13-21 °C, its cultivation extends into winter in milder climates. However, occasional freezing temperatures can significantly reduce yields. Therefore, the development of freezing-tolerant lettuce varieties has become a long-term goal of lettuce breeding programs. Despite its significance, our understanding of freezing tolerance in lettuce remains limited. Plants have evolved a coping mechanism against freezing, known as cold acclimation, whereby they can increase freezing tolerance when pre-exposed to low nonfreezing temperatures. The CBF pathway is well-known for its central role in cold acclimation. Previously, we identified 14 CBF genes in lettuce and discovered that one of them, LsCBF7, had a loss-of-function mutation. In this study, we uncovered that accessions from colder regions carried the wild-type allele of LsCBF7 and this allele likely contributed to increased freezing tolerance, with 14% of the lettuce population carrying this allele. Interestingly, in wild lettuce (L. serriola) that is considered a progenitor of cultivated lettuce, this wild-type allele was much more common, with a frequency of 90%. This finding suggests that this wild-type allele may have undergone negative selection during the domestication or breeding of lettuce. Our data strongly indicate that this allele could be linked to early bolting, an undesirable trait in lettuce, which may have driven the negative selection. While this wild-type allele shows promise for improving freezing tolerance in lettuce, it is crucial to decouple it from the early bolting trait to fully harness its potential in lettuce breeding.

Genome-wide detections for runs of homozygosity and selective signatures reveal novel candidate genes under domestication in chickens.

BACKGROUND: Indigenous chickens were developed through a combination of natural and artificial selection; essentially, changes in genomes led to the formation of these modern breeds via admixture events. However, their confusing genetic backgrounds include a genomic footprint regulating complex traits, which is not conducive to modern animal breeding. RESULTS: To better evaluate the candidate regions under domestication in indigenous chickens, we considered both runs of homozygosity (ROHs) and selective signatures in 13 indigenous chickens. The genomes of Silkie feather chickens presented the highest heterozygosity, whereas the highest inbreeding status and ROH number were found in Luhua chickens. Short ROH (< 1 Mb), were the principal type in all chickens. A total of 291 ROH islands were detected, and QTLdb mapping results indicated that body weight and carcass traits were the most important traits. An ROH on chromosome 2 covering VSTM2A gene was detected in 12 populations. Combined analysis with the Tajima's D index revealed that 18 genes (e.g., VSTM2A, BBOX1, and RYR2) were under selection and covered by ROH islands. Transcriptional analysis results showed that RYR2 and BBOX1 were specifically expressed in the heart and muscle tissue, respectively. CONCLUSION: Based on genome-wide scanning for ROH and selective signatures, we evaluated the genomic characteristics and detected significant candidate genes covered by ROH islands and selective signatures. The findings in this study facilitated the understanding of genetic diversity and provided valuable insights for chicken breeding and conservation strategies.

Archaeobotanical evidence supports indigenous cucurbit long-term use in the Mesoamerican Neotropics.

The squash family (Cucurbitaceae) contains some of the most important crops cultivated worldwide and has played an important ecological, economic, and cultural role for millennia. In the American tropics, squashes were among the first cultivated crop species, but little is known about how their domestication unfolded. Here, we employ direct radiocarbon dating and morphological analyses of desiccated cucurbit seeds, rinds, and stems from El Gigante Rockshelter in Honduras to reconstruct human practices of selection and cultivation of Lagenaria siceraria, Cucurbita pepo, and Cucurbita moschata. Direct radiocarbon dating indicates that humans started using Lagenaria and wild Cucurbita starting ~ 10,950 calendar years before present (cal B.P.), primarily as watertight vessels and possibly as cooking and drinking containers. A rind directly dated to 11,150-10,765 cal B.P. represents the oldest known bottle gourd in the Americas. Domesticated C. moschata subsequently appeared ~ 4035 cal B.P., followed by domesticated C. pepo ~ 2190 cal B.P. associated with increasing evidence for their use as food crops. Multivariate statistical analysis of seed size and shape show that the archaeological C. pepo assemblage exhibits significant variability, representing at least three varieties: one similar to present-day zucchini, another like present-day vegetable marrow, and a native cultivar without modern analogs. Our archaeobotanical data supports the hypothesis that Indigenous cucurbit use started in the Early Holocene, and that agricultural complexity during the Late Holocene involved selective breeding that encouraged crop diversification.

Chorematic modeling to represent dynamics in the quinoa agroecosystems in Peru.

Our research occurred in the Andean region, one of the eight global centers of domestication of plant species grown for agriculture. The shores of Lake Titicaca (located between Peru and Bolivia), at 3800 meters above sea level, are recognized as the center of origin of quinoa (Chenopodium quinoa Willd.). In this region, complex societies have emerged, thanks to the development of water and soil management technologies. They have managed to overcome high mountain territories' extreme and variable climatic conditions. These societies have traditionally protected and preserved quinoa as food for present and future generations through their long-standing knowledge and cultivation practices. The fieldwork occurred in the context of Andean family farming, and our study analyzes nature-society dynamics with a chorematic approach and interviews with local communities. The interest of this work is the transformation of the landscape at the scale of the mountain agroecosystem to understand better the impacts of rural development policies. Chorematic modeling was applied to two periods, before and after 1970, a pivotal year in Peru for agriculture, to show how socio-spatial dynamics in the Andean environment are changing, particularly concerning the evolution of quinoa cultivation. The results show that wild quinoa relatives' distribution is strongly linked to the socio-spatial organization of the agroecosystem. Different species of wild quinoa relatives are maintained by villagers for their multiple foods, medicinal and cultural uses in natural areas, grazed areas, on edge, and also within cultivated fields. However, this management is changing under the pressure of global issues related to the international quinoa market, whose requirements imply reducing the presence of wild relatives in cultivated fields.

Plant genome information facilitates plant functional genomics.

MAIN CONCLUSION: In this review, we give an overview of plant sequencing efforts and how this impacts plant functional genomics research. Plant genome sequence information greatly facilitates the studies of plant biology, functional genomics, evolution of genomes and genes, domestication processes, phylogenetic relationships, among many others. More than two decades of sequencing efforts have boosted the number of available sequenced plant genomes. The first plant genome, of Arabidopsis, was published in the year 2000 and currently, 4604 plant genomes from 1482 plant species have been published. Various large sequence initiatives are running, which are planning to produce tens of thousands of sequenced plant genomes in the near future. In this review, we give an overview on the status of sequenced plant genomes and on the use of genome information in different research areas.

Differential selection of yield and quality traits has shaped genomic signatures of cowpea domestication and improvement.

Cowpeas (tropical legumes) are important in ensuring food and nutritional security in developing countries, especially in sub-Saharan Africa. Herein, we report two high-quality genome assemblies of grain and vegetable cowpeas and we re-sequenced 344 accessions to characterize the genomic variations landscape. We identified 39 loci for ten important agronomic traits and more than 541 potential loci that underwent selection during cowpea domestication and improvement. In particular, the synchronous selections of the pod-shattering loci and their neighboring stress-relevant loci probably led to the enhancement of pod-shattering resistance and the compromise of stress resistance during the domestication from grain to vegetable cowpeas. Moreover, differential selections on multiple loci associated with pod length, grain number per pod, seed weight, pod and seed soluble sugars, and seed crude proteins shaped the yield and quality diversity in cowpeas. Our findings provide genomic insights into cowpea domestication and improvement footprints, enabling further genome-informed cultivar improvement of cowpeas.

A distal enhancer guides the negative selection of toxic glycoalkaloids during tomato domestication.

Steroidal glycoalkaloids (SGAs) are major plant defense metabolites against pests, while they are considered poisonous in food. The genetic basis that guides negative selection of SGAs production during tomato domestication remains poorly understood. Here, we identify a distal enhancer, GAME Enhancer 1 (GE1), as the key regulator of SGAs metabolism in tomato. GE1 recruits MYC2-GAME9 transcriptional complex to regulate the expression of GAME cluster genes via the formation of chromatin loops located in the neighboring DNA region. A naturally occurring GE176 allelic variant is found to be more active in stimulating GAME expression. We show that the weaker GE1 allele has been the main driver for selecting reduced SGAs levels during tomato domestication. Unravelling the "TFs-Enhancer-Promoter" regulatory mechanism operating in SGAs metabolism opens unprecedented prospects for SGAs manipulation in Solanaceae via precision breeding strategies.

Increased DNA methylation contributes to the early ripening of pear fruits during domestication and improvement.

BACKGROUND: DNA methylation is an essential epigenetic modification. However, its contribution to trait changes and diversity in the domestication of perennial fruit trees remains unknown. RESULTS: Here, we investigate the variation in DNA methylation during pear domestication and improvement using whole-genome bisulfite sequencing in 41 pear accessions. Contrary to the significant decrease during rice domestication, we detect a global increase in DNA methylation during pear domestication and improvement. We find this specific increase in pear is significantly correlated with the downregulation of Demeter-like1 (DML1, encoding DNA demethylase) due to human selection. We identify a total of 5591 differentially methylated regions (DMRs). Methylation in the CG and CHG contexts undergoes co-evolution during pear domestication and improvement. DMRs have higher genetic diversity than selection sweep regions, especially in the introns. Approximately 97% of DMRs are not associated with any SNPs, and these DMRs are associated with starch and sucrose metabolism and phenylpropanoid biosynthesis. We also perform correlation analysis between DNA methylation and gene expression. We find genes close to the hypermethylated DMRs that are significantly associated with fruit ripening. We further verify the function of a hyper-DMR-associated gene, CAMTA2, and demonstrate that overexpression of CAMTA2 in tomato and pear callus inhibits fruit ripening. CONCLUSIONS: Our study describes a specific pattern of DNA methylation in the domestication and improvement of a perennial pear tree and suggests that increased DNA methylation plays an essential role in the early ripening of pear fruits.

The Evolution, Expression Patterns, and Domestication Selection Analysis of the Annexin Gene Family in the Barley Pan-Genome.

Plant annexins constitute a conserved protein family that plays crucial roles in regulating plant growth and development, as well as in responses to both biotic and abiotic stresses. In this study, a total of 144 annexin genes were identified in the barley pan-genome, comprising 12 reference genomes, including cultivated barley, landraces, and wild barley. Their chromosomal locations, physical-chemical characteristics, gene structures, conserved domains, and subcellular localizations were systematically analyzed to reveal the certain differences between wild and cultivated populations. Through a cis-acting element analysis, co-expression network, and large-scale transcriptome analysis, their involvement in growth, development, and responses to various stressors was highlighted. It is worth noting that HvMOREXann5 is only expressed in pistils and anthers, indicating its crucial role in reproductive development. Based on the resequencing data from 282 barley accessions worldwide, genetic variations in thefamily were investigated, and the results showed that 5 out of the 12 identified HvMOREXanns were affected by selection pressure. Genetic diversity and haplotype frequency showed notable reductions between wild and domesticated barley, suggesting that a genetic bottleneck occurred on the annexin family during the barley domestication process. Finally, qRT-PCR analysis confirmed the up-regulation of HvMOREXann7 under drought stress, along with significant differences between wild accessions and varieties. This study provides some insights into the genome organization and genetic characteristics of the annexin gene family in barley at the pan-genome level, which will contribute to better understanding its evolution and function in barley and other crops.

A Past Genetic Bottleneck from Argentine Beans and a Selective Sweep Led to the Race Chile of the Common Bean (Phaseolus vulgaris L.).

The domestication process of the common bean gave rise to six different races which come from the two ancestral genetic pools, the Mesoamerican (Durango, Jalisco, and Mesoamerica races) and the Andean (New Granada, Peru, and Chile races). In this study, a collection of 281 common bean landraces from Chile was analyzed using a 12K-SNP microarray. Additionally, 401 accessions representing the rest of the five common bean races were analyzed. A total of 2543 SNPs allowed us to differentiate a genetic group of 165 accessions that corresponds to the race Chile, 90 of which were classified as pure accessions, such as the bean types 'Tórtola', 'Sapito', 'Coscorrón', and 'Frutilla'. Our genetic analysis indicates that the race Chile has a close relationship with accessions from Argentina, suggesting that nomadic ancestral peoples introduced the bean seed to Chile. Previous archaeological and genetic studies support this hypothesis. Additionally, the low genetic diversity (π = 0.053; uHe = 0.53) and the negative value of Tajima' D (D = -1.371) indicate that the race Chile suffered a bottleneck and a selective sweep after its introduction, supporting the hypothesis that a small group of Argentine bean genotypes led to the race Chile. A total of 235 genes were identified within haplotype blocks detected exclusively in the race Chile, most of them involved in signal transduction, supporting the hypothesis that intracellular signaling pathways play a fundamental role in the adaptation of organisms to changes in the environment. To date, our findings are the most complete investigation associated with the origin of the race Chile of common bean.

Resequencing Analyses Revealed Genetic Diversity and Selection Signatures during Rabbit Breeding and Improvement.

Domestication has shaped the diverse characteristics of rabbits, including coat color, fur structure, body size, and various physiological traits. Utilizing whole-genome resequencing (DNBSEQ-T7), we analyzed the genetic diversity, population structure, and genomic selection across 180 rabbits from 17 distinct breeds to uncover the genetic basis of these traits. We conducted whole-genome sequencing on 17 rabbit breeds, identifying 17,430,184 high-quality SNPs and analyzing genomic diversity, patterns of genomic variation, population structure, and selection signatures related to coat color, coat structure, long hair, body size, reproductive capacity, and disease resistance. Through PCA and NJ tree analyses, distinct clusters emerged among Chinese indigenous rabbits, suggesting varied origins and domestication histories. Selective sweep testing pinpointed regions and genes linked to domestication and key morphological and economic traits, including those affecting coat color (TYR, ASIP), structure (LIPH), body size (INSIG2, GLI3), fertility (EDNRA, SRD5A2), heat stress adaptation (PLCB1), and immune response (SEC31A, CD86, LAP3). Our study identified key genomic signatures of selection related to traits such as coat color, fur structure, body size, and fertility; these findings highlight the genetic basis underlying phenotypic diversification in rabbits and have implications for breeding programs aiming to improve productive, reproductive, and adaptive traits. The detected genomic signatures of selection also provide insights into rabbit domestication and can aid conservation efforts for indigenous breeds.

The Effect of Short-Term Artificial Feed Domestication on the Expression of Oxidative-Stress-Related Genes and Antioxidant Capacity in the Liver and Gill Tissues of Mandarin Fish (Siniperca chuatsi).

To investigate whether Mandarin fish developed oxidative stress after being domesticated with artificial feed, we conducted a series of experiments. Oxidative stress is an important factor leading to diseases and aging in the body. The liver integrates functions such as digestion, metabolism, detoxification, coagulation, and immune regulation, while the gills are important respiratory organs that are sensitive to changes in the water environment. Therefore, we used the liver and gills of Mandarin fish as research materials. The aim of this study was to investigate the effects of short-term artificial feed domestication on the expression of oxidative stress genes and the changes in oxidative-stress-related enzyme activity in the liver and gills of Mandarin fish. We divided the Mandarin fish into two groups for treatment. The control group was fed with live bait continuously for 14 days, while the experimental group was fed with half artificial feed and half live bait from 0 to 7 days (T-7 d), followed by solely artificial feed from 7 to 14 days (T-14 d). The experimental results showed that there was no difference in the body weight, length, and standard growth rate of the Mandarin fish between the two groups of treatments; after two treatments, there were differences in the expression of genes related to oxidative stress in the gills (keap1, kappa, gsta, gstt1, gstk1, SOD, and CAT) and in the liver (GPx, keap1, kappa, gsta, gstt1, gr, and SOD). In the liver, GPx activity and the content of MDA were significantly upregulated after 7 days of domestication, while in the gills, SOD activity was significantly upregulated after 7 days of domestication and GPx activity was significantly downregulated after 14 days of domestication. These results suggest that artificial feed domestication is associated with oxidative stress. Moreover, these results provide experimental basic data for increasing the production of aquaculture feed for Mandarin fish.

Tomato domestication rather than subsequent breeding events reduces microbial associations related to phosphorus recovery.

Legacy phosphorus (P) is a reservoir of sparingly available P, and its recovery could enhance sustainable use of nonrenewable mineral fertilizers. Domestication has affected P acquisition, but it is unknown if subsequent breeding efforts, like the Green Revolution (GR), had a similar effect. We examined how domestication and breeding events altered P acquisition by growing wild, traditional (pre-GR), and modern (post-GR) tomato in soil with legacy P but low bioavailable P. Wild tomatoes, particularly accession LA0716 (Solanum pennellii), heavily cultured rhizosphere P solubilizers, suggesting reliance on microbial associations to acquire P. Wild tomato also had a greater abundance of other putatively beneficial bacteria, including those that produce chelating agents and antibiotic compounds. Although wild tomatoes had a high abundance of these P solubilizers, they had lower relative biomass and greater P stress factor than traditional or modern tomato. Compared to wild tomato, domesticated tomato was more tolerant to P deficiency, and both cultivated groups had a similar rhizosphere bacterial community composition. Ultimately, this study suggests that while domestication changed tomato P recovery by reducing microbial associations, subsequent breeding processes have not further impacted microbial P acquisition mechanisms. Selecting microbial P-related traits that diminished with domestication may therefore increase legacy P solubilization.

Genomic variation in weedy and cultivated broomcorn millet accessions uncovers the genetic architecture of agronomic traits.

Large-scale genomic variations are fundamental resources for crop genetics and breeding. Here we sequenced 1,904 genomes of broomcorn millet to an average of 40× sequencing depth and constructed a comprehensive variation map of weedy and cultivated accessions. Being one of the oldest cultivated crops, broomcorn millet has extremely low nucleotide diversity and remarkably rapid decay of linkage disequilibrium. Genome-wide association studies identified 186 loci for 12 agronomic traits. Many causative candidate genes, such as PmGW8 for grain size and PmLG1 for panicle shape, showed strong selection signatures during domestication. Weedy accessions contained many beneficial variations for the grain traits that are largely lost in cultivated accessions. Weedy and cultivated broomcorn millet have adopted different loci controlling flowering time for regional adaptation in parallel. Our study uncovers the unique population genomic features of broomcorn millet and provides an agronomically important resource for cereal crops.

Advances in the evolution research and genetic breeding of peanut.

Peanut is an important cash crop used in oil, food and feed in our country. The rapid development of sequencing technology has promoted the research on the related aspects of peanut genetic breeding. This paper reviews the research progress of peanut origin and evolution, genetic breeding, molecular markers and their applications, genomics, QTL mapping and genome selection techniques. The main problems of molecular genetic breeding in peanut research worldwide include: the narrow genetic resources of cultivated species, unstable genetic transformation and unclear molecular mechanism of important agronomic traits. Considering the severe challenges regarding the supply of edible oil, and the main problems in peanut production, the urgent research directions of peanut are put forward: The de novo domestication and the exploitation of excellent genes from wild resources to improve modern cultivars; Integration of multi-omics data to enhance the importance of big data in peanut genetics and breeding; Cloning the important genes related to peanut agronomic traits and analyzing their fine regulation mechanisms; Precision molecular design breeding and using gene editing technology to accurately improve the key traits of peanut.

The domestication of the larynx: The neural crest connection.

Wolves howl and dogs bark, both are able to produce variants of either vocalization, but we see a distinct difference in usage between wild and domesticate. Other domesticates also show distinct changes to their vocal output: domestic cats retain meows, a distinctly subadult trait in wildcats. Such differences in acoustic output are well-known, but the causal mechanisms remain little-studied. Potential links between domestication and vocal output are intriguing for multiple reasons, and offer a unique opportunity to explore a prominent hypothesis in domestication research: the neural crest/domestication syndrome hypothesis. This hypothesis suggests that in the early stages of domestication, selection for tame individuals decreased neural crest cell (NCCs) proliferation and migration, which led to a downregulation of the sympathetic arousal system, and hence reduced fear and reactive aggression. NCCs are a transitory stem cell population crucial during embryonic development that tie to diverse tissue types and organ systems. One of these neural-crest derived systems is the larynx, the main vocal source in mammals. We argue that this connection between NCCs and the larynx provides a powerful test of the predictions of the neural crest/domestication syndrome hypothesis, discriminating its predictions from those of other current hypotheses concerning domestication.

Cytoplasmic genome contributions to domestication and improvement of modern maize.

BACKGROUND: Studies on maize evolution and domestication are largely limited to the nuclear genomes, and the contribution of cytoplasmic genomes to selection and domestication of modern maize remains elusive. Maize cytoplasmic genomes have been classified into fertile (NA and NB) and cytoplasmic-nuclear male-sterility (CMS-S, CMS-C, and CMS-T) groups, but their contributions to modern maize breeding have not been systematically investigated. RESULTS: Here we report co-selection and convergent evolution between nuclear and cytoplasmic genomes by analyzing whole genome sequencing data of 630 maize accessions modern maize and its relatives, including 24 fully assembled mitochondrial and chloroplast genomes. We show that the NB cytotype is associated with the expansion of modern maize to North America, gradually replaces the fertile NA cytotype probably through unequal division, and predominates in over 90% of modern elite inbred lines. The mode of cytoplasmic evolution is increased nucleotypic diversity among the genes involved in photosynthesis and energy metabolism, which are driven by selection and domestication. Furthermore, genome-wide association study reveals correlation of cytoplasmic nucleotypic variation with key agronomic and reproductive traits accompanied with the diversification of the nuclear genomes. CONCLUSIONS: Our results indicate convergent evolution between cytoplasmic and nuclear genomes during maize domestication and breeding. These new insights into the important roles of mitochondrial and chloroplast genomes in maize domestication and improvement should help select elite inbred lines to improve yield stability and crop resilience of maize hybrids.