The western Mediterranean region provided the founder population of domesticated narrow-leafed lupin

Key message

This study revealed that the western Mediterranean provided the founder population for domesticated narrow-leafed lupin and that genetic diversity decreased significantly during narrow-leafed lupin domestication.

Abstract

The evolutionary history of plants during domestication profoundly shaped the genome structure and genetic diversity of today’s crops. Advances in next-generation sequencing technologies allow unprecedented opportunities to understand genome evolution in minor crops, which constitute the majority of plant domestications. A diverse set of 231 wild and domesticated narrow-leafed lupin (Lupinus angustifolius L.) accessions were subjected to genotyping-by-sequencing using diversity arrays technology. Phylogenetic, genome-wide divergence and linkage disequilibrium analyses were applied to identify the founder population of domesticated narrow-leafed lupin and the genome-wide effect of domestication on its genome. We found wild western Mediterranean population as the founder of domesticated narrow-leafed lupin. Domestication was associated with an almost threefold reduction in genome diversity in domesticated accessions compared to their wild relatives. Selective sweep analysis identified no significant footprints of selection around domestication loci. A genome-wide association study identified single nucleotide polymorphism markers associated with pod dehiscence. This new understanding of the genomic consequences of narrow-leafed lupin domestication along with molecular marker tools developed here will assist plant breeders more effectively access wild genetic diversity for crop improvement.

     

Phenotypic differentiation between wild and domesticated varieties of <Emphasis Type="Italic">Crescentia cujete</Emphasis>L. and culturally relevant uses of their fruits as bowls in the Yucatan Peninsula,Mexico

Background

Selection criteria are important for analyzing domestication of perennial plant species, which experience a selection pressure throughout several human generations. We analyze the preferred morphological characteristics of Crescentia cujete fruits, which are used as bowls by the Maya of Yucatan, according to the uses they are given and the phenotypic consequences of artificial selection between one wild and three domesticated varieties.

Methods

We performed 40 semi-structured interviews in seven communities. We calculated Sutrop’s salience index (S) of five classes of ceremonial and daily life uses, and of each item from the two most salient classes. We sampled 238 bowls at homes of people interviewed and compared their shape, volume and thickness with 139 fruits collected in homegardens and 179 from the wild. Morphology of varieties was assessed in fruit (n?=?114 trees) and vegetative characters (n?=?136 trees). Differences between varieties were evaluated through linear discriminant analysis (LDA).

Results

Use of bowls as containers for the Day of the Dead offerings was the most salient class (S?=?0.489) with chocolate as its most salient beverage (S?=?0.491), followed by consumption of daily beverages (S?=?0.423), especially maize-based pozol (S?=?0.412). The sacred saka’ and balche' are offered in different sized bowls during agricultural and domestic rituals. Roundness was the most relevant character for these uses, as bowls from households showed a strong selection towards round shapes compared with wild and homegarden fruits. Larger fruits from domesticated varieties were also preferred over small wild fruits, although in the household different sizes of the domesticated varieties are useful. LDA separated wild from domesticated trees (p?<?0.001) according to both fruit and vegetative variables, but domesticated varieties were not different among themselves.

Conclusions

The association between C. cujete bowls and traditional beverages in ritual and daily life situations has driven for centuries the selection of preferred fruit morphology in this tree. Selection of fruit roundness and volume has allowed for the differentiation between the wild variety and the three domesticated ones, counteracting gene flow among them. By choosing the best fruits from domesticated varieties propagated in homegardens, the Maya people model the domestication process of this important tree in their culture.

     

Polymorphism pattern at a miniature inverted‐repeat transposable element locus downstream of the domestication gene Teosinte‐branched1 in wild and domesticated pearl millet

