A chromosome‐scale reference genome and genome‐wide genetic variations elucidate adaptation in yak

Yak is an important livestock animal for the people indigenous to the harsh, oxygen‐limited Qinghai‐Tibetan Plateau and Hindu Kush ranges of the Himalayas. The yak genome was sequenced in 2012, but its assembly was fragmented because of the inherent limitations of the Illumina sequencing technology used to analyse it. An accurate and complete reference genome is essential for the study of genetic variations in this species. Long‐read sequences are more complete than their short‐read counterparts and have been successfully applied towards high‐quality genome assembly for various species. In this study, we present a high‐quality chromosome‐scale yak genome assembly (BosGru_PB_v1.0) constructed with long‐read sequencing and chromatin interaction technologies. Compared to an existing yak genome assembly (BosGru_v2.0), BosGru_PB_v1.0 shows substantially improved chromosome sequence continuity, reduced repetitive structure ambiguity, and gene model completeness. To characterize genetic variation in yak, we generated de novo genome assemblies based on Illumina short reads for seven recognized domestic yak breeds in Tibet and Sichuan and one wild yak from Hoh Xil. We compared these eight assemblies to the BosGru_PB_v1.0 genome, obtained a comprehensive map of yak genetic diversity at the whole‐genome level, and identified several protein‐coding genes absent from the BosGru_PB_v1.0 assembly. Despite the genetic bottleneck experienced by wild yak, their diversity was nonetheless higher than that of domestic yak. Here, we identified breed‐specific sequences and genes by whole‐genome alignment, which may facilitate yak breed identification.     

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.     

Accumulation of deleterious mutations in the domestic yak genome

Deleterious mutations play an important functional role, affecting trait phenotypes in ways that decrease the fitness of organisms. Estimating the frequency of occurrence and abundance has been a topic of much interest, especially in crops and livestock. The processes of domestication and breeding allow deleterious mutations to persist at high frequency, and identifying such deleterious mutations is particularly important for breed improvement. Here, we assessed genome‐wide patterns of deleterious variation in 59 domestic and 13 wild yaks using genome resequencing data. Based on the intersection of results given by three methods (provean , polyphen 2 and sift 4g ), we identified 3187 putative deleterious mutation sites affecting 2586 genes in domestic yaks and 2067 affecting 1701 genes in wild yaks. Multiple lines of evidence indicate a significant increase in the load of deleterious mutations in domesticated yaks compared to wild yaks. Private deleterious genes were found to be associated with the perception of smell and detection of chemical stimulus. We also identified 36 genes related to Mendelian genetic diseases involved in sensory perception, skeletal development and the nervous and immune systems. This study not only adds to the understanding of the genetic basis of yak domestication but also provides a rich catalog of variants that will facilitate future breeding‐related research on the yak genome and on other bovid species.     

Structural Variants Selected during Yak Domestication Inferred from Long-Read Whole-Genome Sequencing

Structural variants (SVs) represent an important genetic resource for both natural and artificial selection. Here we present a chromosome-scale reference genome for domestic yak (Bos grunniens) that has longer contigs and scaffolds (N50 44.72 and 114.39 Mb, respectively) than reported for any other ruminant genome. We further obtained long-read resequencing data for 6 wild and 23 domestic yaks and constructed a genetic SV map of 372,220 SVs that covers the geographic range of the yaks. The majority of the SVs contains repetitive sequences and several are in or near genes. By comparing SVs in domestic and wild yaks, we identified genes that are predominantly related to the nervous system, behavior, immunity, and reproduction and may have been targeted by artificial selection during yak domestication. These findings provide new insights in the domestication of animals living at high altitude and highlight the importance of SVs in animal domestication.     

Genetic structure in Red Junglefowl (Gallus gallus) populations: Strong spatial patterns in the wild ancestors of domestic chickens in a core distribution range

