AssemblyQC: A Nextflow pipeline for reproducible reporting of assembly quality.

SUMMARY: Genome assembly projects have grown exponentially due to breakthroughs in sequencing technologies and assembly algorithms. Evaluating the quality of genome assemblies is critical to ensure the reliability of downstream analysis and interpretation. To fulfil this task, we have developed the AssemblyQC pipeline that performs file-format validation, contaminant checking, contiguity measurement, gene- and repeat-space completeness quantification, telomere inspection, taxonomic assignment, synteny alignment, scaffold examination through Hi-C contact-map visualization, and assessments of completeness, consensus quality and phasing through k-mer analysis. It produces a comprehensive HTML report with method descriptions, tables, and visualizations. AVAILABILITY AND IMPLEMENTATION: The pipeline uses Nextflow for workflow orchestration and adheres to the best-practice established by the nf-core community. This pipeline offers a reproducible, scalable, and portable method to assess the quality of genome assemblies-the code is available online at GitHub: https://github.com/Plant-Food-Research-Open/assemblyqc. SUPPLEMENTARY INFORMATION: Pipeline usage documentation, parameter descriptions and example outputs are available on GitHub: https://github.com/Plant-Food-Research-Open/assemblyqc/tree/main/docs. A preview report is also hosted online: https://plant-food-research-open.github.io/assemblyqc.

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.

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.

Assembling large genomes: analysis of the stick insect (Clitarchus hookeri) genome reveals a high repeat content and sex-biased genes associated with reproduction.

BACKGROUND: Stick insects (Phasmatodea) have a high incidence of parthenogenesis and other alternative reproductive strategies, yet the genetic basis of reproduction is poorly understood. Phasmatodea includes nearly 3000 species, yet only the genome of Timema cristinae has been published to date. Clitarchus hookeri is a geographical parthenogenetic stick insect distributed across New Zealand. Sexual reproduction dominates in northern habitats but is replaced by parthenogenesis in the south. Here, we present a de novo genome assembly of a female C. hookeri and use it to detect candidate genes associated with gamete production and development in females and males. We also explore the factors underlying large genome size in stick insects. RESULTS: The C. hookeri genome assembly was 4.2 Gb, similar to the flow cytometry estimate, making it the second largest insect genome sequenced and assembled to date. Like the large genome of Locusta migratoria, the genome of C. hookeri is also highly repetitive and the predicted gene models are much longer than those from most other sequenced insect genomes, largely due to longer introns. Miniature inverted repeat transposable elements (MITEs), absent in the much smaller T. cristinae genome, is the most abundant repeat type in the C. hookeri genome assembly. Mapping RNA-Seq reads from female and male gonadal transcriptomes onto the genome assembly resulted in the identification of 39,940 gene loci, 15.8% and 37.6% of which showed female-biased and male-biased expression, respectively. The genes that were over-expressed in females were mostly associated with molecular transportation, developmental process, oocyte growth and reproductive process; whereas, the male-biased genes were enriched in rhythmic process, molecular transducer activity and synapse. Several genes involved in the juvenile hormone synthesis pathway were also identified. CONCLUSIONS: The evolution of large insect genomes such as L. migratoria and C. hookeri genomes is most likely due to the accumulation of repetitive regions and intron elongation. MITEs contributed significantly to the growth of C. hookeri genome size yet are surprisingly absent from the T. cristinae genome. Sex-biased genes identified from gonadal tissues, including genes involved in juvenile hormone synthesis, provide interesting candidates for the further study of flexible reproduction in stick insects.

Analysis of the genome of the New Zealand giant collembolan (Holacanthella duospinosa) sheds light on hexapod evolution.