Unravelling the mechanisms involved in adaptation to understand plant morphological evolution is a challenging goal. For crop species, identification of molecular causal polymorphisms involved in domestication traits is central to this issue. Pearl millet, a domesticated grass mostly found in semi‐arid areas of Africa and India, is an interesting model to address this topic: the domesticated form shares common derived phenotypes with some other cereals such as a decreased ability to develop basal and axillary branches in comparison with the wild phenotype. Two recent studies have shown that the orthologue of the maize gene Teosinte‐Branched1 in pearl millet (PgTb1) was probably involved in branching evolution during domestication and that a miniature inverted‐repeat transposable element (MITE) of the Tuareg family was inserted in the 3′ untranslated region of PgTb1. For a set of 35 wild and domesticated populations, we compared the polymorphism patterns at this MITE and at microsatellite loci. The Tuareg insertion was nearly absent in the wild populations, whereas a strong longitudinal frequency cline was observed in the domesticated populations. The geographical pattern revealed by neutral microsatellite loci clearly demonstrated that isolation by distance does not account for the existence of this cline. However, comparison of population differentiation at the microsatellite and the MITE loci and analyses of the nucleotide polymorphism pattern in the downstream region of PgTb1 did not show evidence that the cline at the MITE locus has been shaped by selection, suggesting the implication of a neutral process. Alternative hypotheses are discussed.     

An Age-of-Allele Test of Neutrality for Transposable Element Insertions

How natural selection acts to limit the proliferation of transposable elements (TEs) in genomes has been of interest to evolutionary biologists for many years. To describe TE dynamics in populations, previous studies have used models of transposition–selection equilibrium that assume a constant rate of transposition. However, since TE invasions are known to happen in bursts through time, this assumption may not be reasonable. Here we propose a test of neutrality for TE insertions that does not rely on the assumption of a constant transposition rate. We consider the case of TE insertions that have been ascertained from a single haploid reference genome sequence. By conditioning on the age of an individual TE insertion allele (inferred by the number of unique substitutions that have occurred within the particular TE sequence since insertion), we determine the probability distribution of the insertion allele frequency in a population sample under neutrality. Taking models of varying population size into account, we then evaluate predictions of our model against allele frequency data from 190 retrotransposon insertions sampled from North American and African populations of Drosophila melanogaster. Using this nonequilibrium neutral model, we are able to explain ∼80% of the variance in TE insertion allele frequencies based on age alone. Controlling for both nonequilibrium dynamics of transposition and host demography, we provide evidence for negative selection acting against most TEs as well as for positive selection acting on a small subset of TEs. Our work establishes a new framework for the analysis of the evolutionary forces governing large insertion mutations like TEs, gene duplications, or other copy number variants.     

Phylogeny and evolutionary history of the silkworm

The silkworm, Bombyx mori, played an important role in the old Silk Road that connected ancient Asia and Europe. However, to date, there have been few studies of the origins and domestication of this species using molecular methods. In this study, DNA sequences of mitochondrial and nuclear loci were used to infer the phylogeny and evolutionary history of the domesticated silkworm and its relatives. All of the phylogenetic analyses indicated a close relationship between the domesticated silkworm and the Chinese wild silkworm. Domestication was estimated to have occurred about 4100 years ago (ya), and the radiation of the different geographic strains of B. mori about 2000 ya. The Chinese wild silkworm and the Japanese wild silkworm split about 23600 ya. These estimates are in good agreement with the fossil evidence and historical records. In addition, we show that the domesticated silkworm experienced a population expansion around 1000 ya. The divergence times and the population dynamics of silkworms presented in this study will be useful for studies of lepidopteran phylogenetics, in the genetic analysis of domestic animals, and for understanding the spread of human civilizations.     

Insertional variability of four transposable elements and population structure of the midge Chironomus riparius (Diptera)

The dipteran Chironomus riparius is found across the entire Palearctic region; its larvae are among the most abundant macroinvertebrates inhabiting inland waterbodies. Chironomid larvae have been extensively used in ecotoxicological and cytogenetic research, but relatively little is known on the population structure of this species. Transposable elements (TEs) are DNA sequences that are capable of autonomous replication; the number and genomic location of TE insertions varies across individuals; this variability is increasingly being used in population studies. Several TEs had been characterized in Chironomids; this enabled the analysis of insertional variability of four different TEs in six natural populations of C. riparius from Italy, Bulgaria and Russia using a PCR-based method, transposon insertion display (TID). The method allows to obtain dominant markers, similar to AFLP. In all populations, TE insertions showed high individual polymorphism, while median copy numbers of the same TEs did not vary between populations. Analysis of molecular variance (AMOVA) detected significant differentiation between populations for three of the TEs; although no correlation between genetic and geographic distances was found, the corresponding population structures were found to be significantly correlated and indicate a degree of isolation by distance. TEs belonging to different classes have different mechanisms of replication, resulting in different transposition rates of mobilization; the finding of mostly concordant population structuring for three of the TEs indicates that population dynamics contributed significantly in shaping the detected insertional polymorphism.     

