Genome skimming reveals the origin of the Jerusalem Artichoke tuber crop species: neither from Jerusalem nor an artichoke

Jerusalem artichoke: A comprehensive review of nutritional composition, health benefits and emerging trends in food applications

The Jerusalem artichoke (JA), a plantrelated to sunflowers and native to North America, has long been valued for its versatility, especially during periods of food scarcity. This resilient crop serves multiple purposes, functioning as a vegetable, medicinal herb, grazing crop, and even a biofuel source. In recent years, interest in JA has grown, largely due to its high nutritional profile and associated health benefits. This review explores JA's nutritional composition, the benefits of its consumption, and its botanical and agricultural characteristics. Additionally, the various applications of JA in the food industry are discussed, including its use in dairy products, snacks, baked goods, beverages, and functional foods. This review also examines the processing techniques involved in harvesting JA, extracting its valuable components, and incorporating it into food products. Notably, JA is a rich source of fiber and minerals, and incorporating it into food products not only enhances their nutritional value but also improves fermentation processes, lowers the glycemic index, and enhances sensory properties, all while reducing production costs. However, several challenges remain in JA production. These include optimizing growing conditions, addressing high labor costs, developing suitable machinery, determining the ideal harvesting time, increasing pest resistance, identifying suitable packaging materials, and developing sustainable production strategies. These challenges require further research to fully unlock JA's potential as a valuable crop.

Assembly-free reads accurate identification (AFRAID) approach outperforms other methods of DNA barcoding in the walnut family (Juglandaceae)

DNA barcoding has been extensively used for species identification. However, species identification of mixed samples or degraded DNA is limited by current DNA barcoding methods. In this study, we use plant species in Juglandaceae to evaluate an assembly-free reads accurate identification (AFRAID) method of species identification, a novel approach for precise species identification in plants. Specifically, we determined (1) the accuracy of DNA barcoding approaches in delimiting species in Juglandaceae, (2) the minimum size of chloroplast dataset for species discrimination, and (3) minimum amount of next generation sequencing (NGS) data required for species identification. We found that species identification rates were highest when whole chloroplast genomes were used, followed by taxon-specific DNA barcodes, and then universal DNA barcodes. Species identification of 100% was achieved when chloroplast genome sequence coverage reached 20% and the original sequencing data reached 500,000 reads. AFRAID accurately identified species for all samples tested after 500,000 clean reads, with far less computing time than common approaches. These results provide a new approach to accurately identify species, overcoming limitations of traditional DNA barcodes. Our method, which uses next generation sequencing to generate partial chloroplast genomes, reveals that DNA barcode regions are not necessarily fixed, accelerating the process of species identification.

Peering through the hedge: Multiple datasets yield insights into the phylogenetic relationships and incongruences in the tribe Lilieae (Liliaceae)

The increasing use of genome-scale data has significantly facilitated phylogenetic analyses, contributing to the dissection of the underlying evolutionary mechanisms that shape phylogenetic incongruences, such as incomplete lineage sorting (ILS) and hybridization. Lilieae, a prominent member of the Liliaceae family, comprises four genera and approximately 260 species, representing 43% of all species within Liliaceae. They possess high ornamental, medicinal and edible values. Yet, no study has explored the validity of various genome-scale data in phylogenetic analyses within this tribe, nor have potential evolutionary mechanisms underlying its phylogenetic incongruences been investigated. Here, transcriptome, Angiosperms353, plastid and mitochondrial data, were collected from 50 to 93 samples of Lilieae, covering all four recognized genera. Multiple datasets were created and used for phylogenetic analyses based on concatenated and coalescent-based methods. Evolutionary rates of different datasets were calculated, and divergence times were estimated. Various approaches, including coalescence simulation, Quartet Sampling (QS), calculation of concordance factors (gCF and sCF), as well as MSCquartets and reticulate network inference, were carried out to infer the phylogenetic discordances and analyze their underlying mechanisms using a reduced 33-taxon dataset. Despite extensive phylogenetic discordances among gene trees, robust phylogenies were inferred from nuclear and plastid data compared to mitochondrial data, with lower synonymous substitution detected in mitochondrial genes than in nuclear and plastid genes. Significant ILS was detected across the phylogeny of Lilieae, with clear evidence of reticulate evolution identified. Divergence time estimation indicated that most of lineages in Lilieae diverged during a narrow time frame (ranging from 5.0 Ma to 10.0 Ma), consistent with the notion of rapid radiation evolution. Our results suggest that integrating transcriptomic and plastid data can serve as cost-effective and efficient tools for phylogenetic inference and evolutionary analysis within Lilieae, and Angiosperms353 data is also a favorable choice. Mitochondrial data are more suitable for phylogenetic analyses at higher taxonomic levels due to their stronger conservation and lower synonymous substitution rates. Significant phylogenetic incongruences detected in Lilieae were caused by both incomplete lineage sorting (ILS) and reticulate evolution, with hybridization and "ghost introgression" likely prevalent in the evolution of Lilieae species. Our findings provide new insights into the phylogeny of Lilieae, enhancing our understanding of the evolution of species in this tribe.

Characterization and phylogenetic analysis of the complete chloroplast genome of Carpesium lipskyi (Asteraceae, Inuleae)

The species of Carpesium lipskyi C.Winkl. 1998 is an important traditional Chinese medicine in China. In this study, the complete chloroplast (cp) genome of C. lipskyi was determined and analyzed. The result showed that the complete cp genome of C. lipskyi was 151,244 bp in length, consisting of a large single-copy (LSC) region of 82,908 bp, a small single-copy (SSC) region of 18,430 bp, and a pair of inverted repeats (IRs) of 24,953 bp. The overall GC content of the C. lipskyi is 37.68%. The species of C. lipskyi possessed 127 genes, including 83 protein-coding genes, 36 transfer RNA genes, and eight ribosomal RNA genes. The present study found that Inula is sister groups with the closest genetic relationship. The obtained knowledge could provide useful information for future phylogenetic, taxonomic, and evolutionary studies on Inuleae.

Haplotype-resolved chromosome-level genome of hexaploid Jerusalem artichoke provides insights into its origin, evolution, and inulin metabolism

Jerusalem artichoke (Helianthus tuberosus) is a global multifunctional crop. It has wide applications in the food, health, feed, and biofuel industries and in ecological protection; it also serves as a germplasm pool for breeding of the global oil crop common sunflower (Helianthus annuus). However, biological studies of Jerusalem artichoke have been hindered by a lack of genome sequences, and its high polyploidy and large genome size have posed challenges to genome assembly. Here, we report a 21-Gb chromosome-level assembly of the hexaploid Jerusalem artichoke genome, which comprises 17 homologous groups, each with 6 pseudochromosomes. We found multiple large-scale chromosome rearrangements between Jerusalem artichoke and common sunflower, and our results show that the hexaploid genome of Jerusalem artichoke was formed by a hybridization event between a tetraploid and a diploid Helianthus species, followed by chromosome doubling of the hybrid, which occurred approximately 2 million years ago. Moreover, we identified more copies of actively expressed genes involved in inulin metabolism and showed that these genes may still be undergoing loss of function or sub- or neofunctionalization. These genomic resources will promote further biological studies, breeding improvement, and industrial utilization of Helianthus crops.

Genomic evidence for rediploidization and adaptive evolution following the whole-genome triplication

Whole-genome duplication (WGD), or polyploidy, events are widespread and significant in the evolutionary history of angiosperms. However, empirical evidence for rediploidization, the major process where polyploids give rise to diploid descendants, is still lacking at the genomic level. Here we present chromosome-scale genomes of the mangrove tree Sonneratia alba and the related inland plant Lagerstroemia speciosa. Their common ancestor has experienced a whole-genome triplication (WGT) approximately 64 million years ago coinciding with a period of dramatic global climate change. Sonneratia, adapting mangrove habitats, experienced extensive chromosome rearrangements post-WGT. We observe the WGT retentions display sequence and expression divergence, suggesting potential neo- and sub-functionalization. Strong selection acting on three-copy retentions indicates adaptive value in response to new environments. To elucidate the role of ploidy changes in genome evolution, we improve a model of the polyploidization-rediploidization process based on genomic evidence, contributing to the understanding of adaptive evolution during climate change.

Comparative analysis of medicinal plants Scutellaria baicalensis and common adulterants based on chloroplast genome sequencing

BACKGROUND