BACKGROUND: The New Zealand collembolan genus Holacanthella contains the largest species of springtails (Collembola) in the world. Using Illumina technology we have sequenced and assembled a draft genome and transcriptome from Holacanthella duospinosa (Salmon). We have used this annotated assembly to investigate the genetic basis of a range of traits critical to the evolution of the Hexapoda, the phylogenetic position of H. duospinosa and potential horizontal gene transfer events. RESULTS: Our genome assembly was ~375 Mbp in size with a scaffold N50 of ~230 Kbp and sequencing coverage of ~180×. DNA elements, LTRs and simple repeats and LINEs formed the largest components and SINEs were very rare. Phylogenomics (370,877 amino acids) placed H. duospinosa within the Neanuridae. We recovered orthologs of the conserved sex determination genes thought to play a role in sex determination. Analysis of CpG content suggested the absence of DNA methylation, and consistent with this we were unable to detect orthologs of the DNA methyltransferase enzymes. The small subunit rRNA gene contained a possible retrotransposon. The Hox gene complex was broken over two scaffolds. For chemosensory ability, at least 15 and 18 ionotropic glutamate and gustatory receptors were identified, respectively. However, we were unable to identify any odorant receptors or their obligate co-receptor Orco. Twenty-three chitinase-like genes were identified from the assembly. Members of this multigene family may play roles in the digestion of fungal cell walls, a common food source for these saproxylic organisms. We also detected 59 and 96 genes that blasted to bacteria and fungi, respectively, but were located on scaffolds that otherwise contained arthropod genes. CONCLUSIONS: The genome of H. duospinosa contains some unusual features including a Hox complex broken over two scaffolds, in a different manner to other arthropod species, a lack of odorant receptor genes and an apparent lack of environmentally responsive DNA methylation, unlike many other arthropods. Our detection of candidate horizontal gene transfer candidates confirms that this phenomenon is occurring across Collembola. These findings allow us to narrow down the regions of the arthropod phylogeny where key innovations have occurred that have facilitated the evolutionary success of Hexapoda.

A microRNA allele that emerged prior to apple domestication may underlie fruit size evolution.

The molecular genetic mechanisms underlying fruit size remain poorly understood in perennial crops, despite size being an important agronomic trait. Here we show that the expression level of a microRNA gene (miRNA172) influences fruit size in apple. A transposon insertional allele of miRNA172 showing reduced expression associates with large fruit in an apple breeding population, whereas over-expression of miRNA172 in transgenic apple significantly reduces fruit size. The transposon insertional allele was found to be co-located with a major fruit size quantitative trait locus, fixed in cultivated apples and their wild progenitor species with relatively large fruit. This finding supports the view that the selection for large size in apple fruit was initiated prior to apple domestication, likely by large mammals, before being subsequently strengthened by humans, and also helps to explain why signatures of genetic bottlenecks and selective sweeps are normally weaker in perennial crops than in annual crops.

Ascochyta blight disease of pea (Pisum sativum L.): defence-related candidate genes associated with QTL regions and identification of epistatic QTL.

KEY MESSAGE: Advances have been made in our understanding of Ascochyta blight resistance genetics through mapping candidate genes associated with QTL regions and demonstrating the importance of epistatic interactions in determining resistance. Ascochyta blight disease of pea (Pisum sativum L.) is economically significant with worldwide distribution. The causal pathogens are Didymella pinodes, Phoma medicaginis var pinodella and, in South Australia, P. koolunga. This study aimed to identify candidate genes that map to quantitative trait loci (QTL) for Ascochyta blight field disease resistance and to explore the role of epistatic interactions. Candidate genes associated with QTL were identified beginning with 101 defence-related genes from the published literature. Synteny between pea and Medicago truncatula was used to narrow down the candidates for mapping. Fourteen pea candidate sequences were mapped in two QTL mapping populations, A26 × Rovar and A88 × Rovar. QTL peaks, or the intervals containing QTL peaks, for the Asc2.1, Asc4.2, Asc4.3 and Asc7.1 QTL were defined by four of these candidate genes, while another three candidate genes occurred within 1.0 LOD confidence intervals. Epistasis involving QTL × background marker and background marker × background marker interactions contributed to the disease response phenotypes observed in the two mapping populations. For each population, five pairwise interactions exceeded the 5% false discovery rate threshold. Two candidate genes were involved in significant pairwise interactions. Markers in three genomic regions were involved in two or more epistatic interactions. Therefore, this study has identified pea defence-related sequences that are candidates for resistance determination, and that may be useful for marker-assisted selection. The demonstration of epistasis informs breeders that the architecture of this complex quantitative resistance includes epistatic interactions with non-additive effects.

Pheromone receptor evolution in the cryptic leafroller species, Ctenopseustis obliquana and C. herana.

