Generation of an integrated Hieracium genomic and transcriptomic resource enables exploration of small RNA pathways during apomixis initiation.

BACKGROUND: Application of apomixis, or asexual seed formation, in crop breeding would allow rapid fixation of complex traits, economizing improved crop delivery. Identification of apomixis genes is confounded by the polyploid nature, high genome complexity and lack of genomic sequence integration with reproductive tissue transcriptomes in most apomicts. RESULTS: A genomic and transcriptomic resource was developed for Hieracium subgenus Pilosella (Asteraceae) which incorporates characterized sexual, apomictic and mutant apomict plants exhibiting reversion to sexual reproduction. Apomicts develop additional female gametogenic cells that suppress the sexual pathway in ovules. Disrupting small RNA pathways in sexual Arabidopsis also induces extra female gametogenic cells; therefore, the resource was used to examine if changes in small RNA pathways correlate with apomixis initiation. An initial characterization of small RNA pathway genes within Hieracium was undertaken, and ovary-expressed ARGONAUTE genes were identified and cloned. Comparisons of whole ovary transcriptomes from mutant apomicts, relative to the parental apomict, revealed that differentially expressed genes were enriched for processes involved in small RNA biogenesis and chromatin silencing. Small RNA profiles within mutant ovaries did not reveal large-scale alterations in composition or length distributions; however, a small number of differentially expressed, putative small RNA targets were identified. CONCLUSIONS: The established Hieracium resource represents a substantial contribution towards the investigation of early sexual and apomictic female gamete development, and the generation of new candidate genes and markers. Observed changes in small RNA targets and biogenesis pathways within sexual and apomictic ovaries will underlie future functional research into apomixis initiation in Hieracium.

A reference genetic linkage map of apomictic Hieracium species based on expressed markers derived from developing ovule transcripts.

BACKGROUND AND AIMS: Apomixis in plants generates clonal progeny with a maternal genotype through asexual seed formation. Hieracium subgenus Pilosella (Asteraceae) contains polyploid, highly heterozygous apomictic and sexual species. Within apomictic Hieracium, dominant genetic loci independently regulate the qualitative developmental components of apomixis. In H. praealtum, LOSS OF APOMEIOSIS (LOA) enables formation of embryo sacs without meiosis and LOSS OF PARTHENOGENESIS (LOP) enables fertilization-independent seed formation. A locus required for fertilization-independent endosperm formation (AutE) has been identified in H. piloselloides. Additional quantitative loci appear to influence the penetrance of the qualitative loci, although the controlling genes remain unknown. This study aimed to develop the first genetic linkage maps for sexual and apomictic Hieracium species using simple sequence repeat (SSR) markers derived from expressed transcripts within the developing ovaries. METHODS: RNA from microdissected Hieracium ovule cell types and ovaries was sequenced and SSRs were identified. Two different F1 mapping populations were created to overcome difficulties associated with genome complexity and asexual reproduction. SSR markers were analysed within each mapping population to generate draft linkage maps for apomictic and sexual Hieracium species. KEY RESULTS: A collection of 14 684 Hieracium expressed SSR markers were developed and linkage maps were constructed for Hieracium species using a subset of the SSR markers. Both the LOA and LOP loci were successfully assigned to linkage groups; however, AutE could not be mapped using the current populations. Comparisons with lettuce (Lactuca sativa) revealed partial macrosynteny between the two Asteraceae species. CONCLUSIONS: A collection of SSR markers and draft linkage maps were developed for two apomictic and one sexual Hieracium species. These maps will support cloning of controlling genes at LOA and LOP loci in Hieracium and should also assist with identification of quantitative loci that affect the expressivity of apomixis. Future work will focus on mapping AutE using alternative populations.

Genome-wide SNP identification in multiple morphotypes of allohexaploid tall fescue (Festuca arundinacea Schreb).

