Genome-Wide Identification of SSR and SNP Markers Based on Whole-Genome Re-Sequencing of a Thailand Wild Sacred Lotus (Nelumbo nucifera)

Sequence characteristics, genetic diversity and phylogenetic analysis of the Cucurbita ficifolia (Cucurbitaceae) chloroplasts genome

BACKGROUND

Curcubita ficifolia Bouché (Cucurbitaceae) has high value as a food crop and medicinal plant, and also has horticultural value as rootstock for other melon species. China is home to many different cultivars, but the genetic diversity of these resources and the evolutionary relationships among them, as well as the differences between C. ficifolia and other Cucurbita species, remain unclear.

RESULTS

We investigated the chloroplast (cp) genomes of 160 C. ficifolia individuals from 31 populations in Yunnan, a major C. ficifolia production area in China. We found that the cp genome of C. ficifolia is ~151 kb and contains 128 genes, of which 86 are protein coding genes, 34 encode tRNA, and eight encode rRNAs. We also identified 64 SSRs, mainly AT repeats. The cp genome was found to contain a total of 204 SNP and 57 indels, and a total of 21 haplotypes were found in the 160 study individuals. The reverse repeat (IR) region of C. ficifolia contained a few differences compared with this region in the six other Cucurbita species. Sequence difference analysis demonstrated that most of the variable regions were concentrated in the single copy (SC) region. Moreover, the sequences of the coding regions were found to be more similar among species than those of the non-coding regions. The phylogenies reconstructed from the cp genomes of 61 representative species of Cucurbitaceae reflected the currently accepted classification, in which C. ficifolia is sister to the other Cucurbita species, however, different interspecific relationships were found between Cucurbita species.

CONCLUSIONS

These results will be valuable in the classification of C. ficifolia genetic resources and will contribute to our understanding of evolutionary relationships within the genus Cucurbita.

Complete Chloroplast Genomes and the Phylogenetic Analysis of Three Native Species of Paeoniaceae from the Sino-Himalayan Flora Subkingdom

Paeonia delavayi var. lutea, Paeonia delavayi var. angustiloba, and Paeonia ludlowii are Chinese endemics that belong to the Paeoniaceae family and have vital medicinal and ornamental value. It is often difficult to classify Paeoniaceae plants based on their morphological characteristics, and the limited genomic information has strongly hindered molecular evolution and phylogenetic studies of Paeoniaceae. In this study, we sequenced, assembled, and annotated the chloroplast genomes of P. delavayi var. lutea, P. delavayi var. angustiloba, and P. ludlowii. The chloroplast genomes of these strains were comparatively analyzed, and their phylogenetic relationships and divergence times were inferred. These three chloroplast genomes exhibited a typical quadripartite structure and were 152,687-152,759 bp in length. Each genome contains 126-132 genes, including 81-87 protein-coding genes, 37 transfer RNAs, and 8 ribosomal RNAs. In addition, the genomes had 61-64 SSRs, with mononucleotide repeats being the most abundant. The codon bias patterns of the three species tend to use codons ending in A/U. Six regions of high variability were identified (psbK-psbL, trnG-UCC, petN-psbM, psbC, rps8-rpl14, and ycf1) that can be used as DNA molecular markers for phylogenetic and taxonomic analysis. The Ka/Ks ratio indicates positive selection for the rps18 gene associated with self-replication. The phylogenetic analysis of 99 chloroplast genomes from Saxifragales clarified the phylogenetic relationships of Paeoniaceae and revealed that P. delavayi var. lutea, P. delavayi var. angustiloba, and P. ludlowii are monophyletic groups and sisters to P. delavayi. Divergence time estimation revealed two evolutionary divergences of Paeoniaceae species in the early Oligocene and Miocene. Afterward, they underwent rapid adaptive radiation from the Pliocene to the early Pleistocene when P. delavayi var. lutea, P. delavayi var. angustiloba, and P. ludlowii formed. The results of this study enrich the chloroplast genomic information of Paeoniaceae and reveal new insights into the phylogeny of Paeoniaceae.

Comparative Analysis of the Complete Mitochondrial Genomes of Apium graveolens and Apium leptophyllum Provide Insights into Evolution and Phylogeny Relationships