Intraspecific and intraorganismal copy number dynamics of retrotransposons and tandem repeat in <Emphasis Type="Italic">Aegilops speltoides</Emphasis> Tausch (Poaceae,Triticeae)

Transposable elements (TE) and tandem repeats (TR) compose the largest fraction of the plant genome. The abundance and repatterning of repetitive DNA underlie intrapopulation polymorphisms and intraspecific diversification; however, the dynamics of repetitive elements in ontogenesis is not fully understood. Here, we addressed the genotype-specific and tissue-specific abundances and dynamics of the Ty1-copia, Ty3-gypsy, and LINE retrotransposons and species-specific Spelt1 tandem repeat in wild diploid goatgrass, Aegilops speltoides Tausch. Copy numbers of TEs and TR were estimated by real-time quantitative PCR in vegetative and generative tissues in original plants from contrasting allopatric populations and artificial intraspecific hybrids. The results showed that between leaves and somatic spike tissues as well as in progressive microsporogenesis of individual genotypes, the copy numbers of three TEs correlatively oscillated between 2- to 4-fold and the TR copy numbers fluctuated by 18- to 440-fold. Inter-individual and intraorganismal TEs and TR copy number dynamics demonstrate large-scale parallelism with extensive chromosomal rearrangements that were detected using fluorescent in situ hybridization in parental and hybrid genotypes. The data obtained indicate that tissue-specific differences in the abundance and pattern of repetitive sequences emerge during cell proliferation and differentiation in ontogenesis and reflect the reorganization of individual genomes in changing environments, especially in small peripheral population(s) under the influence of rapid climatic changes.     

Domestication history in the Medicago sativa species complex: inferences from nuclear sequence polymorphism

DNA sequence polymorphism carries genealogical information and allows for testing hypotheses on selection and population history, especially through coalescent-based analysis. Understanding the evolutionary forces at work in plant domestication and subsequent selection is of critical importance for the management of genetic resources. In this study, we surveyed DNA sequence diversity at two assumed neutral nuclear loci in the wild-domesticated species complex of alfalfa (Medicago sativa L.). A high level of polymorphism was detected. The domesticated pool contains on average 31% less diversity than the wild pool, but with a high heterogeneity among loci. Coalescent simulations of the domestication process showed that this result cannot be explained by assuming a constant population size but is rather consistent with a demographic bottleneck during domestication. A very low level of divergence was detected between the wild and the domesticated forms as well as between the related subspecies of the M. sativa species complex. However, the originality of the Spanish wild populations, already observed based on mitochondrial DNA polymorphism, was confirmed. These results, together with patterns of intrapopulation polymorphism, suggest that nuclear sequence polymorphism could be a promising tool, complementary to mitochondrial DNA and phenotypic evaluations, to investigate historical demographic and evolutionary processes.     

The dynamics of transposable elements in structured populations

We analyzed the dynamics of transposable elements (TEs) according to Wright's island and continent-island models, assuming that selection tends to counter the deleterious effects of TEs. We showed that migration between host populations has no impact on either the existence or the stability of the TE copy number equilibrium points obtained in the absence of migration. However, if the migration rate is slower than the transposition rate or if selection is weak, then the TE copy numbers in all the populations can be expected to slowly become homogeneous, whereas a heterogeneous TE copy number distribution between populations is maintained if TEs are mobilized in some populations. The mean TE copy number is highly sensitive to the population size, but as a result of migration between populations, it decreases as the sum of the population sizes increases and tends to reach the same value in these populations. We have demonstrated the existence of repulsion between TE insertion sites, which is established by selection and amplified by drift. This repulsion is reduced as much as the migration rate is higher than the recombination rate between the TE insertion sites. Migration and demographic history are therefore strong forces in determining the dynamics of TEs within the genomes and the populations of a species.     

