Genetic mapping of the LOSS OF PARTHENOGENESIS locus in Pilosella piloselloides and the evolution of apomixis in the Lactuceae.

Pilosella piloselloides var. praealta (syn. P. praealta; Hieracium praealtum) is a versatile model used to study gametophytic apomixis. In this system apomixis is controlled by three loci: one that controls the avoidance of meiosis (LOA), one that controls the avoidance of fertilization (LOP) and a third that controls autonomous endosperm formation (AutE). Using a unique polyhaploid mapping approach the LOP locus was mapped to a 654 kb genomic interval syntenic to linkage group 8 of Lactuca sativa. Polyhaploids form through the gametophytic action of a dominant determinant at LOP, so the mapped region represents both a functional and a physical domain for LOP in P. piloselloides. Allele sequence divergence (ASD) analysis of the PARTHENOGENESIS (PAR) gene within the LOP locus revealed that dominant PAR alleles in Pilosella remain highly similar across the genus, whilst the recessive alleles are more divergent. A previous report noted that dominant PAR alleles in both Pilosella and Taraxacum are modified by the presence of a class II transposable element (TE) in the promoter of the gene. This observation was confirmed and further extended to the related genus Hieracium. Sufficient differences were noted in the structure and location of the TE elements to conclude that TE insertional events had occurred independently in the three genera. Measures of allele crossover amongst the polyhaploids revealed that P. piloselloides is an autopolyploid species with tetrasomic inheritance. It was also noted that the dominant determinant of LOP in P. piloselloides could transmit via a diploid gamete (pollen or egg) but not via a haploid gamete. Using this information, a model is presented of how gametophytic apomixis may have evolved in several members of the Lactuceae, a tribe of the Asteraceae.

Recurrent neo-sex chromosome evolution in kiwifruit.

Sex chromosome evolution is thought to be tightly associated with the acquisition and maintenance of sexual dimorphisms. Plant sex chromosomes have evolved independently in many lineages1,2 and can provide a powerful comparative framework to study this. We assembled and annotated genome sequences of three kiwifruit species (genus Actinidia) and uncovered recurrent sex chromosome turnovers in multiple lineages. Specifically, we observed structural evolution of the neo-Y chromosomes, which was driven via rapid bursts of transposable element insertions. Surprisingly, sexual dimorphisms were conserved in the different species studied, despite the fact that the partially sex-linked genes differ between them. Using gene editing in kiwifruit, we demonstrated that one of the two Y-chromosome-encoded sex-determining genes, Shy Girl, shows pleiotropic effects that can explain the conserved sexual dimorphisms. These plant sex chromosomes therefore maintain sexual dimorphisms through the conservation of a single gene, without a process involving interactions between separate sex-determining genes and genes for sexually dimorphic traits.

The Relative Power of Structural Genomic Variation versus SNPs in Explaining the Quantitative Trait Growth in the Marine Teleost Chrysophrys auratus.

BACKGROUND: Genetic diversity provides the basic substrate for evolution. Genetic variation consists of changes ranging from single base pairs (single-nucleotide polymorphisms, or SNPs) to larger-scale structural variants, such as inversions, deletions, and duplications. SNPs have long been used as the general currency for investigations into how genetic diversity fuels evolution. However, structural variants can affect more base pairs in the genome than SNPs and can be responsible for adaptive phenotypes due to their impact on linkage and recombination. In this study, we investigate the first steps needed to explore the genetic basis of an economically important growth trait in the marine teleost finfish Chrysophrys auratus using both SNP and structural variant data. Specifically, we use feature selection methods in machine learning to explore the relative predictive power of both types of genetic variants in explaining growth and discuss the feature selection results of the evaluated methods. METHODS: SNP and structural variant callers were used to generate catalogues of variant data from 32 individual fish at ages 1 and 3 years. Three feature selection algorithms (ReliefF, Chi-square, and a mutual-information-based method) were used to reduce the dataset by selecting the most informative features. Following this selection process, the subset of variants was used as features to classify fish into small, medium, or large size categories using KNN, naïve Bayes, random forest, and logistic regression. The top-scoring features in each feature selection method were subsequently mapped to annotated genomic regions in the zebrafish genome, and a permutation test was conducted to see if the number of mapped regions was greater than when random sampling was applied. RESULTS: Without feature selection, the prediction accuracies ranged from 0 to 0.5 for both structural variants and SNPs. Following feature selection, the prediction accuracy increased only slightly to between 0 and 0.65 for structural variants and between 0 and 0.75 for SNPs. The highest prediction accuracy for the logistic regression was achieved for age 3 fish using SNPs, although generally predictions for age 1 and 3 fish were very similar (ranging from 0-0.65 for both SNPs and structural variants). The Chi-square feature selection of SNP data was the only method that had a significantly higher number of matches to annotated genomic regions of zebrafish than would be explained by chance alone. CONCLUSIONS: Predicting a complex polygenic trait such as growth using data collected from a low number of individuals remains challenging. While we demonstrate that both SNPs and structural variants provide important information to help understand the genetic basis of phenotypic traits such as fish growth, the full complexities that exist within a genome cannot be easily captured by classical machine learning techniques. When using high-dimensional data, feature selection shows some increase in the prediction accuracy of classification models and provides the potential to identify unknown genomic correlates with growth. Our results show that both SNPs and structural variants significantly impact growth, and we therefore recommend that researchers interested in the genotype-phenotype map should strive to go beyond SNPs and incorporate structural variants in their studies as well. We discuss how our machine learning models can be further expanded to serve as a test bed to inform evolutionary studies and the applied management of species.