Red Junglefowl (Gallus gallus) are among the few remaining ancestors of an extant domesticated livestock species, the domestic chicken, that still occur in the wild. Little is known about genetic diversity, population structure, and demography of wild Red Junglefowl in their natural habitats. Extinction threats from habitat loss or genetic alteration from domestic introgression exacerbate further the conservation status of this progenitor species. In a previous study, we reported extraordinary adaptive genetic variation in the MHC B‐locus in wild Red Junglefowl and no evidence of allelic introgression between wild and domestic chickens was observed. In this study, we characterized spatial genetic variation and population structure in naturally occurring populations of Red Junglefowl in their core distribution range in South Central Vietnam. A sample of 212 Red Junglefowl was obtained from geographically and ecologically diverse habitats across an area of 250 × 350 km. We used amplified fragment‐length polymorphism markers obtained from 431 loci to determine whether genetic diversity and population structure varies. We found that Red Junglefowl are widely distributed but form small and isolated populations. Strong spatial genetic patterns occur at both local and regional scales. At local scale, population stratification can be identified to approximately 5 km. At regional scale, we identified distinct populations of Red Junglefowl in the southern lowlands, northern highlands, and eastern coastal portions of the study area. Both local and long‐distance genetic patterns observed in wild Red Junglefowl may reflect the species’ ground‐dwelling and territorial characteristics, including dispersal barriers imposed by the Annamite Mountain Range. Spatially explicit analyses with neutral genetic markers can be highly informative and here elevates the conservation profile of the wild ancestors of domesticated chickens.     

A draft genome sequence of the pulse crop chickpea (Cicer arietinum L.)

Cicer arietinum L. (chickpea) is the third most important food legume crop. We have generated the draft sequence of a desi‐type chickpea genome using next‐generation sequencing platforms, bacterial artificial chromosome end sequences and a genetic map. The 520‐Mb assembly covers 70% of the predicted 740‐Mb genome length, and more than 80% of the gene space. Genome analysis predicts the presence of 27 571 genes and 210 Mb as repeat elements. The gene expression analysis performed using 274 million RNA‐Seq reads identified several tissue‐specific and stress‐responsive genes. Although segmental duplicated blocks are observed, the chickpea genome does not exhibit any indication of recent whole‐genome duplication. Nucleotide diversity analysis provides an assessment of a narrow genetic base within the chickpea cultivars. We have developed a resource for genetic markers by comparing the genome sequences of one wild and three cultivated chickpea genotypes. The draft genome sequence is expected to facilitate genetic enhancement and breeding to develop improved chickpea varieties.     

A comprehensive biomedical variant catalogue based on whole genome sequences of 582 dogs and eight wolves

The domestic dog serves as an excellent model to investigate the genetic basis of disease. More than 400 heritable traits analogous to human diseases have been described in dogs. To further canine medical genetics research, we established the Dog Biomedical Variant Database Consortium (DBVDC) and present a comprehensive list of functionally annotated genome variants that were identified with whole genome sequencing of 582 dogs from 126 breeds and eight wolves. The genomes used in the study have a minimum coverage of 10× and an average coverage of ~24×. In total, we identified 23 133 692 single‐nucleotide variants (SNVs) and 10 048 038 short indels, including 93% undescribed variants. On average, each individual dog genome carried ～4.1 million single‐nucleotide and ~1.4 million short‐indel variants with respect to the reference genome assembly. About 2% of the variants were located in coding regions of annotated genes and loci. Variant effect classification showed 247 141 SNVs and 99 562 short indels having moderate or high impact on 11 267 protein‐coding genes. On average, each genome contained heterozygous loss‐of‐function variants in 30 potentially embryonic lethal genes and 97 genes associated with developmental disorders. More than 50 inherited disorders and traits have been unravelled using the DBVDC variant catalogue, enabling genetic testing for breeding and diagnostics. This resource of annotated variants and their corresponding genotype frequencies constitutes a highly useful tool for the identification of potential variants causative for rare inherited disorders in dogs.     

Whole-Genome Resequencing of Worldwide Wild and Domestic Sheep Elucidates Genetic Diversity,Introgression, and Agronomically Important Loci

