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.     

New chloroplast microsatellite markers suitable for assessing genetic diversity of Lolium perenne and other related grass species

Background and Aims

Lolium perenne (perennial ryegrass) is the most important forage grass species of temperate regions. We have previously released the chloroplast genome sequence of L. perenne ‘Cashel’. Here nine chloroplast microsatellite markers are published, which were designed based on knowledge about genetically variable regions within the L. perenne chloroplast genome. These markers were successfully used for characterizing the genetic diversity in Lolium and different grass species.

Methods

Chloroplast genomes of 14 Poaceae taxa were screened for mononucleotide microsatellite repeat regions and primers designed for their amplification from nine loci. The potential of these markers to assess genetic diversity was evaluated on a set of 16 Irish and 15 European L. perenne ecotypes, nine L. perenne cultivars, other Lolium taxa and other grass species.

Key Results

All analysed Poaceae chloroplast genomes contained more than 200 mononucleotide repeats (chloroplast simple sequence repeats, cpSSRs) of at least 7 bp in length, concentrated mainly in the large single copy region of the genome. Nucleotide composition varied considerably among subfamilies (with Pooideae biased towards poly A repeats). The nine new markers distinguish L. perenne from all non-Lolium taxa. TeaCpSSR28 was able to distinguish between all Lolium species and Lolium multiflorum due to an elongation of an A₈ mononucleotide repeat in L. multiflorum. TeaCpSSR31 detected a considerable degree of microsatellite length variation and single nucleotide polymorphism. TeaCpSSR27 revealed variation within some L. perenne accessions due to a 44-bp indel and was hence readily detected by simple agarose gel electrophoresis. Smaller insertion/deletion events or single nucleotide polymorphisms detected by these new markers could be visualized by polyacrylamide gel electrophoresis or DNA sequencing, respectively.

Conclusions

The new markers are a valuable tool for plant breeding companies, seed testing agencies and the wider scientific community due to their ability to monitor genetic diversity within breeding pools, to trace maternal inheritance and to distinguish closely related species.     

Assembly and annotation of a draft genome sequence for Glycine latifolia,a perennial wild relative of soybean

Glycine latifolia (Benth.) Newell & Hymowitz (2n = 40), one of the 27 wild perennial relatives of soybean, possesses genetic diversity and agronomically favorable traits that are lacking in soybean. Here, we report the 939‐Mb draft genome assembly of G. latifolia (PI 559298) using exclusively linked‐reads sequenced from a single Chromium library. We organized scaffolds into 20 chromosome‐scale pseudomolecules utilizing two genetic maps and the Glycine max (L.) Merr. genome sequence. High copy numbers of putative 91‐bp centromere‐specific tandem repeats were observed in consecutive blocks within predicted pericentromeric regions on several pseudomolecules. No 92‐bp putative centromeric repeats, which are abundant in G. max, were detected in G. latifolia or Glycine tomentella. Annotation of the assembled genome and subsequent filtering yielded a high confidence gene set of 54 475 protein‐coding loci. In comparative analysis with five legume species, genes related to defense responses were significantly overrepresented in Glycine‐specific orthologous gene families. A total of 304 putative nucleotide‐binding site (NBS)‐leucine‐rich‐repeat (LRR) genes were identified in this genome assembly. Different from other legume species, we observed a scarcity of TIR‐NBS‐LRR genes in G. latifolia. The G. latifolia genome was also predicted to contain genes encoding 367 LRR‐receptor‐like kinases, a family of proteins involved in basal defense responses and responses to abiotic stress. The genome sequence and annotation of G. latifolia provides a valuable source of alternative alleles and novel genes to facilitate soybean improvement. This study also highlights the efficacy and cost‐effectiveness of the application of Chromium linked‐reads in diploid plant genome de novo assembly.     

A novel genome‐wide microsatellite resource for species of Eucalyptus with linkage‐to‐physical correspondence on the reference genome sequence

