Comparative transcriptome analysis reveals candidate genes for cold stress response and early flowering in pineapple.

Pineapple originates from tropical regions in South America and is therefore significantly impacted by cold stress. Periodic cold events in the equatorial regions where pineapple is grown may induce early flowering, also known as precocious flowering, resulting in monetary losses due to small fruit size and the need to make multiple passes for harvesting a single field. Currently, pineapple is one of the most important tropical fruits in the world in terms of consumption, and production losses caused by weather can have major impacts on worldwide exportation potential and economics. To further our understanding of and identify mechanisms for low-temperature tolerance in pineapple, and to identify the relationship between low-temperature stress and flowering time, we report here a transcriptomic analysis of two pineapple genotypes in response to low-temperature stress. Using meristem tissue collected from precocious flowering-susceptible MD2 and precocious flowering-tolerant Dole-17, we performed pairwise comparisons and weighted gene co-expression network analysis (WGCNA) to identify cold stress, genotype, and floral organ development-specific modules. Dole-17 had a greater increase in expression of genes that confer cold tolerance. The results suggested that low temperature stress in Dole-17 plants induces transcriptional changes to adapt and maintain homeostasis. Comparative transcriptomic analysis revealed differences in cuticular wax biosynthesis, carbohydrate accumulation, and vernalization-related gene expression between genotypes. Cold stress induced changes in ethylene and abscisic acid-mediated pathways differentially between genotypes, suggesting that MD2 may be more susceptible to hormone-mediated early flowering. The differentially expressed genes and module hub genes identified in this study are potential candidates for engineering cold tolerance in pineapple to develop new varieties capable of maintaining normal reproduction cycles under cold stress. In addition, a total of 461 core genes involved in the development of reproductive tissues in pineapple were also identified in this study. This research provides an important genomic resource for understanding molecular networks underlying cold stress response and how cold stress affects flowering time in pineapple.

An updated phylogeny, biogeography, and PhyloCode-based classification of Cornaceae based on three sets of genomic data.

PREMISE: A major goal of systematic biology is to uncover the evolutionary history of organisms and translate that knowledge into stable classification systems. Here, we integrate three sets of genome-wide data to resolve phylogenetic relationships in Cornaceae (containing only Cornus s.l.), reconstruct the biogeographic history of the clade, and provide a revised classification using the PhyloCode to stabilize names for this taxonomically controversial group. METHODS: We conducted phylogenetic analyses using 312 single-copy nuclear genes and 70 plastid genes from Angiosperms353 Hyb-Seq, plus numerous loci from RAD-Seq. We integrated fossils using morphological data and produced a dated phylogeny for biogeographical analysis. RESULTS: A well-resolved, strongly supported, comprehensive phylogeny was obtained. Biogeographic analyses support an origin and rapid diversification of Cornus into four morphologically distinct major clades in the Northern Hemisphere (with an eastern Asian ancestor) during the late Cretaceous. Dispersal into Africa from eastern Asia likely occurred along the Tethys Seaway during the Paleogene, whereas dispersal into South America likely occurred during the Neogene. Diversification within the northern hemisphere likely involved repeated independent colonization of new areas during the Paleogene and Neogene along the Bering Land Bridge, the North Atlantic Land Bridge, and the Tethys Seaway. Thirteen strongly supported clades were named following rules of the PhyloCode. CONCLUSIONS: Our study provides an example of integrating genomic and morphological data to produce a robust, explicit species phylogeny that includes fossil taxa, which we translate into an updated classification scheme using the PhyloCode to stabilize names.

Discovering variation of secondary metabolite diversity and its relationship with disease resistance in Cornus florida L.

Understanding intraspecific relationships between genetic and functional diversity is a major goal in the field of evolutionary biology and is important for conserving biodiversity. Linking intraspecific molecular patterns of plants to ecological pressures and trait variation remains difficult due to environment-driven plasticity. Next-generation sequencing, untargeted liquid chromatography-mass spectrometry (LC-MS) profiling, and interdisciplinary approaches integrating population genomics, metabolomics, and community ecology permit novel strategies to tackle this problem. We analyzed six natural populations of the disease-threatened Cornus florida L. from distinct ecological regions using genotype-by-sequencing markers and LC-MS-based untargeted metabolite profiling. We tested the hypothesis that higher genetic diversity in C. florida yielded higher chemical diversity and less disease susceptibility (screening hypothesis), and we also determined whether genetically similar subpopulations were similar in chemical composition. Most importantly, we identified metabolites that were associated with candidate loci or were predictive biomarkers of healthy or diseased plants after controlling for environment. Subpopulation clustering patterns based on genetic or chemical distances were largely congruent. While differences in genetic diversity were small among subpopulations, we did observe notable similarities in patterns between subpopulation averages of rarefied-allelic and chemical richness. More specifically, we found that the most abundant compound of a correlated group of putative terpenoid glycosides and derivatives was correlated with tree health when considering chemodiversity. Random forest biomarker and genomewide association tests suggested that this putative iridoid glucoside and other closely associated chemical features were correlated to SNPs under selection.