Domestic sheep and their wild relatives harbor substantial genetic variants that can form the backbone of molecular breeding, but their genome landscapes remain understudied. Here, we present a comprehensive genome resource for wild ovine species, landraces and improved breeds of domestic sheep, comprising high-coverage (∼16.10×) whole genomes of 810 samples from 7 wild species and 158 diverse domestic populations. We detected, in total, ∼121.2 million single nucleotide polymorphisms, ∼61 million of which are novel. Some display significant (P < 0.001) differences in frequency between wild and domestic species, or are private to continent-wide or individual sheep populations. Retained or introgressed wild gene variants in domestic populations have contributed to local adaptation, such as the variation in the HBB associated with plateau adaptation. We identified novel and previously reported targets of selection on morphological and agronomic traits such as stature, horn, tail configuration, and wool fineness. We explored the genetic basis of wool fineness and unveiled a novel mutation (chr25: T7,068,586C) in the 3′-UTR of IRF2BP2 as plausible causal variant for fleece fiber diameter. We reconstructed prehistorical migrations from the Near Eastern domestication center to South-and-Southeast Asia and found two main waves of migrations across the Eurasian Steppe and the Iranian Plateau in the Early and Late Bronze Ages. Our findings refine our understanding of genome variation as shaped by continental migrations, introgression, adaptation, and selection of sheep.     

Genome‐wide survey of artificial mutations induced by ethyl methanesulfonate and gamma rays in tomato

Genome‐wide mutations induced by ethyl methanesulfonate (EMS) and gamma irradiation in the tomato Micro‐Tom genome were identified by a whole‐genome shotgun sequencing analysis to estimate the spectrum and distribution of whole‐genome DNA mutations and the frequency of deleterious mutations. A total of ~370 Gb of paired‐end reads for four EMS‐induced mutants and three gamma‐ray‐irradiated lines as well as a wild‐type line were obtained by next‐generation sequencing technology. Using bioinformatics analyses, we identified 5920 induced single nucleotide variations and insertion/deletion (indel) mutations. The predominant mutations in the EMS mutants were C/G to T/A transitions, while in the gamma‐ray mutants, C/G to T/A transitions, A/T to T/A transversions, A/T to G/C transitions and deletion mutations were equally common. Biases in the base composition flanking mutations differed between the mutagenesis types. Regarding the effects of the mutations on gene function, >90% of the mutations were located in intergenic regions, and only 0.2% were deleterious. In addition, we detected 1 140 687 spontaneous single nucleotide polymorphisms and indel polymorphisms in wild‐type Micro‐Tom lines. We also found copy number variation, deletions and insertions of chromosomal segments in both the mutant and wild‐type lines. The results provide helpful information not only for mutation research, but also for mutant screening methodology with reverse‐genetic approaches.     

野牦牛线粒体基因组序列测定及其系统进化

野牦牛属高寒地区的特有物种,是我国最珍贵的野生动物遗传资源之一,已被列为国家一级重点保护动物。对野牦牛mtDNA进行全序列测定和结构分析,并基于线粒体基因组序列对其系统发生进行了探讨。结果表明:(1)野牦牛线粒体基因组全序列的大小为16 322 bp,整个基因组由37个编码基因和D-loop区组成;22个tRNA基因序列长度为1 524 bp、2个RNA基因序列长度为2 528 bp、13个编码蛋白基因序列长度为11420 bp、D-loop区长度为892 bp。基因组中无间隔序列,基因间排列紧密,基因内无内含子。(2)野牦牛具有较丰富的遗传多样性。(3)分子系统发生关系显示牦牛为牛亚科中的一个独立属,即牦牛属(Poephagus),牦牛属包括家牦牛(Poephagus grunniens)和野牦牛(Poephagus mutus)2个种。野牦牛线粒体基因组全序列的获得和结构解析对研究牦牛的起源、演化和分类,以及野牦牛遗传资源的保护、开发和利用均具有重要的理论和实际意义。     

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.     

Genomewide discovery of DNA polymorphisms in rice cultivars with contrasting drought and salinity stress response and their functional relevance