BACKGROUND: Single nucleotide polymorphisms (SNPs) provide essential tools for the advancement of research in plant genomics, and the development of SNP resources for many species has been accelerated by the capabilities of second-generation sequencing technologies. The current study aimed to develop and use a novel bioinformatic pipeline to generate a comprehensive collection of SNP markers within the agriculturally important pasture grass tall fescue; an outbreeding allopolyploid species displaying three distinct morphotypes: Continental, Mediterranean and rhizomatous. RESULTS: A bioinformatic pipeline was developed that successfully identified SNPs within genotypes from distinct tall fescue morphotypes, following the sequencing of 414 polymerase chain reaction (PCR) - generated amplicons using 454 GS FLX technology. Equivalent amplicon sets were derived from representative genotypes of each morphotype, including six Continental, five Mediterranean and one rhizomatous. A total of 8,584 and 2,292 SNPs were identified with high confidence within the Continental and Mediterranean morphotypes respectively. The success of the bioinformatic approach was demonstrated through validation (at a rate of 70%) of a subset of 141 SNPs using both SNaPshot™ and GoldenGate™ assay chemistries. Furthermore, the quantitative genotyping capability of the GoldenGate™ assay revealed that approximately 30% of the putative SNPs were accessible to co-dominant scoring, despite the hexaploid genome structure. The sub-genome-specific origin of each SNP validated from Continental tall fescue was predicted using a phylogenetic approach based on comparison with orthologous sequences from predicted progenitor species. CONCLUSIONS: Using the appropriate bioinformatic approach, amplicon resequencing based on 454 GS FLX technology is an effective method for the identification of polymorphic SNPs within the genomes of Continental and Mediterranean tall fescue. The GoldenGate™ assay is capable of high-throughput co-dominant SNP allele detection, and minimises the problems associated with SNP genotyping in a polyploid by effectively reducing the complexity to a diploid system. This SNP collection may now be refined and used in applications such as cultivar identification, genetic linkage map construction, genome-wide association studies and genomic selection in tall fescue. The bioinformatic pipeline described here represents an effective general method for SNP discovery within outbreeding allopolyploid species.

Molecular characterisation and interpretation of genetic diversity within globally distributed germplasm collections of tall fescue (Festuca arundinacea Schreb.) and meadow fescue (F. pratensis Huds.).

Allohexaploid tall fescue (Festuca arundinacea Schreb. syn. Lolium arundinaceum [Schreb.] Darbysh.) is an agriculturally important grass cultivated for pasture and turf world-wide. Genetic improvement of tall fescue could benefit from the use of non-domesticated germplasm to diversify breeding populations through the incorporation of novel and superior allele content. However, such potential germplasm must first be characterised, as three major morphotypes (Continental, Mediterranean and rhizomatous) with varying degrees of hybrid interfertility are commonly described within this species. As hexaploid tall fescue is also a member of a polyploid species complex that contains tetraploid, octoploid and decaploid taxa, it is also possible that germplasm collections may have inadvertently sampled some of these sub-species. In this study, 1,040 accessions from the publicly available United States Department of Agriculture tall fescue and meadow fescue germplasm collections were investigated. Sequence of the chloroplast genome-located matK gene and the nuclear ribosomal DNA internal transcribed spacer (rDNA ITS) permitted attribution of accessions to the three previously known morphotypes and also revealed the presence of tall fescue sub-species of varying ploidy levels, as well as other closely related species. The majority of accessions were, however, identified as Continental hexaploid tall fescue. Analysis using 34 simple sequence repeat markers was able to further investigate the level of genetic diversity within each hexaploid tall fescue morphotype group. At least two genetically distinct sub-groups of Continental hexaploid tall fescue were identified which are probably associated with palaeogeographic range expansion of this morphotype. This work has comprehensively characterised a large and complex germplasm collection and has identified genetically diverse accessions which may potentially contribute valuable alleles at agronomic loci for tall fescue cultivar improvement programs.

Identification of QTLs for morphological traits influencing waterlogging tolerance in perennial ryegrass (Lolium perenne L.).