Keystone species in their native ranges, eucalypts, are ecologically and genetically very diverse, growing naturally along extensive latitudinal and altitudinal ranges and variable environments. Besides their ecological importance, eucalypts are also the most widely planted trees for sustainable forestry in the world. We report the development of a novel collection of 535 microsatellites for species of Eucalyptus, 494 designed from ESTs and 41 from genomic libraries. A selected subset of 223 was evaluated for individual identification, parentage testing, and ancestral information content in the two most extensively studied species, Eucalyptus grandis and Eucalyptus globulus. Microsatellites showed high transferability and overlapping allele size range, suggesting they have arisen still in their common ancestor and confirming the extensive genome conservation between these two species. A consensus linkage map with 437 microsatellites, the most comprehensive microsatellite‐only genetic map for Eucalyptus, was built by assembling segregation data from three mapping populations and anchored to the Eucalyptus genome. An overall colinearity between recombination‐based and physical positioning of 84% of the mapped microsatellites was observed, with some ordering discrepancies and sporadic locus duplications, consistent with the recently described whole genome duplication events in Eucalyptus. The linkage map covered 95.2% of the 605.8‐Mbp assembled genome sequence, placing one microsatellite every 1.55 Mbp on average, and an overall estimate of physical to recombination distance of 618 kbp/cM. The genetic parameters estimates together with linkage and physical position data for this large set of microsatellites should assist marker choice for genome‐wide population genetics and comparative mapping in Eucalyptus.     

The effect of handling on the yield of two populations of Lolium perenne selected for differences in mature leaf respiration rate

Two populations of Lolium perenne L. S23 (perennial ryegrass), selected for differences in mature leaf dark respiration, were used in a non-destructive indexing system for individual plants, to determine growth parameters. Population GL66, selected for high respiratory rates and low yield, responded strongly to the indexing treatment, when grown at low plant density. Dry weights of all plant parts decreased strongly, as did dry matter percentages of the leaf blades. At high density this population demonstrated the same trend, but additionally allocation to the shoot increased. In contrast, GL72, selected for low respiratory rates and a high yield, responded only at a high plant density. It is argued that there might be a relation between the dissimilar response of the two populations to mechanical influences and the presence of the genotypes of the low-yielding population in the parent variety. The results also emphasize that non-destructive growth analyses can only be used when their effects on the plants are known.     

Whole‐exome targeted sequencing of the uncharacterized pine genome

The large genome size of many species hinders the development and application of genomic tools to study them. For instance, loblolly pine (Pinus taeda L.), an ecologically and economically important conifer, has a large and yet uncharacterized genome of 21.7 Gbp. To characterize the pine genome, we performed exome capture and sequencing of 14 729 genes derived from an assembly of expressed sequence tags. Efficiency of sequence capture was evaluated and shown to be similar across samples with increasing levels of complexity, including haploid cDNA, haploid genomic DNA and diploid genomic DNA. However, this efficiency was severely reduced for probes that overlapped multiple exons, presumably because intron sequences hindered probe:exon hybridizations. Such regions could not be entirely avoided during probe design, because of the lack of a reference sequence. To improve the throughput and reduce the cost of sequence capture, a method to multiplex the analysis of up to eight samples was developed. Sequence data showed that multiplexed capture was reproducible among 24 haploid samples, and can be applied for high‐throughput analysis of targeted genes in large populations. Captured sequences were de novo assembled, resulting in 11 396 expanded and annotated gene models, significantly improving the knowledge about the pine gene space. Interspecific capture was also evaluated with over 98% of all probes designed from P. taeda that were efficient in sequence capture, were also suitable for analysis of the related species Pinus elliottii Engelm.     

Isolation and characterisation of a sucrose: sucrose 1-fructosyltransferase gene from perennial ryegrass (Lolium perenne)

A sucrose: sucrose 1-fructosyltransferase (1-SST) gene and cDNA (Lp 1-SST) from perennial ryegrass (Lolium perenne) were isolated. The Lp 1-SST gene was fully sequenced and shown to contain three exons and two introns. Nucleotide sequence analysis of the 4824 bp Lp 1-SST genomic sequence revealed 1618 bp of 5' UTR and an open reading frame of 1962 bp encoding a protein of 653 amino acids. Lp 1-SST is 95% identical to the tall fescue 1-SST and contains plant fructosyltransferase functional domains. Lp 1-SST corresponds to a single copy gene in perennial ryegrass, and is expressed in young leaf bases and mature leaf sheaths. The recombinant Lp 1-SST protein from corresponding cDNA expression in Pichia pastoris showed 1-SST activity.     

A fungal endophyte of a palatable grass affects preference of large herbivores

Temperate grasses frequently acquire resistance to herbivores through a symbiosis with epichloid fungi that produces alkaloids of variable deterrent effects. However, in those cases without apparent endophyte negative effects on domestic herbivores, it is not clear whether plant consumption or preference is affected or not. We performed three experiments with 1‐year‐old steers (Bos taurus, Aberdeen Angus) and the annual grass Lolium multiflorum, infected or not by Epichloë occultans to evaluate preference and to identify the underlying tolerance mechanisms. The first experiment evaluated steer preference for L. multiflorum cultivated in plots with three endophyte infection frequencies (low, medium and high), and investigated the role of canopy structure and plant nutritional traits on preference. The second experiment evaluated preference for chopped grass, offered in individual trays with contrasting infection frequencies (low and high), to discard possible effects associated with canopy structure and to focus on nutritional traits. The third experiment was performed with a tray + basket design that separated visual and olfactory stimuli from nutritional traits. High endophyte infection frequencies reduced consistently animal preference, even after short (~10 min) feeding events. However, we did not find significant evidence of plant structural, nutritional, visual or olfactory traits. Our results discarded several potential mechanisms; therefore, the dissuasive effect of fungal endophytes on animal consumption might be related to other mechanisms, including, likely, alkaloids and changes on grass metabolome.     