The genus Apium, belonging to the family Apiaceae, comprises roughly 20 species. Only two species, Apium graveolens and Apium leptophyllum, are available in China and are both rich in nutrients and have favorable medicinal properties. However, the lack of genomic data has severely constrained the study of genetics and evolution in Apium plants. In this study, Illumina NovaSeq 6000 and Nanopore sequencing platforms were employed to identify the mitochondrial genomes of A. graveolens and A. leptophyllum. The complete lengths of the mitochondrial genomes of A. graveolens and A. leptophyllum were 263,017 bp and 260,164 bp, respectively, and contained 39 and 36 protein-coding genes, five and six rRNA genes, and 19 and 20 tRNA genes. Consistent with most angiosperms, both A. graveolens and A. leptophyllum showed a preference for codons encoding leucine (Leu). In the mitochondrial genome of A. graveolens, 335 SSRs were detected, which is higher than the 196 SSRs found in the mitochondrial genome of A. leptophyllum. Studies have shown that the most common RNA editing type is C-to-U, but, in our study, both A. graveolens and A. leptophyllum exhibited the U-C editing type. Furthermore, the transfer of the mitochondrial genomes of A. graveolens and A. leptophyllum into the chloroplast genomes revealed homologous sequences, accounting for 8.14% and 4.89% of the mitochondrial genome, respectively. Lastly, in comparing the mitochondrial genomes of 29 species, it was found that A. graveolens, A. leptophyllum, and Daucus carota form a sister group with a support rate of 100%. Overall, this investigation furnishes extensive insights into the mitochondrial genomes of A. graveolens and A. leptophyllum, thereby enhancing comprehension of the traits and evolutionary patterns within the Apium genus. Additionally, it offers supplementary data for evolutionary and comparative genomic analyses of other species within the Apiaceae family.

Identification of QTLs and Putative Candidate Genes for Plant Architecture of Lotus Revealed by Regional Association Mapping

The lotus (Nelumbo Adans.) is one of the most economically relevant ornamental aquatic plants. Plant architecture (PA) is an important trait for lotus classification, cultivation, breeding, and applications. However, the underlying genetic and molecular basis controlling PA remains poorly understood. In this study, an association study for PA-related traits was performed with 93 genome-wide microsatellite markers (simple sequence repeat, SSR) and 51 insertion-deletion (InDel) markers derived from the candidate regions using a panel of 293 lotus accessions. Phenotypic data analysis of the five PA-related traits revealed a wide normal distribution and high heritability from 2013 to 2016, which indicated that lotus PA-related traits are highly polygenic traits. The population structure (Q-matrix) and the relative kinships (K-matrix) of the association panels were analyzed using 93 SSR markers. The mixed linear model (MLM) taking Q-matrix and K-matrix into account was used to estimate the association between markers and the traits. A total of 26 markers and 65 marker-trait associations were identified by considering associations with p < 0.001 and Q < 0.05. Based on the significant markers, two QTLs on Chromosome 1 were identified, and two candidate genes were preliminarily determined. The results of our study provided useful information for the lotus breeding aiming at different PA phenotypes using a molecular-assisted selection (MAS) method and also laid the foundation for the illustration of the molecular mechanism underlying the major QTL and key markers associated with lotus PA.

Studies on Lotus Genomics and the Contribution to Its Breeding

Lotus (Nelumbo nucifera), under the Nelumbonaceae family, is one of the relict plants possessing important scientific research and economic values. Because of this, much attention has been paid to this species on both its biology and breeding among the scientific community. In the last decade, the genome of lotus has been sequenced, and several high-quality genome assemblies are available, which have significantly facilitated functional genomics studies in lotus. Meanwhile, re-sequencing of the natural and genetic populations along with different levels of omics studies have not only helped to classify the germplasm resources but also to identify the domestication of selected regions and genes controlling different horticultural traits. This review summarizes the latest progress of all these studies on lotus and discusses their potential application in lotus breeding.

Comprehensive analysis of chloroplast genome of Albizia julibrissin Durazz. (Leguminosae sp.)

