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.     

Identification of loci controlling adaptation in Chinese soya bean landraces via a combination of conventional and bioclimatic GWAS

Landraces often contain genetic diversity that has been lost in modern cultivars, including alleles that confer enhanced local adaptation. To comprehensively identify loci associated with adaptive traits in soya bean landraces, for example flowering time, a population of 1938 diverse landraces and 97 accessions of the wild progenitor of cultivated soya bean, Glycine soja was genotyped using tGBS^®. Based on 99 085 high‐quality SNPs, landraces were classified into three sub‐populations which exhibit geographical genetic differentiation. Clustering was inferred from STRUCTURE, principal component analyses and neighbour‐joining tree analyses. Using phenotypic data collected at two locations separated by 10 degrees of latitude, 17 trait‐associated SNPs (TASs) for flowering time were identified, including a stable locus Chr12:5914898 and previously undetected candidate QTL/genes for flowering time in the vicinity of the previously cloned flowering genes, E1 and E2. Using passport data associated with the collection sites of the landraces, 27 SNPs associated with adaptation to three bioclimatic variables (temperature, daylength, and precipitation) were identified. A series of candidate flowering genes were detected within linkage disequilibrium (LD) blocks surrounding 12 bioclimatic TASs. Nine of these TASs exhibit significant differences in flowering time between alleles within one or more of the three individual sub‐populations. Signals of selection during domestication and/or subsequent landrace diversification and adaptation were detected at 38 of the 44 flowering and bioclimatic TASs. Hence, this study lays the groundwork to begin breeding for novel environments predicted to arise following global climate change.     

Whole‐genome resequencing uncovers molecular signatures of natural and sexual selection in wild bighorn sheep

The identification of genes influencing fitness is central to our understanding of the genetic basis of adaptation and how it shapes phenotypic variation in wild populations. Here, we used whole‐genome resequencing of wild Rocky Mountain bighorn sheep (Ovis canadensis) to >50‐fold coverage to identify 2.8 million single nucleotide polymorphisms (SNPs) and genomic regions bearing signatures of directional selection (i.e. selective sweeps). A comparison of SNP diversity between the X chromosome and the autosomes indicated that bighorn males had a dramatically reduced long‐term effective population size compared to females. This probably reflects a long history of intense sexual selection mediated by male–male competition for mates. Selective sweep scans based on heterozygosity and nucleotide diversity revealed evidence for a selective sweep shared across multiple populations at RXFP2, a gene that strongly affects horn size in domestic ungulates. The massive horns carried by bighorn rams appear to have evolved in part via strong positive selection at RXFP2. We identified evidence for selection within individual populations at genes affecting early body growth and cellular response to hypoxia; however, these must be interpreted more cautiously as genetic drift is strong within local populations and may have caused false positives. These results represent a rare example of strong genomic signatures of selection identified at genes with known function in wild populations of a nonmodel species. Our results also showcase the value of reference genome assemblies from agricultural or model species for studies of the genomic basis of adaptation in closely related wild taxa.     

Immunogenetic heterogeneity in a widespread ungulate: the European roe deer (Capreolus capreolus)

Understanding how immune genetic variation is shaped by selective and neutral processes in wild populations is of prime importance in both evolutionary biology and epidemiology. The European roe deer (Capreolus capreolus) has considerably expanded its distribution range these last decades, notably by colonizing agricultural landscapes. This range shift is likely to have led to bottlenecks and increased roe deer exposure to a new range of pathogens that until recently predominantly infected humans and domestic fauna. We therefore investigated the historical and contemporary forces that have shaped variability in a panel of genes involved in innate and acquired immunity in roe deer, including Mhc‐Drb and genes encoding cytokines or toll‐like receptors (TLRs). Together, our results suggest that genetic drift is the main contemporary evolutionary force shaping immunogenetic variation within populations. However, in contrast to the classical view, we found that some innate immune genes involved in micropathogen recognition (e.g. Tlrs) continue to evolve dynamically in roe deer in response to pathogen‐mediated positive selection. Most studied Tlrs (Tlr2, Tlr4 and Tlr5) had similarly high levels of amino acid diversity in the three studied populations including one recently established in southwestern France that showed a clear signature of genetic bottleneck. Tlr2 implicated in the recognition of Gram‐positive bacteria in domestic ungulates, showed strong evidence of balancing selection. The high immunogenetic variation revealed here implies that roe deer are able to cope with a wide spectrum of pathogens and to respond rapidly to emerging infectious diseases.     

Genome‐wide analysis reveals artificial selection on coat colour and reproductive traits in Chinese domestic pigs

