Moulded by humans: The domestication of blue‐veined cheese fungi

It is hard to imagine a world without food‐associated microbes. The production of bread, wine, beer, salami, coffee, chocolate, cheese and many other foods and beverages all rely on specific microbes. In cheese, myriad microbial species collaborate to yield the complex organoleptic properties that are appreciated by millions of people worldwide. In the early days of cheese making, these complex communities emerged spontaneously from the natural flora associated with the raw materials, the equipment, the production environment or craftsmen involved in the production process. However, in some cases, the microbes shifted their natural habitat to the new cheese‐associated environment. The most obvious cause of this is backslopping, where part of a fermented product is used to inoculate the next batch. In addition, some microbes may simply adhere to the tools used in the production process. These microbial communities gradually adapted to the novel man‐made niches, a process referred to as “domestication.” Domestication is associated with specific genomic and phenotypic changes and ultimately leads to lineages that are genetically and phenotypically distinct from their wild ancestors. In this issue of Molecular Ecology, Dumas et al. have investigated a prime example of cheese‐associated microbes, the fungus Penicillium roqueforti. The authors identified several hallmarks of domestication in the genome and phenome of this species, allowing them to hypothesize about the origin of blue‐veined cheese fungi domestication, and the specific evolutionary processes involved in adaptation to the cheese matrix.     

Largely unlinked gene sets targeted by selection for domestication syndrome phenotypes in maize and sorghum

The domestication of diverse grain crops from wild grasses was a result of artificial selection for a suite of overlapping traits producing changes referred to in aggregate as ‘domestication syndrome’. Parallel phenotypic change can be accomplished by either selection on orthologous genes or selection on non‐orthologous genes with parallel phenotypic effects. To determine how often artificial selection for domestication traits in the grasses targeted orthologous genes, we employed resequencing data from wild and domesticated accessions of Zea (maize) and Sorghum (sorghum). Many ‘classic’ domestication genes identified through quantitative trait locus mapping in populations resulting from wild/domesticated crosses indeed show signatures of parallel selection in both maize and sorghum. However, the overall number of genes showing signatures of parallel selection in both species is not significantly different from that expected by chance. This suggests that while a small number of genes will extremely large phenotypic effects have been targeted repeatedly by artificial selection during domestication, the optimization part of domestication targeted small and largely non‐overlapping subsets of all possible genes which could produce equivalent phenotypic alterations.     

In the shadows of snow leopards and the Himalayas: density and habitat selection of blue sheep in Manang,Nepal

There is a growing agreement that conservation needs to be proactive and pay increased attention to common species and to the threats they face. The blue sheep (Pseudois nayaur) plays a key ecological role in sensitive high‐altitude ecosystems of Central Asia and is among the main prey species for the globally vulnerable snow leopard (Panthera uncia). As the blue sheep has been increasingly exposed to human pressures, it is vital to estimate its population dynamics, protect the key populations, identify important habitats, and secure a balance between conservation and local livelihoods. We conducted a study in Manang, Annapurna Conservation Area (Nepal), to survey blue sheep on 60 transects in spring (127.9 km) and 61 transects in autumn (134.7 km) of 2019, estimate their minimum densities from total counts, compare these densities with previous estimates, and assess blue sheep habitat selection by the application of generalized additive models (GAMs). Total counts yielded minimum density estimates of 6.0–7.7 and 6.9–7.8 individuals/km² in spring and autumn, respectively, which are relatively high compared to other areas. Elevation and, to a lesser extent, land cover indicated by the normalized difference vegetation index (NDVI) strongly affected habitat selection by blue sheep, whereas the effects of anthropogenic variables were insignificant. Animals were found mainly in habitats associated with grasslands and shrublands at elevations between 4,200 and 4,700 m. We show that the blue sheep population size in Manang has been largely maintained over the past three decades, indicating the success of the integrated conservation and development efforts in this area. Considering a strong dependence of snow leopards on blue sheep, these findings give hope for the long‐term conservation of this big cat in Manang. We suggest that long‐term population monitoring and a better understanding of blue sheep–livestock interactions are crucial to maintain healthy populations of blue sheep and, as a consequence, of snow leopards.     