A PARTHENOGENESIS allele from apomictic dandelion can induce egg cell division without fertilization in lettuce.

Apomixis, the clonal formation of seeds, is a rare yet widely distributed trait in flowering plants. We have isolated the PARTHENOGENESIS (PAR) gene from apomictic dandelion that triggers embryo development in unfertilized egg cells. PAR encodes a K2-2 zinc finger, EAR-domain protein. Unlike the recessive sexual alleles, the dominant PAR allele is expressed in egg cells and has a miniature inverted-repeat transposable element (MITE) transposon insertion in the promoter. The MITE-containing promoter can invoke a homologous gene from sexual lettuce to complement dandelion LOSS OF PARTHENOGENESIS mutants. A similar MITE is also present in the promoter of the PAR gene in apomictic forms of hawkweed, suggesting a case of parallel evolution. Heterologous expression of dandelion PAR in lettuce egg cells induced haploid embryo-like structures in the absence of fertilization. Sexual PAR alleles are expressed in pollen, suggesting that the gene product releases a block on embryogenesis after fertilization in sexual species while in apomictic species PAR expression triggers embryogenesis in the absence of fertilization.

The genome of New Zealand trevally (Carangidae: Pseudocaranx georgianus) uncovers a XY sex determination locus.

BACKGROUND: The genetic control of sex determination in teleost species is poorly understood. This is partly because of the diversity of mechanisms that determine sex in this large group of vertebrates, including constitutive genes linked to sex chromosomes, polygenic constitutive mechanisms, environmental factors, hermaphroditism, and unisexuality. Here we use a de novo genome assembly of New Zealand silver trevally (Pseudocaranx georgianus) together with sex-specific whole genome sequencing data to detect sexually divergent genomic regions, identify candidate genes and develop molecular makers. RESULTS: The de novo assembly of an unsexed trevally (Trevally_v1) resulted in a final assembly of 579.4 Mb in length, with a N50 of 25.2 Mb. Of the assembled scaffolds, 24 were of chromosome scale, ranging from 11 to 31 Mb in length. A total of 28,416 genes were annotated after 12.8 % of the assembly was masked with repetitive elements. Whole genome re-sequencing of 13 wild sexed trevally (seven males and six females) identified two sexually divergent regions located on two scaffolds, including a 6 kb region at the proximal end of chromosome 21. Blast analyses revealed similarity between one region and the aromatase genes cyp19 (a1a/b) (E-value < 1.00E-25, identity > 78.8 %). Males contained higher numbers of heterozygous variants in both regions, while females showed regions of very low read-depth, indicative of male-specificity of this genomic region. Molecular markers were developed and subsequently tested on 96 histologically-sexed fish (42 males and 54 females). Three markers amplified in absolute correspondence with sex (positive in males, negative in females). CONCLUSIONS: The higher number of heterozygous variants in males combined with the absence of these regions in females support a XY sex-determination model, indicating that the trevally_v1 genome assembly was developed from a male specimen. This sex system contrasts with the ZW sex-determination model documented in closely related carangid species. Our results indicate a sex-determining function of a cyp19a1a-like gene, suggesting the molecular pathway of sex determination is somewhat conserved in this family. The genomic resources developed here will facilitate future comparative work, and enable improved insights into the varied sex determination pathways in teleosts. The sex marker developed in this study will be a valuable resource for aquaculture selective breeding programmes, and for determining sex ratios in wild populations.