The Albizia julibrissin chloroplasts have a classical chloroplast genome structure, containing 93 coding genes and 34 non-coding genes. Our research provides basic data for plant phylogenetic evolutionary studies. There is limited genomic information available for the important Chinese herb Albizia julibrissin Durazz. In this study, we constructed the chloroplast (Cp) genome of A. julibrissin. The length of the assembled Cp genome was 175,922 bp consisting of four conserved regions: a 5145 bp small single-copy (SSC) region, a 91,323 bp large single-copy (LSC) region, and two identical length-inverted repeat (IR) regions (39,725 bp). This Cp genome included 34 non-coding RNAs and 93 unique genes, the former contains 30 transfer and 4 ribosomal RNA genes. Gene annotation indicated some of the coding genes (82) in the A. julibrissin Cp genome classified in the Leguminosae family, with some to other related families (11). The results show that low GC content (36.9%) and codon bias towards A- or T-terminal codons may affect the frequency of gene codon usage. The sequence analysis identified 30 forward, 18 palindrome, and 1 reverse repeat > 30 bp length, and 149 simple sequence repeats (SSR). Fifty-five RNA editing sites in the Cp of A. julibrissin were predicted, most of which are C-to-U conversions. Analysis of the reverse repeat expansion or contraction and divergence area between several species, including A. julibrissin, was performed. The phylogenetic tree revealed that A. julibrissin was most closely related to Albizia odoratissima and Albizia bracteata, followed by Samanea saman, forming an evolutionary branch with Mimosa pudica and Leucaena trichandra. The research results are helpful for breeding and genetic improvement of A. julibrissin, and also provide valuable information for understanding the evolution of this plant.

Characterization of Genomic Variation from Lotus (Nelumbo Adans.) Mutants with Wide and Narrow Tepals

Compared with rose, chrysanthemum, and water lily, the absence of short-wide and long-narrow tepals of ornamental lotus (Nelumbo Adans.) limits the commercial value of flowers. In this study, the genomes of two groups of lotus mutants with wide-short and narrow-long tepals were resequenced to uncover the genomic variation and candidate genes associated with tepal shape. In group NL (short for N. lutea, containing two mutants and one control of N. lutea), 716,656 single nucleotide polymorphisms (SNPs) and 221,688 insertion-deletion mutations (Indels) were obtained, while 639,953 SNPs and 134,6118 Indels were obtained in group WSH (short for ‘Weishan Hong’, containing one mutant and two controls of N. nucifera ‘Weishan Hong’). Only a small proportion of these SNPs and Indels was mapped to exonic regions of genome: 1.92% and 0.47%, respectively, in the NL group, and 1.66% and 0.48%, respectively, in the WSH group. Gene Ontology (GO) analysis showed that out of 4890 (NL group) and 1272 (WSH group) annotated variant genes, 125 and 62 genes were enriched (Q < 0.05), respectively. Additionally, in the Kyoto Encyclopedia of Genes and Genomes (KEGG) database, 104 genes (NL group) and 35 genes (WSH group) were selected (p < 0.05). Finally, there were 306 candidate genes that were sieved to determine the development of tepal shape in lotus plants. It will be an essential reference for future identification of tepal-shaped control genes in lotus plants. This is the first comprehensive report of genomic variation controlling tepal shape in lotus, and the mutants in this study are promising materials for breeding novel lotus cultivars with special tepals.

Genome-wide analysis of AP2/ERF superfamily in lotus (Nelumbo nucifera) and the association between NnADAP and rhizome morphology

Background

The AP2/ERF family is widely present in plants and plays a crucial regulatory role in plant growth and development. As an essential aquatic horticultural model plant, lotus has an increasingly prominent economic and research value.

Results

We have identified and analysed the AP2/ERF gene family in the lotus. Initially, 121 AP2/ERF family genes were identified. By analysing their gene distribution and protein structure, and their expression patterns during the development of lotus rhizome, combined with previous studies, we obtained an SNP (megascaffold_20:3578539) associated with lotus rhizome phenotype. This SNP was in the NnADAP gene of the AP2 subfamily, and the changes in SNP (C/T) caused amino acid conversion (proline/leucine). We constructed a population of 95 lotus varieties for SNP verification. Through population typing experiments, we found that the group with SNP CC had significantly larger lotus rhizome and higher soluble sugar content among the population.