The role of Bh4 in parallel evolution of hull colour in domesticated and weedy rice

The two independent domestication events in the genus Oryza that led to African and Asian rice offer an extremely useful system for studying the genetic basis of parallel evolution. This system is also characterized by parallel de‐domestication events, with two genetically distinct weedy rice biotypes in the US derived from the Asian domesticate. One important trait that has been altered by rice domestication and de‐domestication is hull colour. The wild progenitors of the two cultivated rice species have predominantly black‐coloured hulls, as does one of the two U.S. weed biotypes; both cultivated species and one of the US weedy biotypes are characterized by straw‐coloured hulls. Using Black hull 4 (Bh4) as a hull colour candidate gene, we examined DNA sequence variation at this locus to study the parallel evolution of hull colour variation in the domesticated and weedy rice system. We find that independent Bh4‐coding mutations have arisen in African and Asian rice that are correlated with the straw hull phenotype, suggesting that the same gene is responsible for parallel trait evolution. For the U.S. weeds, Bh4 haplotype sequences support current hypotheses on the phylogenetic relationship between the two biotypes and domesticated Asian rice; straw hull weeds are most similar to indica crops, and black hull weeds are most similar to aus crops. Tests for selection indicate that Asian crops and straw hull weeds deviate from neutrality at this gene, suggesting possible selection on Bh4 during both rice domestication and de‐domestication.     

Population genomics reveals a fine‐scale recombination landscape for genetic improvement of cotton

Recombination breaks up ancestral linkage disequilibrium, creates combinations of alleles, affects the efficiency of natural selection, and plays a major role in crop domestication and improvement. However, there is little knowledge regarding the variation in the population‐scaled recombination rate in cotton. We constructed recombination maps and characterized the difference in the genomic landscape of the population‐scaled recombination rate between Gossypium hirsutum and G. arboreum and sub‐genomes based on the 381 sequenced G. hirsutum and 215 G. arboreum accessions. Comparative genomics identified large structural variations and syntenic genes in the recombination regions, suggesting that recombination was related to structural variation and occurred preferentially in the distal chromosomal regions. Correlation analysis indicated that recombination was only slightly affected by geographical distribution and breeding period. A genome‐wide association study (GWAS) was performed with 15 agronomic traits using 267 cotton accessions and identified 163 quantitative trait loci (QTL) and an important candidate gene (Ghir_COL2) for early maturity traits. Comparative analysis of recombination and a GWAS revealed that the QTL of fibre quality traits tended to be more common in high‐recombination regions than were those of yield and early maturity traits. These results provide insights into the population‐scaled recombination landscape, suggesting that recombination contributed to the domestication and improvement of cotton, which provides a useful reference for studying recombination in other species.     

Global investigation of the co‐evolution of MIRNA genes and microRNA targets during soybean domestication

Although the selection of coding genes during plant domestication has been well studied, the evolution of MIRNA genes (MIRs) and the interaction between microRNAs (miRNAs) and their targets in this process are poorly understood. Here, we present a genome‐wide survey of the selection of MIRs and miRNA targets during soybean domestication and improvement. Our results suggest that, overall, MIRs have higher evolutionary rates than miRNA targets. Nonetheless, they do demonstrate certain similar evolutionary patterns during soybean domestication: MIRs and miRNA targets with high expression and duplication status, and with greater numbers of partners, exhibit lower nucleotide divergence than their counterparts without these characteristics, suggesting that expression level, duplication status, and miRNA–target interaction are essential for evolution of MIRs and miRNA targets. Further investigation revealed that miRNA–target pairs that are subjected to strong purifying selection have greater similarities than those that exhibited genetic diversity. Moreover, mediated by domestication and improvement, the similarities of a large number of miRNA–target pairs in cultivated soybean populations were increased compared to those in wild soybeans, whereas a small number of miRNA–target pairs exhibited decreased similarity, which may be associated with the adoption of particular domestication traits. Taken together, our results shed light on the co‐evolution of MIRs and miRNA targets during soybean domestication.     

