大豆基因组解析与重要农艺性状基因克隆研究进展

自2010年大豆(Glycine max)基因组论文在Nature杂志上正式发表以来, 以中国为主的多国科学家已对众多不同地理来源的野生大豆、农家品种及栽培品种进行了基因组重测序, 并通过比较基因组学等技术手段, 揭示了大豆基因组在进化与驯化过程中的基本规律。利用大豆基因组信息, 科学家已成功克隆了调控大豆生育期、生长(结荚)习性、种子形态性状与组分、结瘤与养分利用效率、生物及非生物胁迫抗(耐)性等重要农艺性状的功能基因, 并初步揭示了其相关的作用机理, 取得了一系列突破性进展, 为今后进一步精细解析大豆基因组与相关性状的基因调控网络奠定了坚实的基础, 同时为通过分子(模块)设计育种来培育新品种创造了必要条件。     

Highly variable patterns of linkage disequilibrium in multiple soybean populations

Prospects for utilizing whole-genome association analysis in autogamous plant populations appear promising due to the reported high levels of linkage disequilibrium (LD). To determine the optimal strategies for implementing association analysis in soybean (Glycine max L. Merr.), we analyzed the structure of LD in three regions of the genome varying in length from 336 to 574 kb. This analysis was conducted in four distinct groups of soybean germplasm: 26 accessions of the wild ancestor of soybean (Glycine soja Seib. et Zucc.); 52 Asian G. max Landraces, the immediate results of domestication from G. soja; 17 Asian Landrace introductions that became the ancestors of North American (N. Am.) cultivars, and 25 Elite Cultivars from N. Am. In G. soja, LD did not extend past 100 kb; however, in the three cultivated G. max groups, LD extended from 90 to 574 kb, likely due to the impacts of domestication and increased self-fertilization. The three genomic regions were highly variable relative to the extent of LD within the three cultivated soybean populations. G. soja appears to be ideal for fine mapping of genes, but due to the highly variable levels of LD in the Landraces and the Elite Cultivars, whole-genome association analysis in soybean may be more difficult than first anticipated.     

Population Structure and Domestication Revealed by High-Depth Resequencing of Korean Cultivated and Wild Soybean Genomes

Despite the importance of soybean as a major crop, genome-wide variation and evolution of cultivated soybeans are largely unknown. Here, we catalogued genome variation in an annual soybean population by high-depth resequencing of 10 cultivated and 6 wild accessions and obtained 3.87 million high-quality single-nucleotide polymorphisms (SNPs) after excluding the sites with missing data in any accession. Nuclear genome phylogeny supported a single origin for the cultivated soybeans. We identified 10-fold longer linkage disequilibrium (LD) in the wild soybean relative to wild maize and rice. Despite the small population size, the long LD and large SNP data allowed us to identify 206 candidate domestication regions with significantly lower diversity in the cultivated, but not in the wild, soybeans. Some of the genes in these candidate regions were associated with soybean homologues of canonical domestication genes. However, several examples, which are likely specific to soybean or eudicot crop plants, were also observed. Consequently, the variation data identified in this study should be valuable for breeding and for identifying agronomically important genes in soybeans. However, the long LD of wild soybeans may hinder pinpointing causal gene(s) in the candidate regions.     

Whole-genome sequence diversity and association analysis of 198 soybean accessions in mini-core collections

We performed whole-genome Illumina resequencing of 198 accessions to examine the genetic diversity and facilitate the use of soybean genetic resources and identified 10 million single nucleotide polymorphisms and 2.8 million small indels. Furthermore, PacBio resequencing of 10 accessions was performed, and a total of 2,033 structure variants were identified. Genetic diversity and structure analysis congregated the 198 accessions into three subgroups (Primitive, World, and Japan) and showed the possibility of a long and relatively isolated history of cultivated soybean in Japan. Additionally, the skewed regional distribution of variants in the genome, such as higher structural variations on the R gene clusters in the Japan group, suggested the possibility of selective sweeps during domestication or breeding. A genome-wide association study identified both known and novel causal variants on the genes controlling the flowering period. Novel candidate causal variants were also found on genes related to the seed coat colour by aligning together with Illumina and PacBio reads. The genomic sequences and variants obtained in this study have immense potential to provide information for soybean breeding and genetic studies that may uncover novel alleles or genes involved in agronomically important traits.     

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.     

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.     

