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.     

Detecting selection signatures in three Iranian sheep breeds

The objective of genome mapping is to achieve valuable insight into the connection between gene variants (genotype) and observed traits (phenotype). Part of that objective is to understand the selective forces that have operated on a population. Finding links between genotype–phenotype changes makes it possible to identify selective sweeps by patterns of genetic variation and linkage disequilibrium. Based on Illumina 50KSNP chip data, two approaches, XP‐EHH (cross‐population extend haplotype homozygosity) and F_ST (fixation index), were carried out in this research to identify selective sweeps in the genome of three Iranian local sheep breeds: Baluchi (n = 86), Lori‐Bakhtiari (n = 45) and Zel (n = 45). Using both methods, 93 candidate genomic regions were identified as harboring putative selective sweeps. Bioinformatics analysis of the genomic regions showed that signatures of selection related to multiple candidate genes, such as HOXB9, HOXB13, ACAN, NPR2, TRIL, AOX1, CSF2, GHR, TNS2, SPAG8, HINT2, ALS2, AAAS, RARG, SYCP2, CAV1, PPP1R3D, PLA2G7, TTLL7 and C20orf10, that play a role in skeletal system and tail, sugar and energy metabolisms, growth, reproduction, immune and nervous system traits. Our findings indicated diverse genomic selection during the domestication of Iranian sheep breeds.     

CRISPR/Cas9‐mediated targeted mutagenesis of GmFT2a delays flowering time in soya bean

Flowering is an indication of the transition from vegetative growth to reproductive growth and has considerable effects on the life cycle of soya bean (Glycine max). In this study, we employed the CRISPR/Cas9 system to specifically induce targeted mutagenesis of GmFT2a, an integrator in the photoperiod flowering pathway in soya bean. The soya bean cultivar Jack was transformed with three sgRNA/Cas9 vectors targeting different sites of endogenous GmFT2a via Agrobacterium tumefaciens‐mediated transformation. Site‐directed mutations were observed at all targeted sites by DNA sequencing analysis. T1‐generation soya bean plants homozygous for null alleles of GmFT2a frameshift mutated by a 1‐bp insertion or short deletion exhibited late flowering under natural conditions (summer) in Beijing, China (N39°58′, E116°20′). We also found that the targeted mutagenesis was stably heritable in the following T2 generation, and the homozygous GmFT2a mutants exhibited late flowering under both long‐day and short‐day conditions. We identified some ‘transgene‐clean’ soya bean plants that were homozygous for null alleles of endogenous GmFT2a and without any transgenic element from the T1 and T2 generations. These ‘transgene‐clean’ mutants of GmFT2a may provide materials for more in‐depth research of GmFT2a functions and the molecular mechanism of photoperiod responses in soya bean. They will also contribute to soya bean breeding and regional introduction.     

Genetic diversity analysis of French goat populations reveals selective sweeps involved in their differentiation

After domestication 11 000 years ago in Asia Minor, the goat followed human migration to Europe and Asia. It was then introduced in Africa and is now raised all over the world. In this study, we exploited a dataset composed of 54 000 SNPs (Illumina goat DNA chip) to analyze the genetic diversity of 223 individuals belonging to eight French breeds (Alpine, Angora, Corse, Fossés, Poitevine, Provençale, Pyrénées and Saanen). Analyses carried out included individual‐based approaches (principal component analysis and population structure) and population‐based approaches (phylogenetic tree constructions). The results of the genetic diversity analyses revealed that French breeds are clearly differentiated, in particular, the Angora breed that originates from south west Asia. The Provençale breed shows a very original genetic pattern that could be the result of ancient admixture. Then, selection signatures were detected by identifying regions of outlying genetic differentiation between populations. Five genomic regions were detected under selection on chromosomes 5, 6, 11, 13 and 20, revealing mainly soft selective sweeps and a few hard selective sweeps and highlighting candidate genes that had been selected for during the evolutionary history of these breeds. Among them, two coat coloration genes (ADAMTS20 and ASIP) and one gene related to milk composition (CSN1S1) were involved.     