Geographical delimitation of a partial selective sweep in African Drosophila melanogaster

Positive selection leaves characteristic footprints on DNA variation but detecting such patterns is challenging as the age, the intensity and the mode of selection as well as demography and evolutionary parameters (mutation and recombination rates) all play roles and these are difficult to disentangle. We recorded nucleotide variation in a sample of isogenic chromosomes from a western African population of Drosophila melanogaster at a locus (Fbp2) for which a partial selective sweep had previously been reported. We compared this locus to four other genes from the same chromosomes and from a European and an East African population. Then, we assessed Fbp2 variation in a sample of 370 chromosomes covering a comprehensive geographic sampling of 16 African localities. The signature of selection was tested while accounting for the demographic history of the populations. We found a significant signal of selection in two West African localities including Ivory Coast. Variation at Fpb2 would thus represent a case of an ongoing selective sweep in the range of this species. A weaker, nonsignificant, signal of selection was, however, apparent in some other populations, thus leaving open several possibilities: (i) the selective sweep originated in Ivory Coast and has spread to the rest of the continent; (ii) several African populations report the signature of a selective event having occurred in an ancestral population; (iii) this genome region is subject to independent selective events in African populations; and (iv) A neutral scenario with population subdivision and local bottleneck cannot be fully excluded to explain the molecular patterns observed in some populations.     

Mitochondrial DNA variation of Nigerian domestic helmeted guinea fowl

We analyzed genetic diversity of 215 mitochondrial DNA (mtDNA) D‐loop sequences from seven populations of domesticated helmeted guinea fowl (Numida meleagris) in Nigeria and compared that with results of samples collected in Kenya (n = 4) and China (n = 22). In total, 241 sequences were assigned to 22 distinct haplotypes. Haplotype diversity in Nigeria was 0.693 ± 0.022. The network grouped most matrilines into two main haplogroups: A and B. There was an absence of a geographic signal, and two haplotypes dominated across all locations with the exception of the Kebbi population in the northwest of the country; AMOVA also confirmed this observation (F_ST = 0.035). The low genetic diversity may be a result of recent domestication, whereas the lack of maternal genetic structure likely suggests the extensive genetic intermixing within the country. Additionally, the differentiation of the Kebbi population may be due to a certain demographic history and/or artificial selection that shaped its haplotype profile. The current data do not permit us to make further conclusions; therefore, more research evidence from genetics and archaeology is still required.     

Analysis of ROH patterns in the Noriker horse breed reveals signatures of selection for coat color and body size

Overlapping runs of homozygosity (ROH islands) shared by the majority of a population are hypothesized to be the result of selection around a target locus. In this study we investigated the impact of selection for coat color within the Noriker horse on autozygosity and ROH patterns. We analyzed overlapping homozygous regions (ROH islands) for gene content in fragments shared by more than 50% of horses. Long‐term assortative mating of chestnut horses and the small effective population size of leopard spotted and tobiano horses resulted in higher mean genome‐wide ROH coverage (S_ROH) within the range of 237.4–284.2 Mb, whereas for bay, black and roan horses, where rotation mating is commonly applied, lower autozygosity (S_ROH from 176.4–180.0 Mb) was determined. We identified seven common ROH islands considering all Noriker horses from our dataset. Specific islands were documented for chestnut, leopard spotted, roan and bay horses. The ROH islands contained, among others, genes associated with body size (ZFAT, LASP1 and LCORL/NCAPG), coat color (MC1R in chestnut and the factor PATN1 in leopard spotted horses) and morphogenesis (HOXB cluster in all color strains except leopard spotted horses). This study demonstrates that within a closed population sharing the same founders and ancestors, selection on a single phenotypic trait, in this case coat color, can result in genetic fragmentation affecting levels of autozygosity and distribution of ROH islands and enclosed gene content.     

Genetic diversity and genomic signatures of selection among cattle breeds from Siberia,eastern and northern Europe