Scutellaria baicalensis Georgi has been extensively used as a medicinal herb in China for over 2000 years. They may be intentionally or inadvertently substituted or blended with comparable species in the local market, threatening clinical medication safety. Molecular markers are effective tools to prevent misidentification and eliminate doping and falsification among Scutellaria plants. This study screened four highly variable regions to identify Scutellaria and its adulterants. In addition, a phylogenetic analysis was performed using the complete cp genome combined with published Scutellaria species samples. Moreover, a comparative analysis of the cp genomes was conducted to investigate the cp genome evolution of S. baicalensis.

RESULTS

The complete cp genome of five species of Scutellaria was sequenced for the first time, and four previously published Scutellaria species were re-sequenced. They all exhibited a conserved quadripartite structure in their cp genomes, including two distinct regions, namely a small and large single copy region, respectively, and two inverted repeats encompassing the majority of ribosomal RNA genes. Furthermore, the nine species exhibited high conservation from aspects of the genome structure, codon usage, repeat sequences, and gene content. Four highly variable regions (matK-rps16, ndhC-trnV-UAC, psbE-petL, and rps16-trnQ-UUG) may function as potential molecular markers for differentiating S. baicalensis from its adulterants. Additionally, the monophyly of Scutellaria was ascertained and could be reclassified into two subgenera, subgenus Anaspis and subgenus Scutellaria, as evidenced by the phylogenetic analyses on sequences of cp genome and shared protein-coding sequences. According to the molecular clock analysis, it has been inferred that the divergence of Scutellaria occurred at approximately 4.0 Mya during the Pliocene Epoch.

CONCLUSION

Our study provides an invaluable theoretical basis for further Scutellaria species identification, phylogenetics, and evolution analysis.

Origin of CMS-PET1 cytotype in cultivated sunflower: A new insight

The vast majority of commercial sunflower hybrids worldwide are produced using cytoplasmic male sterility (CMS) of the PET1 type, resulting from the interspecific hybridization of Helianthus petiolaris with Helianthus annuus. Due to the fact that CMS-PET1 was not previously detected in wild sunflower, it was believed that this cytotype could arise during interspecific hybridization and is specific solely for cultivated sunflower. In this study, the open reading frame, orfH522, associated with the CMS-PET1 phenotype, was revealed for the first time in the 3'-flanking region of the mitochondrial atpA gene in wild H. annuus. An analysis of whole genome data from 1089 accessions showed that the frequency of occurrence of CMS-orfH522 in wild H. annuus populations is 3.58%, while in wild H. petiolaris populations, it is 1.26%. In general, the analysis demonstrated that PET1-CMS is a natural cytotype of H. annuus, and the appearance of the CMS phenotype in cultivated sunflowers is associated with the loss of stabilizing nuclear genes of fertility restorers, which occurred during interspecific hybridization. These data can explain the patterns of differential cytoplasmic and nuclear introgression occurring in wild sunflower and are useful for further evolutionary studies.

Comparative and phylogenetic analyses of nine complete chloroplast genomes of Orchidaceae

The orchid family has 200,000 species and 700 genera, and it is found worldwide in the tropics and subtropics. In China, there are 1247 species and subspecies of orchids belonging to the Orchidaceae family. Orchidaceae is one of the most diverse plant families in the world, known for their lush look, remarkable ecological tolerance, and capability for reproduction. It has significant decorative and therapeutic value. In terms of evolution, the orchid family is one of the more complicated groups, but up until now, little has been known about its affinities. This study examined the properties of 19 chloroplast (cp) genomes, of which 11 had previously been published and nine had only recently been revealed. Following that, topics such as analysis of selection pressure, codon usage, amino acid frequencies, repeated sequences, and reverse repeat contraction and expansion are covered. The Orchidaceae share similar cp chromosomal characteristics, and we have conducted a preliminary analysis of their evolutionary connections. The cp genome of this family has a typical tepartite structure and a high degree of consistency across species. Platanthera urceolata with more tandem repeats of the cp genome. Similar cp chromosomal traits can be seen in the orchidaceae. Galearis roborowskyi, Neottianthe cucullata, Neottianthe monophylla, Platanthera urceolata and Ponerorchis compacta are the closest cousins, according to phylogenetic study.

Phylogenomics of Aralia sect. Aralia (Araliaceae): Signals of hybridization and insights into its species delimitations and intercontinental biogeography

Genome-scale data have significantly increased the number of informative characters for phylogenetic analyses and recent studies have also revealed widespread phylogenomic discordance in many plant lineages. Aralia sect. Aralia is a small plant lineage (14 spp.) of the ginseng family Araliaceae with a disjunct distribution between eastern Asia (11 spp.) and North America (3 spp.). We herein employ sequences of hundreds of nuclear loci and the complete plastomes using targeted sequence capture and genome skimming to reconstruct the phylogenetic and biogeographic history of this section. We detected substantial conflicts among nuclear genes, yet different analytical strategies generated largely congruent topologies from the nuclear data. Significant cytonuclear discordance was detected, especially concerning the positions of the three North American species. The phylogenomic results support two intercontinental disjunctions: (1) Aralia californica of western North America is sister to the eastern Asian clade consisting of A. cordata and A. continentalis in the nuclear tree, and (2) the eastern North American A. racemosa forms a clade with A. bicrenata from southwestern North America, and the North American A. racemosa - A. bicrenata clade is then sister to the eastern Asian clade consisting of A. glabra (Japan), A. fargesii (C China), and A. apioides and A. atropurpurea (the Hengduan Mountains). Aralia cordata is supported to be disjunctly distributed in Japan, Taiwan, the Ulleung island of Korea, and in Central, Southwest and South China, and Aralia continentalis is redefined with a narrower distribution in Northeast China, eastern Russia and peninsular Korea.

Learning to tango with four (or more): the molecular basis of adaptation to polyploid meiosis

Polyploidy, which arises from genome duplication, has occurred throughout the history of eukaryotes, though it is especially common in plants. The resulting increased size, heterozygosity, and complexity of the genome can be an evolutionary opportunity, facilitating diversification, adaptation and the evolution of functional novelty. On the other hand, when they first arise, polyploids face a number of challenges, one of the biggest being the meiotic pairing, recombination and segregation of the suddenly more than two copies of each chromosome, which can limit their fertility. Both for developing polyploidy as a crop improvement tool (which holds great promise due to the high and lasting multi-stress resilience of polyploids), as well as for our basic understanding of meiosis and plant evolution, we need to know both the specific nature of the challenges polyploids face, as well as how they can be overcome in evolution. In recent years there has been a dramatic uptick in our understanding of the molecular basis of polyploid adaptations to meiotic challenges, and that is the focus of this review.

The complete mitochondrial genome of Camellia nitidissima (Theaceae)

The mitochondrial genome of Camellia nitidissima was sequenced by Illumina and Pacbio sequencing. The results of sequences showed that a total length was 949,915 bp, and the GC content was 45.7% in assembled mitochondrial genome of C. nitidissima. 71 unigenes had been found, including 36 coding proteins and 35 non-coding proteins. Subsequently, the phylogenetic tree was built on 24 plants with the maximum-likelihood method, which had high bootstrap value and fited to the angiosperm phylogeny group classification (APG IV). The study's findings unravel the taxonomic status of C. nitidissima and benefit the evolution study.

Broad diversity in monoterpene-sesquiterpene balance across wild sunflowers: Implications of leaf and floral volatiles for biotic interactions

PREMISE

As plant lineages diversify across environmental gradients, species are predicted to encounter divergent biotic pressures. This study investigated the evolution of volatile secondary metabolism across species of Helianthus.

METHODS

Leaves and petals of 40 species of wild Helianthus were analyzed via gas chromatography-mass spectrometry to determine volatile secondary metabolite profiles.

RESULTS

Across all species, 500 compounds were identified; 40% were sesquiterpenes, 18% monoterpenes, 3% diterpenes, 4% fatty acid derivatives, and 35% other compounds such as phenolics and small organic molecules. Qualitatively, annuals and species from more arid western climates had leaf compositions with a higher proportion of total monoterpenes, while erect perennials and species from more mesic eastern habitats contained a higher proportion of total sesquiterpenes. Among species, mass-based leaf monoterpene and sesquiterpene abundance were identified as largely orthogonal axes of variation by principal component analysis. Profiles for leaves were not strongly correlated with those of petals.

CONCLUSIONS

Volatile metabolites were highly diverse among wild Helianthus, indicating the value of this genus as a model system and rich genetic resource. The independence of leaf and petal volatile profiles indicates a low level of phenotypic integration between vegetative and reproductive structures, implying vegetative defense and reproductive defense or pollinator attraction functions mediated by terpene profiles in these two organs can evolve without major trade-offs. The major biosynthetic pathways for the major terpenes in wild Helianthus are already well described, providing a road map to deeper inquiry into the drivers of this diversity.