Clade classification of monolignol biosynthesis gene family members reveals target genes to decrease lignin in Lolium perenne

In monocots, lignin content has a strong impact on the digestibility of the cell wall fraction. Engineering lignin biosynthesis requires a profound knowledge of the role of paralogues in the multigene families that constitute the monolignol biosynthesis pathway. We applied a bioinformatics approach for genome‐wide identification of candidate genes in Lolium perenne that are likely to be involved in the biosynthesis of monolignols. More specifically, we performed functional subtyping of phylogenetic clades in four multigene families: 4CL, COMT, CAD and CCR. Essential residues were considered for functional clade delineation within these families. This classification was complemented with previously published experimental evidence on gene expression, gene function and enzymatic activity in closely related crops and model species. This allowed us to assign functions to novel identified L. perenne genes, and to assess functional redundancy among paralogues. We found that two 4CL paralogues, two COMT paralogues, three CCR paralogues and one CAD gene are prime targets for genetic studies to engineer developmentally regulated lignin in this species. Based on the delineation of sequence conservation between paralogues and a first analysis of allelic diversity, we discuss possibilities to further study the roles of these paralogues in lignin biosynthesis, including expression analysis, reverse genetics and forward genetics, such as association mapping. We propose criteria to prioritise paralogues within multigene families and certain SNPs within these genes for developing genotyping assays or increasing power in association mapping studies. Although L. perenne was the target of the analyses presented here, this functional subtyping of phylogenetic clades represents a valuable tool for studies investigating monolignol biosynthesis genes in other monocot species.     

Overexpression of ubiquitin‐like LpHUB1 gene confers drought tolerance in perennial ryegrass

HUB1, also known as Ubl5, is a member of the subfamily of ubiquitin‐like post‐translational modifiers. HUB1 exerts its role by conjugating with protein targets. The function of this protein has not been studied in plants. A HUB1 gene, LpHUB1, was identified from serial analysis of gene expression data and cloned from perennial ryegrass. The expression of this gene was reported previously to be elevated in pastures during the summer and by drought stress in climate‐controlled growth chambers. Here, pasture‐type and turf‐type transgenic perennial ryegrass plants overexpressing LpHUB1 showed improved drought tolerance, as evidenced by improved turf quality, maintenance of turgor and increased growth. Additional analyses revealed that the transgenic plants generally displayed higher relative water content, leaf water potential, and chlorophyll content and increased photosynthetic rate when subjected to drought stress. These results suggest HUB1 may play an important role in the tolerance of perennial ryegrass to abiotic stresses.     

Genotyping by RAD sequencing enables mapping of fatty acid composition traits in perennial ryegrass (Lolium perenne (L.))

Perennial ryegrass (Lolium perenne L.) is the most important forage crop in temperate livestock agriculture. Its nutritional quality has significant impact on the quality of meat and milk for human consumption. Evidence suggests that higher energy content in forage can assist in reducing greenhouse gas emissions from ruminants. Increasing the fatty acid content (especially α‐linolenic acid, an omega‐3 fatty acid) may thus contribute to better forage, but little is known about the genetic basis of variation for this trait. To this end, quantitative trait loci (QTLs) were identified associated with major fatty acid content in perennial ryegrass using a population derived from a cross between the heterozygous and outbreeding high‐sugar grass variety AberMagic and an older variety, Aurora. A genetic map with 434 restriction‐associated DNA (RAD) and SSR markers was generated. Significant QTLs for the content of palmitic (C16:0) on linkage groups (LGs) 2 and 7; stearic (C18:0) on LGs 3, 4 and 7; linoleic (C18:2n‐6) on LGs 2 and 5; and α‐linolenic acids (C18:3n‐3) on LG 1 were identified. Two candidate genes (a lipase and a beta‐ketoacyl CoA synthase), both associated with C16:0, and separately with C18:2n‐6 and C18:0 contents, were identified. The physical positions of these genes in rice and their genetic positions in perennial ryegrass were consistent with established syntenic relationships between these two species. Validation of these associations is required, but the utility of RAD markers for rapid generation of genetic maps and QTL analysis has been demonstrated for fatty acid composition in a global forage crop.     