Pigs from Asia and Europe were independently domesticated from c. 9000 years ago. During this period, strong artificial selection has led to dramatic phenotypic changes in domestic pigs. However, the genetic basis underlying these morphological and behavioural adaptations is relatively unknown, particularly for indigenous Chinese pigs. Here, we performed a genome‐wide analysis to screen 196 regions with selective sweep signals in Tongcheng pigs, which are a typical indigenous Chinese breed. Genes located in these regions have been found to be involved in lipid metabolism, melanocyte differentiation, neural development and other biological processes, which coincide with the evolutionary phenotypic changes in this breed. A synonymous substitution, c.669T>C, in ESR1, which colocalizes with a major quantitative trait locus for litter size, shows extreme differences in allele frequency between Tongcheng pigs and wild boars. Notably, the variant C allele in this locus exhibits high allele frequency in most Chinese populations, suggesting a consequence of positive selection. Five genes (PRM1, PRM2, TNP2, GPR149 and JMJD1C) related to reproductive traits were found to have high haplotype similarity in Chinese breeds. Two selected genes, MITF and EDNRB, are implied to shape the two‐end black colour trait in Tongcheng pig. Subsequent SNP microarray studies of five Chinese white‐spotted breeds displayed a concordant signature at both loci, suggesting that these two genes are responsible for colour variations in Chinese breeds. Utilizing massively parallel sequencing, we characterized the candidate sites that adapt to artificial and environmental selections during the Chinese pig domestication. This study provides fundamental proof for further research on the evolutionary adaptation of Chinese pigs.     

Social effects of territorial neighbours on the timing of spring breeding in North American red squirrels

Organisms can affect one another's phenotypes when they socially interact. Indirect genetic effects occur when an individual's phenotype is affected by genes expressed in another individual. These heritable effects can enhance or reduce adaptive potential, thereby accelerating or reversing evolutionary change. Quantifying these social effects is therefore crucial for our understanding of evolution, yet estimates of indirect genetic effects in wild animals are limited to dyadic interactions. We estimated indirect phenotypic and genetic effects, and their covariance with direct effects, for the date of spring breeding in North American red squirrels (Tamiasciurus hudsonicus) living in an array of territories of varying spatial proximity. Additionally, we estimated indirect effects and the strength of selection at low and high population densities. Social effects of neighbours on the date of spring breeding were different from zero at high population densities but not at low population densities. Indirect phenotypic effects accounted for a larger amount of variation in the date of breeding than differences attributable to the among‐individual variance, suggesting social interactions are important for determining breeding dates. The genetic component to these indirect effects was however not statistically significant. We therefore showcase a powerful and flexible method that will allow researchers working in organisms with a range of social systems to estimate indirect phenotypic and genetic effects, and demonstrate the degree to which social interactions can influence phenotypes, even in a solitary species.     

Genome‐wide association study of six quality traits reveals the association of the TaRPP13L1 gene with flour colour in Chinese bread wheat

Flour colour, kernel hardness, grain protein content and wet gluten content are important quality properties that determine end use in bread wheat. Here, a wheat 90K genotyping assay was used for a genome‐wide association study (GWAS) of the six quality‐related traits in Chinese wheat cultivars in eight environments over four years. A total of 846 significant single nucleotide polymorphisms (SNPs) were identified, explaining approximately 30% of the phenotypic variation on average, and 103 multienvironment‐significant SNPs were detected in more than four environments. Quantitative trait loci (QTL) mapping in the biparent population confirmed some important SNP loci. Moreover, it was determined that some important genes were associated with the six quality traits, including some known functional genes and annotated unknown functional genes. Of the annotated unknown functional genes, it was verified that TaRPP13L1 was associated with flour colour. Wheat cultivars or lines with TaRPP13L1‐B1a showed extremely significantly higher flour redness and lower yellowness than those with TaRPP13L1‐B1b in the Chinese wheat natural population and the doubled haploid (DH) population. Two tetraploid wheat lines with premature stop codons of the TaRPP13L1 gene mutagenized by ethyl methanesulfonate (EMS) showed extremely significantly higher flour redness and lower yellowness than wild type. Our data suggest that the TaRPP13L1 gene plays an important role in modulating wheat flour colour. This study provides useful information for further dissection of the genetic basis of flour colour and also provides valuable genes or genetic loci for marker‐assisted selection to improve the process of breeding quality wheat in China.     

Genetic variants and underlying mechanisms influencing variance heterogeneity in maize