Mixed Infection and Natural Spread of ‘Candidatus Phytoplasma asteris’ and Mungbean Yellow Mosaic India Virus Affecting Soya Bean Crop in India

Suspected phytoplasma and virus‐like symptoms of little leaf, yellow mosaic and witches’ broom were recorded on soya bean and two weed species (Digitaria sanguinalis and Parthenium hysterophorus), at experimental fields of Indian Agricultural Research Institute, New Delhi, India, in August–September 2013. The phytoplasma aetiology was confirmed in symptomatic soya bean and both the weed species by direct and nested PCR assays with phytoplasma‐specific universal primer pairs (P1/P6 and R16F2n/R16R2n). One major leafhopper species viz. Empoasca motti Pruthi feeding on symptomatic soya bean plants was also found phytoplasma positive in nested PCR assays. Sequencing BLASTn search analysis and phylogenetic analysis revealed that 16Sr DNA sequences of phytoplasma isolates of soya bean, weeds and leafhoppers had 99% sequence identity among themselves and were related to strains of ‘Candidatus Phytoplasma asteris’. PCR assays with Mungbean yellow mosaic India virus (MYMIV) coat‐protein‐specific primers yielded an amplicon of approximately 770 bp both from symptomatic soya bean and from whiteflies (Bemisia tabaci) feeding on soya bean, confirmed the presence of MYMIV in soya bean and whitefly. Hence, this study suggested the mixed infection of MYMIV and ‘Ca. P. asteris’ with soya bean yellow leaf and witches’ broom syndrome. The two weed species (D. sanguinalis and P. hysterophorus) were recorded as putative alternative hosts for ‘Ca. P. asteris’ soya bean Indian strain. However, the leafhopper E. motti was recorded as putative vector for the identified soya bean phytoplasma isolate, and the whitefly (B. tabaci) was identified as vector of MYMIV which belonged to Asia‐II‐1 genotype.     

Multiple origins of the determinate growth habit in domesticated common bean (Phaseolus vulgaris)

Background and Aims

The actual number of domestications of a crop is one of the key questions in domestication studies. Answers to this question have generally been based on relationships between wild progenitors and domesticated descendants determined with anonymous molecular markers. In this study, this question was investigated by determining the number of instances a domestication phenotype had been selected in a crop species. One of the traits that appeared during domestication of common bean (Phaseolus vulgaris) is determinacy, in which stems end with a terminal inflorescence. It has been shown earlier that a homologue of the arabidopsis TFL1 gene – PvTFL1y – controls determinacy in a naturally occurring variation of common bean.

Methods

Sequence variation was analysed for PvTFL1y in a sample of 46 wild and domesticated accessions that included determinate and indeterminate accessions.

Key Results

Indeterminate types – wild and domesticated – showed only synonymous nucleotide substitutions. Determinate types – observed only among domesticated accessions – showed, in addition to synonymous substitutions, non-synonymous substitutions, indels, a putative intron-splicing failure, a retrotransposon insertion and a deletion of the entire locus. The retrotransposon insertion was observed in 70 % of determinate cultivars, in the Americas and elsewhere. Other determinate mutants had a more restricted distribution in the Americas only, either in the Andean or in the Mesoamerican gene pool of common bean.

Conclusions

Although each of the determinacy haplotypes probably does not represent distinct domestication events, they are consistent with the multiple (seven) domestication pattern in the genus Phaseolus. The predominance of determinacy in the Andean gene pool may reflect domestication of common bean prior to maize introduction in the Andes.     

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.     

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.     

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.     

Evolution and selection of Rhg1, a copy‐number variant nematode‐resistance locus