Molecular Characterisation of a Supergene Conditioning Super-High Vitamin C in Kiwifruit Hybrids.

During analysis of kiwifruit derived from hybrids between the high vitamin C (ascorbic acid; AsA) species Actinidia eriantha and A. chinensis, we observed bimodal segregation of fruit AsA concentration suggesting major gene segregation. To test this hypothesis, we performed whole-genome sequencing on pools of hybrid genotypes with either high or low AsA fruit. Pool-GWAS (genome-wide association study) revealed a single Quantitative Trait Locus (QTL) spanning more than 5 Mbp on chromosome 26, which we denote as qAsA26.1. A co-dominant PCR marker was used to validate this association in four diploid (A. chinensis × A. eriantha) × A. chinensis backcross families, showing that the A. eriantha allele at this locus increases fruit AsA levels by 250 mg/100 g fresh weight. Inspection of genome composition and recombination in other A. chinensis genetic maps confirmed that the qAsA26.1 region bears hallmarks of suppressed recombination. The molecular fingerprint of this locus was examined in leaves of backcross validation families by RNA sequencing (RNASEQ). This confirmed strong allelic expression bias across this region as well as differential expression of transcripts on other chromosomes. This evidence suggests that the region harbouring qAsA26.1 constitutes a supergene, which may condition multiple pleiotropic effects on metabolism.

The genomic pool of standing structural variation outnumbers single nucleotide polymorphism by threefold in the marine teleost Chrysophrys auratus.

