Characterizing Floral Symmetry in the Core Goodeniaceae with Geometric Morphometrics

Shape asymmetry - what's new?

Studies of shape asymmetry have become increasingly abundant as the methods of geometric morphometrics have gained widespread use. Most of these studies have focussed on fluctuating asymmetry and have largely obtained similar results as more traditional analyses of asymmetry in distance measurements, but several notable differences have also emerged. A key difference is that shape analyses provide information on the patterns, not just the amount of variation, and therefore tend to be more sensitive. Such analyses have shown that apparently symmetric structures in animals consistently show directional asymmetry for shape, but not for size. Furthermore, the long-standing prediction that phenotypic plasticity in response to environmental heterogeneity can contribute to fluctuating asymmetry has been confirmed for the first time for the shape of flower parts (but not for size). Finally, shape analyses in structures with complex symmetry, such as many flowers, can distinguish multiple types of directional asymmetry, generated by distinct direction-giving factors, which combine to the single component observable in bilaterally symmetric structures. While analyses of shape asymmetry are broadly compatible with traditional analyses of asymmetry, they incorporate more detailed morphological information, particularly for structures with complex symmetry, and therefore can reveal subtle biological effects that would otherwise not be apparent. This makes them a promising tool for a wide range of studies in the basic and applied life sciences.

Directional asymmetry and direction-giving factors: Lessons from flowers with complex symmetry

Directional asymmetry is a systematic difference between the left and right sides for structures with bilateral symmetry or a systematic differentiation among repeated parts for complex symmetry. This study explores factors that produce directional asymmetry in the flower of Iris pumila, a structure with complex symmetry that makes it possible to investigate multiple such factors simultaneously. The shapes and sizes of three types of floral organs, the falls, standards, and style branches, were quantified using the methods of geometric morphometrics. For each flower, this study recorded the compass orientations of floral organs as well as their anatomical orientations relative to the two spathes subtending each flower. To characterize directional asymmetry at the whole-flower level, differences in the average sizes and shapes according to compass orientation and relative orientation were computed, and the left-right asymmetry was also evaluated for each individual organ. No size or shape differences within flowers were found in relation to anatomical position; this may relate to the terminal position of flowers in Iris pumila, suggesting that there may be no adaxial-abaxial polarity, which is very prominent in many other taxa. There was clear directional asymmetry of shape in relation to compass orientation, presumably driven by a consistent environmental gradient such as solar irradiance. There was also clear directional asymmetry between left and right halves of every floral organ, most likely related to the arrangement of organs in the bud. These findings indicate that different factors are acting to produce directional asymmetry at different levels. In conventional analyses not recording flower orientations, these effects would be impossible to disentangle from each other and would probably be included as part of fluctuating asymmetry.

The symmetry spectrum in a hybridising, tropical group of rhododendrons

Many diverse plant clades possess bilaterally symmetrical flowers and specialised pollination syndromes, suggesting that these traits may promote diversification. We examined the evolution of diverse floral morphologies in a species-rich tropical radiation of Rhododendron. We used restriction-site associated DNA sequencing on 114 taxa from Rhododendron sect. Schistanthe to reconstruct phylogenetic relationships and examine hybridisation. We then captured and quantified floral variation using geometric morphometric analyses, which we interpreted in a phylogenetic context. We uncovered phylogenetic conflict and uncertainty caused by introgression within and between clades. Morphometric analyses revealed flower symmetry to be a morphological continuum without clear transitions between radial and bilateral symmetry. Tropical Rhododendron species that began diversifying into New Guinea c. 6 million years ago expanded into novel floral morphological space. Our results showed that the evolution of tropical Rhododendron is characterised by recent speciation, recurrent hybridisation and the origin of floral novelty. Floral variation evolved via changes to multiple components of the corolla that are only recognised in geometric morphometrics with both front and side views of flowers.

The concluding chapter: recircumscription of Goodenia (Goodeniaceae) to include four allied genera with an updated infrageneric classification