Inference of domestication history and differentiation between early‐ and late‐flowering varieties in pearl millet

Pearl millet (Pennisetum glaucum) is a staple crop in Sahelian Africa. Farmers usually grow varieties with different cycle lengths and complementary functions in Sahelian agrosystems. Both the level of genetic differentiation of these varieties and the domestication history of pearl millet have been poorly studied. We investigated the neutral genetic diversity and population genetic structure of early‐ and late‐flowering domesticated and wild pearl millet populations using 18 microsatellite loci and 8 nucleotide sequences. Strikingly, early‐ and late‐flowering domesticated varieties were not differentiated over their whole distribution area, despite a clear difference in their isolation‐by‐distance pattern. Conversely, our data brought evidence for two well‐differentiated genetic pools in wild pearl millet, allowing us to test scenarios with different numbers and origins of domestication using approximate Bayesian computation (ABC). The ABC analysis showed the likely existence of asymmetric migration between wild and domesticated populations. The model choice procedure indicated that a single domestication from the eastern wild populations was the more likely scenario to explain the polymorphism patterns observed in cultivated pearl millet.     

Synergistic epistasis of the deleterious effects of transposable elements

The replicative nature and generally deleterious effects of transposable elements (TEs) raise an outstanding question about how TE copy number is stably contained in host populations. Classic theoretical analyses predict that, when the decline in fitness due to each additional TE insertion is greater than linear, or when there is synergistic epistasis, selection against TEs can result in a stable equilibrium of TE copy number. While several mechanisms are predicted to yield synergistic deleterious effects of TEs, we lack empirical investigations of the presence of such epistatic interactions. Purifying selection with synergistic epistasis generates repulsion linkage between deleterious alleles. We investigated this population genetic signal in the likely ancestral Drosophila melanogaster population and found evidence supporting the presence of synergistic epistasis among TE insertions, especially TEs expected to exert large fitness impacts. Even though synergistic epistasis of TEs has been predicted to arise through ectopic recombination and TE-mediated epigenetic silencing mechanisms, we only found mixed support for the associated predictions. We observed signals of synergistic epistasis for a large number of TE families, which is consistent with the expectation that such epistatic interaction mainly happens among copies of the same family. Curiously, significant repulsion linkage was also found among TE insertions from different families, suggesting the possibility that synergism of TEs’ deleterious fitness effects could arise above the family level and through mechanisms similar to those of simple mutations. Our findings set the stage for investigating the prevalence and importance of epistatic interactions in the evolutionary dynamics of TEs.     

Evidence of selection at melanin synthesis pathway loci during silkworm domestication

The domesticated silkworm (Bombyx mori) was domesticated from wild silkworm (Bombyx mandarina) more than 5,000 years ago. During domestication, body color between B. mandarina and B. mori changed dramatically. However, the molecular mechanism of the silkworm body color transition is not known. In the present study, we examined within- and between-species nucleotide diversity for eight silkworm melanin synthesis pathway genes, which play a key role in cuticular pigmentation of insects. Our results showed that the genetic diversity of B. mori was significantly lower than that of B. mandarina and 40.7% of the genetic diversity of wild silkworm was lost in domesticated silkworm. We also examined whether position effect exists among melanin synthesis pathway genes in B. mandarina and B. mori. We found that the upstream genes have significantly lower levels of genetic diversity than the downstream genes, supporting a functional constraint hypothesis (FCH) of metabolic pathway, that is, upstream enzymes are under greater selective constraint than downstream enzymes because upstream enzymes participate in biosynthesis of a number of metabolites. We also investigated whether some of the melanin synthesis pathway genes experienced selection during domestication. Neutrality test, coalescent simulation, as well as network and phylogenetic analyses showed that tyrosine hydroxylase (TH) gene was a domestication locus. Sequence analysis further suggested that a putative expression enhancer (Abd-B-binding site) in the intron of TH gene might be disrupted during domestication. TH is the rate-limiting enzyme of melanin synthesis pathway in insects. Real-time polymerase chain reaction assay did show that the relative expression levels of TH gene in B. mori were significantly lower than that in B. mandarina at three different developmental stages, which is consistent with light body color of domesticated silkworm relative to wild silkworm. Therefore, we speculated that expression change of TH gene may contribute to the body color transition from B. mandarina to B. mori. Our results emphasize the exceptional role of gene expression regulation in morphological transition of domesticated animals.     