Recent studies have highlighted an important role of structural variation (SV) in ecological and evolutionary processes, but few have studied nonmodel species in the wild. As part of our long-term research programme on the nonmodel teleost fish Australasian snapper (Chrysophrys auratus), we aim to build one of the first catalogues of genomic variants (SNPs and indels, and deletions, duplications and inversions) in fishes and evaluate overlap of genomic variants with regions under putative selection (Tajima's D and π), and coding sequences (genes). For this, we analysed six males and six females from three locations in New Zealand and generated a high-resolution genomic variation catalogue. We characterized 20,385 SVs and found they intersected with almost a third of all annotated genes. Together with small indels, SVs account for three times more variation in the genome in terms of bases affected compared to SNPs. We found that a sizeable portion of detected SVs was in the upper and lower genomic regions of Tajima's D and π, indicating that some of these have an effect on the phenotype. Together, these results shed light on the often neglected genomic variation that is produced by SVs and highlights the need to go beyond the mere measure of SNPs when investigating evolutionary processes, such as species diversification and adaptation.

Comparative transcriptome analysis of the wild-type model apomict Hieracium praealtum and its loss of parthenogenesis (lop) mutant.

BACKGROUND: Asexual seed formation (apomixis) has been observed in diverse plant families but is rare in crop plants. The generation of apomictic crops would revolutionize agriculture, as clonal seed production provides a low cost and efficient way to produce hybrid seed. Hieracium (Asteraceae) is a model system for studying the molecular components of gametophytic apomixis (asexual seed reproduction). RESULTS: In this study, a reference transcriptome was produced from apomictic Hieracium undergoing the key apomictic events of apomeiosis, parthenogenesis and autonomous endosperm development. In addition, transcriptome sequences from pre-pollination and post-pollination stages were generated from a loss of parthenogenesis (lop) mutant accession that exhibits loss of parthenogenesis and autonomous endosperm development. The transcriptome is composed of 147,632 contigs, 50% of which were annotated with orthologous genes and their probable function. The transcriptome was used to identify transcripts differentially expressed during apomictic and pollination dependent (lop) seed development. Gene Ontology enrichment analysis of differentially expressed transcripts showed that an important difference between apomictic and pollination dependent seed development was the expression of genes relating to epigenetic gene regulation. Genes that mark key developmental stages, i.e. aposporous embryo sac development and seed development, were also identified through their enhanced expression at those stages. CONCLUSION: The production of a comprehensive floral reference transcriptome for Hieracium provides a valuable resource for research into the molecular basis of apomixis and the identification of the genes underlying the LOP locus.

A manually annotated Actinidia chinensis var. chinensis (kiwifruit) genome highlights the challenges associated with draft genomes and gene prediction in plants.

BACKGROUND: Most published genome sequences are drafts, and most are dominated by computational gene prediction. Draft genomes typically incorporate considerable sequence data that are not assigned to chromosomes, and predicted genes without quality confidence measures. The current Actinidia chinensis (kiwifruit) 'Hongyang' draft genome has 164 Mb of sequences unassigned to pseudo-chromosomes, and omissions have been identified in the gene models. RESULTS: A second genome of an A. chinensis (genotype Red5) was fully sequenced. This new sequence resulted in a 554.0 Mb assembly with all but 6 Mb assigned to pseudo-chromosomes. Pseudo-chromosomal comparisons showed a considerable number of translocation events have occurred following a whole genome duplication (WGD) event some consistent with centromeric Robertsonian-like translocations. RNA sequencing data from 12 tissues and ab initio analysis informed a genome-wide manual annotation, using the WebApollo tool. In total, 33,044 gene loci represented by 33,123 isoforms were identified, named and tagged for quality of evidential support. Of these 3114 (9.4%) were identical to a protein within 'Hongyang' The Kiwifruit Information Resource (KIR v2). Some proportion of the differences will be varietal polymorphisms. However, as most computationally predicted Red5 models required manual re-annotation this proportion is expected to be small. The quality of the new gene models was tested by fully sequencing 550 cloned 'Hort16A' cDNAs and comparing with the predicted protein models for Red5 and both the original 'Hongyang' assembly and the revised annotation from KIR v2. Only 48.9% and 63.5% of the cDNAs had a match with 90% identity or better to the original and revised 'Hongyang' annotation, respectively, compared with 90.9% to the Red5 models. CONCLUSIONS: Our study highlights the need to take a cautious approach to draft genomes and computationally predicted genes. Our use of the manual annotation tool WebApollo facilitated manual checking and correction of gene models enabling improvement of computational prediction. This utility was especially relevant for certain types of gene families such as the EXPANSIN like genes. Finally, this high quality gene set will supply the kiwifruit and general plant community with a new tool for genomics and other comparative analysis.

Seeds of doubt: Mendel's choice of Hieracium to study inheritance, a case of right plant, wrong trait.

KEY MESSAGE: In this review, we explore Gregor Mendel's hybridization experiments with Hieracium , update current knowledge on apomictic reproduction and describe approaches now being used to develop true-breeding hybrid crops. From our perspective, it is easy to conclude that Gregor Mendel's work on pea was insightful, but his peers clearly did not regard it as being either very convincing or of much importance. One apparent criticism was that his findings only applied to pea. We know from a letter he wrote to Carl von Nägeli, a leading botanist, that he believed he needed to "verify, with other plants, the results obtained with Pisum". For this purpose, Mendel adopted Hieracium subgenus Pilosella, a phenotypically diverse taxon under botanical study at the time. What Mendel could not have known, however, is that the majority of these plants are not sexual plants like pea, but instead are facultatively apomictic. In these forms, the majority of seed arises asexually, and such progeny are, therefore, clones of the maternal parent. Mendel obtained very few hybrids in his Hieracium crosses, yet we calculate that he probably emasculated in excess of 5000 Hieracium florets to even obtain the numbers he did. Despite that effort, he was perplexed by the results, and they ultimately led him to conclude that "the hybrids of Hieracium show a behaviour exactly opposite to those of Pisum". Apomixis is now a topic of intense research interest, and in an ironic twist of history, Hieracium subgenus Pilosella has been developed as a molecular model to study this trait. In this paper, we explore further Mendel's hybridization experiments with Hieracium, update current knowledge on apomictic reproduction and describe approaches now being used to develop true-breeding hybrid crops.

The sweet potato IbMYB1 gene as a potential visible marker for sweet potato intragenic vector system.

MYB transcription factors play important roles in transcriptional regulation of many secondary metabolites including anthocyanins. We cloned the R2R3-MYB type IbMYB1 complementary DNAs from the purple-fleshed sweet potato (Ipomoea batatas L. cv Sinzami) and investigated the expression patterns of IbMYB1 gene with IbMYB1a and IbMYB1b splice variants in leaf and root tissues of various sweet potato cultivars by reverse transcription-polymerase chain reaction. The transcripts of IbMYB1 were predominantly expressed in the purple-fleshed storage roots and they were also detectable in the leaf tissues accumulating anthocyanin pigments. In addition, transcript levels of IbMYB1 gene were up-regulated by treatment with methyl jasmonate or salicylic acid in leaf and root tissues of cv. White Star. To set up the intragenic vector system in sweet potato, we first evaluated the utilization of the IbMYB1 gene as a visible selectable marker. The IbMYB1a was transiently expressed in tobacco leaves under the control of a constitutive cauliflower mosaic virus 35S promoter, a root-specific and sucrose-inducible sporamin promoter, and an oxidative stress-inducible sweet potato anionic peroxidase2 promoter. We also showed that overexpression of IbMYB1a induced massive anthocyanin pigmentation in tobacco leaves and up-regulated the transcript levels of the structural genes in anthocyanin biosynthetic pathway. Furthermore, high-performance liquid chromatography analysis revealed that the expression of IbMYB1a led to production of cyanidin as a major core molecule of anthocyanidins in tobacco leaves. These results suggest that the IbMYB1 gene can be applicable to a visible marker for sweet potato transformation with intragenic vectors, as well as the production of anthocyanin as important nutritive value in other plant species.

Deletion mapping of genetic regions associated with apomixis in Hieracium.

Although apomixis has been quoted as a technology with the potential to deliver benefits similar in scale to those achieved with the Green Revolution, very little is currently known of the genetic mechanisms that control this trait in plants. To address this issue, we developed Hieracium, a genus of daisies native to Eurasia and North America, as a genetic model to study apomixis. In a molecular mapping study, we defined the number of genetic loci involved in apomixis, and we explored dominance and linkage relationships between these loci. To avoid difficulties often encountered with inheritance studies of apomicts, we based our mapping effort on the use of deletion mutagenesis, coupled with amplified fragment length polymorphism (AFLP) as a genomic fingerprinting tool. The results indicate that apomixis in Hieracium caespitosum is controlled at two principal loci, one of which regulates events associated with the avoidance of meiosis (apomeiosis) and the other, an unlinked locus that controls events associated with the avoidance of fertilization (parthenogenesis). AFLP bands identified as central to both loci were isolated, sequenced, and used to develop sequence-characterized amplified region (SCAR) markers. The validity of the AFLP markers was verified by using a segregating population generated by hybridization. The validity of the SCAR markers was verified by their pattern of presence/absence in specific mutants. The mutants, markers, and genetic data derived from this work are now being used to isolate genes controlling apomixis in this system.

Agrobacterium tumefaciens-mediated transformation of leek (Allium porrum) and garlic (Allium sativum).

Transgenic leek (Allium porrum) and garlic (Allium sativum) plants have been recovered by the selective culturing of immature leek and garlic embryos via Agrobacterium-mediated transformation using a method similar to that described by Eady et al. (Plant Cell Rep 19:376-381, 2000) for onion transformation. This method involved the use of a binary vector containing the m-gfp-ER reporter gene and nptII selectable marker, and followed the protocol developed previously for the transformation of onions with only minor modifications pertaining to the post-transformation selection procedure which was simplified to have just a single selection regime. Transgenic cultures were selected for their ability to express the m-gfp-ER reporter gene and grown in the presence of geneticin (20 mg/l). The presence of transgenes in the genome of the plants was confirmed using TAIL-PCR and Southern analysis. This is the first report of leek and "true seed" garlic transformation. It now makes possible the integration of useful agronomic and quality traits into these crops.

A unique set of 11,008 onion expressed sequence tags reveals expressed sequence and genomic differences between the monocot orders Asparagales and Poales.

Enormous genomic resources have been developed for plants in the monocot order Poales; however, it is not clear how representative the Poales are for the monocots as a whole. The Asparagales are a monophyletic order sister to the lineage carrying the Poales and possess economically important plants such as asparagus, garlic, and onion. To assess the genomic differences between the Asparagales and Poales, we generated 11,008 unique ESTs from a normalized cDNA library of onion. Sequence analyses of these ESTs revealed microsatellite markers, single nucleotide polymorphisms, and homologs of transposable elements. Mean nucleotide similarity between rice and the Asparagales was 78% across coding regions. Expressed sequence and genomic comparisons revealed strong differences between the Asparagales and Poales for codon usage and mean GC content, GC distribution, and relative GC content at each codon position, indicating that genomic characteristics are not uniform across the monocots. The Asparagales were more similar to eudicots than to the Poales for these genomic characteristics.