Close scrutiny of Goodenia (Goodeniaceae) and allied genera in the 'Core Goodeniaceae' over recent years has clarified our understanding of this captivating group. While expanded sampling, sequencing of multiple regions, and a genome skimming reinforced backbone clearly supported Goodenia s.l. as monophyletic and distinct from Scaevola and Coopernookia, there appears to be no synapomorphic characters that uniquely characterise this morphologically diverse clade. Within Goodenia s.l., there is strong support from nuclear, chloroplast and mitochondrial data for three major clades (Goodenia Clades A, B and C) and various subclades, which lead to earlier suggestions for the possible recognition of these as distinct genera. Through ongoing work, it has become evident that this is impractical, as conflict remains within the most recently diverged Clade C, likely due to recent radiation and incomplete lineage sorting. In light of this, it is proposed that a combination of morphological characters is used to circumscribe an expanded Goodenia that now includes Velleia, Verreauxia, Selliera and Pentaptilon, and an updated infrageneric classification is proposed to accommodate monophyletic subclades. A total of twenty-five new combinations, three reinstatements, and seven new names are published herein including Goodenia subg. Monochila sect. Monochila subsect. Infracta K.A.Sheph. subsect. nov. Also, a type is designated for Goodenia subg. Porphyranthus sect. Ebracteolatae (K.Krause) K.A.Sheph. comb. et stat. nov., and lectotypes or secondstep lectotypes are designated for a further three names.

A Binary-Based Matrix Model for Malus Corolla Symmetry and Its Variational Significance

Floral symmetry (corolla symmetry) has important biological significance in plant genetics and evolution. However, it is often multi-dimensional and difficult to quantify. Here, we constructed a multi-dimensional data matrix [X Y Z] by extracting three qualitative variables with binary properties (X: corolla regularity of interval and coplanarity; Y: petal regularity of shape and size; Z: petal local regularity of curling and wrinkle) from different dimensions of petals (overall to individual, and then to the local): all petals (corolla), individual petals, and local areas of petals. To quantitatively express the degree of Malus corolla symmetry, these variables were then combined with weight assignments (X: 2² > Y: 2¹ > Z: 2⁰) based on their contributions to the corolla symmetry and the algorithm rule of converting binary to decimal values, which facilitated the unification of qualitative and quantitative analyses. Our results revealed significant reduction in degrees of Malus corolla symmetry along the direction of local to overall. Species showed higher degree of corolla symmetry than cultivars; however, taxa with stronger corolla symmetry might not necessarily be species. These findings provide new insights into the circumscription of Malus controversial species. The matrix model should be reference for future evaluation of angiosperm flower symmetry (lack of corolla fusion).

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.

FlowerMorphology: fully automatic flower morphometry software

MAIN CONCLUSION

The software FlowerMorphology is designed for automatic morphometry of actinomorphic flowers. The novel complex parameters of flowers calculated by FlowerMorphology allowed us to quantitatively characterize a polyploid series of tobacco. Morphological differences of plants representing closely related lineages or mutants are mostly quantitative. Very often, there are only very fine variations in plant morphology. Therefore, accurate and high-throughput methods are needed for their quantification. In addition, new characteristics are necessary for reliable detection of subtle changes in morphology. FlowerMorphology is an all-in-one software package to automatically image and analyze five-petal actinomorphic flowers of the dicotyledonous plants. Sixteen directly measured parameters and ten calculated complex parameters of a flower allow us to characterize variations with high accuracy. The program was developed for the needs of automatic characterization of Nicotiana tabacum flowers, but is applicable to many other plants with five-petal actinomorphic flowers and can be adopted for flowers of other merosity. A genetically similar polyploid series of N. tabacum plants was used to investigate differences in flower morphology. For the first time, we could quantify the dependence between ploidy and size and form of the tobacco flowers. We found that the radius of inner petal incisions shows a persistent positive correlation with the chromosome number. In contrast, a commonly used parameter-radius of outer corolla-does not discriminate 2n and 4n plants. Other parameters show that polyploidy leads to significant aberrations in flower symmetry and are also positively correlated with chromosome number. Executables of FlowerMorphology, source code, documentation, and examples are available at the program website: https://github.com/Deyneko/FlowerMorphology .

