No abstract

Rhododendron delavayi Franch. is globally famous as an ornamental plant. Its distribution in southwest China covers several different habitats and environments. However, not much research had been conducted on Rhododendron spp. at the molecular level, which hinders understanding of its evolution, speciation, and synthesis of secondary metabolites, as well as its wide adaptability to different environments. Here, we report the genome assembly and gene annotation of R. delavayi var. delavayi (the second genome sequenced in the Ericaceae), which will facilitate the study of the family. The genome assembly will have further applications in genome-assisted cultivar breeding. The final size of the assembled R. delavayi var. delavayi genome (695.09 Mb) was close to the 697.94 Mb, estimated by k-mer analysis. A total of 336.83 gigabases (Gb) of raw Illumina HiSeq 2000 reads were generated from 9 libraries (with insert sizes ranging from 170 bp to 40 kb), achieving a raw sequencing depth of ×482.6. After quality filtering, 246.06 Gb of clean reads were obtained, giving ×352.55 coverage depth. Assembly using Platanus gave a total scaffold length of 695.09 Mb, with a contig N50 of 61.8 kb and a scaffold N50 of 637.83 kb. Gene prediction resulted in the annotation of 32 938 protein-coding genes. The genome completeness was evaluated by CEGMA and BUSCO and reached 95.97% and 92.8%, respectively. The gene annotation completeness was also evaluated by CEGMA and BUSCO and reached 97.01% and 87.4%, respectively. Genome annotation revealed that 51.77% of the R. delavayi genome is composed of transposable elements, and 37.48% of long terminal repeat elements (LTRs). The de novo assembled genome of R. delavayi var. delavayi (hereinafter referred to as R. delavayi) is the second genomic resource of the family Ericaceae and will provide a valuable resource for research on future comparative genomic studies in Rhododendron species. The availability of the R. delavayi genome sequence will hopefully provide a tool for scientists to tackle open questions regarding molecular mechanisms underlying environmental interactions in the genus Rhododendron, more accurately understand the evolutionary processes and systematics of the genus, facilitate the identification of genes encoding pharmaceutically important compounds, and accelerate molecular breeding to release elite varieties.

No abstract

Although natural hybridization is common in Rhododendron, it is rare in R. cyanocarpum and is apparently not accompanied by backcrossing toward R. cyanocarpum. Hence, there is no immediate risk of genetic swamping, unless habitat disturbance increases and changes the patterns of hybridization. Our study is the first to report a plant hybrid zone dominated by F(2) hybrids. This pattern might contribute to species barrier maintenance.

Hybridization is strongly asymmetric, with P. bulleyana the only possible mother and all detected hybrids being backcrosses in this direction. Partial ecological isolation and inhibition of heterospecific pollen, and possibly complete barriers to F1 formation on P. beesiana, may be enough to make F1 formation very rare in these species. Therefore, with no F1 detected, this hybrid zone may have a finite life span as successive generations become more similar to P. bulleyana.

Summary With advanced sequencing technology, dozens of complex polyploid plant genomes have been characterized. However, for many polyploid species, their diploid ancestors are unknown or extinct, making it impossible to unravel the subgenomes and genome evolution directly. We developed a novel subgenome‐phasing algorithm, SubPhaser, specifically designed for a neoallopolyploid or a homoploid hybrid. SubPhaser first searches for the subgenome‐specific sequence (k‐mer), then assigns homoeologous chromosomes into subgenomes, and further provides tools to annotate and investigate specific sequences. SubPhaser works well on neoallopolyploids and homoploid hybrids containing subgenome‐specific sequences like wheat, but fails on autopolyploids lacking subgenome‐specific sequences like alfalfa, indicating that SubPhaser can phase neoallopolyploid/homoploid hybrids with high accuracy, sensitivity and performance. This highly accurate, highly sensitive, ancestral data free chromosome phasing algorithm, SubPhaser, offers significant application value for subgenome phasing in neoallopolyploids and homoploid hybrids, and for the subsequent exploration of genome evolution and related genetic/epigenetic mechanisms.

BackgroundNatural hybridization in plants is universal and plays an important role in evolution. Based on morphology it has been presumed that hybridization occurred in the genus Buddleja, though genetic studies confirming this assumption have not been conducted to date. The two species B. crispa and B. officinalis overlap in their distributions over a wide range in South-West China, and we aimed to provide genetic evidence for ongoing hybridization in this study.ResultsWe investigated the occurrence of hybrids between the two species at the southern-most edge of the distribution of B. crispa using five nuclear loci and pollination experiments. The genetic data suggest substantial differentiation between the two species as species-specific alleles are separated by at least 7–28 mutations. The natural hybrids found were nearly all F1s (21 of 23), but backcrosses were detected, and some individuals, morphologically indistinguishable from the parental species, showed introgression. Pollen viability test shows that the percentage of viable pollen grains was 50 ± 4 % for B. crispa, and 81 ± 2 % for B. officinalis. This difference is highly significant (t = 7.382, p < 0.0001). Hand cross-pollination experiments showed that B. crispa is not successful as pollen-parent, but B. officinalis is able to pollinate B. crispa to produce viable hybrid seed. Inter-specific seed-set is low (8 seeds per fruit, as opposed to about 65 for intra-specific pollinations), suggesting post-zygotic reproductive barriers. In addition, one of the reference populations also suggests a history of introgression at other localities.ConclusionsThe occurrence of morphologically intermediate individuals between B. crispa and B. officinalis at Xishan Mountain is unequivocally linked to hybridization and almost all examined individuals of the putative hybrids were likely F1s. Despite pollination experiments indicating higher chances for introgression into B. officinalis (hybrids only produced viable seed when crossed with B. officinalis), observed introgression was asymmetrical into B. crispa. This could be due to seeds produced by hybrids not contributing to seedlings, or other factors favoring the establishment of backcrosses towards B. crispa. However, further research will be needed to confirm these observations, as the small number of plants used for the pollination experiments could have introduced an artifact, for example if used individuals were more or less compatible than the species average, and also the small number of loci used could convey a picture of introgression that is not representative for the whole genome.Electronic supplementary materialThe online version of this article (doi:10.1186/s12870-015-0539-9) contains supplementary material, which is available to authorized users.

No abstract