Perennial ryegrass is a globally cultivated obligate outbreeding diploid species (2n = 2x = 14) which is subjected to periods of waterlogging stress due to flood irrigation during winter and the lead-up to summer. Reduction of oxygen supply to root systems due to waterlogging produces consequent deleterious effects on plant performance. Framework genetic maps for a large-scale genetic mapping family [F1(NA(x) × AU6)] were constructed containing 91 simple sequence repeat and 24 single nucleotide polymorphism genetic markers. Genetic trait dissection using both control and waterlogging treatments was performed in the glasshouse, a total of 143 maximally recombinant genotypes being selected from the overall sib-ship and replicated threefold in the trial. Analysis was performed for nine quantitative morphological traits measured 8 weeks after stress treatments were applied. A total of 37 quantitative trait loci (QTLs) were identified; 19 on the NA(x) parental genetic map, and 18 on the AU6 parental genetic map. Regions of particular interest were identified on linkage groups (LGs) 4 and 3 of the respective maps, which have been targeted for further analysis by selection of critical recombinants. This first study of genetic control of waterlogging tolerance in ryegrasses has important implications for breeding improvement of abiotic stress adaptation.

Evolutionary history of tall fescue morphotypes inferred from molecular phylogenetics of the Lolium-Festuca species complex.

BACKGROUND: The agriculturally important pasture grass tall fescue (Festuca arundinacea Schreb. syn. Lolium arundinaceum (Schreb.) Darbysh.) is an outbreeding allohexaploid, that may be more accurately described as a species complex consisting of three major (Continental, Mediterranean and rhizomatous) morphotypes. Observation of hybrid infertility in some crossing combinations between morphotypes suggests the possibility of independent origins from different diploid progenitors. This study aims to clarify the evolutionary relationships between each tall fescue morphotype through phylogenetic analysis using two low-copy nuclear genes (encoding plastid acetyl-CoA carboxylase [Acc1] and centroradialis [CEN]), the nuclear ribosomal DNA internal transcribed spacer (rDNA ITS) and the chloroplast DNA (cpDNA) genome-located matK gene. Other taxa within the closely related Lolium-Festuca species complex were also included in the study, to increase understanding of evolutionary processes in a taxonomic group characterised by multiple inter-specific hybridisation events. RESULTS: Putative homoeologous sequences from both nuclear genes were obtained from each polyploid species and compared to counterparts from 15 diploid taxa. Phylogenetic reconstruction confirmed F. pratensis and F. arundinacea var. glaucescens as probable progenitors to Continental tall fescue, and these species are also likely to be ancestral to the rhizomatous morphotype. However, these two morphotypes are sufficiently distinct to be located in separate clades based on the ITS-derived data set. All four of the generated data sets suggest independent evolution of the Mediterranean and Continental morphotypes, with minimal affinity between cognate sequence haplotypes. No obvious candidate progenitor species for Mediterranean tall fescues were identified, and only two putative sub-genome-specific haplotypes were identified for this morphotype. CONCLUSIONS: This study describes the first phylogenetic analysis of the Festuca genus to include representatives of each tall fescue morphotype, and to use low copy nuclear gene-derived sequences to identify putative progenitors of the polyploid species. The demonstration of distinct tall fescue lineages has implications for both taxonomy and molecular breeding strategies, and may facilitate the generation of morphotype and/or sub-genome-specific molecular markers.

Comparison of homoeolocus organisation in paired BAC clones from white clover (Trifolium repens L.) and microcolinearity with model legume species.

