Barley: From Brittle to Stable Harvest.

Selection and domestication of plants with genes that prevent grains from shattering in cereals was essential for human civilization's transition to agriculture-based societies. In this issue, Pourkheirandish et al. show that domestication of barley required evolution of a molecular system distinct from other grains, such as rice and maize, and reveal that present-day cultivars derive from two ancient domestication centers.

The mosaic oat genome gives insights into a uniquely healthy cereal crop.

Cultivated oat (Avena sativa L.) is an allohexaploid (AACCDD, 2n = 6x = 42) thought to have been domesticated more than 3,000 years ago while growing as a weed in wheat, emmer and barley fields in Anatolia1,2. Oat has a low carbon footprint, substantial health benefits and the potential to replace animal-based food products. However, the lack of a fully annotated reference genome has hampered efforts to deconvolute its complex evolutionary history and functional gene dynamics. Here we present a high-quality reference genome of A. sativa and close relatives of its diploid (Avena longiglumis, AA, 2n = 14) and tetraploid (Avena insularis, CCDD, 2n = 4x = 28) progenitors. We reveal the mosaic structure of the oat genome, trace large-scale genomic reorganizations in the polyploidization history of oat and illustrate a breeding barrier associated with the genome architecture of oat. We showcase detailed analyses of gene families implicated in human health and nutrition, which adds to the evidence supporting oat safety in gluten-free diets, and we perform mapping-by-sequencing of an agronomic trait related to water-use efficiency. This resource for the Avena genus will help to leverage knowledge from other cereal genomes, improve understanding of basic oat biology and accelerate genomics-assisted breeding and reanalysis of quantitative trait studies.

The barley pan-genome reveals the hidden legacy of mutation breeding.

Genetic diversity is key to crop improvement. Owing to pervasive genomic structural variation, a single reference genome assembly cannot capture the full complement of sequence diversity of a crop species (known as the 'pan-genome'1). Multiple high-quality sequence assemblies are an indispensable component of a pan-genome infrastructure. Barley (Hordeum vulgare L.) is an important cereal crop with a long history of cultivation that is adapted to a wide range of agro-climatic conditions2. Here we report the construction of chromosome-scale sequence assemblies for the genotypes of 20 varieties of barley-comprising landraces, cultivars and a wild barley-that were selected as representatives of global barley diversity. We catalogued genomic presence/absence variants and explored the use of structural variants for quantitative genetic analysis through whole-genome shotgun sequencing of 300 gene bank accessions. We discovered abundant large inversion polymorphisms and analysed in detail two inversions that are frequently found in current elite barley germplasm; one is probably the product of mutation breeding and the other is tightly linked to a locus that is involved in the expansion of geographical range. This first-generation barley pan-genome makes previously hidden genetic variation accessible to genetic studies and breeding.

Multiple wheat genomes reveal global variation in modern breeding.

Advances in genomics have expedited the improvement of several agriculturally important crops but similar efforts in wheat (Triticum spp.) have been more challenging. This is largely owing to the size and complexity of the wheat genome1, and the lack of genome-assembly data for multiple wheat lines2,3. Here we generated ten chromosome pseudomolecule and five scaffold assemblies of hexaploid wheat to explore the genomic diversity among wheat lines from global breeding programs. Comparative analysis revealed extensive structural rearrangements, introgressions from wild relatives and differences in gene content resulting from complex breeding histories aimed at improving adaptation to diverse environments, grain yield and quality, and resistance to stresses4,5. We provide examples outlining the utility of these genomes, including a detailed multi-genome-derived nucleotide-binding leucine-rich repeat protein repertoire involved in disease resistance and the characterization of Sm16, a gene associated with insect resistance. These genome assemblies will provide a basis for functional gene discovery and breeding to deliver the next generation of modern wheat cultivars.

Natural alleles of the abscisic acid catabolism gene ZmAbh4 modulate water use efficiency and carbon isotope discrimination in maize.