Traditional genetic studies focus on identifying genetic variants associated with the mean difference in a quantitative trait. Because genetic variants also influence phenotypic variation via heterogeneity, we conducted a variance‐heterogeneity genome‐wide association study to examine the contribution of variance heterogeneity to oil‐related quantitative traits. We identified 79 unique variance‐controlling single nucleotide polymorphisms (vSNPs) from the sequences of 77 candidate variance‐heterogeneity genes for 21 oil‐related traits using the Levene test (P < 1.0 × 10^?5). About 30% of the candidate genes encode enzymes that work in lipid metabolic pathways, most of which define clear expression variance quantitative trait loci. Of the vSNPs specifically associated with the genetic variance heterogeneity of oil concentration, 89% can be explained by additional linked mean‐effects genetic variants. Furthermore, we demonstrated that gene × gene interactions play important roles in the formation of variance heterogeneity for fatty acid compositional traits. The interaction pattern was validated for one gene pair (GRMZM2G035341 and GRMZM2G152328) using yeast two‐hybrid and bimolecular fluorescent complementation analyses. Our findings have implications for uncovering the genetic basis of hidden additive genetic effects and epistatic interaction effects, and we indicate opportunities to stabilize efficient breeding and selection of high‐oil maize (Zea mays L.).     

Rabbits – their domestication and molecular genetics of hair coat development and quality

The European rabbit (Oryctolagus cuniculus) is the only representative of its genus living in present‐day Europe and North Africa, and all domestic rabbits are descendants of this one species, which is native to the Iberian Peninsula. There are over 300 breeds of rabbits that differ in size, coat color, length of ears and type of fur. Rabbits are bred for various reasons, such as for laboratory animals and a source of meat, wool and fur, as well as for pets and exhibition animals. The hair coat is a important economic trait of rabbits. Its development and quality are influenced by various factors, both environmental and genetic. The genetic mechanisms underlying its development have not been thoroughly researched. The aim of this review is to discuss the domestication of rabbits and the different aspects of rabbit genetics. A brief review of the properties of rabbit hair coat, hair coat development and hair cycle will be provided, followed by discussion of the factors regulating hair coat development, molecular control of hair coat development and the role of non‐coding RNAs in the regulation of gene expression in the hair follicles of rabbits. Information about genetic regulation of pathways could provide useful tools for improving hair coat quality and be of practical use in rabbit breeding.     

The Ficus erecta genome aids Ceratocystis canker resistance breeding in common fig (F. carica)

Ficus erecta, a wild relative of the common fig (F. carica), is a donor of Ceratocystis canker resistance in fig breeding programmes. Interspecific hybridization followed by recurrent backcrossing is an effective method to transfer the resistance trait from wild to cultivated fig. However, this process is time consuming and labour intensive for trees, especially for gynodioecious plants such as fig. In this study, genome resources were developed for F. erecta to facilitate fig breeding programmes. The genome sequence of F. erecta was determined using single‐molecule real‐time sequencing technology. The resultant assembly spanned 331.6 Mb with 538 contigs and an N50 length of 1.9 Mb, from which 51 806 high‐confidence genes were predicted. Pseudomolecule sequences corresponding to the chromosomes of F. erecta were established with a genetic map based on single nucleotide polymorphisms from double‐digest restriction‐site‐associated DNA sequencing. Subsequent linkage analysis and whole‐genome resequencing identified a candidate gene for the Ceratocystis canker resistance trait. Genome‐wide genotyping analysis enabled the selection of female lines that possessed resistance and effective elimination of the donor genome from the progeny. The genome resources provided in this study will accelerate and enhance disease‐resistance breeding programmes in fig.     

Genome‐wide association and epistasis studies unravel the genetic architecture of sudden death syndrome resistance in soybean

Soybean [Glycine max (L.) Merr.] is an economically important crop that is grown worldwide. Sudden death syndrome (SDS), caused by Fusarium virguliforme, is one of the top yield‐limiting diseases in soybean. However, the genetic basis of SDS resistance, especially with respect to epistatic interactions, is still unclear. To better understand the genetic architecture of soybean SDS resistance, genome‐wide association and epistasis studies were performed using a population of 214 germplasm accessions and 31 914 SNPs from the SoySNP50K Illumina Infinium BeadChip. Twelve loci and 12 SNP–SNP interactions associated with SDS resistance were identified at various time points after inoculation. These additive and epistatic loci together explained 24–52% of the phenotypic variance. Disease‐resistant, pathogenesis‐related and chitin‐ and wound‐responsive genes were identified in the proximity of peak SNPs, including stress‐induced receptor‐like kinase gene 1 (SIK1), which is pinpointed by a trait‐associated SNP and encodes a leucine‐rich repeat‐containing protein. We report that the proportion of phenotypic variance explained by identified loci may be considerably improved by taking epistatic effects into account. This study shows the necessity of considering epistatic effects in soybean SDS resistance breeding using marker‐assisted and genomic selection approaches. Based on our findings, we propose a model for soybean root defense against the SDS pathogen. Our results facilitate identification of the molecular mechanism underlying SDS resistance in soybean, and provide a genetic basis for improvement of soybean SDS resistance through breeding strategies based on additive and epistatic effects.     

