Chromosome-Scale Assembly and Annotation of the Macadamia Genome (Macadamia integrifolia HAES 741)

Comparative Transcriptome Analysis Reveals Key Functions of MiMYB Gene Family in Macadamia Nut Pericarp Formation

Macadamia nuts are one of the most important economic food items in the world. Pericarp thickness and flavonoid composition are the key quality traits of Macadamia nuts, but the underlying mechanism of pericarp formation is still unknown. In this study, three varieties with significantly different pericarp thicknesses, namely, A38, Guire No.1, and HAES 900, at the same stage of maturity, were used for transcriptome analysis, and the results showed that there were significant differences in their gene expression profile. A total of 3837 new genes were discovered, of which 1532 were functionally annotated. The GO, COG, and KEGG analysis showed that the main categories in which there were significant differences were flavonoid biosynthesis, phenylpropanoid biosynthesis, and the cutin, suberine, and wax biosynthesis pathways. Furthermore, 63 MiMYB transcription factors were identified, and 56 R2R3-MYB transcription factors were clustered into different subgroups compared with those in Arabidopsis R2R3-MYB. Among them, the S4, S6, and S7 subgroups were involved in flavonoid biosynthesis and pericarp formation. A total of 14 MiMYBs' gene expression were verified by RT-qPCR analysis. These results provide fundamental knowledge of the pericarp formation regulatory mechanism in macadamia nuts.

Effects of Substrate Composition on the Growth Traits of Grafted Seedling in Macadamia (Macadamia integrifolia) Nuts

Macadamia nut plantings in China are expanding year by year. In order to breed and promote superior varieties, this study analyzed the effects of different rootstocks and scions on the survival rate of grafted seedlings, and then selected the best substrate composition for plant growth. The results showed that the survival rate of the HAES788 variety as rootstock and Guire No. 1 as scion was the highest, reaching 96%. The optimal grafting time in December was better than that in March. Furthermore, among 16 substrate formulations, T12, T13, T15, and T16 had advantages of agglomerated soil and more well-developed root systems compared to the CK made of loess. The plant height, stem diameter, leaf length, leaf width, and dry weight of the aboveground and underground parts of the grafted seedlings planted in these substrate formulations were significantly higher than those plants planted in the CK. In addition, the substrate formulations T12, T13, T15, and T16 significantly improved the organic matter, total nitrogen, and total potassium content of the substrate soils, but little improvement was observed for total phosphorus content after 13 months. Overall, macadamia grafting times are best in December, with HAES788 and Guire No. 1 being the best rootstock and scion. The optimal substrate formulations are T12, T13, T15, and T16. This study provides a solid foundation for the production of high-quality macadamia plants.

Genome-Wide Characterization of PEBP Gene Family and Functional Analysis of TERMINAL FLOWER 1 Homologs in Macadamia integrifolia

Edible Macadamia is one of the most important commercial nut trees cultivated in many countries, but its large tree size and long juvenile period pose barriers to commercial cultivation. The short domestication period and well-annotated genome of Macadamia integrifolia create great opportunities to breed commercial varieties with superior traits. Recent studies have shown that members of the phosphatidylethanolamine binding protein (PEBP) family play pivotal roles in regulating plant architecture and flowering time in various plants. In this study, thirteen members of MiPEBP were identified in the genome of M. integrifolia, and they are highly similarity in both motif and gene structure. A phylogenetic analysis divided the MiPEBP genes into three subfamilies: MFT-like, FT-like and TFL1-like. We subsequently identified two TERMINAL FLOWER 1 homologues from the TFL1-like subfamily, MiTFL1 and MiTFL1-like, both of which were highly expressed in stems and vegetative shoots, while MiTFL1-like was highly expressed in young leaves and early flowers. A subcellular location analysis revealed that both MiTFL1 and MiTFL1-like are localized in the cytoplasm and nucleus. The ectopic expression of MiTFL1 can rescue the early-flowering and terminal-flower phenotypes in the tfl1-14 mutant of Arabidopsis thaliana, and it indicates the conserved functions in controlling the inflorescence architecture and flowering time. This study will provide insight into the isolation of PEBP family members and the key targets for breeding M. integrifolia with improved traits in plant architecture and flowering time.

Macadamia germplasm and genomic database (MacadamiaGGD): A comprehensive platform for germplasm innovation and functional genomics in Macadamia