360度群体遗传变异扫描——大豆泛基因组研究

大豆(Glycine max)是重要的油料和蛋白作物, 其丰富的遗传变异为生物学性状挖掘和育种改良提供了重要的资源基础。然而, 单个基因组信息无法全面揭示种质资源的遗传变异, 泛基因组研究为解决这一不足提供了新方案。近日, 中国科学院遗传与发育生物学研究所田志喜和梁承志研究团队从2 898份大豆种质中选取26份代表性材料, 并整合已有的3个基因组, 构建了包含野生和栽培大豆的泛基因组和图基因组(graph-based genome), 鉴定了整个群体的绝大多数结构变异数据集, 确定了大豆种质的核心、非必需和个体特异的基因集。利用这些数据系统地揭示了生育期位点E3的等位基因变异和基因融合事件、种皮颜色基因I的单体型和演化关系以及结构变异对铁离子转运基因表达和地区适应性选择的影响。该研究为作物基因组学研究提供了一个新的模式, 同时将加速推动大豆遗传变异的鉴定、性状解析和种质创新。     

Genomic analyses provide comprehensive insights into the domestication of bast fiber crop ramie (Boehmeria nivea)

Ramie (Boehmeria nivea) is an economically important natural fiber-producing crop that has been cultivated for thousands of years in China; however, the evolution of this crop remains largely unknown. Here, we report a ramie domestication analysis based on genome assembly and resequencing of cultivated and wild accessions. Two chromosome-level genomes representing wild and cultivated ramie were assembled de novo. Numerous structural variations between two assemblies, together with the genetic variations from population resequencing, constituted a comprehensive genomic variation map for ramie. Domestication analysis identified 71 high-confidence selective sweeps comprising 320 predicted genes, and 29 genes from sweeps were associated with fiber growth in the expression. In addition, we identified seven genetic loci associated with the fiber yield trait in the segregated population derived from the crossing of two assembled accessions, and two of which showed an overlap with the selective sweeps. These findings indicated that bast fiber traits were focused on during the domestication history of ramie. This study sheds light on the domestication of ramie and provides a valuable resource for biological and breeding studies of this important crop.     

Genome-Wide Study of Structural Variants in Bovine Holstein,Montbéliarde and Normande Dairy Breeds

High-throughput sequencing technologies have offered in recent years new opportunities to study genome variations. These studies have mostly focused on single nucleotide polymorphisms, small insertions or deletions and on copy number variants. Other structural variants, such as large insertions or deletions, tandem duplications, translocations, and inversions are less well-studied, despite that some have an important impact on phenotypes. In the present study, we performed a large-scale survey of structural variants in cattle. We report the identification of 6,426 putative structural variants in cattle extracted from whole-genome sequence data of 62 bulls representing the three major French dairy breeds. These genomic variants affect DNA segments greater than 50 base pairs and correspond to deletions, inversions and tandem duplications. Out of these, we identified a total of 547 deletions and 410 tandem duplications which could potentially code for CNVs. Experimental validation was carried out on 331 structural variants using a novel high-throughput genotyping method. Out of these, 255 structural variants (77%) generated good quality genotypes and 191 (75%) of them were validated. Gene content analyses in structural variant regions revealed 941 large deletions removing completely one or several genes, including 10 single-copy genes. In addition, some of the structural variants are located within quantitative trait loci for dairy traits. This study is a pan-genome assessment of genomic variations in cattle and may provide a new glimpse into the bovine genome architecture. Our results may also help to study the effects of structural variants on gene expression and consequently their effect on certain phenotypes of interest.     