Conclusions

In conclusion, we speculate that the alteration of the SNP in the NnADAP can affect the size and sugar content of the lotus rhizome.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07473-w.

The Latest Studies on Lotus (Nelumbo nucifera)-an Emerging Horticultural Model Plant

Lotus (Nelumbo nucifera) is a perennial aquatic basal eudicot belonging to a small family Nelumbonaceace, which contains only one genus with two species. It is an important horticultural plant, with its uses ranging from ornamental, nutritional to medicinal values, and has been widely used, especially in Southeast Asia. Recently, the lotus obtained a lot of attention from the scientific community. An increasing number of research papers focusing on it have been published, which have shed light on the mysteries of this species. Here, we comprehensively reviewed the latest advancement of studies on the lotus, including phylogeny, genomics and the molecular mechanisms underlying its unique properties, its economic important traits, and so on. Meanwhile, current limitations in the research of the lotus were addressed, and the potential prospective were proposed as well. We believe that the lotus will be an important model plant in horticulture with the generation of germplasm suitable for laboratory operation and the establishment of a regeneration and transformation system.

Origin and parental genome characterization of the allotetraploid Stylosanthes scabra Vogel (Papilionoideae, Leguminosae), an important legume pasture crop

BACKGROUNDS AND AIMS

The genus Stylosanthes includes nitrogen-fixing and drought-tolerant species of considerable economic importance for perennial pasture, green manure and land recovery. Stylosanthes scabra is adapted to variable soil conditions, being cultivated to improve pastures and soils worldwide. Previous studies have proposed S. scabra as an allotetraploid species (2n = 40) with a putative diploid A genome progenitor S. hamata or S. seabrana (2n = 20) and the B genome progenitor S. viscosa (2n = 20). We aimed to provide conclusive evidence for the origin of S. scabra.

METHODS

We performed fluorescence in situ hybridization (FISH) and genomic in situ hybridization (GISH) experiments and Illumina paired-end sequencing of S. scabra, S. hamata and S. viscosa genomic DNA, to assemble and compare complete ribosomal DNA (rDNA) units and chloroplast genomes. Plastome- and genome-wide single nucleotide variation detection was also performed.

KEY RESULTS

GISH and phylogenetic analyses of plastid DNA and rDNA sequences support that S. scabra is an allotetraploid formed from 0.63 to 0.52 million years ago (Mya), from progenitors with a similar genome structure to the maternal donor S. hamata and the paternal donor S. viscosa. FISH revealed a non-additive number of 35S rDNA sites in S. scabra compared with its progenitors, indicating the loss of one locus from A genome origin. In S. scabra, most 5S rDNA units were similar to S. viscosa, while one 5S rDNA site of reduced size most probably came from an A genome species as revealed by GISH and in silico analysis.

CONCLUSIONS

Our approach combined whole-plastome and rDNA assembly with additional cytogenetic analysis to shed light successfully on the allotetraploid origin of S. scabra. We propose a Middle Pleistocene origin for S. scabra involving species with maternal A and paternal B genomes. Our data also suggest that variation found in rDNA units in S. scabra and its progenitors reveals differences that can be explained by homogenization, deletion and amplification processes that have occurred since its origin.

Detection of Highly Differentiated Genomic Regions Between Lotus (Nelumbo nucifera Gaertn.) With Contrasting Plant Architecture and Their Functional Relevance to Plant Architecture