Genome sequence of M6, a diploid inbred clone of the high‐glycoalkaloid‐producing tuber‐bearing potato species Solanum chacoense,reveals residual heterozygosity

Cultivated potato (Solanum tuberosum L.) is a highly heterozygous autotetraploid that presents challenges in genome analyses and breeding. Wild potato species serve as a resource for the introgression of important agronomic traits into cultivated potato. One key species is Solanum chacoense and the diploid, inbred clone M6, which is self‐compatible and has desirable tuber market quality and disease resistance traits. Sequencing and assembly of the genome of the M6 clone of S. chacoense generated an assembly of 825 767 562 bp in 8260 scaffolds with an N50 scaffold size of 713 602 bp. Pseudomolecule construction anchored 508 Mb of the genome assembly into 12 chromosomes. Genome annotation yielded 49 124 high‐confidence gene models representing 37 740 genes. Comparative analyses of the M6 genome with six other Solanaceae species revealed a core set of 158 367 Solanaceae genes and 1897 genes unique to three potato species. Analysis of single nucleotide polymorphisms across the M6 genome revealed enhanced residual heterozygosity on chromosomes 4, 8 and 9 relative to the other chromosomes. Access to the M6 genome provides a resource for identification of key genes for important agronomic traits and aids in genome‐enabled development of inbred diploid potatoes with the potential to accelerate potato breeding.     

The differential contributions of herkogamy and dichogamy as mechanisms of avoiding self‐interference in four self‐incompatible Epimedium species

Self‐interference is one of the most important selective forces in shaping floral evolution. Herkogamy and dichogamy both can achieve reductions in the extent of self‐interference, but they may have different roles in minimizing self‐interference in a single species. We used four self‐incompatible Epimedium species to explore the roles of herkogamy and dichogamy in avoiding self‐interference and to test the hypothesis that herkogamy and dichogamy may be separated and become selected preferentially in the taxa. Two species (E. franchetii and E. mikinorii) expressed strong herkogamy and weak protogyny (adichogamy), whereas another two species (E. sutchuenense and E. leptorrhizum) expressed slight herkogamy and partial protandry. Field investigations indicated that there was no physical self‐interference between male function and female function regarding pollen removal and pollen deposition in all species. Self‐pollination (autonomous or facilitated) was greater in species with slight herkogamy than in those with strong herkogamy. Artificial pollination treatments revealed that self‐pollination could reduce outcrossed female fertility in all species, and we found evidence that self‐interference reduced seed set in E. sutchuenense and E. leptorrhizum in the field, but not in E. franchetii and E. mikinorii. These results indicate that well‐developed herkogamy is more effective compared with dichogamy in avoiding self‐interference in the four species. In genus Epimedium, herkogamy instead of dichogamy should be selected preferentially and evolved as an effective mechanism for avoiding self‐interference and might not need to evolve linked with dichogamy.     

Haplotype‐based genotyping‐by‐sequencing in oat genome research

In a de novo genotyping‐by‐sequencing (GBS) analysis of short, 64‐base tag‐level haplotypes in 4657 accessions of cultivated oat, we discovered 164741 tag‐level (TL) genetic variants containing 241224 SNPs. From this, the marker density of an oat consensus map was increased by the addition of more than 70000 loci. The mapped TL genotypes of a 635‐line diversity panel were used to infer chromosome‐level (CL) haplotype maps. These maps revealed differences in the number and size of haplotype blocks, as well as differences in haplotype diversity between chromosomes and subsets of the diversity panel. We then explored potential benefits of SNP vs. TL vs. CL GBS variants for mapping, high‐resolution genome analysis and genomic selection in oats. A combined genome‐wide association study (GWAS) of heading date from multiple locations using both TL haplotypes and individual SNP markers identified 184 significant associations. A comparative GWAS using TL haplotypes, CL haplotype blocks and their combinations demonstrated the superiority of using TL haplotype markers. Using a principal component‐based genome‐wide scan, genomic regions containing signatures of selection were identified. These regions may contain genes that are responsible for the local adaptation of oats to Northern American conditions. Genomic selection for heading date using TL haplotypes or SNP markers gave comparable and promising prediction accuracies of up to r = 0.74. Genomic selection carried out in an independent calibration and test population for heading date gave promising prediction accuracies that ranged between r = 0.42 and 0.67. In conclusion, TL haplotype GBS‐derived markers facilitate genome analysis and genomic selection in oat.