Altering plant water use efficiency (WUE) is a promising approach for achieving sustainable crop production in changing climate scenarios. Here, we show that WUE can be tuned by alleles of a single gene discovered in elite maize (Zea mays) breeding material. Genetic dissection of a genomic region affecting WUE led to the identification of the gene ZmAbh4 as causative for the effect. CRISPR/Cas9-mediated ZmAbh4 inactivation increased WUE without growth reductions in well-watered conditions. ZmAbh4 encodes an enzyme that hydroxylates the phytohormone abscisic acid (ABA) and initiates its catabolism. Stomatal conductance is regulated by ABA and emerged as a major link between variation in WUE and discrimination against the heavy carbon isotope (Δ13C) during photosynthesis in the C4 crop maize. Changes in Δ13C persisted in kernel material, which offers an easy-to-screen proxy for WUE. Our results establish a direct physiological and genetic link between WUE and Δ13C through a single gene with potential applications in maize breeding.

A chromosome conformation capture ordered sequence of the barley genome.

Cereal grasses of the Triticeae tribe have been the major food source in temperate regions since the dawn of agriculture. Their large genomes are characterized by a high content of repetitive elements and large pericentromeric regions that are virtually devoid of meiotic recombination. Here we present a high-quality reference genome assembly for barley (Hordeum vulgare L.). We use chromosome conformation capture mapping to derive the linear order of sequences across the pericentromeric space and to investigate the spatial organization of chromatin in the nucleus at megabase resolution. The composition of genes and repetitive elements differs between distal and proximal regions. Gene family analyses reveal lineage-specific duplications of genes involved in the transport of nutrients to developing seeds and the mobilization of carbohydrates in grains. We demonstrate the importance of the barley reference sequence for breeding by inspecting the genomic partitioning of sequence variation in modern elite germplasm, highlighting regions vulnerable to genetic erosion.

De Novo Genome Assembly of the Japanese Wheat Cultivar Norin 61 Highlights Functional Variation in Flowering Time and Fusarium-Resistant Genes in East Asian Genotypes.

Bread wheat is a major crop that has long been the focus of basic and breeding research. Assembly of its genome has been difficult because of its large size and allohexaploid nature (AABBDD genome). Following the first reported assembly of the genome of the experimental strain Chinese Spring (CS), the 10+ Wheat Genomes Project was launched to produce multiple assemblies of worldwide modern cultivars. The only Asian cultivar in the project is Norin 61, a representative Japanese cultivar adapted to grow across a broad latitudinal range, mostly characterized by a wet climate and a short growing season. Here, we characterize the key aspects of its chromosome-scale genome assembly spanning 15 Gb with a raw scaffold N50 of 22 Mb. Analysis of the repetitive elements identified chromosomal regions unique to Norin 61 that encompass a tandem array of the pathogenesis-related 13 family. We report novel copy-number variations in the B homeolog of the florigen gene FT1/VRN3, pseudogenization of its D homeolog and the association of its A homeologous alleles with the spring/winter growth habit. Furthermore, the Norin 61 genome carries typical East Asian functional variants different from CS, ranging from a single nucleotide to multi-Mb scale. Examples of such variation are the Fhb1 locus, which confers Fusarium head-blight resistance, Ppd-D1a, which confers early flowering, Glu-D1f for Asian noodle quality and Rht-D1b, which introduced semi-dwarfism during the green revolution. The adoption of Norin 61 as a reference assembly for functional and evolutionary studies will enable comprehensive characterization of the underexploited Asian bread wheat diversity.

An improved assembly and annotation of the allohexaploid wheat genome identifies complete families of agronomic genes and provides genomic evidence for chromosomal translocations.

Advances in genome sequencing and assembly technologies are generating many high-quality genome sequences, but assemblies of large, repeat-rich polyploid genomes, such as that of bread wheat, remain fragmented and incomplete. We have generated a new wheat whole-genome shotgun sequence assembly using a combination of optimized data types and an assembly algorithm designed to deal with large and complex genomes. The new assembly represents >78% of the genome with a scaffold N50 of 88.8 kb that has a high fidelity to the input data. Our new annotation combines strand-specific Illumina RNA-seq and Pacific Biosciences (PacBio) full-length cDNAs to identify 104,091 high-confidence protein-coding genes and 10,156 noncoding RNA genes. We confirmed three known and identified one novel genome rearrangements. Our approach enables the rapid and scalable assembly of wheat genomes, the identification of structural variants, and the definition of complete gene models, all powerful resources for trait analysis and breeding of this key global crop.

Under fire-simultaneous volatilome and transcriptome analysis unravels fine-scale responses of tansy chemotypes to dual herbivore attack.