The lotus (Nelumbo nucifera Gaertn.) is one of the most economically and ornamentally important perennial aquatic plants. Plant architecture is an important trait for lotus classification, cultivation, breeding, and applications. In this study, traits representing plant architecture were measured in 390 lotus germplasms for 3 years. According to the phenotypic distribution, 21 large architecture (LA) and 22 small architecture (SA) germplasms exhibiting extreme phenotypes were selected as representatives of plant architecture. Microscopy analyses revealed that LA lotuses possessed far more vertical cells and longer cell lengths than SA lotuses, and there was a closer linear relationship between vertical cell number and plant architecture than cell length and plant architecture. Furthermore, based on whole genome re-sequencing data from 10 LA and 10 SA lotus germplasms, fixation index (FST) genome scan identified 11.02 Mb of genomic regions that were highly differentiated between the LA and SA lotus groups. Chi-square test revealed that 17,154 single nucleotide polymorphisms (SNPs) and 1,554 insertions and deletions (InDels) showed distinct allelic distribution between the LA and SA lotus groups within these regions. A total of 126 variants with distinct allelic distribution in the highly differentiated region were predicted to cause amino acid changes in 60 genes. Among the 41 genes with functional annotation, the expression patterns of six genes involved in cell division and cell wall construction were confirmed using quantitative reverse-transcription PCR (qRT-PCR). In addition, 34 plant architecture-associated InDel markers were developed and verified in the remaining 11 LA and 12 SA lotus plant architecture representatives. This study identified promising functional markers and candidates for molecular breeding and will facilitate further elucidation of the genetic mechanisms underlying plant architecture in the lotus.

Improving Nelumbo nucifera genome assemblies using high-resolution genetic maps and BioNano genome mapping reveals ancient chromosome rearrangements

Genetic and physical maps are powerful tools to anchor fragmented draft genome assemblies generated from next-generation sequencing. Currently, two draft assemblies of Nelumbo nucifera, the genomes of 'China Antique' and 'Chinese Tai-zi', have been released. However, there is presently no information on how the sequences are assembled into chromosomes in N. nucifera. The lack of physical maps and inadequate resolution of available genetic maps hindered the assembly of N. nucifera chromosomes. Here, a linkage map of N. nucifera containing 2371 bin markers [217 577 single nucleotide polymorphisms (SNPs)] was constructed using restriction-site associated DNA sequencing data of 181 F₂ individuals and validated by adding 197 simple sequence repeat (SSR) markers. Additionally, a BioNano optical map covering 86.20% of the 'Chinese Tai-zi' genome was constructed. The draft assembly of 'Chinese Tai-zi' was improved based on the BioNano optical map, showing an increase of the scaffold N50 from 0.989 to 1.48 Mb. Using a combination of multiple maps, 97.9% of the scaffolds in the 'Chinese Tai-zi' draft assembly and 97.6% of the scaffolds in the 'China Antique' draft assembly were anchored into pseudo-chromosomes, and the centromere regions along the pseudo-chromosomes were identified. An evolutionary scenario was proposed to reach the modern N. nucifera karyotype from the seven ancestral eudicot chromosomes. The present study provides the highest-resolution linkage map, the optical map and chromosome level genome assemblies for N. nucifera, which are valuable for the breeding and cultivation of N. nucifera and future studies of comparative and evolutionary genomics in angiosperms.

Whole genome re-sequencing reveals evolutionary patterns of sacred lotus (Nelumbo nucifera)

Sacred lotus (Nelumbo nucifera or lotus) is an important aquatic plant in horticulture and ecosystems. As a foundation for exploring genomic variation and evolution among different germplasms, we re-sequenced 19 individuals from three cultivated temperate lotus subgroups (rhizome, seed and flower lotus), one wild temperate lotus subgroup (wild lotus), one tropical lotus group (Thai lotus) and an outgroup (Nelumbo lutea). Through genetic diversity and polymorphism analysis by non-missing SNP sites widely distributed in the whole genome, we confirmed that wild and Thai lotus exhibited greater differentiation with a higher genomic diversity compared to cultivated lotus. Rhizome lotus had the lowest genomic diversity and a closer relationship to wild lotus, whereas the genomes of seed and flower lotus were admixed. Genes in energy metabolism process and plant immunity evolved rapidly in lotus, reflecting local adaptation. We established that candidate genes in genomic regions with significant differentiation associated with temperate and tropical lotus divergence always exhibited highly divergent expression pattern. Together, this study comprehensive and credible interpretates important patterns of genetic diversity and relationships, gene evolution, and genomic signature from ecotypic differentiation of sacred lotus.

The Complete Chloroplast Genome Sequences of the Medicinal Plant Forsythia suspensa (Oleaceae)