Domestication in the near eastern region had a major impact on the gene pool of humpless taurine cattle (Bos taurus). As a result of subsequent natural and artificial selection, hundreds of different breeds have evolved, displaying a broad range of phenotypic traits. Here, 10 Eurasian B. taurus breeds from different biogeographic and production conditions, which exhibit different demographic histories and have been under artificial selection at various intensities, were investigated using the Illumina BovineSNP50 panel to understand their genetic diversity and population structure. In addition, we scanned genomes from eight breeds for signatures of diversifying selection. Our population structure analysis indicated six distinct breed groups, the most divergent being the Yakutian cattle from Siberia. Selection signals were shared (experimental P‐value < 0.01) with more than four breeds on chromosomes 6, 7, 13, 16 and 22. The strongest selection signals in the Yakutian cattle were found on chromosomes 7 and 21, where a miRNA gene and genes related to immune system processes are respectively located. In general, genomic regions indicating selection overlapped with known QTL associated with milk production (e.g. on chromosome 19), reproduction (e.g. on chromosome 24) and meat quality (e.g. on chromosome 7). The selection map created in this study shows that native cattle breeds and their genetic resources represent unique material for future breeding.     

Domestication Process of the Goat Revealed by an Analysis of the Nearly Complete Mitochondrial Protein-Encoding Genes

Goats (Capra hircus) are one of the oldest domesticated species, and they are kept all over the world as an essential resource for meat, milk, and fiber. Although recent archeological and molecular biological studies suggested that they originated in West Asia, their domestication processes such as the timing of population expansion and the dynamics of their selection pressures are little known. With the aim of addressing these issues, the nearly complete mitochondrial protein-encoding genes were determined from East, Southeast, and South Asian populations. Our coalescent time estimations suggest that the timing of their major population expansions was in the Late Pleistocene and significantly predates the beginning of their domestication in the Neolithic era (≈10,000 years ago). The ω (ratio of non-synonymous rate/synonymous substitution rate) for each lineage was also estimated. We found that the ω of the globally distributed haplogroup A which is inherited by more than 90% of goats examined, turned out to be extremely low, suggesting that they are under severe selection pressure probably due to their large population size. Conversely, the ω of the Asian-specific haplogroup B inherited by about 5% of goats was relatively high. Although recent molecular studies suggest that domestication of animals may tend to relax selective constraints, the opposite pattern observed in our goat mitochondrial genome data indicates the process of domestication is more complex than may be presently appreciated and cannot be explained only by a simple relaxation model.     

A draft genome of sweet cherry (Prunus avium L.) reveals genome‐wide and local effects of domestication

Sweet cherry (Prunus avium L.) trees are both economically important fruit crops but also important components of natural forest ecosystems in Europe, Asia and Africa. Wild and domesticated trees currently coexist in the same geographic areas with important questions arising on their historical relationships. Little is known about the effects of the domestication process on the evolution of the sweet cherry genome. We assembled and annotated the genome of the cultivated variety “Big Star*” and assessed the genetic diversity among 97 sweet cherry accessions representing three different stages in the domestication and breeding process (wild trees, landraces and modern varieties). The genetic diversity analysis revealed significant genome‐wide losses of variation among the three stages and supports a clear distinction between wild and domesticated trees, with only limited gene flow being detected between wild trees and domesticated landraces. We identified 11 domestication sweeps and five breeding sweeps covering, respectively, 11.0 and 2.4 Mb of the P. avium genome. A considerable fraction of the domestication sweeps overlaps with those detected in the related species, Prunus persica (peach), indicating that artificial selection during domestication may have acted independently on the same regions and genes in the two species. We detected 104 candidate genes in sweep regions involved in different processes, such as the determination of fruit texture, the regulation of flowering and fruit ripening and the resistance to pathogens. The signatures of selection identified will enable future evolutionary studies and provide a valuable resource for genetic improvement and conservation programs in sweet cherry.     

Comparative population genomics identified genomic regions and candidate genes associated with fruit domestication traits in peach