Production and cytology of Brassica autoallohexaploids with two and four copies of two subgenomes

Different digenomic Brassica autoallohexaploids were produced from the crosses of three allotetraploids and ancestral diploids and characterized for the cytological behavior of two subgenomes with two and four copies. Interspecific hybridization and allopolyploidization present an important pathway for plant evolution and breeding. In this study, different types of digenomic autoallohexaploids with two or four copies of two subgenomes (AAAACC, AACCCC, AAAABB, BBBBCC, BBCCCC) were synthesized by the crosses between three Brassica allotetraploids and their diploid progenitors and the chromosome doubling, and their meiotic behaviors were analyzed by fluorescence in situ hybridization (FISH). These autoallohexaploids showed some variations in pollen fertility and seed-sets and produced both euploid and aneuploid progenies with some chromosomes lost. Two subgenomes in these autoallohexaploids showed some aberrant pairings and segregations, and the degrees of meiotic regularity were negatively associated with the genome affinities. The chromosomes of the subgenome with four copies formed few quadrivalents with the average number < 2, and mainly paired as bivalents, and majority of the chromosomes from the subgenome with two copies gave the expected bivalents. The different extents of the equal and unequal segregations corresponded to the chromosome pairings. The development and cytological investigation of these autoallohexaploids provide not only the new germplasm for genetic research and breeding but also the new clues for the genome behavior and interplay between these subgenomes with different copies.

Allele-aware chromosome-level genome assembly of Artemisia annua reveals the correlation between ADS expansion and artemisinin yield

Artemisia annua is the major natural source of artemisinin, an anti-malarial medicine commonly used worldwide. Here, we present chromosome-level haploid maps for two A. annua strains with different artemisinin contents to explore the relationships between genomic organization and artemisinin production. High-fidelity sequencing, optical mapping, and chromatin conformation capture sequencing were used to assemble the heterogeneous and repetitive genome and resolve the haplotypes of A. annua. Approximately 50,000 genes were annotated for each haplotype genome, and a triplication event that occurred approximately 58.12 million years ago was examined for the first time in this species. A total of 3,903,467-5,193,414 variants (SNPs, indels, and structural variants) were identified in the 1.5-Gb genome during pairwise comparison between haplotypes, consistent with the high heterozygosity of this species. Genomic analyses revealed a correlation between artemisinin concents and the copy number of amorpha-4,11-diene synthase genes. This correlation was further confirmed by resequencing of 36 A. annua samples with varied artemisinin contents. Circular consensus sequencing of transcripts facilitated the detection of paralog expression. Collectively, our study provides chromosome-level allele-aware genome assemblies for two A. annua strains and new insights into the biosynthesis of artemisinin and its regulation, which will contribute to conquering malaria worldwide.

Comparative analysis of complete chloroplast genome sequences of five endangered species and new insights into phylogenetic relationships of Paris

Paris L. genus has been a precious traditional herb for more than 2000 years in China. However, due to overexploitation and habitat destruction, Paris is threatened by extinction. Similar morphological features cause the classification of Paris species in dispute. The chloroplast (cp) genome approach has been used to investigate the evolution of Paris. However, some studies confirm that the cp genome may result in misleading relationships because of the length variation, gaps/indels deletion, and incorrect models of sequence evolution in concatenated datasets. Therefore, there is a high demand for a reconstructed phylogenetic relationship and developed genetic markers to conserve these species. Recent studies have demonstrated that the protein-coding genes could provide a better phylogenetic relationship in the phylogenetic investigation. In this study, the complete cp genomes of five species were characterized, and the length of five cp genomes ranges from 162,927 bp to 165,267 bp, covering 89 protein-coding genes, 38 tRNA, and eight rRNA. The analysis of the repeat sequences, codon usage, RNA-editing sites, and comparison of cp genomes shared a high degree of conservation. Based on the protein-coding genes, the phylogenetic tree confirmed Paris's position in the order Melanthiaceae, providing maximum support for a sister relationship of the subgenera Paris sensu strict (Paris s.s.) with the Daiswa and Trillium. In addition, the molecular clock showed that subgenus Paris was inferred to have occurred at about 52.81 Mya, whereas subgenus Daiswa has originated at 24.56 Mya, which was consistent with the phylogenetic investigation. This study provided a valuable insight into the evolutionary dynamics of cp genome structure in the family Melanthiaceae, and it also contributes to the bioprospecting and conservation of Paris species.

Ultrasound affects physical and chemical properties of Jerusalem artichoke and chicory inulin

Jerusalem artichoke (Helianthus tuberosus L.) and Chicory (Cichorium intybus L.) have a heterogeneous collection of fructose polymers, known as inulin. This study was aimed to explore the effects of ultrasound (US) and autoclave (AC) on inulin physico-chemical properties as well as investigate structural characterizations and relationships with inulin physico-chemical properties. More specifically, Jerusalem artichoke powder (JA, 69.99% inulin in dry basis), purified inulin from Jerusalem artichoke (PJAI) and chicory inulin (CI) were studied to determine the effects of both treatments on reducing sugar contents, degree of polymerization (DP), water-holding capacity (WHC) and particle size. US (90 W, 20 KHZ) treatments had increased reducing sugar content up to 12.27% for PJAI, 10.86% for JA powder and 2.18% for CI. HPLC analysis showed that the DP of inulin decreased for PJAI after 2 min US treatment. WHC analysis showed that both treatments did not have significant effects (p > .05) on WHC for JA powder. This study suggests that US can be a preferable treatment for reducing the DP of inulin from JA for designing variety of food formulations. PRACTICAL APPLICATIONS: Ultrasound treatments could result in more inulin breaking down into reducing sugars, and in the decrease of inulin DP. This research suggested that the DP of inulin might be a very important factor in ultrasound treatment for their affect in the absorption of energy from ultrasound. Therefore, ultrasound can be a desirable treatment for changing the degree of polymerization of inulin from JA for designing different food products. Future studies need to investigate the relationship between the viscosity of inulin solution and the de-polymerization of inulin caused by ultrasound treatment.

Organelle Phylogenomics and Extensive Conflicting Phylogenetic Signals in the Monocot Order Poales

The Poales is one of the largest orders of flowering plants with significant economic and ecological values. Reconstructing the phylogeny of the Poales is important for understanding its evolutionary history that forms the basis for biological studies. However, due to sparse taxon sampling and limited molecular data, previous studies have resulted in a variety of contradictory topologies. In particular, there are three nodes surrounded by incongruence: the phylogenetic ambiguity near the root of the Poales tree, the sister family of Poaceae, and the delimitation of the xyrid clade. We conducted a comprehensive sampling and reconstructed the phylogenetic tree using plastid and mitochondrial genomic data from 91 to 66 taxa, respectively, representing all the 16 families of Poales. Our analyses support the finding of Bromeliaceae and Typhaceae as the earliest diverging groups within the Poales while having phylogenetic relationships with the polytomy. The clade of Ecdeiocoleaceae and Joinvilleaceae is recovered as the sister group of Poaceae. The three families, Mayacaceae, Eriocaulaceae, and Xyridaceae, of the xyrid assembly diverged successively along the backbone of the Poales phylogeny, and thus this assembly is paraphyletic. Surprisingly, we find substantial phylogenetic conflicts within the plastid genomes of the Poales, as well as among the plastid, mitochondrial, and nuclear data. These conflicts suggest that the Poales could have a complicated evolutionary history, such as rapid radiation and polyploidy, particularly allopolyploidy through hybridization. In sum, our study presents a new perspicacity into the complex phylogenetic relationships and the underlying phylogenetic conflicts within the Poales.

The Insights into Mitochondrial Genomes of Sunflowers

The significant difference in the mtDNA size and structure with simultaneous slow evolving genes makes the mitochondrial genome paradoxical among all three DNA carriers in the plant cell. Such features make mitochondrial genome investigations of particular interest. The genus Helianthus is a diverse taxonomic group, including at least two economically valuable species-common sunflower (H. annuus) and Jerusalem artichoke (H. tuberosus). The successful investigation of the sunflower nuclear genome provided insights into some genomics aspects and significantly intensified sunflower genetic studies. However, the investigations of organelles' genetic information in Helianthus, especially devoted to mitochondrial genomics, are presented by limited studies. Using NGS sequencing, we assembled the complete mitochondrial genomes for H. occidentalis (281,175 bp) and H. tuberosus (281,287 bp) in the current investigation. Besides the master circle chromosome, in the case of H. tuberosus, the 1361 bp circular plasmid was identified. The mitochondrial gene content was found to be identical for both sunflower species, counting 32 protein-coding genes, 3 rRNA, 23 tRNA genes, and 18 ORFs. The comparative analysis between perennial sunflowers revealed common and polymorphic SSR and SNPs. Comparison of perennial sunflowers with H. annuus allowed us to establish similar rearrangements in mitogenomes, which have possibly been inherited from a common ancestor after the divergence of annual and perennial sunflower species. It is notable that H. occidentalis and H. tuberosus mitogenomes are much more similar to H. strumosus than H. grosseserratus.