栽培大豆和野生大豆线粒体基因组密码子使用偏性的比较分析

为分析栽培大豆和野生大豆线粒体基因组的密码子使用特征差异,该文以其线粒体基因组编码序列为研究对象,比较其密码子偏性形成的影响因素和演化过程。结果表明:(1)栽培大豆和野生大豆线粒体基因组编码区的GC含量分别为44.56%和44.58%,说明栽培大豆和野生大豆线粒体编码基因均富含A/T碱基。(2)栽培大豆和野生大豆线粒体基因组密码子第1位、第2位GC含量平均值与第3位GC含量的相关性均呈极显著水平,说明突变在其密码子偏性形成中的作用不可忽略; PR2-plot分析显示,在同义密码子第3位碱基的使用频率上,嘌呤低于嘧啶; Nc-plot分析中Nc比值位于-0.1～0.2区间的基因数占总基因数的95%以上;突变和选择等多重因素共同作用影响了大豆线粒体基因组编码序列密码子使用偏性的形成。(3)有20、21个密码子分别被确定为栽培大豆和野生大豆线粒体基因组编码序列的最优密码子,其中除丝氨酸TCC密码子外均以A或T结尾。综上结果认为,栽培大豆线粒体密码子偏性的形成受选择的影响要高于野生大豆,这可能是栽培大豆由野生大豆经长期人工栽培驯化的结果。     

Differential responses of molecular mechanisms and physiochemical characters in wild and cultivated soybeans against invasion by the pathogenic Fusarium oxysporum Schltdl

Cultivated soybean (Glycine max) was derived from the wild soybean (Glycine soja), which has genetic resources that can be critically important for improving plant stress resistance. However, little information is available pertaining to the molecular and physiochemical comparison between the cultivated and wild soybeans in response to the pathogenic Fusarium oxysporum Schltdl. In this study, we first used comparative phenotypic and paraffin section analyses to indicate that wild soybean is indeed more resistant to F. oxysporum than cultivated soybean. Genome‐wide RNA‐sequencing approach was then used to elucidate the genetic mechanisms underlying the differential physiological and biochemical responses of the cultivated soybean, and its relative, to F. oxysporum. A greater number of genes related to cell wall synthesis and hormone metabolism were significantly altered in wild soybean than in cultivated soybean under F. oxysporum infection. Accordingly, a higher accumulation of lignins was observed in wild soybean than cultivated soybean under F. oxysporum infection. Collectively, these results indicated that secondary metabolites and plant hormones may play a vital role in differentiating the response between cultivated and wild soybeans against the pathogen. These important findings may provide future direction to breeding programs to improve resistance to F. oxysporum in the elite soybean cultivars by taking advantage of the genetic resources within wild soybean germplasm.     

A Comprehensive Analysis of the Cupin Gene Family in Soybean (Glycine max)

Cupin superfamily of proteins, including germin and germin-like proteins (GLPs) from higher plants, is known to play crucial roles in plant development and defense. To date, no systematic analysis has been conducted in soybean (Glycine max) incorporating genome organization, gene structure, expression compendium. In this study, 69 putative Cupin genes were identified from the whole-genome of soybean, which were non-randomly distributed on 17 of the 20 chromosomes. These Gmcupin proteins were phylogenetically clustered into ten distinct subgroups among which the gene structures were highly conserved. Eighteen pairs (52.2%) of duplicate paralogous genes were preferentially retained in duplicated regions of the soybean genome. The distributions of GmCupin genes implied that long segmental duplications contributed significantly to the expansion of the GmCupin gene family. According to the RNA-seq data analysis, most of the Gmcupins were differentially expressed in tissue-specific expression pattern and the expression of some duplicate genes were partially redundant while others showed functional diversity, suggesting the Gmcupins have been retained by substantial subfunctionalization during soybean evolutionary processes. Selective analysis based on single nucleotide polymorphisms (SNPs) in cultivated and wild soybeans revealed sixteen Gmcupins had selected site(s), with all SNPs in Gmcupin10.3 and Gmcupin07.2 genes were selected sites, which implied these genes may have undergone strong selection effects during soybean domestication. Taken together, our results contribute to the functional characterization of Gmcupin genes in soybean.     