BACKGROUND: White clover (Trifolium repens L.) is an outbreeding allotetraploid species and an important forage legume in temperate grassland agriculture. Comparison of sub-genome architecture and study of nucleotide sequence diversity within allopolyploids provides insight into evolutionary divergence mechanisms, and is also necessary for the development of whole-genome sequencing strategies. This study aimed to evaluate the degree of divergence between the O and P' sub-genomes of white clover through sequencing of BAC clones containing paired homoeoloci. The microsyntenic relationships between the genomes of white clover and the model legumes Lotus japonicus and Medicago truncatula as well as Arabidopsis thaliana were also characterised. RESULTS: A total of four paired homoeologous BACs were selected and sequenced to generate 173 kb of overlapping sequence between the O and P' sub-genomes. Equivalent gene content was generally observed, apart from small-scale deletions, in contrast to conservation of intergenic sequences, which varied between the four selected regions. Measurement of the number of synonymous substitutions between homoeologous genes led to estimation of a 4.2 million year divergence time between the two sub-genomes. Microsynteny was observed between the genomes of white clover and L. japonicus for all four targeted regions, but corresponding M. truncatula genomic regions were only identified for two BAC pairs. CONCLUSIONS: This study describes the first analysis of sub-genome structural conservation across selected genomic regions in white clover. Although the high levels of sequence conservation between the O and P' sub-genomes would complicate efforts for whole genome sequence assembly, the conserved microsynteny with model legume genomes, especially that of L. japonicus, will be highly valuable for the future of white clover genomics and molecular breeding.

Identification of genetic factors influencing salt stress tolerance in white clover (Trifolium repens L.) by QTL analysis.

Allotetraploid (2n = 4x = 32) white clover (Trifolium repens L.) is the most commonly cultivated legume component of temperate pastures, sown in swards with a companion grass species. Genetic control of growth performance of white clover on saline land is highly important for dairy industries, due to increasing soil salinity problems. The objective of this study was to identify quantitative trait loci (QTLs) for salinity tolerance in terms of vegetative growth under stress. Two parental genetic maps consisting of 213 and 159 marker loci and spanning 1,973.0 and 1,837.6 cM, respectively, were constructed using simple sequence repeat (SSR) and single nucleotide polymorphism (SNP) markers from a two-way pseudo-test cross F(1) population derived from pair-crossing of the Haifa(2) and LCL(2) genotypes. A total of 8 unique genomic regions on 8 linkage groups (LGs) of the Haifa(2) parental map and 6 unique regions on 5 LGs in the LCL(2) parental map were associated with plant growth under salt stress and relative growth under stress, as compared to control conditions. The results of this study indicate that salt tolerance in white clover is controlled by multiple QTLs, some at common locations, but each of limited magnitude. Location of these QTLs provides the genetic basis and potential for pyramiding of salt tolerance genes in breeding improvement.

A detached leaf assay for testing transient gene expression and gene editing in cowpea (Vigna unguiculata [L.] Walp.).

BACKGROUND: The legume cowpea (Vigna unguiculata L.) is extensively grown in sub-Saharan Africa. Cowpea, like many legumes has proved recalcitrant to plant transformation. A rapid transient leaf assay was developed for testing gene expression and editing constructs prior to stable cowpea transformation, to accelerate cowpea and legume crop improvement. RESULTS: Attempts to develop a transient protoplast system for cowpea were unsuccessful. Leaflets from plants 3-4 weeks post-germination were age selected to establish a rapid Agrobacterium (Agro) infiltration-mediated transient system for efficacy testing of gene expression and CRISPR/Cas9 gene editing constructs. In planta, Agro-infiltration of leaflets with fluorescent expression constructs, resulted in necrosis. By contrast, Agro-infiltration of detached leaflets with an Arabidopsis (At) ubiquitin3 promoter:ZsGreen construct, followed by culture on solid nutrient medium resulted in fluorescence in over 48% of leaf cells. Expression efficiency was leaf age-dependent. Three cowpea meiosis genes were identified for CRISPR/Cas9 gene-editing, with the forward aim of meiosis-knock out for asexual seed induction in cowpea. Constructs were designed and tested containing candidate gene-specific guide RNAs, expressed using either the cowpea or Arabidopsis U6 promoters with Cas9 expression directed by either the Arabidopsis 40S ribosomal protein or parsley ubiquitin4-2 promoters. Leaflets were infiltrated with test gene-editing constructs and analytical methods developed to identify gene-specific mutations. A construct that produced mutations predicted to induce functional knockout of in the VuSPO11-1 meiosis gene was tested for efficacy in primary transgenic cowpea plants using a previously established stable transformation protocol. Vuspo11-1 mutants were identified, that cytologically phenocopied spo11-1 mutants previously characterized in Arabidopsis, and rice. Importantly, a biallelic male and female sterile mutant was identified in primary transgenics, exhibiting the expected defects in 100% of examined male and female meiocytes. CONCLUSION: The transient, detached cowpea leaf assay, and supporting analytical methods developed, provide a rapid and reproducible means for testing gene expression constructs, and constructs for inducing mutagenesis in genes involved in both vegetative and reproductive developmental programs. The method and tested editing constructs and components have potential application for a range of crop legumes.