Forsythia suspensa is an important medicinal plant and traditionally applied for the treatment of inflammation, pyrexia, gonorrhea, diabetes, and so on. However, there is limited sequence and genomic information available for F. suspensa. Here, we produced the complete chloroplast genomes of F. suspensa using Illumina sequencing technology. F. suspensa is the first sequenced member within the genus Forsythia (Oleaceae). The gene order and organization of the chloroplast genome of F. suspensa are similar to other Oleaceae chloroplast genomes. The F. suspensa chloroplast genome is 156,404 bp in length, exhibits a conserved quadripartite structure with a large single-copy (LSC; 87,159 bp) region, and a small single-copy (SSC; 17,811 bp) region interspersed between inverted repeat (IRa/b; 25,717 bp) regions. A total of 114 unique genes were annotated, including 80 protein-coding genes, 30 tRNA, and four rRNA. The low GC content (37.8%) and codon usage bias for A- or T-ending codons may largely affect gene codon usage. Sequence analysis identified a total of 26 forward repeats, 23 palindrome repeats with lengths >30 bp (identity > 90%), and 54 simple sequence repeats (SSRs) with an average rate of 0.35 SSRs/kb. We predicted 52 RNA editing sites in the chloroplast of F. suspensa, all for C-to-U transitions. IR expansion or contraction and the divergent regions were analyzed among several species including the reported F. suspensa in this study. Phylogenetic analysis based on whole-plastome revealed that F. suspensa, as a member of the Oleaceae family, diverged relatively early from Lamiales. This study will contribute to strengthening medicinal resource conservation, molecular phylogenetic, and genetic engineering research investigations of this species.

The evolution of plant microRNAs: insights from a basal eudicot sacred lotus

microRNAs (miRNAs) are important noncoding small RNAs that regulate mRNAs in eukaryotes. However, under which circumstances different miRNAs/miRNA families exhibit different evolutionary trajectories in plants remains unclear. In this study, we sequenced the small RNAs and degradome from a basal eudicot, sacred lotus (Nelumbo nucifera or lotus), to identify miRNAs and their targets. Combining with public miRNAs, we predicted 57 pre-eudicot miRNA families from different evolutionary stages. We found that miRNA families featuring older age, higher copy and target number tend to show lower propensity for miRNA family loss (PGL) and stronger signature of purifying selection during divergence of temperate and tropical lotus. Further analyses of lotus genome revealed that there is an association between loss of miRNA families in descendent plants and in duplicated genomes. Gene dosage balance is crucial in maintaining those preferentially retained MIRNA duplicates by imposing stronger purifying selection. However, these factors and selection influencing miRNA family evolution are not applicable to the putative MIRNA-likes. Additionally, the MIRNAs participating in lotus pollen-pistil interaction, a conserved process in angiosperms, also have a strong signature of purifying selection. Functionally, sequence divergence in MIRNAs escalates expression divergence of their target genes between temperate and tropical lotus during rhizome and leaf growth. Overall, our study unravels several important factors and selection that determine the miRNA family distribution in plants and duplicated genomes, and provides evidence for functional impact of MIRNA sequence evolution.

The Complete Chloroplast Genome Sequences of the Medicinal Plant Pogostemon cablin

Pogostemon cablin, the natural source of patchouli alcohol, is an important herb in the Lamiaceae family. Here, we present the entire chloroplast genome of P. cablin. This genome, with 38.24% GC content, is 152,460 bp in length. The genome presents a typical quadripartite structure with two inverted repeats (each 25,417 bp in length), separated by one small and one large single-copy region (17,652 and 83,974 bp in length, respectively). The chloroplast genome encodes 127 genes, of which 107 genes are single-copy, including 79 protein-coding genes, four rRNA genes, and 24 tRNA genes. The genome structure, GC content, and codon usage of this chloroplast genome are similar to those of other species in the family, except that it encodes less protein-coding genes and tRNA genes. Phylogenetic analysis reveals that P. cablin diverged from the Scutellarioideae clade about 29.45 million years ago (Mya). Furthermore, most of the simple sequence repeats (SSRs) are short polyadenine or polythymine repeats that contribute to high AT content in the chloroplast genome. Complete sequences and annotation of P. cablin chloroplast genome will facilitate phylogenic, population and genetic engineering research investigations involving this particular species.