Comparative analysis of chloroplast genomes of seven perennial Helianthus species

Sequencing and a comparative analysis of the complete chloroplast genomes of seven perennial Helianthus species were carried out. The chloroplast genomes have a typical quadripartite structure, including large and small single regions and a pair of inverted repeats. Genome sizes were between 151,152 bp and 151,289 bp. The genome of H. grosseserratus was the smallest, while that of H. microcephalus was the largest. The size variation of the chloroplast genomes is substantially due to the change in the length of simple sequence repeats (SSRs) in non-coding regions. An analysis of these SSRs revealed 35 polymorphic loci (average PIC value > 0.5) that can be used to examine ecological and evolutionary processes in wild Helianthus species. Eight divergence hotspots, including five intergenic regions (petN-psbM, clpP intron, rps3-rpl16, ndhD-ccsA, and ndhF-rpl32) and three gene regions (rbcL, ycf1, and ndhF) were also identified in Helianthus chloroplast genomes. The evolutionary selection pressure analysis revealed a strong purifying selection. Only the rbcL gene experienced efficiency of positive selection at the annual/perennial transitions. The inverted repeat (IR)/single copy (SC) boundaries were identical in all of these (Helianthus) species. In general, the comparison of the genomes revealed low levels of sequence variability (Pi = 0.00051). This indicates that the chloroplast genomes of the studied perennial species of Helianthus, in addition to purifying selection, are closely related and have a recent divergence time.

Host shifting and host sharing in a genus of specialist flies diversifying alongside their sunflower hosts

Congeneric parasites are unlikely to specialize on the same tissues of the same host species, likely because of strong multifarious selection against niche overlap. Exceptions where >1 congeneric species use the same tissues reveal important insights into ecological factors underlying the origins and maintenance of diversity. Larvae of sunflower maggot flies in the genus Strauzia feed on plants in the family Asteraceae. Although Strauzia tend to be host specialists, some species specialize on the same hosts. To resolve the origins of host sharing among these specialist flies, we used reduced representation genomic sequencing to infer the first multilocus phylogeny of genus Strauzia. Our results show that Helianthus tuberosus and Helianthus grosseserratus each host three different Strauzia species and that the flies co-occurring on a host are not one another's closest relatives. Though this pattern implies that host sharing is most likely the result of host shifts, these may not all be host shifts in the conventional sense of an insect moving onto an entirely new plant. Many hosts of Strauzia belong to a clade of perennial sunflowers that arose 1-2 MYA and are noted for frequent introgression and hybrid speciation events. Our divergence time estimates for all of the Helianthus-associated Strauzia are within this same time window (<1 MYA), suggesting that rapid and recent adaptive introgression and speciation in Helianthus may have instigated the diversification of Strauzia, with some flies converging upon a single plant host after their respective ancestral host plants hybridized to form a new sunflower species.

Genome skimming and exploration of DNA barcodes for Taiwan endemic cypresses

Cypresses are characterized by their longevity and valuable timber. In Taiwan, two endemic cypress species, Chamaecyparis formosensis and C. obtusa var. formosana, are threatened by prevalent illegal logging. A DNA barcode system is urgently needed for reforestation and conservation of these two cypresses. In this study, both plastomes and 35S rDNAs from 16, 10, and 6 individuals of C. formosensis, C. obtusa var. formosana, and C. obtusa var. obtusa were sequenced, respectively. We show that the loss of plastid trnT-GGU readily distinguishes C. formosensis from its congeneric species. We demonstrate that entire sequences of plastomes or 35S rDNAs are capable of correctly identifying cypress species and varieties, suggesting that they are effective super-barcodes. We also discover three short hypervariable loci (i.e., 3′ETS, ITS1, and trnH-psbA) that are promising barcodes for identifying cypress species and varieties. Moreover, nine species-specific indels of > 100 bp were detected in the cypress plastomes. These indels, together with the three aforementioned short barcodes, constitute an alternative and powerful barcode system crucial for identifying specimens that are fragmentary or contain degraded/poor DNA. Our sequenced data and barcode systems not only enrich the genetic reference for cypresses, but also contribute to future reforestation, conservation, and forensic investigations.

Comparative Analyses of the Chloroplast Genomes of Patchouli Plants and Their Relatives in Pogostemon (Lamiaceae)

Pogostemon Desf., the largest genus of the tribe Pogostemoneae (Lamiaceae), consists of ca. 80 species distributed mainly from South and Southeast Asia to China. The genus contains many patchouli plants, which are of great economic importance but taxonomically difficult. Therefore, it is necessary to characterize more chloroplast (cp) genomes for infrageneric phylogeny analyses and species identification of Pogostemon, especially for patchouli plants. In this study, we newly generated four cp genomes for three patchouli plants (i.e., Pogostemon plectranthoides Desf., P. septentrionalis C. Y. Wu et Y. C. Huang, and two cultivars of P. cablin (Blanoco) Benth.). Comparison of all samples (including online available cp genomes of P. yatabeanus (Makino) Press and P. stellatus (Lour.) Kuntze) suggested that Pogostemon cp genomes are highly conserved in terms of genome size and gene content, with a typical quadripartite circle structure. Interspecific divergence of cp genomes has been maintained at a relatively low level, though seven divergence hotspot regions were identified by stepwise window analysis. The nucleotide diversity (Pi) value was correlated significantly with gap proportion (indels), but significantly negative with GC content. Our phylogenetic analyses based on 80 protein-coding genes yielded high-resolution backbone topologies for the Lamiaceae and Pogostemon. For the overall mean substitution rates, the synonymous (d_S) and nonsynonymous (d_N) substitution rate values of protein-coding genes varied approximately threefold, while the d_N values among different functional gene groups showed a wider variation range. Overall, the cp genomes of Pogostemon will be useful for phylogenetic reconstruction, species delimitation and identification in the future.

The Investigation of Perennial Sunflower Species (Helianthus L.) Mitochondrial Genomes

The genus Helianthus is a diverse taxonomic group with approximately 50 species. Most sunflower genomic investigations are devoted to economically valuable species, e.g., H. annuus, while other Helianthus species, especially perennial, are predominantly a blind spot. In the current study, we have assembled the complete mitogenomes of two perennial species: H. grosseserratus (273,543 bp) and H. strumosus (281,055 bp). We analyzed their sequences and gene profiles in comparison to the available complete mitogenomes of H. annuus. Except for sdh4 and trnA-UGC, both perennial sunflower species had the same gene content and almost identical protein-coding sequences when compared with each other and with annual sunflowers (H. annuus). Common mitochondrial open reading frames (ORFs) (orf117, orf139, and orf334) in sunflowers and unique ORFs for H. grosseserratus (orf633) and H. strumosus (orf126, orf184, orf207) were identified. The maintenance of plastid-derived coding sequences in the mitogenomes of both annual and perennial sunflowers and the low frequency of nonsynonymous mutations point at an extremely low variability of mitochondrial DNA (mtDNA) coding sequences in the Helianthus genus.

Beyond DNA barcoding: The unrealized potential of genome skim data in sample identification

Genetic tools are increasingly used to identify and discriminate between species. One key transition in this process was the recognition of the potential of the ca 658bp fragment of the organelle cytochrome c oxidase I (COI) as a barcode region, which revolutionized animal bioidentification and lead, among others, to the instigation of the Barcode of Life Database (BOLD), containing currently barcodes from >7.9 million specimens. Following this discovery, suggestions for other organellar regions and markers, and the primers with which to amplify them, have been continuously proposed. Most recently, the field has taken the leap from PCR-based generation of DNA references into shotgun sequencing-based "genome skimming" alternatives, with the ultimate goal of assembling organellar reference genomes. Unfortunately, in genome skimming approaches, much of the nuclear genome (as much as 99% of the sequence data) is discarded, which is not only wasteful, but can also limit the power of discrimination at, or below, the species level. Here, we advocate that the full shotgun sequence data can be used to assign an identity (that we term for convenience its "DNA-mark") for both voucher and query samples, without requiring any computationally intensive pretreatment (e.g. assembly) of reads. We argue that if reference databases are populated with such "DNA-marks," it will enable future DNA-based taxonomic identification to complement, or even replace PCR of barcodes with genome skimming, and we discuss how such methodology ultimately could enable identification to population, or even individual, level.