Assembled genomic and tissue-specific transcriptomic data resources for two genetically distinct lines of Cowpea ( Vigna unguiculata (L.) Walp).

Cowpea ( Vigna unguiculata (L.) Walp) is an important legume crop for food security in areas of low-input and smallholder farming throughout Africa and Asia. Genetic improvements are required to increase yield and resilience to biotic and abiotic stress and to enhance cowpea crop performance. An integrated cowpea genomic and gene expression data resource has the potential to greatly accelerate breeding and the delivery of novel genetic traits for cowpea. Extensive genomic resources for cowpea have been absent from the public domain; however, a recent early release reference genome for IT97K-499-35 ( Vigna unguiculata v1.0, NSF, UCR, USAID, DOE-JGI, http://phytozome.jgi.doe.gov/) has now been established in a collaboration between the Joint Genome Institute (JGI) and University California (UC) Riverside. Here we release supporting genomic and transcriptomic data for two transformable cowpea varieties, IT97K-499-35 and IT86D-1010. The transcriptome resource includes six tissue-specific datasets for each variety, with particular emphasis on reproductive tissues that extend and support the V. unguiculata v1.0 reference. Annotations have been included in our resource to allow direct mapping to the v1.0 cowpea reference. The resource described here is supported by downloadable raw and assembled sequence data.

A Comparison of In Vitro and In Vivo Asexual Embryogenesis.

In plants, embryogenesis generally occurs through the sexual process of double fertilization, which involves a haploid sperm cell fusing with a haploid egg cell to ultimately give rise to a diploid embryo. Embryogenesis can also occur asexually in the absence of fertilization, both in vitro and in vivo. Somatic or gametic cells are able to differentiate into embryos in vitro following the application of plant growth regulators or stress treatments. Asexual embryogenesis also occurs naturally in some plant species in vivo, from either ovule cells as part of a process defined as apomixis, or from somatic leaf tissue in other species. In both in vitro and in vivo asexual embryogenesis, the embryo precursor cells must attain an embryogenic fate without the act of fertilization. This review compares the processes of in vitro and in vivo asexual embryogenesis including what is known regarding the genetic and epigenetic regulation of each process, and considers how the precursor cells are able to change fate and adopt an embryogenic pathway.

The genetic control of apomixis: asexual seed formation.

Apomixis (asexual seed formation) is the result of a plant gaining the ability to bypass the most fundamental aspects of sexual reproduction: meiosis and fertilization. Without the need for male fertilization, the resulting seed germinates a plant that develops as a maternal clone. This dramatic shift in reproductive process has been documented in many flowering plant species, although no major seed crops have been shown to be capable of apomixis. The ability to generate maternal clones and therefore rapidly fix desirable genotypes in crop species could accelerate agricultural breeding strategies. The potential of apomixis as a next-generation breeding technology has contributed to increasing interest in the mechanisms controlling apomixis. In this review, we discuss the progress made toward understanding the genetic and molecular control of apomixis. Research is currently focused on two fronts. One aims to identify and characterize genes causing apomixis in apomictic species that have been developed as model species. The other aims to engineer or switch the sexual seed formation pathway in non-apomictic species, to one that mimics apomixis. Here we describe the major apomictic mechanisms and update knowledge concerning the loci that control them, in addition to presenting candidate genes that may be used as tools for switching the sexual pathway to an apomictic mode of reproduction in crops.

Nucleotide diversity of vernalization and flowering-time-related genes in a germplasm collection of meadow fescue (Festuca pratensis Huds. syn. Lolium pratense (Huds.) Darbysh.).

In plant species, control of flowering time is an important factor for adaptation to local natural environments. The Vrn1 , CO , FT1 and CK2α genes are key components in the flowering-specific signaling pathway of grass species. Meadow fescue is an agronomically important forage grass species, which is naturally distributed across Europe and Western Asia. In this study, meadow fescue flowering-time-related genes were resequenced to assess nucleotide diversity in European and Western Asian subpopulations. Identified sequence polymorphisms were then converted into PCR-based molecular genetic markers, and a meadow fescue germplasm collection was genotyped to investigate global allelic variation. Lower nucleotide diversities were observed for the Vrn1 and CO orthologs, while relatively higher values were observed for the FT1 and casein kinase II α-subunit (CK2α) orthologs. The nucleotide diversity for FT1 orthologs in the Western Asian subpopulation was significantly higher than those of the European subpopulation. Similarly, significant differences in nucleotide diversity for the remaining genes were observed between several combinations of subpopulation. The global allele distribution pattern was consistent with observed level of nucleotide diversity. These results suggested that the degree of purifying selection acting on the genes differs according to geographical location. As previously shown for model plant species, functional specificities of flowering-time-related genes may also vary according to environmental conditions.

Plastome Sequence Determination and Comparative Analysis for Members of the Lolium-Festuca Grass Species Complex.

Chloroplast genome sequences are of broad significance in plant biology, due to frequent use in molecular phylogenetics, comparative genomics, population genetics, and genetic modification studies. The present study used a second-generation sequencing approach to determine and assemble the plastid genomes (plastomes) of four representatives from the agriculturally important Lolium-Festuca species complex of pasture grasses (Lolium multiflorum, Festuca pratensis, Festuca altissima, and Festuca ovina). Total cellular DNA was extracted from either roots or leaves, was sequenced, and the output was filtered for plastome-related reads. A comparison between sources revealed fewer plastome-related reads from root-derived template but an increase in incidental bacterium-derived sequences. Plastome assembly and annotation indicated high levels of sequence identity and a conserved organization and gene content between species. However, frequent deletions within the F. ovina plastome appeared to contribute to a smaller plastid genome size. Comparative analysis with complete plastome sequences from other members of the Poaceae confirmed conservation of most grass-specific features. Detailed analysis of the rbcL-psaI intergenic region, however, revealed a "hot-spot" of variation characterized by independent deletion events. The evolutionary implications of this observation are discussed. The complete plastome sequences are anticipated to provide the basis for potential organelle-specific genetic modification of pasture grasses.

Identification of homologous, homoeologous and paralogous sequence variants in an outbreeding allopolyploid species based on comparison with progenitor taxa.

The combination of homologous, homoeologous and paralogous classes of sequence variation presents major challenges for SNP discovery in outbreeding allopolyploid species. Previous in vitro gene-associated SNP discovery studies in the allotetraploid forage legume white clover (Trifolium repens L.) were vulnerable to such effects, leading to prohibitive levels of attrition during SNP validation. Identification of T. occidentale and T. pallescens as the putative diploid progenitors of white clover has permitted discrimination of the different sequence variant categories. Amplicons from selected abiotic stress tolerance-related genes were obtained using mapping family parents and individuals from each diploid species. Following cloning, progenitor comparison allowed tentative assignment of individual haplotypes to one or other sub-genome, as well as to gene copies within sub-genomes. A high degree of coincidence and identity between SNPs and HSVs was observed. Close similarity was observed between the genome of T. occidentale and one white clover sub-genome, but the affinity between T. pallescens and the other sub-genome was weaker, suggesting that a currently uncharacterised taxon may be the true second progenitor. Selected validated SNPs were attributed to individual sub-genomes by assignment to and naming of homoeologous linkage groups, providing the basis for improved genetic trait-dissection studies. The approach described in this study is broadly applicable to a range of allopolyploid taxa of equivocal ancestry.