A novel quantitative trait locus on chromosome A9 controlling oleic acid content in Brassica napus

One of the most important goals in the breeding of oilseed crops, including Brassica napus, is to improve the quality of edible vegetable oil, which is mainly determined by the seed fatty acid composition, particularly the C18:1 content. Previous studies have indicated that the C18:1 content is a polygenic trait, and no stable quantitative trait loci (QTLs) except for FAD2 have been reported. By performing a GWAS using 375 low erucic acid B. napus accessions genotyped with the Brassica 60K SNP array and constructing a high‐density SNP‐based genetic map of a 150 DH population, we identified a novel QTL on the A9 chromosome. The novel locus could explain 11.25%, 5.72% and 6.29% of phenotypic variation during three consecutive seasons and increased the C18:1 content by approximately 3%–5%. By fine mapping and gene expression analysis, we found three potential candidate genes and verified the fatty acids in a homologous gene mutant of Arabidopsis. A metal ion‐binding protein was found to be the most likely candidate gene in the region. Thus, the C18:1 content can be further increased to about 80% with this novel locus together with FAD2 mutant allele without compromise of agronomic performance. A closely linked marker, BnA129, for this novel QTL (OLEA9) was developed so that we can effectively identify materials with high C18:1 content at an early growth stage by marker‐assisted selection. Our results may also provide new insight for understanding the complex genetic mechanism of fatty acid metabolism.     

Selection on an extreme weapon in the frog‐legged leaf beetle (Sagra femorata)

Biologists have been fascinated with the extreme products of sexual selection for decades. However, relatively few studies have characterized patterns of selection acting on ornaments and weapons in the wild. Here, we measure selection on a wild population of weapon‐bearing beetles (frog‐legged leaf beetles: Sagra femorata) for two consecutive breeding seasons. We consider variation in both weapon size (hind leg length) and in relative weapon size (deviations from the population average scaling relationship between hind leg length and body size), and provide evidence for directional selection on weapon size per se and stabilizing selection on a particular scaling relationship in this population. We suggest that whenever growth in body size is sensitive to external circumstance such as nutrition, then considering deviations from population‐level scaling relationships will better reflect patterns of selection relevant to evolution of the ornament or weapon than will variation in trait size per se. This is because trait‐size versus body‐size scaling relationships approximate underlying developmental reaction norms relating trait growth with body condition in these species. Heightened condition‐sensitive expression is a hallmark of the exaggerated ornaments and weapons favored by sexual selection, yet this plasticity is rarely reflected in the way we think about—and measure—selection acting on these structures in the wild.     

Evolutionary stasis despite selection on a heritable trait in an invasive zooplankton

Invasive species are one of the greatest threats to ecosystems, and there is evidence that evolution plays an important role in the success or failure of invasions. Yet, few studies have measured natural selection and evolutionary responses to selection in invasive species, particularly invasive animals. We quantified the strength of natural selection on the defensive morphology (distal spine) of an invasive zooplankton, Bythotrephes longimanus, in Lake Michigan across multiple months during three growing seasons. We used multiple lines of evidence, including historic and contemporary wild‐captured individuals and palaeoecology of retrieved spines, to assess phenotypic change in distal spine length since invasion. We found evidence of temporally variable selection, with selection for decreased distal spine length early in the growing season and selection for increased distal spine length later in the season. This trend in natural selection is consistent with seasonal changes in the relative strength of non‐gape‐limited and gape‐limited fish predation. Yet, despite net selection for increased distal spine length and a known genetic basis for distal spine length, we observed little evidence of an evolutionary response to selection. Multiple factors likely limit an evolutionary response to selection, including genetic correlations, trade‐offs between components of fitness, and phenotypic plasticity.     

Chromosome-level wild Hevea brasiliensis genome provides new tools for genomic-assisted breeding and valuable loci to elevate rubber yield