The soybean cyst nematode (SCN) resistance locus Rhg1 is a tandem repeat of a 31.2 kb unit of the soybean genome. Each 31.2‐kb unit contains four genes. One allele of Rhg1, Rhg1‐b, is responsible for protecting most US soybean production from SCN. Whole‐genome sequencing was performed, and PCR assays were developed to investigate allelic variation in sequence and copy number of the Rhg1 locus across a population of soybean germplasm accessions. Four distinct sequences of the 31.2‐kb repeat unit were identified, and some Rhg1 alleles carry up to three different types of repeat unit. The total number of copies of the repeat varies from 1 to 10 per haploid genome. Both copy number and sequence of the repeat correlate with the resistance phenotype, and the Rhg1 locus shows strong signatures of selection. Significant linkage disequilibrium in the genome outside the boundaries of the repeat allowed the Rhg1 genotype to be inferred using high‐density single nucleotide polymorphism genotyping of 15 996 accessions. Over 860 germplasm accessions were found likely to possess Rhg1 alleles. The regions surrounding the repeat show indications of non‐neutral evolution and high genetic variability in populations from different geographic locations, but without evidence of fixation of the resistant genotype. A compelling explanation of these results is that balancing selection is in operation at Rhg1.     

Recent Selective Sweeps in North American Drosophila melanogaster Show Signatures of Soft Sweeps

Adaptation from standing genetic variation or recurrent de novo mutation in large populations should commonly generate soft rather than hard selective sweeps. In contrast to a hard selective sweep, in which a single adaptive haplotype rises to high population frequency, in a soft selective sweep multiple adaptive haplotypes sweep through the population simultaneously, producing distinct patterns of genetic variation in the vicinity of the adaptive site. Current statistical methods were expressly designed to detect hard sweeps and most lack power to detect soft sweeps. This is particularly unfortunate for the study of adaptation in species such as Drosophila melanogaster, where all three confirmed cases of recent adaptation resulted in soft selective sweeps and where there is evidence that the effective population size relevant for recent and strong adaptation is large enough to generate soft sweeps even when adaptation requires mutation at a specific single site at a locus. Here, we develop a statistical test based on a measure of haplotype homozygosity (H12) that is capable of detecting both hard and soft sweeps with similar power. We use H12 to identify multiple genomic regions that have undergone recent and strong adaptation in a large population sample of fully sequenced Drosophila melanogaster strains from the Drosophila Genetic Reference Panel (DGRP). Visual inspection of the top 50 candidates reveals that in all cases multiple haplotypes are present at high frequencies, consistent with signatures of soft sweeps. We further develop a second haplotype homozygosity statistic (H2/H1) that, in combination with H12, is capable of differentiating hard from soft sweeps. Surprisingly, we find that the H12 and H2/H1 values for all top 50 peaks are much more easily generated by soft rather than hard sweeps. We discuss the implications of these results for the study of adaptation in Drosophila and in species with large census population sizes.     

Rapid Diagnosis of Soya Bean Root Rot Caused by Fusarium culmorum Using a Loop‐Mediated Isothermal Amplification Assay

We report a rapid diagnosis of soya bean (Glycine max L.) root rot caused by Fusarium culmorum, using a loop‐mediated isothermal amplification (LAMP) assay. We used the CYP51C gene sequence to design LAMP assay primers specific for F. culmorum. The LAMP assay amplified the target gene efficiently in 60 min at 63°C. The sensitivity of the assay was 100 pg/μl of genomic DNA. Among the tested soya bean pathogens, a positive colour (sky blue) was only observed in the presence of F. culmorum with the addition of hydroxynaphthol blue (HNB) dye prior to amplification, whereas other species isolates showed no colour change. Suspected diseased soya bean samples collected in the field from Jiangsu, Shandong and Anhui provinces and Beijing were diagnosed successfully using the LAMP assay reported here. This study provides a new and readily available method for rapid diagnosis of soya bean root rot caused by F. culmorum.     

A microRNA allele that emerged prior to apple domestication may underlie fruit size evolution