How new mate recognition systems evolve when changes are required in both the male and female components remains a conundrum. Here, we investigated the molecular basis of pheromone reception in two species of tortricid (leafroller) moth, Ctenopseustis obliquana and C. herana. Male C. obliquana are attracted to a 90:10 blend of (Z)-8-tetradecenyl acetate (Z8-14:OAc) and (Z)-5-tetradecenyl acetate (Z5-14:OAc), whereas C. herana males are attracted to Z5-14:OAc alone. We used a transcriptome sequencing approach from adult male and female antennae to identify 47 olfactory receptors (ORs) from each species and assessed their expression levels in male and female antennae using RNA-Seq counting and quantitative RT-PCR. Three male-biased and one female-biased OR were identified in C. obliquana by quantitative RT-PCR, and four male-biased and one female-biased receptor in C. herana. The male-biased receptors, CoblOR7, CoblOR30, CherOR7, CherOR30, CherOR1a and CherOR1b were tested for their ability to respond to sex pheromone components in a HEK293 cell calcium assay. CoblOR7 and CherOR7 responded to Z8-14:OAc, however, no receptor for Z5-14:OAc was identified. In addition to Z8-14:OAc, CherOR7 also responded to Z7-14:OAc, indicating that this receptor may be under relaxed constraint. Of the 29 amino acid differences between CoblOR7 and CherOR7, significantly more are located in the third and the sixth transmembrane domain regions. Overall, these findings are consistent with studies revealing the presence of neurons tuned to both Z8-14:OAc and Z5-14:OAc in both species, but that for C. herana males, the ability to detect Z8-14:OAc is currently not required.

Staying green postharvest: how three mutations in the Arabidopsis chlorophyll b reductase gene NYC1 delay degreening by distinct mechanisms.

Stresses such as energy deprivation, wounding and water-supply disruption often contribute to rapid deterioration of harvested tissues. To uncover the genetic regulation behind such stresses, a simple assessment system was used to detect senescence mutants in conjunction with two rapid mapping techniques to identify the causal mutations. To demonstrate the power of this approach, immature inflorescences of Arabidopsis plants that contained ethyl methanesulfonate-induced lesions were detached and screened for altered timing of dark-induced senescence. Numerous mutant lines displaying accelerated or delayed timing of senescence relative to wild type were discovered. The underlying mutations in three of these were identified using High Resolution Melting analysis to map to a chromosomal arm followed by a whole-genome sequencing-based mapping method, termed 'Needle in the K-Stack', to identify the causal lesions. All three mutations were single base pair changes and occurred in the same gene, NON-YELLOW COLORING1 (NYC1), a chlorophyll b reductase of the short-chain dehydrogenase/reductase (SDR) superfamily. This was consistent with the mutants preferentially retaining chlorophyll b, although substantial amounts of chlorophyll b were still lost. The single base pair mutations disrupted NYC1 function by three distinct mechanisms, one by producing a termination codon, the second by interfering with correct intron splicing and the third by replacing a highly conserved proline with a non-equivalent serine residue. This non-synonymous amino acid change, which occurred in the NADPH binding domain of NYC1, is the first example of such a mutation in an SDR protein inhibiting a physiological response in plants.

Genomes of 'Candidatus Liberibacter solanacearum' Haplotype A from New Zealand and the United States Suggest Significant Genome Plasticity in the Species.

'Candidatus Liberibacter solanacearum' contains two solanaceous crop-infecting haplotypes, A and B. Two haplotype A draft genomes were assembled and compared with ZC1 (haplotype B), revealing inversion and relocation genomic rearrangements, numerous single-nucleotide polymorphisms, and differences in phage-related regions. Differences in prophage location and sequence were seen both within and between haplotype comparisons. OrthoMCL and BLAST analyses identified 46 putative coding sequences present in haplotype A that were not present in haplotype B. Thirty-eight of these loci were not found in sequences from other Liberibacter spp. Quantitative polymerase chain reaction (qPCR) assays designed to amplify sequences from 15 of these loci were screened against a panel of 'Ca. L. solanacearum'-positive samples to investigate genetic diversity. Seven of the assays demonstrated within-haplotype diversity; five failed to amplify loci in at least one haplotype A sample while three assays produced amplicons from some haplotype B samples. Eight of the loci assays showed consistent A-B differentiation. Differences in genome arrangements, prophage, and qPCR results suggesting locus diversity within the haplotypes provide more evidence for genetic complexity in this emerging bacterial species.

An ancient duplication of apple MYB transcription factors is responsible for novel red fruit-flesh phenotypes.