Next‐generation sequencing technologies provide opportunities to understand the genetic basis of phenotypic differences, such as abiotic stress response, even in the closely related cultivars via identification of large number of DNA polymorphisms. We performed whole‐genome resequencing of three rice cultivars with contrasting responses to drought and salinity stress (sensitive IR64, drought‐tolerant Nagina 22 and salinity‐tolerant Pokkali). More than 356 million 90‐bp paired‐end reads were generated, which provided about 85% coverage of the rice genome. Applying stringent parameters, we identified a total of 1 784 583 nonredundant single‐nucleotide polymorphisms (SNPs) and 154 275 InDels between reference (Nipponbare) and the three resequenced cultivars. We detected 401 683 and 662 509 SNPs between IR64 and Pokkali, and IR64 and N22 cultivars, respectively. The distribution of DNA polymorphisms was found to be uneven across and within the rice chromosomes. One‐fourth of the SNPs and InDels were detected in genic regions, and about 3.5% of the total SNPs resulted in nonsynonymous changes. Large‐effect SNPs and InDels, which affect the integrity of the encoded protein, were also identified. Further, we identified DNA polymorphisms present in the differentially expressed genes within the known quantitative trait loci. Among these, a total of 548 SNPs in 232 genes, located in the conserved functional domains, were identified. The data presented in this study provide functional markers and promising target genes for salinity and drought tolerance and present a valuable resource for high‐throughput genotyping and molecular breeding for abiotic stress traits in rice.     

Towards a whole‐genome sequence for rye (Secale cereale L.)

We report on a whole‐genome draft sequence of rye (Secale cereale L.). Rye is a diploid Triticeae species closely related to wheat and barley, and an important crop for food and feed in Central and Eastern Europe. Through whole‐genome shotgun sequencing of the 7.9‐Gbp genome of the winter rye inbred line Lo7 we obtained a de novo assembly represented by 1.29 million scaffolds covering a total length of 2.8 Gbp. Our reference sequence represents nearly the entire low‐copy portion of the rye genome. This genome assembly was used to predict 27 784 rye gene models based on homology to sequenced grass genomes. Through resequencing of 10 rye inbred lines and one accession of the wild relative S. vavilovii, we discovered more than 90 million single nucleotide variants and short insertions/deletions in the rye genome. From these variants, we developed the high‐density Rye600k genotyping array with 600 843 markers, which enabled anchoring the sequence contigs along a high‐density genetic map and establishing a synteny‐based virtual gene order. Genotyping data were used to characterize the diversity of rye breeding pools and genetic resources, and to obtain a genome‐wide map of selection signals differentiating the divergent gene pools. This rye whole‐genome sequence closes a gap in Triticeae genome research, and will be highly valuable for comparative genomics, functional studies and genome‐based breeding in rye.     

Development,validation and genetic analysis of a large soybean SNP genotyping array

Cultivated soybean (Glycine max) suffers from a narrow germplasm relative to other crop species, probably because of under‐use of wild soybean (Glycine soja) as a breeding resource. Use of a single nucleotide polymorphism (SNP) genotyping array is a promising method for dissecting cultivated and wild germplasms to identify important adaptive genes through high‐density genetic mapping and genome‐wide association studies. Here we describe a large soybean SNP array for use in diversity analyses, linkage mapping and genome‐wide association analyses. More than four million high‐quality SNPs identified from high‐depth genome re‐sequencing of 16 soybean accessions and low‐depth genome re‐sequencing of 31 soybean accessions were used to select 180 961 SNPs for creation of the Axiom^® SoyaSNP array. Validation analysis for a set of 222 diverse soybean lines showed that 170 223 markers were of good quality for genotyping. Phylogenetic and allele frequency analyses of the validation set data indicated that accessions showing an intermediate morphology between cultivated and wild soybeans collected in Korea were natural hybrids. More than 90 unanchored scaffolds in the current soybean reference sequence were assigned to chromosomes using this array. Finally, dense average spacing and preferential distribution of the SNPs in gene‐rich chromosomal regions suggest that this array may be suitable for genome‐wide association studies of soybean germplasm. Taken together, these results suggest that use of this array may be a powerful method for soybean genetic analyses relating to many aspects of soybean breeding.     

Novel Y‐chromosome short tandem repeats in Sus scrofa and their variation in European wild boar and domestic pig populations