Unraveling the Chloroplast Genomes of Two Prosopis Species to Identify Its Genomic Information, Comparative Analyses and Phylogenetic Relationship

Genus Prosopis (family Fabaceae) are shrubby trees, native to arid and semi-arid regions of Asia, Africa, and America and known for nitrogen fixation. Here, we have sequenced the complete chloroplast (cp) genomes of two Prosopis species (P. juliflora and P. cineraria) and compared them with previously sequenced P. glandulosa, Adenanthera microsperma, and Parkia javanica belonging to the same family. The complete genome sequences of Prosopis species and related species ranged from 159,389 bp (A. microsperma) to 163,677 bp (P. cineraria). The overall GC contents of the genomes were almost the similar (35.9–36.6%). The P. juliflora and P. cineraria genomes encoded 132 and 131 genes, respectively, whereas both the species comprised of 85 protein-coding genes higher than other compared species. About 140, 134, and 129 repeats were identified in P. juliflora, P. cineraria and P. glandulosa cp genomes, respectively. Similarly, the maximum number of simple sequence repeats were determined in P. juliflora (88), P. cineraria (84), and P. glandulosa (78). Moreover, complete cp genome comparison determined a high degree of sequence similarity among P. juliflora, P. cineraria, and P. glandulosa, however some divergence in the intergenic spacers of A. microsperma and Parkia javanica were observed. The phylogenetic analysis showed that P. juliflora is closer to P. cineraria than P. glandulosa.

Comparative analyses of chloroplast genomes from 13 Lagerstroemia (Lythraceae) species: identification of highly divergent regions and inference of phylogenetic relationships

Seven divergence hotspots as plastid markers for DNA barcoding was selected, and the phylogeny of 13 Lagerstroemia species based on the cp genome data was reconstructed within Myrtales. The Lagerstroemia species used in this study originated in China and have high economic and ecological value. The shared interspecific morphological characteristics and intraspecific morphological variation resulting from hybridization among Lagerstroemia taxa have made resolving their classification problems and phylogenetic relationships difficult. Systematic comparative genomic analysis has been shown to resolve phylogenetic relationships. We sequenced and annotated 6 Lagerstroemia cp genomes (Lagerstroemia excelsa, Lagerstroemia limii, Lagerstroemia siamica, Lagerstroemia tomentosa, Lagerstroemia venusta, and Lagerstroemia calyculata) for the first time and combined them with previously published genomes for Lagerstroemia species. Bioinformatics was used to analyse the 13 cp genomes in terms of gene structure and organization, codon usage, contraction and expansion of inverted repeat regions, repeat structure, divergence hotspots, species pairwise Ka/Ks ratios and phylogenetic relationships. The length varied between 152,049 bp in Lagerstroemia subcostata and 152,521 bp in L. venusta. We selected seven divergence hotspots in the cp genomes that had the potential to act as plastid markers to distinguish Lagerstroemia species. The phylogenetic relationships within Myrtales inferred from the cp genomes of 13 Lagerstroemia species and 27 other Myrtales species were highly supported, which illustrated several novel relationships within Myrtales. Taken together, our results provide comprehensive chloroplast genomic resources, which can be used further for species identification and molecular breeding of Lagerstroemia species.

Whole-Transcriptome Analysis Unveils the Synchronized Activities of Genes for Fructans in Developing Tubers of the Jerusalem Artichoke

Helianthus tuberosus L., known as the Jerusalem artichoke, is a hexaploid plant species, adapted to low-nutrient soils, that accumulates high levels of inulin in its tubers. Inulin is a fructose-based polysaccharide used either as dietary fiber or for the production of bioethanol. Key enzymes involved in inulin biosynthesis are well known. However, the gene networks underpinning tuber development and inulin accumulation in H. tuberous remain elusive. To fill this gap, we selected 6,365 expressed sequence tags (ESTs) from an H. tuberosus library to set up a microarray platform and record their expression across three tuber developmental stages, when rhizomes start enlarging (T0), at maximum tuber elongation rate (T3), and at tuber physiological maturity (Tm), in "VR" and "K8-HS142"clones. The former was selected as an early tuberizing and the latter as a late-tuberizing clone. We quantified inulin and starch levels, and qRT-PCR confirmed the expression of critical genes accounting for inulin biosynthesis. The microarray analysis revealed that the differences in morphological and physiological traits between tubers of the two clones are genetically determined since T0 and that is relatively low the number of differentially expressed ESTs across the stages shared between the clones (93). The expression of ESTs for sucrose:sucrose 1-fructosyltransferase (1-SST) and fructan:fructan 1-fructosyltransferase (1-FFT), the two critical genes for fructans polymerization, resulted to be temporarily synchronized and mirror the progress of inulin accumulation and stretching. The expression of ESTs for starch biosynthesis was insignificant throughout the developmental stages of the clones in line with the negligible level of starch into their mature tubers, where inulin was the dominant polysaccharide. Overall, our study disclosed candidate genes underpinning the development and storage of carbohydrates in the tubers of two H. tuberosus clones. A model according to which the steady-state levels of 1-SST and 1-FFT transcripts are developmentally controlled and might represent a limiting factor for inulin accumulation has been provided. Our finding may have significant repercussions for breeding clones with improved levels of inulin for food and chemical industry.

Identification of miR390-TAS3-ARF pathway in response to salt stress in Helianthus tuberosus L

MicroRNA390 (miR390), an ancient and highly conserved miRNA family in land plants, plays multiple roles in plant growth, development and stress responses. In this study, we isolated and identified MIR390, miR390, TAS3a/b/c, tasiARF-1/2/3 (trans-acting small interfering RNAs influencing Auxin Response Factors) and ARF2/3/4 in Jerusalem artichoke (Helianthus tuberosus L.). Treatment with 100 mM NaCl induced expression of miR390, increased cleavage of TAS3, produced high levels of tasiARFs, and subsequently enhanced cleavage of ARF3/4, which was most likely associated with salt tolerance of the plants. In contrast, treatment with 300 mM NaCl inhibited expression of miR390, attenuated cleavage of TAS3, produced a small amount of tasiARFs, and reduced cleavage of ARF3/4. We proposed that ARF2, one of the targets of tasiARFs, induced under salinity was likely to play an active role in salt tolerance of Jerusalem artichoke. The study of the miR390-TAS3-ARF model in Jerusalem artichoke may broaden our understanding of salt tolerance mechanisms, and provides a theoretical support for further genetic identification and breeding crops with increased tolerance to salt stress.

Skim-Sequencing Reveals the Likely Origin of the Enigmatic Endangered Sunflower Helianthus schweinitzii

Resolving the origin of endangered taxa is an essential component of conservation. This information can be used to guide efforts of bolstering genetic diversity, and also enables species recovery and future evolutionary studies. Here, we used low-coverage whole genome sequencing to clarify the origin of Helianthus schweinitzii, an endangered tetraploid sunflower that is endemic to the Piedmont Plateau in the eastern United States. We surveyed four accessions representing four populations of H. schweinitzii and 38 accessions of six purported parental species. Using de novo approaches, we assembled 87,004 bp of the chloroplast genome and 6770 bp of the nuclear 35S rDNA. Phylogenetic reconstructions based on the chloroplast genome revealed no reciprocal monophyly of taxa. In contrast, nuclear rDNA data strongly supported the currently accepted sections of the genus Helianthus. Information from combined cpDNA and rDNA provided evidence that H. schweinitzii is likely an allo-tetraploid that formed as a result of hybridization between the diploids Helianthus giganteus and Helianthus microcephalus.

Plastome phylogenomics, biogeography, and clade diversification of Paris (Melanthiaceae)

BACKGROUND

Paris (Melanthiaceae) is an economically important but taxonomically difficult genus, which is unique in angiosperms because some species have extremely large nuclear genomes. Phylogenetic relationships within Paris have long been controversial. Based on complete plastomes and nuclear ribosomal DNA (nrDNA) sequences, this study aims to reconstruct a robust phylogenetic tree and explore historical biogeography and clade diversification in the genus.

RESULTS