Crop evolution is a long‐term process involving selection by natural evolutionary forces and anthropogenic influences; however, the genetic mechanisms underlying the domestication and improvement of fruit crops have not been well studied to date. Here, we performed a population structure analysis in peach (Prunus persica) based on the genome‐wide resequencing of 418 accessions and confirmed the presence of an obvious domestication event during evolution. We identified 132 and 106 selective sweeps associated with domestication and improvement, respectively. Analysis of their tissue‐specific expression patterns indicated that the up‐regulation of selection genes during domestication occurred mostly in fruit and seeds as opposed to other organs. However, during the improvement stage, more up‐regulated selection genes were identified in leaves and seeds than in the other organs. Genome‐wide association studies (GWAS) using 4.24 million single nucleotide polymorphisms (SNPs) revealed 171 loci associated with 26 fruit domestication traits. Among these loci, three candidate genes were highly associated with fruit weight and the sorbitol and catechin content in fruit. We demonstrated that as the allele frequency of the SNPs associated with high polyphenol composition decreased during peach evolution, alleles associated with high sugar content increased significantly. This indicates that there is genetic potential for the breeding of more nutritious fruit with enhanced bioactive polyphenols without disturbing a harmonious sugar and acid balance by crossing with wild species. This study also describes the development of the genomic resources necessary for evolutionary research in peach and provides the large‐scale characterization of key agronomic traits in this crop species.     

Population genomic,olfactory, dietary,and gut microbiota analyses demonstrate the unique evolutionary trajectory of feral pigs

Domestication is an intriguing evolutionary process. Many domestic populations are subjected to strong human-mediated selection, and when some individuals return to the wild, they are again subjected to selective forces associated with new environments. Generally, these feral populations evolve into something different from their wild predecessors and their members typically possess a combination of both wild and human selected traits. Feralisation can manifest in different forms on a spectrum from a wild to a domestic phenotype. This depends on how the rewilded domesticated populations can readapt to natural environments based on how much potential and flexibility the ancestral genome retains after its domestication signature. Whether feralisation leads to the evolution of new traits that do not exist in the wild or to convergence with wild forms, however, remains unclear. To address this question, we performed population genomic, olfactory, dietary, and gut microbiota analyses on different populations of Sus scrofa (wild boar, hybrid, feral and several domestic pig breeds). Porcine single nucleotide polymorphisms (SNPs) analysis shows that the feral population represents a cluster distinctly separate from all others. Its members display signatures of past artificial selection, as demonstrated by values of F_ST in specific regions of the genome and bottleneck signature, such as the number and length of runs of homozygosity. Generalised F_ST values, reacquired olfactory abilities, diet, and gut microbiota variation show current responses to natural selection. Our results suggest that feral pigs are an independent evolutionary unit which can persist so long as levels of human intervention remain unchanged.     

Detection of signatures of selective sweeps in the Blonde d'Aquitaine cattle breed

Identifying recent positive selection signatures in domesticated animals could provide information on genome response to strong directional selection from domestication and artificial selection and therefore could help in identifying mutations responsible for improved traits. We used genotyping data generated using Illumina's BovineSNP50 Genotyping BeadChips to identify selection signatures in the Blonde d'Aquitaine breed, a well‐muscled French beef breed. For this purpose, we employed a hidden Markov model‐based test, which detects selection by studying local variations in the allele frequency spectrum along the genome, within a single population. Three regions containing selective sweeps were identified. Annotation of genes located within these regions revealed interesting candidate genes. For example, myostatin (also known as GDF8), a known muscle growth factor inhibitor, is located within the selection signature region found on chromosome 2. In addition, we have identified chromosomal regions that show some evidence of selection within QTL regions for economically important traits. The results of this study could help to better understand the mechanisms related to the selection of the Blonde d'Aquitaine breed.     

Signatures of soft sweeps across the Dt1 locus underlying determinate growth habit in soya bean [Glycine max (L.) Merr.]