Characterization of the sequence and expression pattern of LFY homologues from dogwood species (Cornus) with divergent inflorescence architectures.

BACKGROUND AND AIMS: LFY homologues encode transcription factors that regulate the transition from vegetative to reproductive growth in flowering plants and have been shown to control inflorescence patterning in model species. This study investigated the expression patterns of LFY homologues within the diverse inflorescence types (head-like, umbel-like and inflorescences with elongated internodes) in closely related lineages in the dogwood genus (Cornus s.l.). The study sought to determine whether LFY homologues in Cornus species are expressed during floral and inflorescence development and if the pattern of expression is consistent with a function in regulating floral development and inflorescence architectures in the genus. METHODS: Total RNAs were extracted using the CTAB method and the first-strand cDNA was synthesized using the SuperScript III first-strand synthesis system kit (Invitrogen). Expression of CorLFY was investigated by RT-PCR and RNA in situ hybridization. Phylogenetic analyses were conducted using the maximum likelihood methods implemented in RAxML-HPC v7.2.8. KEY RESULTS: cDNA clones of LFY homologues (designated CorLFY) were isolated from six Cornus species bearing different types of inflorescence. CorLFY cDNAs were predicted to encode proteins of approximately 375 amino acids. The detection of CorLFY expression patterns using in situ RNA hybridization demonstrated the expression of CorLFY within the inflorescence meristems, inflorescence branch meristems, floral meristems and developing floral organ primordia. PCR analyses for cDNA libraries derived from reverse transcription of total RNAs showed that CorLFY was also expressed during the late-stage development of flowers and inflorescences, as well as in bracts and developing leaves. Consistent differences in the CorLFY expression patterns were not detected among the distinct inflorescence types. CONCLUSIONS: The results suggest a role for CorLFY genes during floral and inflorescence development in dogwoods. However, the failure to detect expression differences between the inflorescence types in the Cornus species analysed suggests that the evolutionary shift between major inflorescence types in the genus is not controlled by dramatic alterations in the levels of CorLFY gene transcript accumulation. However, due to spatial, temporal and quantitative limitations of the expression data, it cannot be ruled out that subtle differences in the level or location of CorLFY transcripts may underlie the different inflorescence architectures that are observed across these species. Alternatively, differences in CorLFY protein function or the expression or function of other regulators (e.g. TFL1 and UFO homologues) may support the divergent developmental trajectories.

Complete cpDNA genome sequence of Smilax china and phylogenetic placement of Liliales--influences of gene partitions and taxon sampling.

The complete nucleotide sequence of the chloroplast genome (cpDNA) of Smilax china L. (Smilacaceae) is reported. It is the first complete cp genome sequence in Liliales. Genomic analyses were conducted to examine the rate and pattern of cpDNA genome evolution in Smilax relative to other major lineages of monocots. The cpDNA genomic sequences were combined with those available for Lilium to evaluate the phylogenetic position of Liliales and to investigate the influence of taxon sampling, gene sampling, gene function, natural selection, and substitution rate on phylogenetic inference in monocots. Phylogenetic analyses using sequence data of gene groups partitioned according to gene function, selection force, and total substitution rate demonstrated evident impacts of these factors on phylogenetic inference of monocots and the placement of Liliales, suggesting potential evolutionary convergence or adaptation of some cpDNA genes in monocots. Our study also demonstrated that reduced taxon sampling reduced the bootstrap support for the placement of Liliales in the cpDNA phylogenomic analysis. Analyses of sequences of 77 protein genes with some missing data and sequences of 81 genes (all protein genes plus the rRNA genes) support a sister relationship of Liliales to the commelinids-Asparagales clade, consistent with the APG III system. Analyses of 63 cpDNA protein genes for 32 taxa with few missing data, however, support a sister relationship of Liliales (represented by Smilax and Lilium) to Dioscoreales-Pandanales. Topology tests indicated that these two alignments do not significantly differ given any of these three cpDNA genomic sequence data sets. Furthermore, we found no saturation effect of the data, suggesting that the cpDNA genomic sequence data used in the study are appropriate for monocot phylogenetic study and long-branch attraction is unlikely to be the cause to explain the result of two well-supported, conflict placements of Liliales. Further analyses using sufficient nuclear data remain necessary to evaluate these two phylogenetic hypotheses regarding the position of Liliales and to address the causes of signal conflict among genes and partitions.