Genome-wide identification and evolution of <Emphasis Type="Italic">TC1</Emphasis>/<Emphasis Type="Italic">Mariner</Emphasis> in the silkworm (<Emphasis Type="Italic">Bombyx mori</Emphasis>) genome

TC1/Mariner transposons belong to class II transposable elements (TEs) that use DNA-mediated “cut and paste” mechanism to transpose, and they have been identified in almost all organisms. Although silkworm (Bombyx mori) has a large amount of TC1/Mariner elements, the genome wide information of this superfamily in the silkworm is unknown. In this study, we have identified 2670 TC1/Mariner (Bmmar) elements in the silkworm genome. All the TEs were classified into 22 families by means of fgclust, a tool of repetitive sequence classification, seven of which was first reported in this study. Phylogenetic and structure analyses based on the catalytic domain (DDxD/E) of transposase sequences indicated that all members of TC1/Mariner were grouped into five subgroups: Mariner, Tc1, maT, DD40D and DD41D/E. Of these five subgroups, maT rather than Mariner possessed most members of TC1/Mariner (51.23%) in the silkworm genome. In particular, phylogenetic analysis and structure analysis revealed that Bmmar15 (DD40D) formed a new basal subgroup of TC1/Mariner element in insects, which was referred to as bmori. Furthermore, we concluded that DD40D appeared to intermediate between mariner and Tc1. Finally, we estimated the insertion time for each copy of TC1/Mariner in the silkworm and found that most of members were dramatically amplified during a period from 0 to 1 mya. Moreover, the detailed functional data analysis showed that Bmmar1, Bmmar6 and Bmmar9 had EST evidence and intact transposases. These implied that TC1/Mariner might have potential transpositional activity. In conclusion, this study provides some new insights into the landscape, origin and evolution of TC1/Mariner in the insect genomes.     

Identifying Domesticated and Wild Kañawa (<Emphasis Type="Italic">Chenopodium pallidicaule</Emphasis>) in the Archeobotanical Record of the Lake Titicaca Basin of the Andes

Here, we present new data on how to identify both the wild and domesticated forms of kañawa (Chenopodium pallidicaule Aellen) in the archeological record of Andean South America using characteristics of their seed morphology. The ability to identify both the domesticated and wild forms of kañawa is an essential step in advancing our understanding of the processes of its domestication, diversification, and the role it has played in past food systems throughout the Andes.     

Molecular Evidence and the Origin and Development of the Domesticated Sunflower (<Emphasis Type="Italic">Helianthus annum</Emphasis>, Asteraceae)

The domesticated sunflower,Helianthus annuus, is an important economic crop, yet molecular data regarding its evolution are limited. Here we review morphological, geographical, archaeological, and molecular evidence pertaining to its origin and development. New isozyme and chloroplast DNA (cpDNA) evidence is also presented.Morphological, geographical, and archaeological evidence has led to the hypothesis that the domesticated sunflower was derived from a wild/weedy form ofH. annuus possibly in the Midwest. Molecular evidence was concordant with this hypothesis. A high degree of enzymatic and cpDNA sequence similarity was observed between wild and domesticatedH. annuus, and domesticatedH. annuus contained a subset of the alleles and cpDNAs found in wildH. annuus. The extensive polymorphism in the wild plants and the virtual monomorphism in cultivated lines for both isozyme and cpDNA phenotypes further suggest a single origin of the domesticated sunflower from a very limited gene pool. In addition, Native American varieties of the domesticated sunflower were genetically more variable than other cultivated lines, possibly indicating that they gave rise to the other cultivated stocks. Molecular evidence did not, however, allow conclusions as to the exact geographic origin of the domesticated sunflower.     