All 29 species currently recognized in Paris were sampled. Whole plastomes and nrDNA sequences were generated by the genome skimming approach. Phylogenetic relationships were reconstructed using the maximum likelihood and Bayesian inference methods. Based on the phylogenetic framework and molecular dating, biogeographic scenarios and historical diversification of Paris were explored. Significant conflicts between plastid and nuclear datasets were identified, and the plastome tree is highly congruent with past interpretations of the morphology. Ancestral area reconstruction indicated that Paris may have originated in northeastern Asia and northern China, and has experienced multiple dispersal and vicariance events during its diversification. The rate of clade diversification has sharply accelerated since the Miocene/Pliocene boundary.

CONCLUSIONS

Our results provide important insights for clarifying some of the long-standing taxonomic debates in Paris. Cytonuclear discordance may have been caused by ancient and recent hybridizations in the genus. The climatic and geological changes since the late Miocene, such as the intensification of Asian monsoon and the rapid uplift of Qinghai-Tibet Plateau, as well as the climatic fluctuations during the Pleistocene, played essential roles in driving range expansion and radiative diversification in Paris. Our findings challenge the theoretical prediction that large genome sizes may limit speciation.

Whole Plastome Sequencing Within Silene Section Psammophilae Reveals Mainland Hybridization and Divergence With the Balearic Island Populations

Reconstructing the phylogenetic relationships within Caryophyllaceae tribe Sileneae has been obscured by hybridization and incomplete lineage sorting. Silene is the largest genus in the Caryophyllaceae, and unraveling its evolutionary history has been particularly challenging. In order to infer the phylogenetic relationships among the five species in Silene section Psammophilae, we have performed a genome skimming approach to acquire the complete plastid genome (cpDNA), nuclear ribosomal cistron (nrDNA), and partial mitochondrial genome (mtDNA). We have included 26 populations, representing the range of each species' distribution. This section includes five morphologically similar species endemic to the Iberian Peninsula and Balearic Islands (Ibiza and Formentera), yet some of them occupy distinct edaphic habitats (e.g. maritime sands, calcareous sandstones). In addition to phylogeographic analyses, genetic structuring using the chloroplast data set was inferred with Discriminant Analysis of Principal Components (DAPC), analyses of molecular variance (AMOVA), and a partial Mantel test. Reference-guided assembly of 50 bp single-end and 250 bp paired-end Illumina reads produced the nearly complete cpDNA genome (154 kbp), partial mtDNA genome (from 81 to 114 kbp), and the nrDNA cistron (6.4 kbp). Selected variable regions of the cpDNA and mtDNA assemblies were confirmed by Sanger sequencing. Phylogenetic analyses of the mainland populations reveal incongruence among the three genomes. None of the three data sets produced relationships consistent with taxonomy or geography. In contrast, Silene cambessedesii, present in the Balearic Islands, is the only species that forms a strongly supported monophyletic clade in the cpDNA genome and is strongly differentiated with respect to the remaining taxa of the Iberian Peninsula. These results contrast with those obtained for mainland populations. Across the entire analysis, only one well-supported mainland clade of Silene littorea and Silene stockenii emerges from the southern region of the Iberian Peninsula. DAPC and AMOVA results suggest the absence of genetic structure among mainland populations of Silene section Psammophilae, whereas partial Mantel test discarded spatial correlation of genetic differentiation. The widespread incongruence between morphology-based taxonomic boundaries and phylogeography suggests a history of interspecific hybridization, in which only a substantial geographic barrier, like isolation by the Mediterranean Sea, was sufficient to create and maintain species boundaries in Silene section Psammophilae.

Plastome of Quercus xanthoclada and comparison of genomic diversity amongst selected Quercus species using genome skimming

The genus Quercus L. contains several of the most economically important species for timber production in the Northern Hemisphere. It was one of the first genera described, but genetic diversity at a global scale within and amongst oak species remains unclear, despite numerous regional or species-specific assessments. To evaluate global plastid diversity in oaks, we sequenced the complete chloroplast of Quercus xanthoclada and compared its sequence with those available from other main taxonomic groups in Quercus. We quantify genomic divergence amongst oaks and performed a sliding window analysis to detect the most variable regions amongst members of the various clades, as well as divergent regions occurring in specific pairs of species. We identified private and shared SNPs amongst oaks species and sections and stress the need for a large global assessment of genetic diversity in this economically and ecologically important genus.

Identification and analysis of RNA editing sites in chloroplast transcripts of kenaf (Hibiscus cannabinus L.)

RNA editing is one of the post-transcriptional modification processes and can lead to changes in sequencing and functioning of corresponding proteins and genetic information. To reveal the composition and characteristic of RNA editing of kenaf chloroplast genome, the RNA editing sites in kenaf chloroplast were predicted and identified using bioinformatics and RT-PCR analysis. The prediction results showed a total of 48 editing sites distributed in 22 genes, all of them were C to U conversion leading to amino acid changes. Further analysis of the position of RNA editing sites revealed that except 11 editing sites located at the first codon base, the other editing sites were found at the second codon base. Then four genes were randomly selected to validate the editing sites. Results showed that it was accurate to study the chloroplast RNA editing sites by bioinformatics method accompanied with cloning sequencing. Furthermore, the protein secondary structure and transmembrane domain of ndhD and atpA that had undergone gene editing also changed after editing. This implied that proteins with structural changes may have an impact on kenaf growth. Meanwhile, the differential editing site was found in chloroplast transcripts in kenaf CMS line and its maintainer line, indicating that chloroplast RNA editing could be associated with kenaf CMS. Therefore, the present study laid a foundation to further reveal the biological functioning of chloroplast RNA editing in CMS and its maintainer lines in kenaf.

The complete chloroplast genome of the Jerusalem artichoke (Helianthus tuberosus L.) and an adaptive evolutionary analysis of the ycf2 gene

Jerusalem artichoke (Helianthus tuberosus L.) is widely cultivated in Northwest China, and it has become an emerging economic crop that is rapidly developing. Because of its elevated inulin content and high resistance, it is widely used in functional food, inulin processing, feed, and ecological management. In this study, Illumina sequencing technology was utilized to assemble and annotate the complete chloroplast genome sequences of Jerusalem artichoke. The total length was 151,431 bp, including four conserved regions: A pair of reverse repeat regions (IRa 24,568 bp and IRb 24,603 bp), a large single-copy region (83,981 bp), and a small single-copy region (18,279 bp). The genome had a total of 115 genes, with 19 present in the reverse direction in the IR region. A total of 36 simple sequence repeats (SSRs) were identified in the coding and non-coding regions, most of which were biased toward A/T bases. A total of 32 SSRs were distributed in the non-coding regions. A comparative analysis of the chloroplast genome sequence of the Jerusalem artichoke and other species of the composite family revealed that the chloroplast genome sequences of plants of the composite family were highly conserved. Differences were observed in 24 gene loci in the coding region, with the degree of differentiation of the ycf2 gene being the most obvious. A phylogenetic analysis showed that H. petiolaris subsp. fallax had the closest relationship with Jerusalem artichoke, both members of the Helianthus genus. Selective locus detection of the ycf2 gene in eight species of the composite family was performed to explore adaptive evolution traits of the ycf2 gene in Jerusalem artichoke. The results show that there are significant and extremely significant positive selection sites at the 1239N and 1518R loci, respectively, indicating that the ycf2 gene has been subject to adaptive evolution. Insights from our assessment of the complete chloroplast genome sequences of Jerusalem artichoke will aid in the in-depth study of the evolutionary relationship of the composite family and provide significant sequencing information for the genetic improvement of Jerusalem artichoke.

Chloroplast genome analysis of box-ironbark Eucalyptus

Eucalyptus L'Hérit. (Myrtaceae) is a taxonomically complex and highly speciose genus that dominates much of Australia's woody vegetation. However, very little information is available about the molecular biology and chloroplast diversity of certain groups, such as Eucalyptus section Adnataria, which is found in many woodland habitats of eastern Australia. We report four new complete chloroplast genomes of Eucalyptus, including three genomes from species previously lacking any chloroplast reference sequences. Plastomes of E. albens, E. conica, E. crebra and E. melliodora assembled using a de novo approach were shown to be largely identical to each other, and similar in size and structure to previously published chloroplast genomes from Eucalyptus. A total of 132 genes (114 single-copy genes and 18 duplicated genes in the IR regions) were identified, and shown to be highly conserved in terms of gene order, content and organization. Slightly higher divergence in the intergenic spacers was identified through comparative genomic analyses. Chloroplast sequences of 35 additional individuals representing 12 species were assembled using a reference guided approach. Rates of nucleotide substitution varied among the protein coding genes, with 17 genes under possible positive selection, and 29 invariant genes. Phylogenetic analysis of either the whole reconstructed plastome sequences or the individual genes revealed extreme discordance with expected species boundaries or higher-level relationships. Plastome relationships were better predicted by geography than by nuclear DNA or taxonomic relationships, suggesting a substantial influence of gene flow over and above the effects of incomplete lineage sorting. These results provide resources for future research and valuable insights into the prevalence of interspecific gene flow among Eucalyptus species.