The molecular genetic mechanisms underlying fruit size remain poorly understood in perennial crops, despite size being an important agronomic trait. Here we show that the expression level of a microRNA gene (miRNA172) influences fruit size in apple. A transposon insertional allele of miRNA172 showing reduced expression associates with large fruit in an apple breeding population, whereas over‐expression of miRNA172 in transgenic apple significantly reduces fruit size. The transposon insertional allele was found to be co‐located with a major fruit size quantitative trait locus, fixed in cultivated apples and their wild progenitor species with relatively large fruit. This finding supports the view that the selection for large size in apple fruit was initiated prior to apple domestication, likely by large mammals, before being subsequently strengthened by humans, and also helps to explain why signatures of genetic bottlenecks and selective sweeps are normally weaker in perennial crops than in annual crops.     

Construction and comparison of three reference‐quality genome assemblies for soybean

We report reference‐quality genome assemblies and annotations for two accessions of soybean (Glycine max) and for one accession of Glycine soja, the closest wild relative of G. max. The G. max assemblies provided are for widely used US cultivars: the northern line Williams 82 (Wm82) and the southern line Lee. The Wm82 assembly improves the prior published assembly, and the Lee and G. soja assemblies are new for these accessions. Comparisons among the three accessions show generally high structural conservation, but nucleotide difference of 1.7 single‐nucleotide polymorphisms (snps) per kb between Wm82 and Lee, and 4.7 snps per kb between these lines and G. soja. snp distributions and comparisons with genotypes of the Lee and Wm82 parents highlight patterns of introgression and haplotype structure. Comparisons against the US germplasm collection show placement of the sequenced accessions relative to global soybean diversity. Analysis of a pan‐gene collection shows generally high conservation, with variation occurring primarily in genomically clustered gene families. We found approximately 40–42 inversions per chromosome between either Lee or Wm82v4 and G. soja, and approximately 32 inversions per chromosome between Wm82 and Lee. We also investigated five domestication loci. For each locus, we found two different alleles with functional differences between G. soja and the two domesticated accessions. The genome assemblies for multiple cultivated accessions and for the closest wild ancestor of soybean provides a valuable set of resources for identifying causal variants that underlie traits for the domestication and improvement of soybean, serving as a basis for future research and crop improvement efforts for this important crop species.     

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.     

Increased copy number of gibberellin 2-oxidase 8 genes reduced trailing growth and shoot length during soybean domestication

Copy number variations (CNVs) play important roles in crop domestication. However, there is only very limited information on the involvement of CNVs in soybean domestication. Trailing growth and long shoots are soybean adaptations for natural habitats but cause lodging that hampers yield in cultivation. Previous studies have focused on Dt1/2 affecting the indeterminate/determinate growth habit, whereas the possible role of the gibberellin pathway remained unclear. In the present study, quantitative trait locus (QTL) mapping of a recombinant inbred population of 460 lines revealed a trailing-growth-and-shoot-length QTL. A CNV region within this QTL was identified, featuring the apical bud-expressed gibberellin 2-oxidase 8A/B, the copy numbers of which were positively correlated with expression levels and negatively with trailing growth and shoot length, and their effects were demonstrated by transgenic soybean and Arabidopsis thaliana. Based on the fixation index, this CNV region underwent intense selection during the initial domestication process.     

A haplotype method detects diverse scenarios of local adaptation from genomic sequence variation