Chromosomal Distribution of Transposable Elements in Drosophila Melanogaster: Test of the Ectopic Recombination Model for Maintenance of Insertion Site Number

Data of insertion site localization and site occupancy frequency of P, hobo, I, copia, mdg1, mdg3, 412, 297, and roo transposable elements (TEs) on the polytene chromosomes of Drosophila melanogaster were extracted from the literature. We show that TE insertion site number per chromosomal division was significantly correlated with the amount of DNA. The insertion site number weighted by DNA content was not correlated with recombination rate for all TEs except hobo, for which a positive correlation was detected. No global tendency emerged in the relationship between TE site occupancy frequency, weighted by DNA content, and recombination rate; a strong negative correlation was, however, found for the 3L arm. A possible dominant deleterious effect of chromosomal rearrangements due to recombination between TE insertions is thus not the main factor explaining the dynamics of TEs, since this hypothesis implies a negative relationship between recombination rate and both TE insertion site number and site occupancy frequency. The alternative hypothesis of selection against deleterious effects of insertional mutations is discussed.     

Tagging the signatures of domestication in common bean (Phaseolus vulgaris) by means of pooled DNA samples

BACKGROUND AND AIMS: The main aim of this study was to use an amplified fragment length polymorphism (AFLP)-based, large-scale screening of the whole genome of Phaseolus vulgaris to determine the effects of selection on the structure of the genetic diversity in wild and domesticated populations. METHODS: Using pooled DNA samples, seven each of wild and domesticated populations of P. vulgaris were studied using 2506 AFLP markers (on average, one every 250 kb). About 10 % of the markers were also analysed on individual genotypes and were used to infer allelic frequencies empirically from bulk data. In both data sets, tests were made to determine the departure from neutral expectation for each marker using an F(ST)-based method. KEY RESULTS: The most important outcome is that a large fraction of the genome of the common bean (16 %; P < 0.01) appears to have been subjected to effects of selection during domestication. Markers obtained in individual genotypes were also mapped and classified according to their proximities to known genes and quantitative trait loci (QTLs) of the domestication syndrome. Most of the markers that were found to be potentially under the effects of selection were located in the proximity of previously mapped genes and QTLs related to the domestication syndrome. CONCLUSIONS: Overall, the results indicate that in P. vulgaris a large portion of the genome appears to have been subjected to the effects of selection, probably because of linkage to the loci selected during domestication. As most of the markers that are under the effects of selection are linked to known loci related to the domestication syndrome, it is concluded that population genomics approaches are very efficient in detecting QTLs. A method based on bulk DNA samples is presented that is effective in pre-screening for a large number of markers to determine selection signatures.     

Evolution of the population structure of Venturia inaequalis,the apple scab fungus,associated with the domestication of its host

Evaluating the impact of plant domestication on the population structure of the associated pathogens provides an opportunity to increase our understanding of how and why diseases emerge. Here, we investigated the evolution of the population structure of the apple scab fungus Venturia inaequalis in response to the domestication of its host. Inferences were drawn from multilocus microsatellite data obtained from samples collected on (i) the Central Asian Malus sieversii, the main progenitor of apple, (ii) the European crabapple, Malus sylvestris, a secondary progenitor of apple, and (iii) the cultivated apple, Malus×domestica, in orchards from Europe and Central Asia. Using clustering methods, we identified three distinct populations: (i) a large European population on domesticated and wild apples, (ii) a large Central Asian population on domesticated and wild apples in urban and agricultural areas, and (iii) a more geographically restricted population in M. sieversii forests growing in the eastern mountains of Kazakhstan. Unique allele richness and divergence time estimates supported a host‐tracking co‐evolutionary scenario in which this latter population represents a relict of the ancestral populations from which current populations found in human‐managed habitats were derived. Our analyses indicated that the domestication of apple induced a significant change in the genetic differentiation of populations of V. inaequalis in its centre of origin, but had little impact on its population dynamics and mating system. We discuss how the structure of the apple‐based agrosystem may have restricted changes in the population structure of the fungus in response to the domestication of its host.