As an important nut crop species, macadamia continues to gain increased amounts of attention worldwide. Nevertheless, with the vast increase in macadamia omic data, it is becoming difficult for researchers to effectively process and utilize the information. In this work, we developed the first integrated germplasm and genomic database for macadamia (MacadamiaGGD), which includes five genomes of four species; three chloroplast and mitochondrial genomes; genome annotations; transcriptomic data for three macadamia varieties, germplasm data for four species and 262 main varieties; nine genetic linkage maps; and 35 single-nucleotide polymorphisms (SNPs). The database serves as a valuable collection of simple sequence repeat (SSR) markers, including both markers that are based on macadamia genomic sequences and developed in this study and markers developed previously. MacadamiaGGD is also integrated with multiple bioinformatic tools, such as search, JBrowse, BLAST, primer designer, sequence fetch, enrichment analysis, multiple sequence alignment, genome alignment, and gene homology annotation, which allows users to conveniently analyze their data of interest. MacadamiaGGD is freely available online (http://MacadamiaGGD.net). We believe that the database and additional information of the SSR markers can help scientists better understand the genomic sequence information of macadamia and further facilitate molecular breeding efforts of this species.

Paleopolyploidies and Genomic Fractionation in Major Eudicot Clades

Eudicots account for ~75% of living angiosperms, containing important food and energy crops. Recently, high-quality genome sequences of several eudicots including Aquilegia coerulea and Nelumbo nucifera have become available, providing an opportunity to investigate the early evolutionary characteristics of eudicots. We performed genomic hierarchical and event-related alignments to infer homology within and between representative species of eudicots. The results provide strong evidence for multiple independent polyploidization events during the early diversification of eudicots, three of which are likely to be allopolyploids: The core eudicot-common hexaploidy (ECH), Nelumbo-specific tetraploidy (NST), and Ranunculales-common tetraploidy (RCT). Using different genomes as references, we constructed genomic alignment to list the orthologous and paralogous genes produced by polyploidization and speciation. This could provide a fundamental framework for studying other eudicot genomes and gene(s) evolution. Further, we revealed significantly divergent evolutionary rates among these species. By performing evolutionary rate correction, we dated RCT to be ~118-134 million years ago (Mya), after Ranunculales diverged with core eudicots at ~123-139 Mya. Moreover, we characterized genomic fractionation resulting from gene loss and retention after polyploidizations. Notably, we revealed a high degree of divergence between subgenomes. In particular, synonymous nucleotide substitutions at synonymous sites (Ks) and phylogenomic analyses implied that A. coerulea might provide the subgenome(s) for the gamma-hexaploid hybridization.

Horticultural innovation by viral-induced gene regulation of carotenogenesis

Multipartite viral vectors provide a simple, inexpensive and effective biotechnological tool to transiently manipulate (i.e. reduce or increase) gene expression in planta and characterise the function of genetic traits. The development of virus-induced gene regulation (VIGR) systems usually involve the targeted silencing or overexpression of genes involved in pigment biosynthesis or degradation in plastids, thereby providing rapid visual assessment of success in establishing RNA- or DNA-based VIGR systems in planta. Carotenoids pigments provide plant tissues with an array of yellow, orange, and pinkish-red colours. VIGR-induced transient manipulation of carotenoid-related gene expression has advanced our understanding of carotenoid biosynthesis, regulation, accumulation and degradation, as well as plastid signalling processes. In this review, we describe mechanisms of VIGR, the importance of carotenoids as visual markers of technology development, and knowledge gained through manipulating carotenogenesis in model plants as well as horticultural crops not always amenable to transgenic approaches. We outline how VIGR can be utilised in plants to fast-track the characterisation of gene function(s), accelerate fruit tree breeding programs, edit genomes, and biofortify plant products enriched in carotenoid micronutrients for horticultural innovation.

Chromosome-level de novo genome assembly of Telopea speciosissima (New South Wales waratah) using long-reads, linked-reads and Hi-C

Telopea speciosissima, the New South Wales waratah, is an Australian endemic woody shrub in the family Proteaceae. Waratahs have great potential as a model clade to better understand processes of speciation, introgression and adaptation, and are significant from a horticultural perspective. Here, we report the first chromosome-level genome for T. speciosissima. Combining Oxford Nanopore long-reads, 10x Genomics Chromium linked-reads and Hi-C data, the assembly spans 823 Mb (scaffold N50 of 69.0 Mb) with 97.8% of Embryophyta BUSCOs "Complete". We present a new method in Diploidocus (https://github.com/slimsuite/diploidocus) for classifying, curating and QC-filtering scaffolds, which combines read depths, k-mer frequencies and BUSCO predictions. We also present a new tool, DepthSizer (https://github.com/slimsuite/depthsizer), for genome size estimation from the read depth of single-copy orthologues and estimate the genome size to be approximately 900 Mb. The largest 11 scaffolds contained 94.1% of the assembly, conforming to the expected number of chromosomes (2n = 22). Genome annotation predicted 40,158 protein-coding genes, 351 rRNAs and 728 tRNAs. We investigated CYCLOIDEA (CYC) genes, which have a role in determination of floral symmetry, and confirm the presence of two copies in the genome. Read depth analysis of 180 "Duplicated" BUSCO genes using a new tool, DepthKopy (https://github.com/slimsuite/depthkopy), suggests almost all are real duplications, increasing confidence in the annotation and highlighting a possible need to revise the BUSCO set for this lineage. The chromosome-level T. speciosissima reference genome (Tspe_v1) provides an important new genomic resource of Proteaceae to support the conservation of flora in Australia and further afield.

Signatures of selection in recently domesticated macadamia

Macadamia is a high value nut crop that is recently domesticated, ideal for testing the effect of artificial selection. Here, we sequence the genome of Hawaiian cultivar ‘Kau’ and assemble into 794 Mb in 14 pseudo-chromosomes with 37,728 genes. Genome analysis reveals a whole-genome duplication event, occurred 46.8 million years ago. Gene expansions occurred in gene families involves in fatty acid biosynthesis. Gene duplication of MADS-Box transcription factors in proanthocyanidin biosynthesis are relevant for seed coat development. Genome re-sequencing of 112 accessions reveals the origin of Hawaiian cultivars from Mount Bauple in southeast Queensland in Australia. Selective sweeps are detected in macadamia cultivars, including genes involved in fatty acid biosynthesis, seed coat development, and heat stress response. Such strong effects of artificial selection in few generations reveals the genomic basis for ‘one-step operation’ for clonal crop domestication. The knowledge gained could accelerate domestication of new crops from wild species.

Macadamia is a recently domesticated nut crop. Here, the authors report the genome assembly of Hawaiian cultivar ‘Kau’ and conduct population genomic analyses to reveal the origin of Hawaiian cultivars and the genomic basis for one-step operation for the clonal crop domestication.

The genome of the endangered Macadamia jansenii displays little diversity but represents an important genetic resource for plant breeding

Macadamia, a recently domesticated expanding nut crop in the tropical and subtropical regions of the world, is one of the most economically important genera in the diverse and widely adapted Proteaceae family. All four species of Macadamia are rare in the wild with the most recently discovered, M. jansenii, being endangered. The M. jansenii genome has been used as a model for testing sequencing methods using a wide range of long read sequencing techniques. Here, we report a chromosome level genome assembly, generated using a combination of Pacific Biosciences sequencing and Hi-C, comprising 14 pseudo-molecules, with a N50 of 52 Mb and a total genome assembly size of 758 Mb of which 56% is repetitive. Completeness assessment revealed that the assembly covered -97.1% of the conserved single copy genes. Annotation predicted 31,591 protein coding genes and allowed the characterization of genes encoding biosynthesis of cyanogenic glycosides, fatty acid metabolism, and anti-microbial proteins. Re-sequencing of seven other genotypes confirmed low diversity and low heterozygosity within this endangered species. Important morphological characteristics of this species such as small tree size and high kernel recovery suggest that M. jansenii is an important source of these commercial traits for breeding. As a member of a small group of families that are sister to the core eudicots, this high-quality genome also provides a key resource for evolutionary and comparative genomics studies.

The genome of the endangered Macadamia jansenii displays little diversity but represents an important genetic resource for plant breeding

Macadamia, a recently domesticated expanding nut crop in the tropical and subtropical regions of the world, is one of the most economically important genera in the diverse and widely adapted Proteaceae family. All four species of Macadamia are rare in the wild with the most recently discovered, M. jansenii, being endangered. The M. jansenii genome has been used as a model for testing sequencing methods using a wide range of long read sequencing techniques. Here, we report a chromosome level genome assembly, generated using a combination of Pacific Biosciences sequencing and Hi-C, comprising 14 pseudo-molecules, with a N50 of 52 Mb and a total genome assembly size of 758 Mb of which 56% is repetitive. Completeness assessment revealed that the assembly covered -97.1% of the conserved single copy genes. Annotation predicted 31,591 protein coding genes and allowed the characterization of genes encoding biosynthesis of cyanogenic glycosides, fatty acid metabolism, and anti-microbial proteins. Re-sequencing of seven other genotypes confirmed low diversity and low heterozygosity within this endangered species. Important morphological characteristics of this species such as small tree size and high kernel recovery suggest that M. jansenii is an important source of these commercial traits for breeding. As a member of a small group of families that are sister to the core eudicots, this high-quality genome also provides a key resource for evolutionary and comparative genomics studies.

Advances in genomics and genome editing for breeding next generation of fruit and nut crops

Fruit tree crops are an essential part of the food production systems and are key to achieve food and nutrition security. Genetic improvement of fruit trees by conventional breeding has been slow due to the long juvenile phase. Advancements in genomics and molecular biology have paved the way for devising novel genetic improvement tools like genome editing, which can accelerate the breeding of these perennial crops to a great extent. In this article, advancements in genomics of fruit trees covering genome sequencing, transcriptome sequencing, genome editing technologies (GET), CRISPR-Cas system based genome editing, potential applications of CRISPR-Cas9 in fruit tree crops improvement, the factors influencing the CRISPR-Cas editing efficiency and the challenges for CRISPR-Cas9 applications in fruit tree crops improvement are reviewed. Besides, base editing, a recently emerging more precise editing system, and the future perspectives of genome editing in the improvement of fruit and nut crops are covered.

An exploration of assembly strategies and quality metrics on the accuracy of the rewarewa (Knightia excelsa) genome

We used long read sequencing data generated from Knightia excelsa, a nectar-producing Proteaceae tree endemic to Aotearoa (New Zealand), to explore how sequencing data type, volume and workflows can impact final assembly accuracy and chromosome reconstruction. Establishing a high-quality genome for this species has specific cultural importance to Māori and commercial importance to honey producers in Aotearoa. Assemblies were produced by five long read assemblers using data subsampled based on read lengths, two polishing strategies and two Hi-C mapping methods. Our results from subsampling the data by read length showed that each assembler tested performed differently depending on the coverage and the read length of the data. Subsampling highlighted that input data with longer read lengths but perhaps lower coverage constructed more contiguous, kmers and gene-complete assemblies than short read length input data with higher coverage. The final genome assembly was constructed into 14 pseudochromosomes using an initial flye long read assembly, a racon/medaka/pilon combined polishing strategy, salsa2 and allhic scaffolding, juicebox curation, and Macadamia linkage map validation. We highlighted the importance of developing assembly workflows based on the volume and read length of sequencing data and established a robust set of quality metrics for generating high-quality assemblies. Scaffolding analyses highlighted that problems found in the initial assemblies could not be resolved accurately by Hi-C data and that assembly scaffolding was more successful when the underlying contig assembly was of higher accuracy. These findings provide insight into how quality assessment tools can be implemented throughout genome assembly pipelines to inform the de novo reconstruction of a high-quality genome assembly for nonmodel organisms.

Improvements in the sequencing and assembly of plant genomes

Advances in DNA sequencing have made it easier to sequence and assemble plant genomes. Here, we extend an earlier study, and compare recent methods for long read sequencing and assembly. Updated Oxford Nanopore Technology software improved assemblies. Using more accurate sequences produced by repeated sequencing of the same molecule (Pacific Biosciences HiFi) resulted in less fragmented assembly of sequencing reads. Using data for increased genome coverage resulted in longer contigs, but reduced total assembly length and improved genome completeness. The original model species, Macadamia jansenii, was also compared with three other Macadamia species, as well as avocado (Persea americana) and jojoba (Simmondsia chinensis). In these angiosperms, increasing sequence data volumes caused a linear increase in contig size, decreased assembly length and further improved already high completeness. Differences in genome size and sequence complexity influenced the success of assembly. Advances in long read sequencing technology continue to improve plant genome sequencing and assembly. However, results were improved by greater genome coverage, with the amount needed to achieve a particular level of assembly being species dependent.

Genomic selection and genetic gain for nut yield in an Australian macadamia breeding population

Background

Improving yield prediction and selection efficiency is critical for tree breeding. This is vital for macadamia trees with the time from crossing to production of new cultivars being almost a quarter of a century. Genomic selection (GS) is a useful tool in plant breeding, particularly with perennial trees, contributing to an increased rate of genetic gain and reducing the length of the breeding cycle. We investigated the potential of using GS methods to increase genetic gain and accelerate selection efficiency in the Australian macadamia breeding program with comparison to traditional breeding methods. This study evaluated the prediction accuracy of GS in a macadamia breeding population of 295 full-sib progeny from 32 families (29 parents, reciprocals combined), along with a subset of parents. Historical yield data for tree ages 5 to 8years were used in the study, along with a set of 4113 SNP markers. The traits of focus were average nut yield from tree ages 5 to 8years and yield stability, measured as the standard deviation of yield over these 4 years. GBLUP GS models were used to obtain genomic estimated breeding values for each genotype, with a five-fold cross-validation method and two techniques: prediction across related populations and prediction across unrelated populations.

Results

Narrow-sense heritability of yield and yield stability was low (h²=0.30 and 0.04, respectively). Prediction accuracy for yield was 0.57 for predictions across related populations and 0.14 when predicted across unrelated populations. Accuracy of prediction of yield stability was high (r=0.79) for predictions across related populations. Predicted genetic gain of yield using GS in related populations was 474g/year, more than double that of traditional breeding methods (226g/year), due to the halving of generation length from 8 to 4years.

Conclusions

The results of this study indicate that the incorporation of GS for yield into the Australian macadamia breeding program may accelerate genetic gain due to reduction in generation length, though the cost of genotyping appears to be a constraint at present.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07694-z.

Comparison of long-read methods for sequencing and assembly of a plant genome

BACKGROUND

Sequencing technologies have advanced to the point where it is possible to generate high-accuracy, haplotype-resolved, chromosome-scale assemblies. Several long-read sequencing technologies are available, and a growing number of algorithms have been developed to assemble the reads generated by those technologies. When starting a new genome project, it is therefore challenging to select the most cost-effective sequencing technology, as well as the most appropriate software for assembly and polishing. It is thus important to benchmark different approaches applied to the same sample.

RESULTS

Here, we report a comparison of 3 long-read sequencing technologies applied to the de novo assembly of a plant genome, Macadamia jansenii. We have generated sequencing data using Pacific Biosciences (Sequel I), Oxford Nanopore Technologies (PromethION), and BGI (single-tube Long Fragment Read) technologies for the same sample. Several assemblers were benchmarked in the assembly of Pacific Biosciences and Nanopore reads. Results obtained from combining long-read technologies or short-read and long-read technologies are also presented. The assemblies were compared for contiguity, base accuracy, and completeness, as well as sequencing costs and DNA material requirements.

CONCLUSIONS

The 3 long-read technologies produced highly contiguous and complete genome assemblies of M. jansenii. At the time of sequencing, the cost associated with each method was significantly different, but continuous improvements in technologies have resulted in greater accuracy, increased throughput, and reduced costs. We propose updating this comparison regularly with reports on significant iterations of the sequencing technologies.

Comparison of long-read methods for sequencing and assembly of a plant genome

BACKGROUND

Sequencing technologies have advanced to the point where it is possible to generate high-accuracy, haplotype-resolved, chromosome-scale assemblies. Several long-read sequencing technologies are available, and a growing number of algorithms have been developed to assemble the reads generated by those technologies. When starting a new genome project, it is therefore challenging to select the most cost-effective sequencing technology, as well as the most appropriate software for assembly and polishing. It is thus important to benchmark different approaches applied to the same sample.

RESULTS

Here, we report a comparison of 3 long-read sequencing technologies applied to the de novo assembly of a plant genome, Macadamia jansenii. We have generated sequencing data using Pacific Biosciences (Sequel I), Oxford Nanopore Technologies (PromethION), and BGI (single-tube Long Fragment Read) technologies for the same sample. Several assemblers were benchmarked in the assembly of Pacific Biosciences and Nanopore reads. Results obtained from combining long-read technologies or short-read and long-read technologies are also presented. The assemblies were compared for contiguity, base accuracy, and completeness, as well as sequencing costs and DNA material requirements.

CONCLUSIONS

The 3 long-read technologies produced highly contiguous and complete genome assemblies of M. jansenii. At the time of sequencing, the cost associated with each method was significantly different, but continuous improvements in technologies have resulted in greater accuracy, increased throughput, and reduced costs. We propose updating this comparison regularly with reports on significant iterations of the sequencing technologies.