The rubber tree (Hevea brasiliensis) is grown in tropical regions and is the major source of natural rubber. Using traditional breeding approaches, the latex yield has increased by sixfold in the last century. However, the underlying genetic basis of rubber yield improvement is largely unknown. Here, we present a high-quality, chromosome-level genome sequence of the wild rubber tree, the first report on selection signatures and a genome-wide association study (GWAS) of its yield traits. Population genomic analysis revealed a moderate population divergence between the Wickham clones and wild accessions. Interestingly, it is suggestive that H. brasiliensis and six relatives of the Hevea genus might belong to the same species. The selective sweep analysis found 361 obvious signatures in the domesticated clones associated with 245 genes. In a 15-year field trial, GWAS identified 155 marker–trait associations with latex yield, in which 326 candidate genes were found. Notably, six genes related to sugar transport and metabolism, and four genes related to ethylene biosynthesis and signalling are associated with latex yield. The homozygote frequencies of the causal nonsynonymous SNPs have been greatly increased under selection, which may have contributed to the fast latex yield improvement during the short domestication history. Our study provides insights into the genetic basis of the latex yield trait and has implications for genomic-assisted breeding by offering valuable resources in this new domesticated crop.     

Animal evolution during domestication: the domesticated fox as a model

We review the evolution of domestic animals, emphasizing the effect of the earliest steps of domestication on its course. Using the first domesticated species, the dog (Canis familiaris), for illustration, we describe the evolutionary peculiarities during the historical domestication, such as the high level and wide range of diversity. We suggest that the process of earliest domestication via unconscious and later conscious selection of human‐defined behavioral traits may accelerate phenotypic variations. The review is based on the results of a long‐term experiment designed to reproduce early mammalian domestication in the silver fox (Vulpes vulpes) selected for tameability or amenability to domestication. We describe changes in behavior, morphology and physiology that appeared in the fox during its selection for tameability, which were similar to those observed in the domestic dog. Based on the data of the fox experiment and survey of relevant data, we discuss the developmental, genetic and possible molecular genetic mechanisms underlying these changes. We ascribe the causative role in evolutionary transformation of domestic animals to the selection for behavior and to the neurospecific regulatory genes it affects.     

Indirect genetic effects in a sex‐limited trait: the case of breeding time in red‐billed gulls

Female reproductive performance can be strongly affected by male care, so that breeding time, a trait expressed only by females, can be seen as one trait determined by both male and female genotypes. Animal model analyses of a 46‐year study of red‐billed gulls (Larus novaehollandiae scopulinus) revealed that laying date was not heritable in females (h² = 0.001 ± 0.030), but significantly so in males (h² = 0.134 ± 0.029). Heritability of breeding time in males probably reflects genetic variability in some other trait such as courtship feeding ability. In line with predictions of evolutionary models incorporating indirect genetic effects, the strong and consistent directional selection for advanced breeding time has not resulted in detectable selection response in males. Our results demonstrate that a female trait is largely determined by genetic characteristics of its mate, and hence, any evolutionary change in red‐billed gull breeding time depends critically on genetic variation in males.     

Pro‐opiomelanocortin gene and melanin‐based colour polymorphism in a reptile

Colour polymorphism is widespread among vertebrates and plays important roles in prey–predator interactions, thermoregulation, social competition, and sexual selection. However, the genetic mechanisms involved in colour variation have been studied mainly in domestic mammals and birds, whereas information on wild animals remains scarce. Interestingly, the pro‐opiomelanocortin gene (POMC) gives rise to melanocortin hormones that trigger melanogenesis (by binding the melanocortin‐1‐receptor; Mc1r) and other physiological and behavioural functions (by binding the melanocortin receptors Mc1‐5rs). Owing to its pleiotropic effect, the POMC gene could therefore account for the numerous covariations between pigmentation and other phenotypic traits. We screened the POMC and Mc1r genes in 107 wild asp vipers (Vipera aspis) that can exhibit four discrete colour morphs (two unpatterned morphs: concolor or melanistic; two patterned morphs: blotched or lined) in a single population. Our study revealed a correlation between a single nucleotide polymorphism situated within the 3′‐untranslated region of the POMC gene and colour variation, whereas Mc1r was not found to be polymorphic. To the best of our knowledge, we disclose for the first time a relationship between a mutation at the POMC gene and coloration in a wild animal, as well as a correlation between a genetic marker and coloration in a snake species. Interestingly, similar mutations within the POMC 3′‐untranslated region are linked to human obesity and alcohol and drug dependence. Combined with our results, this suggests that the 3′‐untranslated region of the POMC gene may play a role in its regulation in distant vertebrates. © 2013 The Linnean Society of London, Biological Journal of the Linnean Society, 2014, 111 , 160–168.