BACKGROUND: Tansy plants (Tanacetum vulgare L.) are known for their high intraspecific chemical variation, especially of volatile organic compounds (VOC) from the terpenoid compound group. These VOCs are closely involved in plant-insect interactions and, when profiled, can be used to classify plants into groups known as chemotypes. Tansy chemotypes have been shown to influence plant-aphid interactions, however, to date no information is available on the response of different tansy chemotypes to simultaneous herbivory by more than one insect species. RESULTS: Using a multi-cuvette system, we investigated the responses of five tansy chemotypes to feeding by sucking and/or chewing herbivores (aphids and caterpillars; Metopeurum fuscoviride Stroyan and Spodoptera littoralis Boisduval). Herbivory by caterpillars following aphid infestation led to a plant chemotype-specific change in the patterns of terpenoids stored in trichome hairs and in VOC emissions. The transcriptomic analysis of a plant chemotype represents the first de novo assembly of a transcriptome in tansy and demonstrates priming effects of aphids on a subsequent herbivory. Overall, we show that the five chemotypes do not react in the same way to the two herbivores. As expected, we found that caterpillar feeding increased VOC emissions, however, a priori aphid infestation only led to a further increase in VOC emissions for some chemotypes. CONCLUSIONS: We were able to show that different chemotypes respond to the double herbivore attack in different ways, and that pre-treatment with aphids had a priming effect on plants when they were subsequently exposed to a chewing herbivore. If neighbouring chemotypes in a field population react differently to herbivory/dual herbivory, this could possibly have effects from the individual level to the group level. Individuals of some chemotypes may respond more efficiently to herbivory stress than others, and in a group environment these "louder" chemotypes may affect the local insect community, including the natural enemies of herbivores, and other neighbouring plants.

Expression Pattern Similarities Support the Prediction of Orthologs Retaining Common Functions after Gene Duplication Events.

The identification of functionally equivalent, orthologous genes (functional orthologs) across genomes is necessary for accurate transfer of experimental knowledge from well-characterized organisms to others. This frequently relies on automated, coding sequence-based approaches such as OrthoMCL, Inparanoid, and KOG, which usually work well for one-to-one homologous states. However, this strategy does not reliably work for plants due to the occurrence of extensive gene/genome duplication. Frequently, for one query gene, multiple orthologous genes are predicted in the other genome, and it is not clear a priori from sequence comparison and similarity which one preserves the ancestral function. We have studied 11 organ-dependent and stress-induced gene expression patterns of 286 Arabidopsis lyrata duplicated gene groups and compared them with the respective Arabidopsis (Arabidopsis thaliana) genes to predict putative expressologs and nonexpressologs based on gene expression similarity. Promoter sequence divergence as an additional tool to substantiate functional orthology only partially overlapped with expressolog classification. By cloning eight A. lyrata homologs and complementing them in the respective four Arabidopsis loss-of-function mutants, we experimentally proved that predicted expressologs are indeed functional orthologs, while nonexpressologs or nonfunctionalized orthologs are not. Our study demonstrates that even a small set of gene expression data in addition to sequence homologies are instrumental in the assignment of functional orthologs in the presence of multiple orthologs.

The Sorghum bicolor genome and the diversification of grasses.

Sorghum, an African grass related to sugar cane and maize, is grown for food, feed, fibre and fuel. We present an initial analysis of the approximately 730-megabase Sorghum bicolor (L.) Moench genome, placing approximately 98% of genes in their chromosomal context using whole-genome shotgun sequence validated by genetic, physical and syntenic information. Genetic recombination is largely confined to about one-third of the sorghum genome with gene order and density similar to those of rice. Retrotransposon accumulation in recombinationally recalcitrant heterochromatin explains the approximately 75% larger genome size of sorghum compared with rice. Although gene and repetitive DNA distributions have been preserved since palaeopolyploidization approximately 70 million years ago, most duplicated gene sets lost one member before the sorghum-rice divergence. Concerted evolution makes one duplicated chromosomal segment appear to be only a few million years old. About 24% of genes are grass-specific and 7% are sorghum-specific. Recent gene and microRNA duplications may contribute to sorghum's drought tolerance.

GABI-DUPLO: a collection of double mutants to overcome genetic redundancy in Arabidopsis thaliana.