Reduced representation genome sequencing reveals patterns of genetic diversity and selection in apple

Identifying DNA sequence variations is a fundamental step towards deciphering the genetic basis of traits of interest.Here,a total of 20 cultivated and 10 wild apples were genotyped using specific-locus amplified fragment sequencing,and 39,635 single nucleotide polymorphisms with no missing genotypes and evenly distributed along the genome were selected to investigate patterns of genome-wide genetic variations between cultivated and wild apples.Overall,wild apples displayed higher levels of genetic diversity than cultivated apples.Linkage disequilibrium(LD) decays were observed quite rapidly in cultivated and wild apples,with an r~2-value below 0.2 at 440 and 280 bp,respectively.Moreover,bidirectional gene flow and different distribution patterns of LD blocks were detected between domesticated and wild apples.Most LD blocks unique to cultivated apples were located within QTL regions controlling fruit quality,thus suggesting that fruit quality had probably undergone selection during apple domestication.The genome of the earliest cultivated apple in China,Nai,was highly similar to that of Malus sieversii,and contained a small portion of genetic material from other wild apple species.This suggested that introgression could have been an important driving force during initial domestication of apple.These findings will facilitate future breeding and genetic dissection of complex traits in apple.     

Proteomic and genetic analysis of glycinin subunits of sixteen soybean genotypes.

We investigated proteomic and genomic profiles of glycinin, a family of major storage proteins in 16 different soybean genotypes consisting of four groups including wild soybean (Glycine soja), unimproved cultivated soybean landraces from Asia (G. max), ancestors of N. American soybean (G. max), and modern soybean (G. max) genotypes. We observed considerable variation in all five glycinin subunits, G1, G2 G3, G4 and G5 using proteomics and genetic analysis. Two-dimensional polyacrylamide gel electrophoresis (2D-PAGE) and mass spectrometry (MS) analysis showed that the wild genotypes had a range of 25-29 glycinin protein spots that included both acidic and basic polypeptides followed by the ancestors with 24-28, modern cultivars with 24-25, and landraces with 17-23 protein spots. Overall, the wild genotypes have a higher number of protein spots when compared to the other three genotypes. Major variation was observed in acidic polypeptides of G3, G4 and G5 compared to G1 and G2, and minor variation was observed in basic polypeptides of all subunits. Our data indicated that there are major variations of glycinin subunits between wild and cultivated genotypes rather than within the same groups. Based on Southern blot DNA analysis, we observed genetic polymorphisms in group I genes (G1, G2, and G3) between and within the four genotype groups, but not in group II genes (G4 and G5). This is the first study reporting the comparative analysis of glycinin in a diverse set of soybean genotypes using combined proteomic and genetic analysis.     

Ten years of AFLP in ecology and evolution: why so few animals?

Researchers in the field of molecular ecology and evolution require versatile and low-cost genetic typing methods. The AFLP (amplified fragment length polymorphism) method was introduced 10 years ago and shows many features that fulfil these requirements. With good quality genomic DNA at hand, it is relatively easy to generate anonymous multilocus DNA profiles in most species and the start-up time before data can be generated is often less than a week. Built-in dynamic, yet simple modifications make it possible to find a protocol suitable to the genome size of the species and to screen thousands of loci in hundreds of individuals for a relatively low cost. Until now, the method has primarily been applied in studies of plants, bacteria and fungi, with a strong bias towards economically important cultivated species and their pests. In this review we identify a number of research areas in the study of wild species of animals where the AFLP method, presently very much underused, should be a very valuable tool. These aspects include classical problems such as studies of population genetic structure and phylogenetic reconstructions, and also new challenges such as finding markers for genes governing adaptations in wild populations and modifications of the protocol that makes it possible to measure expression variation of multiple genes (cDNA-AFLP) and the distribution of DNA methylation. We hope this review will help molecular ecologists to identify when AFLP is likely to be superior to other more established methods, such as microsatellites, SNP (single nucleotide polymorphism) analyses and multigene DNA sequencing.     