Characterization of the Tibet plateau Jerusalem artichoke (Helianthus tuberosus L.) transcriptome by de novo assembly to discover genes associated with fructan synthesis and SSR analysis

BACKGROUND

Jerusalem artichoke (Helianthus tuberosus L.) is a characteristic crop in the Qinghai-Tibet Plateau which has rapidly developed and gained socioeconomic importance in recent years. Fructans are abundant in tubers and represent the foundation for their formation, processing and utilization of yield; and are also widely used in new sugar-based materials, bioenergy processing, ecological management, and functional feed. To identify key genes in the metabolic pathway of fructans in Jerusalem artichoke, high-throughput sequencing was performed using Illumina Hi Seq™ 2500 equipment to construct a transcriptome library.

RESULTS

Qinghai-Tibet Plateau Jerusalem artichoke "Qingyu No.1" was used as the material; roots, stems, leaves, flowers and tubers of Jerusalem artichoke in its flowering stage were mixed into a mosaic of the Jerusalem artichoke transcriptome library, obtaining 63,089 unigenes with an average length of 713.6 bp. Gene annotation through the Nr, Swiss Prot, GO, KOG and KEGG databases revealed 34.95 and 46.91% of these unigenes had similar sequences in the Nr and Swiss Prot databases. The GO classification showed the Jerusalem artichoke unigenes were divided into three ontologies, with a total of 49 functional groups encompassing biological processes, cellular components, and molecular functions. Among them, there were more unigenes involved in the functional groups for cellular processes, metabolic processes, and single-organism processes. 38,999 unigenes were annotated by KOG and divided into 25 categories according to their functions; the most common annotation being general function prediction. A total of 13,878 unigenes (22%) were annotated in the KEGG database, with the largest proportion corresponding to pathways related to carbohydrate metabolism. A total of 12 unigenes were involved in the synthesis and degradation of fructan. Cluster analysis revealed the candidate 12 unigene proteins were dispersed in the 5 major families of proteins involved in fructan synthesis and degradation. The synergistic effect of INV gene is necessary during fructose synthesis and degradation in Jerusalem artichoke tuber development. The sequencing data from the transcriptome of this species can provide a reliable data basis for the identification and assessment of the expression of the members of the INV gene family.A simple sequence repeat (SSR) loci search was performed on the transcriptome data of Jerusalem artichoke, identifying 6635 eligible SSR loci with a large proportion of dinucleotide and trinucleotide repeats, and the most different motifs were repeated 5 times and 6 times. Dinucleotide and trinucleotide repeat motifs were the most frequent, with AG/CT and ACC/GGT repeat motifs accounting for the highest proportion.

CONCLUSIONS

In this study, a database search of the transcriptome of the Jerusalem artichoke from the Qinghai Tibet Plateau was conducted by high throughput sequencing technology to obtain important transcriptional and SSR loci information. This allowed characterization of the overall expression features of the Jerusalem artichoke transcriptome, identifying the key genes involved in metabolism in this species. In turn, this offers a foundation for further research on the regulatory mechanisms of fructan metabolism in Jerusalem artichoke.

Phylogenetic trends and environmental correlates of nuclear genome size variation in Helianthus sunflowers

Flowering plants serve as a powerful model for studying the evolution of nuclear genome size (GS) given the tremendous GS variation that exists both within and across angiosperm lineages. Helianthus sunflowers consist of c. 50 species native to North America that occupy diverse habitats and vary in ploidy level. In the current study, we generated a comprehensive GS database for 49 Helianthus species using flow cytometric approaches. We examined variability across the genus and present a comparative phylogenetic analysis of GS evolution in diploid Helianthus species. Results demonstrated that different clades of diploid Helianthus species showed evolutionary patterns of GS contraction, expansion and relative stasis, with annual diploid species evolving smaller GS with the highest rate of evolution. Phylogenetic comparative analyses of diploids revealed significant negative associations of GS with temperature seasonality and cell production rate, indicating that the evolution of larger GS in Helianthus diploids may be more permissible in habitats with longer growing seasons where selection for more rapid growth may be relaxed. The Helianthus GS database presented here and corresponding analyses of environmental and phenotypic correlates will facilitate ongoing and future research on the ultimate drivers of GS evolution in this well-studied North American plant genus.

An evaluation of alternative explanations for widespread cytonuclear discordance in annual sunflowers (Helianthus)

Cytonuclear discordance is commonly observed in phylogenetic studies, yet few studies have tested whether these patterns reflect incomplete lineage sorting or organellar introgression. Here, we used whole-chloroplast sequence data in combination with over 1000 nuclear single-nucleotide polymorphisms to clarify the extent of cytonuclear discordance in wild annual sunflowers (Helianthus), and to test alternative explanations for such discordance. Our phylogenetic analyses indicate that cytonuclear discordance is widespread within this group, both in terms of the relationships among species and among individuals within species. Simulations of chloroplast evolution show that incomplete lineage sorting cannot explain these patterns in most cases. Instead, most of the observed discordance is better explained by cytoplasmic introgression. Molecular tests of evolution further indicate that selection may have played a role in driving patterns of plastid variation - although additional experimental work is needed to fully evaluate the importance of selection on organellar variants in different parts of the geographic range. Overall, this study represents one of the most comprehensive tests of the drivers of cytonuclear discordance and highlights the potential for gene flow to lead to extensive organellar introgression in hybridizing taxa.

Wild Origins of Macadamia Domestication Identified Through Intraspecific Chloroplast Genome Sequencing

Identifying the geographic origins of crops is important for the conservation and utilization of novel genetic variation. Even so, the origins of many food crops remain elusive. The tree nut crop macadamia has a remarkable domestication history, from subtropical rain forests in Australia through Hawaii to global cultivation all within the last century. The industry is based primarily on Macadamia integrifolia and M. integrifolia-M. tetraphylla hybrid cultivars with Hawaiian cultivars the main contributors to world production. Sequence data from the chloroplast genome assembled using a genome skimming strategy was used to determine population structure among remnant populations of the main progenitor species, M. integrifolia. Phylogenetic analysis of a 506 bp chloroplast SNP alignment from 64 wild and cultivated accessions identified phylogeographic structure and deep divergences between clades providing evidence for historical barriers to seed dispersal. High levels of variation were detected among wild accessions. Most Hawaiian cultivars, however, shared a single chlorotype that was also present at two wild sites at Mooloo and Mt Bauple from the northernmost distribution of the species in south-east Queensland. Our results provide evidence for a maternal genetic bottleneck during early macadamia domestication, and pinpoint the likely source of seed used to develop the Hawaiian cultivars. The extensive variability and structuring of M. integrifolia chloroplast genomic variation detected in this study suggests much unexploited genetic diversity is available for improvement of this recently domesticated crop.

Eriobotrya Belongs to Rhaphiolepis (Maleae, Rosaceae): Evidence From Chloroplast Genome and Nuclear Ribosomal DNA Data

The Eriobotrya-Rhaphiolepis (ER) clade consists of about 46 species distributed in East and Southeast Asia. Although Eriobotrya and Rhaphiolepis have been supported to form a clade, the monophyly of Eriobotrya and Rhaphiolepis at the genus level has never been well tested and their phylogenetic positions in Maleae still remain uncertain. This study aims to reconstruct a robust phylogeny of the ER clade in the framework of Maleae with a broad taxon sampling and clarify the phylogenetic relationship between Eriobotrya and Rhaphiolepis. This study employed sequences of the whole plastome (WP) and entire nuclear ribosomal DNA (nrDNA) repeats assembled from the genome skimming approach and included 83 samples representing 76 species in 32 genera of Rosaceae, especially Maleae. The Maximum Likelihood (ML) and Bayesian Analysis (BI) based on three datasets, i.e., WP, coding sequences of plastome (CDS), and nrDNA, strongly supported the paraphyly of Eriobotrya, within which Rhaphiolepis was nested. Our plastid tree supported the sister relationship between the ER clade and Heteromeles, and the nrDNA tree, however, did not resolve the phylogenetic placement of the ER clade in Maleae. Strong incongruence between the plastid and the nuclear trees is most likely explained by hybridization events, which may have played an important role in the evolutionary history of the ER clade. Molecular, morphological, and geographic evidence all supports the merge of Eriobotrya with Rhaphiolepis, which has the nomenclatural priority. We herein transferred 36 taxa of Eriobotrya to Rhaphiolepis. We also proposed a new name, Rhaphiolepis loquata B.B.Liu & J.Wen, for the economically important loquat, as the specific epithet "japonica" was pre-occupied in Rhaphiolepis.