Determinate growth habit is an agronomically important trait associated with domestication in soya bean. Previous studies have demonstrated that the emergence of determinacy is correlated with artificial selection on four nonsynonymous mutations in the Dt1 gene. To better understand the signatures of the soft sweeps across the Dt1 locus and track the origins of the determinate alleles, we examined patterns of nucleotide variation in Dt1 and the surrounding genomic region of approximately 800 kb. Four local, asymmetrical hard sweeps on four determinate alleles, sized approximately 660, 120, 220 and 150 kb, were identified, which constitute the soft sweeps for the adaptation. These variable‐sized sweeps substantially reflected the strength and timing of selection and indicated that the selection on the alleles had been completed rapidly within half a century. Statistics of EHH, iHS, H12 and H2/H1 based on haplotype data had the power to detect the soft sweeps, revealing distinct signatures of extensive long‐range LD and haplotype homozygosity, and multiple frequent adaptive haplotypes. A haplotype network constructed for Dt1 and a phylogenetic tree based on its extended haplotype block implied independent sources of the adaptive alleles through de novo mutations or rare standing variation in quick succession during the selective phase, strongly supporting multiple origins of the determinacy. We propose that the adaptation of soya bean determinacy is guided by a model of soft sweeps and that this model might be indispensable during crop domestication or evolution.     

Detection of selection signatures of population‐specific genomic regions selected during domestication process in Jinhua pigs

Chinese pigs have been undergoing both natural and artificial selection for thousands of years. Jinhua pigs are of great importance, as they can be a valuable model for exploring the genetic mechanisms linked to meat quality and other traits such as disease resistance, reproduction and production. The purpose of this study was to identify distinctive footprints of selection between Jinhua pigs and other breeds utilizing genome‐wide SNP data. Genotyping by genome reducing and sequencing was implemented in order to perform cross‐population extended haplotype homozygosity to reveal strong signatures of selection for those economically important traits. This work was performed at a 2% genome level, which comprised 152 006 SNPs genotyped in a total of 517 individuals. Population‐specific footprints of selective sweeps were searched for in the genome of Jinhua pigs using six native breeds and three European breeds as reference groups. Several candidate genes associated with meat quality, health and reproduction, such as GH1, CRHR2, TRAF4 and CCK, were found to be overlapping with the significantly positive outliers. Additionally, the results revealed that some genomic regions associated with meat quality, immune response and reproduction in Jinhua pigs have evolved directionally under domestication and subsequent selections. The identified genes and biological pathways in Jinhua pigs showed different selection patterns in comparison with the Chinese and European breeds.     

Did maize domestication and early spread mediate the population genetics of corn leafhopper?

Investigating how crop domestication and early farming mediated crop attributes, distributions, and interactions with antagonists may shed light on today's agricultural pest problems. Crop domestication generally involved artificial selection for traits desirable to early farmers, for example, in creased productivity or yield, and enhanced qualities, though invariably it altered the interactions between crops and insects, and expanded the geographical ranges of crops. Thus, some studies suggest that with crop domestication and spread, insect populations on wild crop ancestors gave rise to pestiferous insect populations on crops. Here, we addressed whether the emergence of corn leafhopper (Dalbulus ma id is) as an agricultural pest may be associated with domestication and early spread of maize (Zea mays mays). We used AFLP markers and mitochondrial COI sequences to assess population genetic structuring and haplotype relationships among corn leafhopper samples from maize and its wild relative Zea diploperennis from multiple locations in Mexico and Argentina. We uncovered seven corn leafhopper haplotypes contained within two haplogroups, one haplogroup containing haplotypes associated with maize and the other containing haplotypes associated with Z. diploperennis in a mountainous habitat. Within the first haplogroup, one haplotype was predominant across Mexican locations, and another across Argentinean locations;both were considered pestiferous. We suggested that the divergence times of the maize-associated haplogroup and of the "pestiferous" haplotypes are correlated with the chronology of maize spread following its domestication. Overall, our results support a hypothesis positing that maize domestication favored corn leafhopper genotypes preadapted for exploiting maize so that they became pestiferous, and that with the geographical expansi on of maize farming, corn leafhopper colonized Z. diploperennis, a host exclusive to secluded habitats that serves as a refuge for archaic corn leafhopper genotypic diversity. Broadly, our results help explain the extents to which crop domestication and early spread may have mediated the emergence of today's agricultural pests.