Owing to duplication events in its progenitor, more than 90% of the genes in the Arabidopsis thaliana genome are members of multigene families. A set of 2108 gene families, each consisting of precisely two unlinked paralogous genes, was identified in the nuclear genome of A. thaliana on the basis of sequence similarity. A systematic method for the creation of double knock-out lines for such gene pairs, designated as DUPLO lines, was established and 200 lines are now publicly available. Their initial phenotypic characterisation led to the identification of seven lines with defects that emerge only in the adult stage. A further six lines display seedling lethality and 23 lines were lethal before germination. Another 14 lines are known to show phenotypes under non-standard conditions or at the molecular level. Knock-out of gene pairs with very similar coding sequences or expression profiles is more likely to produce a mutant phenotype than inactivation of gene pairs with dissimilar profiles or sequences. High coding sequence similarity and highly similar expression profiles are only weakly correlated, implying that promoter and coding regions of these gene pairs display different degrees of diversification.

A powerful tool for genome analysis in maize: development and evaluation of the high density 600 k SNP genotyping array.

BACKGROUND: High density genotyping data are indispensable for genomic analyses of complex traits in animal and crop species. Maize is one of the most important crop plants worldwide, however a high density SNP genotyping array for analysis of its large and highly dynamic genome was not available so far. RESULTS: We developed a high density maize SNP array composed of 616,201 variants (SNPs and small indels). Initially, 57 M variants were discovered by sequencing 30 representative temperate maize lines and then stringently filtered for sequence quality scores and predicted conversion performance on the array resulting in the selection of 1.2 M polymorphic variants assayed on two screening arrays. To identify high-confidence variants, 285 DNA samples from a broad genetic diversity panel of worldwide maize lines including the samples used for sequencing, important founder lines for European maize breeding, hybrids, and proprietary samples with European, US, semi-tropical, and tropical origin were used for experimental validation. We selected 616 k variants according to their performance during validation, support of genotype calls through sequencing data, and physical distribution for further analysis and for the design of the commercially available Affymetrix® Axiom® Maize Genotyping Array. This array is composed of 609,442 SNPs and 6,759 indels. Among these are 116,224 variants in coding regions and 45,655 SNPs of the Illumina® MaizeSNP50 BeadChip for study comparison. In a subset of 45,974 variants, apart from the target SNP additional off-target variants are detected, which show only a minor bias towards intermediate allele frequencies. We performed principal coordinate and admixture analyses to determine the ability of the array to detect and resolve population structure and investigated the extent of LD within a worldwide validation panel. CONCLUSIONS: The high density Affymetrix® Axiom® Maize Genotyping Array is optimized for European and American temperate maize and was developed based on a diverse sample panel by applying stringent quality filter criteria to ensure its suitability for a broad range of applications. With 600 k variants it is the largest currently publically available genotyping array in crop species.

Development of a high density 600K SNP genotyping array for chicken.

BACKGROUND: High density (HD) SNP genotyping arrays are an important tool for genetic analyses of animals and plants. Although the chicken is one of the most important farm animals, no HD array is yet available for high resolution genetic analysis of this species. RESULTS: We report here the development of a 600 K Affymetrix® Axiom® HD genotyping array designed using SNPs segregating in a wide variety of chicken populations. In order to generate a large catalogue of segregating SNPs, we re-sequenced 243 chickens from 24 chicken lines derived from diverse sources (experimental, commercial broiler and layer lines) by pooling 10-15 samples within each line. About 139 million (M) putative SNPs were detected by mapping sequence reads to the new reference genome (Gallus_gallus_4.0) of which ~78 M appeared to be segregating in different lines. Using criteria such as high SNP-quality score, acceptable design scores predicting high conversion performance in the final array and uniformity of distribution across the genome, we selected ~1.8 M SNPs for validation through genotyping on an independent set of samples (n = 282). About 64% of the SNPs were polymorphic with high call rates (>98%), good cluster separation and stable Mendelian inheritance. Polymorphic SNPs were further analysed for their population characteristics and genomic effects. SNPs with extreme breach of Hardy-Weinberg equilibrium (P < 0.00001) were excluded from the panel. The final array, designed on the basis of these analyses, consists of 580,954 SNPs and includes 21,534 coding variants. SNPs were selected to achieve an essentially uniform distribution based on genetic map distance for both broiler and layer lines. Due to a lower extent of LD in broilers compared to layers, as reported in previous studies, the ratio of broiler and layer SNPs in the array was kept as 3:2. The final panel was shown to genotype a wide range of samples including broilers and layers with over 100 K to 450 K informative SNPs per line. A principal component analysis was used to demonstrate the ability of the array to detect the expected population structure which is an important pre-investigation step for many genome-wide analyses. CONCLUSIONS: This Affymetrix® Axiom® array is the first SNP genotyping array for chicken that has been made commercially available to the public as a product. This array is expected to find widespread usage both in research and commercial application such as in genomic selection, genome-wide association studies, selection signature analyses, fine mapping of QTLs and detection of copy number variants.