Anthocyanin accumulation is coordinated in plants by a number of conserved transcription factors. In apple (Malus × domestica), an R2R3 MYB transcription factor has been shown to control fruit flesh and foliage anthocyanin pigmentation (MYB10) and fruit skin color (MYB1). However, the pattern of expression and allelic variation at these loci does not explain all anthocyanin-related apple phenotypes. One such example is an open-pollinated seedling of cv Sangrado that has green foliage and develops red flesh in the fruit cortex late in maturity. We used methods that combine plant breeding, molecular biology, and genomics to identify duplicated MYB transcription factors that could control this phenotype. We then demonstrated that the red-flesh cortex phenotype is associated with enhanced expression of MYB110a, a paralog of MYB10. Functional characterization of MYB110a showed that it was able to up-regulate anthocyanin biosynthesis in tobacco (Nicotiana tabacum). The chromosomal location of MYB110a is consistent with a whole-genome duplication event that occurred during the evolution of apple within the Maloideae family. Both MYB10 and MYB110a have conserved function in some cultivars, but they differ in their expression pattern and response to fruit maturity.

A toolkit for bulk PCR-based marker design from next-generation sequence data: application for development of a framework linkage map in bulb onion (Allium cepa L.).

BACKGROUND: Although modern sequencing technologies permit the ready detection of numerous DNA sequence variants in any organisms, converting such information to PCR-based genetic markers is hampered by a lack of simple, scalable tools. Onion is an example of an under-researched crop with a complex, heterozygous genome where genome-based research has previously been hindered by limited sequence resources and genetic markers. RESULTS: We report the development of generic tools for large-scale web-based PCR-based marker design in the Galaxy bioinformatics framework, and their application for development of next-generation genetics resources in a wide cross of bulb onion (Allium cepa L.). Transcriptome sequence resources were developed for the homozygous doubled-haploid bulb onion line 'CUDH2150' and the genetically distant Indian landrace 'Nasik Red', using 454™ sequencing of normalised cDNA libraries of leaf and shoot. Read mapping of 'Nasik Red' reads onto 'CUDH2150' assemblies revealed 16836 indel and SNP polymorphisms that were mined for portable PCR-based marker development. Tools for detection of restriction polymorphisms and primer set design were developed in BioPython and adapted for use in the Galaxy workflow environment, enabling large-scale and targeted assay design. Using PCR-based markers designed with these tools, a framework genetic linkage map of over 800cM spanning all chromosomes was developed in a subset of 93 F(2) progeny from a very large F(2) family developed from the 'Nasik Red' x 'CUDH2150' inter-cross. The utility of tools and genetic resources developed was tested by designing markers to transcription factor-like polymorphic sequences. Bin mapping these markers using a subset of 10 progeny confirmed the ability to place markers within 10 cM bins, enabling increased efficiency in marker assignment and targeted map refinement. The major genetic loci conditioning red bulb colour (R) and fructan content (Frc) were located on this map by QTL analysis. CONCLUSIONS: The generic tools developed for the Galaxy environment enable rapid development of sets of PCR assays targeting sequence variants identified from Illumina and 454 sequence data. They enable non-specialist users to validate and exploit large volumes of next-generation sequence data using basic equipment.

The Genomics and Population Genomics of the Light Brown Apple Moth, Epiphyas postvittana, an Invasive Tortricid Pest of Horticulture.

The light brown apple moth, Epiphyas postvittana is an invasive, polyphagous pest of horticultural systems around the world. With origins in Australia, the pest has subsequently spread to New Zealand, Hawaii, California and Europe, where it has been found on over 500 plants, including many horticultural crops. We have produced a genomic resource, to understand the biological basis of the polyphagous and invasive nature of this and other lepidopteran pests. The assembled genome sequence encompassed 598 Mb and has an N50 of 301.17 kb, with a BUSCO completion rate of 97.9%. Epiphyas postvittana has 34% of its assembled genome represented as repetitive sequences, with the majority of the known elements made up of longer DNA transposable elements (14.07 Mb) and retrotransposons (LINE 17.83 Mb). Of the 31,389 predicted genes, 28,714 (91.5%) were assigned to 11,438 orthogroups across the Lepidoptera, of which 945 were specific to E. postvittana. Twenty gene families showed significant expansions in E. postvittana, including some likely to have a role in its pest status, such as cytochrome p450s, glutathione-S-transferases and UDP-glucuronosyltransferases. Finally, using a RAD-tag approach, we investigated the population genomics of this pest, looking at its likely patterns of invasion.

Comparative analysis of rosaceous genomes and the reconstruction of a putative ancestral genome for the family.

