Unraveling the Developmental and Genetic Mechanisms Underpinning Floral Architecture in Proteaceae

Genome-wide characterization leading to simple sequence repeat (SSR) markers development in Shorea robusta

Tropical rainforests in Southeast Asia are enriched by multifarious biota dominated by Dipterocarpaceae. In this family, Shorea robusta is an ecologically sensitive and economically important timber species whose genomic diversity and phylogeny remain understudied due to lack of datasets on genetic resources. Smattering availability of molecular markers impedes population genetic studies indicating a necessity to develop genomic databases and species-specific markers in S. robusta. Accordingly, the present study focused on fostering de novo low-depth genome sequencing, identification of reliable microsatellites markers, and their validation in various populations of S. robusta in Uttarakhand Himalayas. With 69.88 million raw reads assembled into 1,97,489 contigs (read mapped to 93.2%) and a genome size of 357.11 Mb (29 × coverage), Illumina paired-end sequencing technology arranged a library of sequence data of ~ 10 gigabases (Gb). From 57,702 microsatellite repeats, a total of 35,049 simple sequence repeat (SSR) primer pairs were developed. Afterward, among randomly selected 60 primer pairs, 50 showed successful amplification and 24 were found as polymorphic. Out of which, nine polymorphic loci were further used for genetic analysis in 16 genotypes each from three different geographical locations of Uttarakhand (India). Prominently, the average number of alleles per locus (Na), observed heterozygosity (Ho), expected heterozygosity (He), and the polymorphism information content (PIC) were recorded as 2.44, 0.324, 0.277 and 0.252, respectively. The accessibility of sequence information and novel SSR markers potentially enriches the current knowledge of the genomic background for S. robusta and to be utilized in various genetic studies in species under tribe Shoreae.

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.

Genome skimming-based simple sequence repeat (SSR) marker discovery and characterization in Grevillea robusta

Proteaceae, a largely southern hemisphere family consisting of 80 genera distributed in Australia and southern Africa as its centres of greatest diversity, also extends well in northern and southern America. Under this family, Grevillea robusta is a fast-growing species got popularity in farm and avenue plantations. Despite the ecological and economic importance, the species has not yet been investigated for its genetic improvement and genome-based studies. Only a few molecular markers are available for the species or its close relatives, which hinders  genomic and population genetics studies. Genetic markers have been intensively applied for the main strategies in breeding programs, especially for the economically important traits. Hence, it is of utmost priority to develop genomic database resources and species-specific markers for studying quantitative genetics in G. robusta. Given this, the present study aimed to develop de novo genome sequencing, robust microsatellites markers, sequence annotation and their validation in different stands of G. robusta in northern India. Library preparation and sequencing were carried out using Illumina paired-end sequencing technology. Approximately, ten gigabases (Gb) sequence data with 70.87 million raw reads assembled into 425,923 contigs (read mapped to 76.48%) comprising 455 Mb genome size (23 × coverage) generated through genome skimming approach. In total, 9421 simple sequence repeat (SSR) primer pairs were successfully designed from 13,335 microsatellite repeats. Afterward, a subset of 161 primer pairs was randomly selected, synthesized and validated. All the tested primers showed successful amplification but only 13 showed polymorphisms. The polymorphic SSRs were further used to estimate the measures of genetic diversity in 12 genotypes each from the states of Punjab, Haryana, Himachal Pradesh and Uttarakhand. Importantly, the average number of alleles (Na), observed heterozygosity (Ho), expected heterozygosity (He), and the polymorphism information content (PIC) were recorded as 2.69, 0.356, 0.557 and 0.388, respectively. The availability of sequence information and newly developed SSR markers could potentially be used in various genetic analyses and improvements through molecular breeding strategies for G. robusta.

SUPPLEMENTARY INFORMATION

The online version contains supplementary material available at 10.1007/s12298-021-01035-w.

Transcriptome Analysis Reveals Putative Target Genes of APETALA3-3 During Early Floral Development in Nigella damascena L