Impact of natural genetic variation on the transcriptome of autotetraploid Arabidopsis thaliana.

Polyploidy, the presence of more than two complete sets of chromosomes in an organism, has significantly shaped the genomes of angiosperms during evolution. Two forms of polyploidy are often considered: allopolyploidy, which originates from interspecies hybrids, and autopolyploidy, which originates from intraspecies genome duplication events. Besides affecting genome organization, polyploidy generates other genetic effects. Synthetic allopolyploid plants exhibit considerable transcriptome alterations, part of which are likely caused by the reunion of previously diverged regulatory hierarchies. In contrast, autopolyploids have relatively uniform genomes, suggesting lower alteration of gene expression. To evaluate the impact of intraspecies genome duplication on the transcriptome, we generated a series of unique Arabidopsis thaliana autotetraploids by using different ecotypes. A. thaliana autotetraploids show transcriptome alterations that strongly depend on their parental genome composition and include changed expression of both new genes and gene groups previously described from allopolyploid Arabidopsis. Alterations in gene expression are stable, nonstochastic, developmentally specific, and associated with changes in DNA methylation. We propose that Arabidopsis possesses an inherent and heritable ability to sense and respond to elevated, yet balanced chromosome numbers. The impact of natural variation on alteration of autotetraploid gene expression stresses its potential importance in the evolution and breeding of plants.

Intracompartmental and intercompartmental transcriptional networks coordinate the expression of genes for organellar functions.

Genes for mitochondrial and chloroplast proteins are distributed between the nuclear and organellar genomes. Organelle biogenesis and metabolism, therefore, require appropriate coordination of gene expression in the different compartments to ensure efficient synthesis of essential multiprotein complexes of mixed genetic origin. Whereas organelle-to-nucleus signaling influences nuclear gene expression at the transcriptional level, organellar gene expression (OGE) is thought to be primarily regulated posttranscriptionally. Here, we show that intracompartmental and intercompartmental transcriptional networks coordinate the expression of genes for organellar functions. Nearly 1,300 ATH1 microarray-based transcriptional profiles of nuclear and organellar genes for mitochondrial and chloroplast proteins in the model plant Arabidopsis (Arabidopsis thaliana) were analyzed. The activity of genes involved in organellar energy production (OEP) or OGE in each of the organelles and in the nucleus is highly coordinated. Intracompartmental networks that link the OEP and OGE gene sets serve to synchronize the expression of nucleus- and organelle-encoded proteins. At a higher regulatory level, coexpression of organellar and nuclear OEP/OGE genes typically modulates chloroplast functions but affects mitochondria only when chloroplast functions are perturbed. Under conditions that induce energy shortage, the intercompartmental coregulation of photosynthesis genes can even override intracompartmental networks. We conclude that dynamic intracompartmental and intercompartmental transcriptional networks for OEP and OGE genes adjust the activity of organelles in response to the cellular energy state and environmental stresses, and we identify candidate cis-elements involved in the transcriptional coregulation of nuclear genes. Regarding the transcriptional regulation of chloroplast genes, novel tentative target genes of σ factors are identified.

Chromosome-scale assembly of barley cv. 'Haruna Nijo' as a resource for barley genetics.