BACKGROUND: Comparative genome mapping studies in Rosaceae have been conducted until now by aligning genetic maps within the same genus, or closely related genera and using a limited number of common markers. The growing body of genomics resources and sequence data for both Prunus and Fragaria permits detailed comparisons between these genera and the recently released Malus × domestica genome sequence. RESULTS: We generated a comparative analysis using 806 molecular markers that are anchored genetically to the Prunus and/or Fragaria reference maps, and physically to the Malus genome sequence. Markers in common for Malus and Prunus, and Malus and Fragaria, respectively were 784 and 148. The correspondence between marker positions was high and conserved syntenic blocks were identified among the three genera in the Rosaceae. We reconstructed a proposed ancestral genome for the Rosaceae. CONCLUSIONS: A genome containing nine chromosomes is the most likely candidate for the ancestral Rosaceae progenitor. The number of chromosomal translocations observed between the three genera investigated was low. However, the number of inversions identified among Malus and Prunus was much higher than any reported genome comparisons in plants, suggesting that small inversions have played an important role in the evolution of these two genera or of the Rosaceae.

The Gillenia trifoliata genome reveals dynamics correlated with growth and reproduction in Rosaceae.

The Rosaceae family has striking phenotypic diversity and high syntenic conservation. Gillenia trifoliata is sister species to the Maleae tribe of apple and ~1000 other species. Gillenia has many putative ancestral features, such as herb/sub-shrub habit, dry fruit-bearing and nine base chromosomes. This coalescence of ancestral characters in a phylogenetically important species, positions Gillenia as a 'rosetta stone' for translational science within Rosaceae. We present genomic and phenological resources to facilitate the use of Gillenia for this purpose. The Gillenia genome is the first fully annotated chromosome-level assembly with an ancestral genome complement (x = 9), and with it we developed an improved model of the Rosaceae ancestral genome. MADS and NAC gene family analyses revealed genome dynamics correlated with growth and reproduction and we demonstrate how Gillenia can be a negative control for studying fleshy fruit development in Rosaceae.

Construction of a high-density genetic map for hexaploid kiwifruit (Actinidia chinensis var. deliciosa) using genotyping by sequencing.

Commercially grown kiwifruit (genus Actinidia) are generally of two sub-species which have a base haploid genome of 29 chromosomes. The yellow-fleshed Actinidia chinensis var. chinensis, is either diploid (2n = 2x = 58) or tetraploid (2n = 4x = 116) and the green-fleshed cultivar A. chinensis var. deliciosa "Hayward," is hexaploid (2n = 6x = 174). Advances in breeding green kiwifruit could be greatly sped up by the use of molecular resources for more efficient and faster selection, for example using marker-assisted selection (MAS). The key genetic marker that has been implemented for MAS in hexaploid kiwifruit is for gender testing. The limited marker-trait association has been reported for other polyploid kiwifruit for fruit and production traits. We have constructed a high-density linkage map for hexaploid green kiwifruit using genotyping-by-sequence (GBS). The linkage map obtained consists of 3686 and 3940 markers organized in 183 and 176 linkage groups for the female and male parents, respectively. Both parental linkage maps are co-linear with the A. chinensis "Red5" reference genome of kiwifruit. The linkage map was then used for quantitative trait locus (QTL) mapping, and successfully identified QTLs for king flower number, fruit number and weight, dry matter accumulation, and storage firmness. These are the first QTLs to be reported and discovered for complex traits in hexaploid kiwifruit.

A gene-rich linkage map in the dioecious species Actinidia chinensis (kiwifruit) reveals putative X/Y sex-determining chromosomes.