The evolution of floral ontogenetic allometry in the Andean genus Caiophora (Loasaceae, subfam. Loasoideae)

The astounding variety of angiosperm flower morphologies has evolved in response to many selective forces. Flower development is highly coordinated and involves developmental associations between size and shape, ontogenetic allometry, which in turn affect the morphology of mature flowers. Although ontogenetic allometries can act as a developmental constraint and may influence adaptive evolution, allometries can evolve themselves and may change rapidly in response to selection. We explored the evolution of ontogenetic allometry in the flowers of 11 species of Loasoideae. Seven species belong to Caiophora, which radiated recently in the central Andes, and contains species that are pollinated by bees, hummingbirds, and small rodents. According to a previous study, the diversification of Caiophora involved departures from simple allometric scaling, but the changes to allometry that enabled flower diversification have not been explored yet. We characterized the ontogenetic allometry of each species with the methods of geometric morphometrics. We studied the evolution of allometries by constructing allometric spaces, in which the allometry of each species is represented by a point and the arrangement of points indicates the relations among allometric trajectories. To examine the history of changes of ontogenetic allometries, we projected the phylogeny into the allometric spaces. Inspection of allometric spaces suggests that ontogenetic variation is limited to a few dominant features. The allometries of the two main functional flower parts under study differ in their evolutionary labilities, and patterns of variation reflect pollination systems, differences in structural organization, and abiotic environmental factors.

Phenotypic plasticity in response to environmental heterogeneity contributes to fluctuating asymmetry in plants: first empirical evidence

Fluctuating asymmetry (FA) is widely used to quantify developmental instability (DI) in ecological and evolutionary studies. It has long been recognized that FA may not exclusively originate from DI for sessile organisms such as plants, because phenotypic plasticity in response to heterogeneities in the environment might also produce FA. This study provides the first empirical evidence for this hypothesis. We reasoned that solar irradiance, which is greater on the southern side than on the northern side of plants growing in the temperate zone of the Northern Hemisphere, would cause systematic morphological differences and asymmetry associated with the orientation of plant parts. We used geometric morphometrics to characterize the size and shape of flower parts in Iris pumila grown in a common garden. The size of floral organs was not significantly affected by orientation. Shape and particularly its asymmetric component differed significantly according to orientation for three different floral parts. Orientation accounted for 10.4% of the total shape asymmetry within flowers in the falls, for 11.4% in the standards and for 2.2% in the style branches. This indicates that phenotypic plasticity in response to a directed environmental factor, most likely solar irradiance, contributes to FA of flowers under natural conditions. That FA partly results from phenotypic plasticity and not just from DI needs to be considered by studies of FA in plants and other sessile organisms.

The unexpected depths of genome-skimming data: A case study examining Goodeniaceae floral symmetry genes

PREMISE OF THE STUDY

The use of genome skimming allows systematists to quickly generate large data sets, particularly of sequences in high abundance (e.g., plastomes); however, researchers may be overlooking data in low abundance that could be used for phylogenetic or evo-devo studies. Here, we present a bioinformatics approach that explores the low-abundance portion of genome-skimming next-generation sequencing libraries in the fan-flowered Goodeniaceae.

METHODS

Twenty-four previously constructed Goodeniaceae genome-skimming Illumina libraries were examined for their utility in mining low-copy nuclear genes involved in floral symmetry, specifically the CYCLOIDEA (CYC)-like genes. De novo assemblies were generated using multiple assemblers, and BLAST searches were performed for CYC1, CYC2, and CYC3 genes.

RESULTS

Overall Trinity, SOAPdenovo-Trans, and SOAPdenovo implementing lower k-mer values uncovered the most data, although no assembler consistently outperformed the others. Using SOAPdenovo-Trans across all 24 data sets, we recovered four CYC-like gene groups (CYC1, CYC2, CYC3A, and CYC3B) from a majority of the species. Alignments of the fragments included the entire coding sequence as well as upstream and downstream regions.

DISCUSSION

Genome-skimming data sets can provide a significant source of low-copy nuclear gene sequence data that may be used for multiple downstream applications.