Y‐chromosome markers are important tools for studying male‐specific gene flow within and between populations, hybridization patterns and kinship. However, their use in non‐human mammals is often hampered by the lack of Y‐specific polymorphic markers. We identified new male‐specific short tandem repeats (STRs) in Sus scrofa using the available genome sequence. We selected four polymorphic loci (5–10 alleles per locus), falling in one duplicated and two single‐copy regions. A total of 32 haplotypes were found by screening 211 individuals from eight wild boar populations across Europe and five domestic pig populations. European wild boar were characterized by significantly higher levels of haplotype diversity compared to European domestic pigs (H_D = 0.904 ± 0.011 and H_D = 0.491 ± 0.077 respectively). Relationships among STR haplotypes were investigated by combining them with single nucleotide polymorphisms at two linked genes (AMELY and UTY) in a network analysis. A differentiation between wild and domestic populations was observed (F_ST = 0.229), with commercial breeds sharing no Y haplotype with the sampled wild boar. Similarly, a certain degree of geographic differentiation was observed across Europe, with a number of local private haplotypes and high diversity in northern populations. The described Y‐chromosome markers can be useful to track male inheritance and gene flow in wild and domestic populations, promising to provide insights into evolutionary and population genetics in Sus scrofa.     

SNP deserts of Asian cultivated rice: genomic regions under domestication

When performing a genome‐wide comparison between indica (93‐11) and japonica (Nipponbare), we find 8% of the genome, which have an extremely low SNP rate (< 1 SNP/kb). Inside these ‘SNP deserts’, experimentally confirmed genes show increased K_a/K_s that indicate adaptive selection. To further elucidate this connection, we survey the level and pattern of genetic variation in both cultivated and wild rice groups, using 155 noncoding regions located within SNP deserts. The results suggest that cultivated rice has greatly reduced genetic variation within SNP deserts as compared to either the nondesert or corresponding genomic regions in wild rice. Consistent with this reduction in genetic variation, we find a biased distribution of derived allele frequency in the cultivated group, indicative of positive selection. Furthermore, over half of the confirmed, domestication‐related genes are found within SNP deserts, also suggesting that SNP deserts are strongly related to domestication, and might be the key sites in the process of domestication.     

Diversity and population structure of northern switchgrass as revealed through exome capture sequencing

Panicum virgatum L. (switchgrass) is a polyploid, perennial grass species that is native to North America, and is being developed as a future biofuel feedstock crop. Switchgrass is present primarily in two ecotypes: a northern upland ecotype, composed of tetraploid and octoploid accessions, and a southern lowland ecotype, composed of primarily tetraploid accessions. We employed high‐coverage exome capture sequencing (~2.4 Tb) to genotype 537 individuals from 45 upland and 21 lowland populations. From these data, we identified ~27 million single‐nucleotide polymorphisms (SNPs), of which 1 590 653 high‐confidence SNPs were used in downstream analyses of diversity within and between the populations. From the 66 populations, we identified five primary population groups within the upland and lowland ecotypes, a result that was further supported through genetic distance analysis. We identified conserved, ecotype‐restricted, non‐synonymous SNPs that are predicted to affect the protein function of CONSTANS (CO) and EARLY HEADING DATE 1 (EHD1), key genes involved in flowering, which may contribute to the phenotypic differences between the two ecotypes. We also identified, relative to the near‐reference Kanlow population, 17 228 genes present in more copies than in the reference genome (up‐CNVs), 112 630 genes present in fewer copies than in the reference genome (down‐CNVs) and 14 430 presence/absence variants (PAVs), affecting a total of 9979 genes, including two upland‐specific CNV clusters. In total, 45 719 genes were affected by an SNP, CNV, or PAV across the panel, providing a firm foundation to identify functional variation associated with phenotypic traits of interest for biofuel feedstock production.     

Copy number variation in Fayoumi and Leghorn chickens analyzed using array comparative genomic hybridization