Identifying genomic targets of population‐specific positive selection is a major goal in several areas of basic and applied biology. However, it is unclear how often such selection should act on new mutations versus standing genetic variation or recurrent mutation, and furthermore, favoured alleles may either become fixed or remain variable in the population. Very few population genetic statistics are sensitive to all of these modes of selection. Here, we introduce and evaluate the Comparative Haplotype Identity statistic (χ_MD), which assesses whether pairwise haplotype sharing at a locus in one population is unusually large compared with another population, relative to genomewide trends. Using simulations that emulate human and Drosophila genetic variation, we find that χ_MD is sensitive to a wide range of selection scenarios, and for some very challenging cases (e.g. partial soft sweeps), it outperforms other two‐population statistics. We also find that, as with F_ST, our haplotype approach has the ability to detect surprisingly ancient selective sweeps. Particularly for the scenarios resembling human variation, we find that χ_MD outperforms other frequency‐ and haplotype‐based statistics for soft and/or partial selective sweeps. Applying χ_MD and other between‐population statistics to published population genomic data from D. melanogaster, we find both shared and unique genes and functional categories identified by each statistic. The broad utility and computational simplicity of χ_MD will make it an especially valuable tool in the search for genes targeted by local adaptation.     

Towards the development of a sustainable soya bean‐based feedstock for aquaculture

Soya bean (Glycine max (L.) Merr.) is sought after for both its oil and protein components. Genetic approaches to add value to either component are ongoing efforts in soya bean breeding and molecular biology programmes. The former is the primary vegetable oil consumed in the world. Hence, its primary usage is in direct human consumption. As a means to increase its utility in feed applications, thereby expanding the market of soya bean coproducts, we investigated the simultaneous displacement of marine ingredients in aquafeeds with soya bean‐based protein and a high Omega‐3 fatty acid soya bean oil, enriched with alpha‐linolenic and stearidonic acids, in both steelhead trout (Oncorhynchus mykiss) and Kampachi (Seriola rivoliana). Communicated herein are aquafeed formulations with major reduction in marine ingredients that translates to more total Omega‐3 fatty acids in harvested flesh. Building off of these findings, subsequent efforts were directed towards a genetic strategy that would translate to a prototype design of an optimal identity‐preserved soya bean‐based feedstock for aquaculture, whereby a multigene stack approach for the targeted synthesis of two value‐added output traits, eicosapentaenoic acid and the ketocarotenoid, astaxanthin, were introduced into the crop. To this end, the systematic introduction of seven transgenic cassettes into soya bean, and the molecular and phenotypic evaluation of the derived novel events are described.     

CRISPR/Cas9 and TALENs generate heritable mutations for genes involved in small RNA processing of Glycine max and Medicago truncatula

Processing of double‐stranded RNA precursors into small RNAs is an essential regulator of gene expression in plant development and stress response. Small RNA processing requires the combined activity of a functionally diverse group of molecular components. However, in most of the plant species, there are insufficient mutant resources to functionally characterize each encoding gene. Here, mutations in loci encoding protein machinery involved in small RNA processing in soya bean and Medicago truncatula were generated using the CRISPR/Cas9 and TAL‐effector nuclease (TALEN) mutagenesis platforms. An efficient CRISPR/Cas9 reagent was used to create a bi‐allelic double mutant for the two soya bean paralogous Double‐stranded RNA‐binding2 (GmDrb2a and GmDrb2b) genes. These mutations, along with a CRISPR/Cas9‐generated mutation of the M. truncatula Hua enhancer1 (MtHen1) gene, were determined to be germ‐line transmissible. Furthermore, TALENs were used to generate a mutation within the soya bean Dicer‐like2 gene. CRISPR/Cas9 mutagenesis of the soya bean Dicer‐like3 gene and the GmHen1a gene was observed in the T₀ generation, but these mutations failed to transmit to the T₁ generation. The irregular transmission of induced mutations and the corresponding transgenes was investigated by whole‐genome sequencing to reveal a spectrum of non‐germ‐line‐targeted mutations and multiple transgene insertion events. Finally, a suite of combinatorial mutant plants were generated by combining the previously reported Gmdcl1a, Gmdcl1b and Gmdcl4b mutants with the Gmdrb2ab double mutant. Altogether, this study demonstrates the synergistic use of different genome engineering platforms to generate a collection of useful mutant plant lines for future study of small RNA processing in legume crops.