Draft Genome Sequence of Eggplant (Solanum melongena L.): the Representative Solanum Species Indigenous to the Old World

Unlike other important Solanaceae crops such as tomato, potato, chili pepper, and tobacco, all of which originated in South America and are cultivated worldwide, eggplant (Solanum melongena L.) is indigenous to the Old World and in this respect it is phylogenetically unique. To broaden our knowledge of the genomic nature of solanaceous plants further, we dissected the eggplant genome and built a draft genome dataset with 33,873 scaffolds termed SME_r2.5.1 that covers 833.1 Mb, ca. 74% of the eggplant genome. Approximately 90% of the gene space was estimated to be covered by SME_r2.5.1 and 85,446 genes were predicted in the genome. Clustering analysis of the predicted genes of eggplant along with the genes of three other solanaceous plants as well as Arabidopsis thaliana revealed that, of the 35,000 clusters generated, 4,018 were exclusively composed of eggplant genes that would perhaps confer eggplant-specific traits. Between eggplant and tomato, 16,573 pairs of genes were deduced to be orthologous, and 9,489 eggplant scaffolds could be mapped onto the tomato genome. Furthermore, 56 conserved synteny blocks were identified between the two species. The detailed comparative analysis of the eggplant and tomato genomes will facilitate our understanding of the genomic architecture of solanaceous plants, which will contribute to cultivation and further utilization of these crops.     

中国野生大豆(Glycine soja)遗传资源主要形态、遗传变异和结构

介绍了我国野生大豆遗传资源主要的形态类型、物种内遗传关系和遗传结构。进化的瓶颈不仅发生在由野生大豆到栽培大豆,也以另一种方式"分化瓶颈"出现于同性状的不同表型类型之间。野生大豆种内种子大小类型是否存在遗传分化?野生和半野生大豆的边界在哪?半野生大豆如何产生的?半野生大豆遗传上密切于栽培种还是野生种?百粒重3～4g的小粒半野生大豆与百粒重8.5g以上的特大粒半野生大豆是否有遗传差异?百粒重8.5g以上的特大粒半野生大豆是否属于栽培大豆?野生大豆的种皮色和种子大小哪个更能反映进化程度?栽培大豆基因是否已经渗入到野生大豆?对这些在学术界常年存在的疑问本文介绍了我们的研究答案。我们认为"真"半野生大豆不存在于现在中国半野生资源收集品中;一些野生大豆中的白花、灰毛、无泥膜性状来源于栽培大豆的基因渗透。     

Setaria genome sequencing: an overview

The genus Setaria includes two important C₄ Panicoid grass species, namely S. italica (cultivated) and S. viridis (weed; wild ancestor), which together represent an appropriate model system for architectural, physiological, evolutionary, and genomic studies in related grasses. It is a diploid, inbreeder, self-fertile annual cereal grass having short life cycle and minimal growth requirements. There close relatedness to biofuel crops like switch grass and napier grass further signifies their importance. Further, foxtail millet is an important food and fodder grain crop grown in arid and semi-arid regions in many parts of the world. Therefore, an increasing interest in these species has led to a gradual accumulation and development of genomic data and genetic resources. Setaria genome sequencing is an outcome of such endeavors. These sequencing efforts uncovered several distinctive attributes of Setaria genome that may help in understanding its physiology, evolution and adaptation. This will not only aid in comparative genomics studies of Setaria and related crops including bioenergy grasses but also help in rapid advancements of genomics information for developing varieties with superior traits either through marker-assisted selection (MAS) or using transgenic approaches in these crops.