Copy number variation refers to regions along chromosomes that harbor a type of structural variation, such as duplications or deletions. Copy number variants (CNVs) play a role in many important traits as well as in genetic diversity. Previous analyses of chickens using array comparative genomic hybridizations or single‐nucleotide polymorphism chip assays have been performed on various breeds and genetic lines to discover CNVs. In this study, we assessed individuals from two highly inbred (inbreeding coefficiency > 99.99%) lines, Leghorn G‐B2 and Fayoumi M15.2, to discover novel CNVs in chickens. These lines have been previously studied for disease resistance, and to our knowledge, this represents the first global assessment of CNVs in the Fayoumi breed. Genomic DNA from individuals was examined using the Agilent chicken 244 K comparative genomic hybridization array and quantitative PCR. We identified a total of 273 CNVs overall, with 112 CNVs being novel and not previously reported. Quantitative PCR using the standard curve method validated a subset of our array data. Through enrichment analysis of genes within CNV regions, we observed multiple chromosomes, terms and pathways that were significantly enriched, largely dealing with the major histocompatibility complex and immune responsiveness. Using an additional round of computational and statistical analysis with a different bioinformatic pipeline, we identified 43 CNVs among these as high‐confidence regions, 14 of which were found to be novel. We further compared and contrasted individuals of the two inbred lines to discover regions that have a significant difference in copy number between lines. A total of 40 regions had significant deletions or duplications between the lines. Gene Ontology analysis of genomic regions containing CNVs between lines also was performed. This between‐line candidate CNV list will be useful in studies with these two unique genetic lines, which may harbor variations that underlie quantitative trait loci for disease resistance and other important traits. Through the global discovery of novel CNVs in chicken, these data also provide resources for further genetic and functional genomics studies.     

Association mapping of morphological traits in wild and captive zebra finches: reliable within,but not between populations

Identifying causal genetic variants underlying heritable phenotypic variation is a long‐standing goal in evolutionary genetics. We previously identified several quantitative trait loci (QTL) for five morphological traits in a captive population of zebra finches (Taeniopygia guttata) by whole‐genome linkage mapping. We here follow up on these studies with the aim to narrow down on the quantitative trait variants (QTN) in one wild and three captive populations. First, we performed an association study using 672 single nucleotide polymorphisms (SNPs) within candidate genes located in the previously identified QTL regions in a sample of 939 wild‐caught zebra finches. Then, we validated the most promising SNP–phenotype associations (n = 25 SNPs) in 5228 birds from four populations. Genotype–phenotype associations were generally weak in the wild population, where linkage disequilibrium (LD) spans only short genomic distances. In contrast, in captive populations, where LD blocks are large, apparent SNP effects on morphological traits (i.e. associations) were highly repeatable with independent data from the same population. Most of those SNPs also showed significant associations with the same trait in other captive populations, but the direction and magnitude of these effects varied among populations. This suggests that the tested SNPs are not the causal QTN but rather physically linked to them, and that LD between SNPs and causal variants differs between populations due to founder effects. While the identification of QTN remains challenging in nonmodel organisms, we illustrate that it is indeed possible to confirm the location and magnitude of QTL in a population with stable linkage between markers and causal variants.     

The draft genome of a wild barley genotype reveals its enrichment in genes related to biotic and abiotic stresses compared to cultivated barley

Wild barley (Hordeum spontaneum) is the progenitor of cultivated barley (Hordeum vulgare) and provides a rich source of genetic variations for barley improvement. Currently, the genome sequences of wild barley and its differences with cultivated barley remain unclear. In this study, we report a high‐quality draft assembly of wild barley accession (AWCS276; henceforth named as WB1), which consists of 4.28 Gb genome and 36 395 high‐confidence protein‐coding genes. BUSCO analysis revealed that the assembly included full lengths of 95.3% of the 956 single‐copy plant genes, illustrating that the gene‐containing regions have been well assembled. By comparing with the genome of the cultivated genotype Morex, it is inferred that the WB1 genome contains more genes involved in resistance and tolerance to biotic and abiotic stresses. The presence of the numerous WB1‐specific genes indicates that, in addition to enhance allele diversity for genes already existing in the cultigen, exploiting the wild barley taxon in breeding should also allow the incorporation of novel genes. Furthermore, high levels of genetic variation in the pericentromeric regions were detected in chromosomes 3H and 5H between the wild and cultivated genotypes, which may be the results of domestication. This H. spontaneum draft genome assembly will help to accelerate wild barley research and be an invaluable resource for barley improvement and comparative genomics research.