BACKGROUND: The genus Actinidia (kiwifruit) consists of woody, scrambling vines, native to China, and only recently propagated as a commercial crop. All species described are dioecious, but the genetic mechanism for sex-determination is unknown, as is the genetic basis for many of the cluster of characteristics making up the unique fruit. It is, however, an important crop in the New Zealand economy, and a classical breeding program would benefit greatly by knowledge of the trait alleles carried by both female and male parents. The application of marker assisted selection (MAS) in seedling populations would also aid the accurate and efficient development of novel fruit types for the market. RESULTS: Gene-rich female, male and consensus linkage maps of the diploid species A. chinensis have been constructed with 644 microsatellite markers. The maps consist of twenty-nine linkage groups corresponding to the haploid number n = 29. We found that sex-linked sequence characterized amplified region (SCAR) markers and the 'Flower-sex' phenotype consistently mapped to a single linkage group, in a subtelomeric region, in a section of inconsistent marker order. The region also contained markers of expressed genes, some of unknown function. Recombination, assessed by allelic distribution and marker order stability, was, in the remainder of the linkage group, in accordance with other linkage groups. Fully informative markers to other genes in this linkage group identified the comparative linkage group in the female map, where recombination ratios determining marker order were similar to the autosomes. CONCLUSION: We have created genetic linkage maps that define the 29 linkage groups of the haploid genome, and have revealed the position and extent of the sex-determining locus in A. chinensis. As all Actinidia species are dioecious, we suggest that the sex-determining loci of other Actinidia species will be similar to that region defined in our maps. As the extent of the non-recombining region is limited, our result supports the suggestion that the subtelomeric region of an autosome is in the early stages of developing the characteristics of a sex chromosome. The maps provide a reference of genetic information in Actinidia for use in genetic analysis and breeding programs.

Development of a set of SNP markers present in expressed genes of the apple.

Molecular markers associated with gene coding regions are useful tools for bridging functional and structural genomics. Due to their high abundance in plant genomes, single nucleotide polymorphisms (SNPs) are present within virtually all genomic regions, including most coding sequences. The objective of this study was to develop a set of SNPs for the apple by taking advantage of the wealth of genomics resources available for the apple, including a large collection of expressed sequenced tags (ESTs). Using bioinformatics tools, a search for SNPs within an EST database of approximately 350,000 sequences developed from a variety of apple accessions was conducted. This resulted in the identification of a total of 71,482 putative SNPs. As the apple genome is reported to be an ancient polyploid, attempts were made to verify whether those SNPs detected in silico were attributable either to allelic polymorphisms or to gene duplication or paralogous or homeologous sequence variations. To this end, a set of 464 PCR primer pairs was designed, PCR was amplified using two subsets of plants, and the PCR products were sequenced. The SNPs retrieved from these sequences were then mapped onto apple genetic maps, including a newly constructed map of a Royal Gala x A689-24 cross and a Malling 9 x Robusta 5, map using a bin mapping strategy. The SNP genotyping was performed using the high-resolution melting (HRM) technique. A total of 93 new markers containing 210 coding SNPs were successfully mapped. This new set of SNP markers for the apple offers new opportunities for understanding the genetic control of important horticultural traits using quantitative trait loci (QTL) or linkage disequilibrium analysis. These also serve as useful markers for aligning physical and genetic maps, and as potential transferable markers across the Rosaceae family.

Common Variants of the Plant microRNA-168a Exhibit Differing Silencing Efficacy for Human Low-Density Lipoprotein Receptor Adaptor Protein 1 (LDLRAP1).

BACKGROUND: The discovery that a plant microRNA (miRNAs) from rice (Oryza sativa miR168a) can modify post-transcriptional expression of the mammalian. Low-Density Lipoprotein Receptor Adaptor Protein 1 (LDLRAP1) gene highlights the potential for cross-kingdom miRNAmRNA interactions. OBJECTIVE: To investigate whether common variants of the conserved miR168a family have the capability for similar cross-kingdom regulatory functions, we selected sequences from three dietary plant sources: rice (Oryza sativa), tomato (Solanum lycopersicum), apple (Malus domestica) and compared their ability to regulate human LDLRAP1 expression. METHODS: Target prediction software intaRNA and RNAhybrid were used to analyze and calculate the energy and alignment score between the miR168a variants and human LDLRAP1 mRNA. An in vitro cell-based Dual-Luciferase® Reporter Assay (pmirGLO, Promega), was then used to validate the miRNA-mRNA interaction experimentally. RESULTS: Computational analyses revealed that a single nucleotide difference at position 14 (from the 5' end of the miRNA) creates a G:U wobble in the miRNA-mRNA duplex formed by tomato and apple miR168a variants. This G:U wobble had only a small effect on the free energy score (-33.8-34.7 kcal/mol). However, despite reasonable hybridization energy scores (<-20 kcal/mol) for all miR168a variants, only the rice miR168a variant lacking a G:U wobble significantly reduced LDLRAP1 transcript expression by 25.8 + 7.3% (p<0.05), as measured by relative luciferase activity. CONCLUSION: In summary, single nucleotide differences at key positions can have a marked influence on regulatory function despite similar predicted energy scores and miRNA-mRNA duplex structures.

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.