The unexpected depths of genome-skimming data: A case study examining Goodeniaceae floral symmetry genes

PREMISE OF THE STUDY

The use of genome skimming allows systematists to quickly generate large data sets, particularly of sequences in high abundance (e.g., plastomes); however, researchers may be overlooking data in low abundance that could be used for phylogenetic or evo-devo studies. Here, we present a bioinformatics approach that explores the low-abundance portion of genome-skimming next-generation sequencing libraries in the fan-flowered Goodeniaceae.

METHODS

Twenty-four previously constructed Goodeniaceae genome-skimming Illumina libraries were examined for their utility in mining low-copy nuclear genes involved in floral symmetry, specifically the CYCLOIDEA (CYC)-like genes. De novo assemblies were generated using multiple assemblers, and BLAST searches were performed for CYC1, CYC2, and CYC3 genes.

RESULTS

Overall Trinity, SOAPdenovo-Trans, and SOAPdenovo implementing lower k-mer values uncovered the most data, although no assembler consistently outperformed the others. Using SOAPdenovo-Trans across all 24 data sets, we recovered four CYC-like gene groups (CYC1, CYC2, CYC3A, and CYC3B) from a majority of the species. Alignments of the fragments included the entire coding sequence as well as upstream and downstream regions.

DISCUSSION

Genome-skimming data sets can provide a significant source of low-copy nuclear gene sequence data that may be used for multiple downstream applications.

A potential pitfall in studies of biological shape: Does size matter?

The number of published studies using geometric morphometrics (GM) for analysing biological shape has increased steadily since the beginning of the 1990s, covering multiple research areas such as ecology, evolution, development, taxonomy and palaeontology. Unfortunately, we have observed that many published studies using GM do not evaluate the potential allometric effects of size on shape, which normally require consideration or assessment. This might lead to misinterpretations and flawed conclusions in certain cases, especially when size effects explain a large part of the shape variation. We assessed, for the first time and in a systematic manner, how often published studies that have applied GM consider the potential effects of allometry on shape. We reviewed the 300 most recent published papers that used GM for studying biological shape. We also estimated how much of the shape variation was explained by allometric effects in the reviewed papers. More than one-third (38%) of the reviewed studies did not consider the allometric component of shape variation. In studies where the allometric component was taken into account, it was significant in 88% of the cases, explaining up to 87.3% of total shape variation. We believe that one reason that may cause the observed results is a misunderstanding of the process that superimposes landmark configurations, i.e. the Generalized Procrustes Analysis, which removes isometric effects of size on shape, but not allometric effects. Allometry can be a crucial component of shape variation. We urge authors to address, and report, size effects in studies of biological shape. However, we do not propose to always remove size effects, but rather to evaluate the research question with and without the allometric component of shape variation. This approach can certainly provide a thorough understanding of how much size contributes to the observed shaped variation.

A robust, high-throughput method for computing maize ear, cob, and kernel attributes automatically from images

Grain yield of the maize plant depends on the sizes, shapes, and numbers of ears and the kernels they bear. An automated pipeline that can measure these components of yield from easily-obtained digital images is needed to advance our understanding of this globally important crop. Here we present three custom algorithms designed to compute such yield components automatically from digital images acquired by a low-cost platform. One algorithm determines the average space each kernel occupies along the cob axis using a sliding-window Fourier transform analysis of image intensity features. A second counts individual kernels removed from ears, including those in clusters. A third measures each kernel's major and minor axis after a Bayesian analysis of contour points identifies the kernel tip. Dimensionless ear and kernel shape traits that may interrelate yield components are measured by principal components analysis of contour point sets. Increased objectivity and speed compared to typical manual methods are achieved without loss of accuracy as evidenced by high correlations with ground truth measurements and simulated data. Millimeter-scale differences among ear, cob, and kernel traits that ranged more than 2.5-fold across a diverse group of inbred maize lines were resolved. This system for measuring maize ear, cob, and kernel attributes is being used by multiple research groups as an automated Web service running on community high-throughput computing and distributed data storage infrastructure. Users may create their own workflow using the source code that is staged for download on a public repository.