Ancient DNA reveals the timing and persistence of organellar genetic bottlenecks over 3,000 years of sunflower domestication and improvement

Here, we report a comprehensive paleogenomic study of archaeological and ethnographic sunflower remains that provides significant new insights into the process of domestication of this important crop. DNA from both ancient and historic contexts yielded high proportions of endogenous DNA, and although archaeological DNA was found to be highly degraded, it still provided sufficient coverage to analyze genetic changes over time. Shotgun sequencing data from specimens from the Eden's Bluff archaeological site in Arkansas yielded organellar DNA sequence from specimens up to 3,100 years old. Their sequences match those of modern cultivated sunflowers and are consistent with an early domestication bottleneck in this species. Our findings also suggest that recent breeding of sunflowers has led to a loss of genetic diversity that was present only a century ago in Native American landraces. These breeding episodes also left a profound signature on the mitochondrial and plastid haplotypes in cultivars, as two types were intentionally introduced from other Helianthus species for crop improvement. These findings gained from ancient and historic sunflower specimens underscore how future in-depth gene-based analyses can advance our understanding of the pace and targets of selection during the domestication of sunflower and other crop species.

Assembly of a Complete Mitogenome of Chrysanthemum nankingense Using Oxford Nanopore Long Reads and the Diversity and Evolution of Asteraceae Mitogenomes

Diversity in structure and organization is one of the main features of angiosperm mitochondrial genomes (mitogenomes). The ultra-long reads of Oxford Nanopore Technology (ONT) provide an opportunity to obtain a complete mitogenome and investigate the structural variation in unprecedented detail. In this study, we compared mitogenome assembly methods using Illumina and/or ONT sequencing data and obtained the complete mitogenome (208 kb) of Chrysanthemum nankingense based on the hybrid assembly method. The mitogenome encoded 19 transfer RNA genes, three ribosomal RNA genes, and 34 protein-coding genes with 21 group II introns disrupting eight intron-contained genes. A total of seven medium repeats were related to homologous recombination at different frequencies as supported by the long ONT reads. Subsequently, we investigated the variations in gene content and constitution of 28 near-complete mitogenomes from Asteraceae. A total of six protein-coding genes were missing in all Asteraceae mitogenomes, while four other genes were not detected in some lineages. The core fragments (~88 kb) of the Asteraceae mitogenomes had a higher GC content (~46.7%) than the variable and specific fragments. The phylogenetic topology based on the core fragments of the Asteraceae mitogenomes was highly consistent with the topologies obtained from the corresponding plastid datasets. Our results highlighted the advantages of the complete assembly of the C. nankingense mitogenome and the investigation of its structural variation based on ONT sequencing data. Moreover, the method based on local collinear blocks of the mitogenomes could achieve the alignment of highly rearrangeable and variable plant mitogenomes as well as construct a robust phylogenetic topology.

Comparative Analyses of Chloroplast Genomes of Cucurbitaceae Species: Lights into Selective Pressures and Phylogenetic Relationships

Cucurbitaceae is the fourth most important economic plant family with creeping herbaceous species mainly distributed in tropical and subtropical regions. Here, we described and compared the complete chloroplast genome sequences of ten representative species from Cucurbitaceae. The lengths of the ten complete chloroplast genomes ranged from 155,293 bp (C. sativus) to 158,844 bp (M. charantia), and they shared the most common genomic features. 618 repeats of three categories and 813 microsatellites were found. Sequence divergence analysis showed that the coding and IR regions were highly conserved. Three protein-coding genes (accD, clpP, and matK) were under selection and their coding proteins often have functions in chloroplast protein synthesis, gene transcription, energy transformation, and plant development. An unconventional translation initiation codon of psbL gene was found and provided evidence for RNA editing. Applying BI and ML methods, phylogenetic analysis strongly supported the position of Gomphogyne, Hemsleya, and Gynostemma as the relatively original lineage in Cucurbitaceae. This study suggested that the complete chloroplast genome sequences were useful for phylogenetic studies. It would also determine potential molecular markers and candidate DNA barcodes for coming studies and enrich the valuable complete chloroplast genome resources of Cucurbitaceae.

Neo-Domestication of an Interspecific Tetraploid Helianthus annuus × Helianthus tuberous Population That Segregates for Perennial Habit

Perennial agriculture has been proposed as an option to improve the sustainability of cropping systems, by increasing the efficiency of resource use, while also providing ecosystem services. Neo-domestication, the contemporary domestication of plants that have not previously been used in agriculture, can be used to generate new crops for these systems. Here we explore the potential of a tetraploid (2n = 4x = 68) interspecific hybrid sunflower as a perennial oilseed for use in multifunctional agricultural systems. A population of this novel tetraploid was obtained from crosses between the annual diploid oilseed crop Helianthus annuus (2n = 2x = 34) and the perennial hexaploid tuber crop Helianthus tuberosus (2n = 6x = 102). We selected for classic domestication syndrome traits for three generations. Substantial phenotypic gains were made, in some cases approaching 320%. We also analyzed the genetic basis of tuber production (i.e., perenniality), with the goal of obtaining molecular markers that could be used to facilitate future breeding in this system. Results from quantitative trait locus (QTL) mapping suggest that tuber production has an oligogenic genetic basis. Overall, this study indicates that substantial gains towards domestication goals can be achieved over contemporary time scales.

Sequencing and Analysis of Chrysanthemum carinatum Schousb and Kalimeris indica. The Complete Chloroplast Genomes Reveal Two Inversions and rbcL as Barcoding of the Vegetable

Chrysanthemum carinatum Schousb and Kalimeris indica are widely distributed edible vegetables and the sources of the Chinese medicine Asteraceae. The complete chloroplast (cp) genome of Asteraceae usually occurs in the inversions of two regions. Hence, the cp genome sequences and structures of Asteraceae species are crucial for the cp genome genetic diversity and evolutionary studies. Hence, in this paper, we have sequenced and analyzed for the first time the cp genome size of C. carinatum Schousb and K. indica, which are 149,752 bp and 152,885 bp, with a pair of inverted repeats (IRs) (24,523 bp and 25,003) separated by a large single copy (LSC) region (82,290 bp and 84,610) and a small single copy (SSC) region (18,416 bp and 18,269), respectively. In total, 79 protein-coding genes, 30 distinct transfer RNA (tRNA) genes, four distinct rRNA genes and two pseudogenes were found not only in C. carinatum Schousb but also in the K. indica cp genome. Fifty-two (52) and fifty-nine (59) repeats, and seventy (70) and ninety (90) simple sequence repeats (SSRs) were found in the C. carinatum Schousb and K. indica cp genomes, respectively. Codon usage analysis showed that leucine, isoleucine, and serine are the most frequent amino acids and that the UAA stop codon was the significantly favorite stop codon in both cp genomes. The two inversions, the LSC region ranging from trnC-GCA to trnG-UCC and the whole SSC region were found in both of them. The complete cp genome comparison with other Asteraceae species showed that the coding area is more conservative than the non-coding area. The phylogenetic analysis revealed that the rbcL gene is a good barcoding marker for identifying different vegetables. These results give an insight into the identification, the barcoding, and the understanding of the evolutionary model of the Asteraceae cp genome.

Employing hypothesis testing and data from multiple genomic compartments to resolve recalcitrant backbone nodes in Goodenia s.l. (Goodeniaceae)

Goodeniaceae is a primarily Australian flowering plant family with a complex taxonomy and evolutionary history. Previous phylogenetic analyses have successfully resolved the backbone topology of the largest clade in the family, Goodenia s.l., but have failed to clarify relationships within the species-rich and enigmatic Goodenia clade C, a prerequisite for taxonomic revision of the group. We used genome skimming to retrieve sequences for chloroplast, mitochondrial, and nuclear markers for 24 taxa representing Goodenia s.l., with a particular focus on Goodenia clade C. We performed extensive hypothesis tests to explore incongruence in clade C and evaluate statistical support for clades within this group, using datasets from all three genomic compartments. The mitochondrial dataset is comparable to the chloroplast dataset in providing resolution within Goodenia clade C, though backbone support values within this clade remain low. The hypothesis tests provided an additional, complementary means of evaluating support for clades. We propose that the major subclades of Goodenia clade C (C1-C3 + Verreauxia) are the result of a rapid radiation, and each represents a distinct lineage.