Even though petals are homoplastic structures, their identity consistently involves genes of the APETALA3 (AP3) lineage. However, the extent to which the networks downstream of AP3 are conserved in species with petals of different evolutionary origins is unknown. In Ranunculaceae, the specificity of the AP3-III lineage offers a great opportunity to identify the petal gene regulatory network in a comparative framework. Using a transcriptomic approach, we investigated putative target genes of the AP3-III ortholog NdAP3-3 in Nigella damascena at early developmental stages when petal identity is determined, and we compared our data with that from selected eudicot species. We generated a de novo reference transcriptome to carry out a differential gene expression analysis between the wild-type and mutant NdAP3-3 genotypes differing by the presence vs. absence of petals at early stages of floral development. Among the 1,620 genes that were significantly differentially expressed between the two genotypes, functional annotation suggested a large involvement of nuclear activities, including regulation of transcription, and enrichment in processes linked to cell proliferation. Comparing with Arabidopsis data, we found that highly conserved genes between the two species are enriched in homologs of direct targets of the AtAP3 protein. Integrating AP3-3 binding site data from another Ranunculaceae species, Aquilegia coerulea, allowed us to identify a set of 18 putative target genes that were conserved between the three species. Our results suggest that, despite the independent evolutionary origin of petals in core eudicots and Ranunculaceae, a small conserved set of genes determines petal identity and early development in these taxa.

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

Macadamia integrifolia is a representative of the large basal eudicot family Proteaceae and the main progenitor species of the Australian native nut crop macadamia. Since its commercialisation in Hawaii fewer than 100 years ago, global production has expanded rapidly. However, genomic resources are limited in comparison to other horticultural crops. The first draft assembly of M. integrifolia had good coverage of the functional gene space but its high fragmentation has restricted its use in comparative genomics and association studies. Here we have generated an improved assembly of cultivar HAES 741 (4,094 scaffolds, 745 Mb, N50 413 kb) using a combination of Illumina paired and PacBio long read sequences. Scaffolds were anchored to 14 pseudo-chromosomes using seven genetic linkage maps. This assembly has improved contiguity and coverage, with >120 Gb of additional sequence. Following annotation, 34,274 protein-coding genes were predicted, representing 90% of the expected gene content. Our results indicate that the macadamia genome is repetitive and heterozygous. The total repeat content was 55% and genome-wide heterozygosity, estimated by read mapping, was 0.98% or an average of one SNP per 102 bp. This is the first chromosome-scale genome assembly for macadamia and the Proteaceae. It is expected to be a valuable resource for breeding, gene discovery, conservation and evolutionary genomics.

Evolution of Class II TCP genes in perianth bearing Piperales and their contribution to the bilateral calyx in Aristolochia

Controlled spatiotemporal cell division and expansion are responsible for floral bilateral symmetry. Genetic studies have pointed to class II TCP genes as major regulators of cell division and floral patterning in model core eudicots. Here we study their evolution in perianth-bearing Piperales and their expression in Aristolochia, a rare occurrence of bilateral perianth outside eudicots and monocots. The evolution of class II TCP genes reveals single-copy CYCLOIDEA-like genes and three paralogs of CINCINNATA (CIN) in early diverging angiosperms. All class II TCP genes have independently duplicated in Aristolochia subgenus Siphisia. Also CIN2 genes duplicated before the diversification of Saruma and Asarum. Sequence analysis shows that CIN1 and CIN3 share motifs with Cyclin proteins and CIN2 genes have lost the miRNA319a binding site. Expression analyses of all paralogs of class II TCP genes in Aristolochia fimbriata point to a role of CYC and CIN genes in maintaining differential perianth expansion during mid- and late flower developmental stages by promoting cell division in the distal and ventral portion of the limb. It is likely that class II TCP genes also contribute to cell division in the leaf, the gynoecium and the ovules in A. fimbriata.