Cultivated barley (Hordeum vulgare ssp. vulgare) is used for food, animal feed, and alcoholic beverages and is widely grown in temperate regions. Both barley and its wild progenitor (H. vulgare ssp. spontaneum) have large 5.1-Gb genomes. High-quality chromosome-scale assemblies for several representative barley genotypes, both wild and domesticated, have been constructed recently to populate the nascent barley pan-genome infrastructure. Here, we release a chromosome-scale assembly of the Japanese elite malting barley cultivar 'Haruna Nijo' using a similar methodology as in the barley pan-genome project. The 4.28-Gb assembly had a scaffold N50 size of 18.9 Mb. The assembly showed high collinearity with the barley reference genome 'Morex' cultivar, with some inversions. The pseudomolecule assembly was characterized using transcript evidence of gene projection derived from the reference genome and de novo gene annotation achieved using published full-length cDNA sequences and RNA-Seq data for 'Haruna Nijo'. We found good concordance between our whole-genome assembly and the publicly available BAC clone sequence of 'Haruna Nijo'. Interesting phenotypes have since been identified in Haruna Nijo; its genome sequence assembly will facilitate the identification of the underlying genes.

The Barley and Wheat Pan-Genomes.

To unlock the genetic potential in crops, multi-genome comparisons are an essential tool. Decreasing costs and improved sequencing technologies have democratized plant genome sequencing and led to a vast increase in the amount of available reference sequences on the one hand and enabled the assembly of even the largest and most complex and repetitive crops genomes such as wheat and barley. These developments have led to the era of pan-genomics in recent years. Pan-genome projects enable the definition of the core and dispensable genome for various crop species as well as the analysis of structural and functional variation and hence offer unprecedented opportunities for exploring and utilizing the genetic basis of natural variation in crops. Comparing, analyzing, and visualizing these multiple reference genomes and their diversity requires powerful and specialized computational strategies and tools.

A reference-guided TILLING by amplicon-sequencing platform supports forward and reverse genetics in barley.

Barley is a diploid species with a genome smaller than those of other members of the Triticeae tribe, making it an attractive model for genetic studies in Triticeae crops. The recent development of barley genomics has created a need for a high-throughput platform to identify genetically uniform mutants for gene function investigations. In this study, we report an ethyl methanesulfonate (EMS)-mutagenized population consisting of 8525 M3 lines in the barley landrace "Hatiexi" (HTX), which we complement with a high-quality de novo assembly of a reference genome for this genotype. The mutation rate within the population ranged from 1.51 to 4.09 mutations per megabase, depending on the treatment dosage of EMS and the mutation discrimination platform used for genotype analysis. We implemented a three-dimensional DNA pooling strategy combined with multiplexed amplicon sequencing to create a highly efficient and cost-effective TILLING (targeting induced locus lesion in genomes) platform in barley. Mutations were successfully identified from 72 mixed amplicons within a DNA pool containing 64 individual mutants and from 56 mixed amplicons within a pool containing 144 individuals. We discovered abundant allelic mutants for dozens of genes, including the barley Green Revolution contributor gene Brassinosteroid insensitive 1 (BRI1). As a proof of concept, we rapidly determined the causal gene responsible for a chlorotic mutant by following the MutMap strategy, demonstrating the value of this resource to support forward and reverse genetic studies in barley.

Identification and genetic analysis of the APOSPORY locus in Hypericum perforatum L.

The introduction of apomixis - seed formation without fertilization - into crop plants is a long-held goal of breeding research, since it would allow for the ready fixation of heterozygosity. The genetic basis of apomixis, whether of the aposporous or the diplosporous type, is still only poorly understood. Hypericum perforatum (St John's wort), a plant with a small genome and a short generation time, can be aposporous and/or parthenogenetic, and so represents an interesting model dicot for apomixis research. Here we describe a genetic analysis which first defined and then isolated a locus (designated HAPPY for Hypericum APOSPORY) associated with apospory. Amplified fragment length polymorphism (AFLP) profiling was used to generate a cleaved amplified polymorphic sequence (CAPS) marker for HAPPY which co-segregated with apospory but not with parthenogenesis, showing that these two components of apomixis are independently controlled. Apospory was inherited as a dominant simplex gene at the tetraploid level. Part of the HAPPY sequence is homologous to the Arabidopsis thaliana gene ARI7 encoding the ring finger protein ARIADNE7. This protein is predicted to be involved in various regulatory processes, including ubiquitin-mediated protein degradation. While the aposporous and sexual alleles of the HAPPY component HpARI were co-expressed in many parts of the plant, the gene product of the apomict's allele is truncated. Cloning HpARI represents the first step towards the full characterization of HAPPY and the elucidation of the molecular mechanisms underlying apomixis in H. perforatum.