Diploid Species Research Articles

Polyploid Nicotiana section Suaveolentes originated by hybridization of two ancestral Nicotiana clades.

Nicotiana section Suaveolentes is an almost all-Australian clade of allopolyploid tobacco species that emerged through hybridization between diploid relatives of the genus. In this study, we aimed to assess the phylogenetic relationship of the Suaveolentes section with several Nicotiana diploid species based on both plastidial and nuclear genes. The Nicotiana plastome-based phylogenetic analysis representing 47 newly re-built plastid genomes suggested that an ancestor of N. section Noctiflorae is the most likely maternal donor of the Suaveolentes clade. Nevertheless, we found clear evidence of plastid recombination with an ancestor from the Sylvestres clade. We analyzed 411 maximum likelihood-based phylogenetic trees from a set of conserved nuclear diploid single copy gene families following an approach that assessed the genomic origin of each homeolog. We found that Nicotiana section Suaveolentes is monophyletic with contributions from the sections Alatae, Sylvestres, Petunioides and Noctiflorae. The dating of the divergence between these sections indicates that the Suaveolentes hybridization predates the split between Alatae/Sylvestres, and Noctiflorae/Petunioides. We propose that Nicotiana section Suaveolentes arose from the hybridization of two ancestral species from which the Noctiflorae/Petunioides and Alatae/Sylvestres sections are derived, with Noctiflorae the maternal parent. This study is a good example in which the use of genome wide data provided additional evidence about the origin of a complex polyploid clade.

Transcript profiling of hairy vetch (Vicia villosa Roth) identified interesting genes for seed dormancy.

Hairy vetch, a diploid annual legume species, has a robust growth habit, high biomass yield, and winter hardy characteristics. Seed hardness is a major constraint for growing hairy vetch commercially. Hard seeded cultivars are valuable as forages, whereas soft seeded and shatter resistant cultivars have advantages for their use as a cover crop. Transcript analysis of hairy vetch was performed to understand the genetic mechanisms associated with important hairy vetch traits. RNA was extracted from leaves, flowers, immature pods, seed coats, and cotyledons of contrasting soft and hard seeded "AU Merit" plants. A range of 31.22-79.18Gb RNA sequence data per tissue sample were generated with estimated coverage of 1040-2639×. RNA sequence assembly and mapping of the contigs against the Medicago truncatula (V4.0) genome identified 76,422 gene transcripts. A total of 24,254 transcripts were constitutively expressed in hairy vetch tissues. Key genes, such as KNOX4 (a class II KNOTTED-like homeobox KNOXII gene), qHs1 (endo-1,4-β-glucanase), GmHs1-1 (calcineurin-like metallophosphoesterase), chitinase, shatterproof 1 and 2 (SHP1, SHP2), shatter resistant 1-5 (SHAT1-5)(NAC transcription factor), PDH1 (prephenate dehydrogenase 1), and pectin methylesterases with a potential role in seed hardness and pod shattering, were further explored based on genes involved in seed hardness from other species to query the hairy vetch transcriptome data. Identification of interesting candidate genes in hairy vetch can facilitate the development of improved cultivars with desirable seed characteristics for use as a forage and as a cover crop.

First report of chromosome and karyological analysis of Gekko nutaphandi (Gekkonidae, Squamata) from Thailand: Neo-diploid chromosome number in genus Gekko

The karyotypes of red-eyed Gecko are not reported yet. Herein, we describe the karyotypes of red-eyed Gecko (Gekko nutaphandi Bauer, Sumontha & Pauwels, 2008) from Thailand. Gecko chromosome preparations were directly conducted from bone marrow and testis. Chromosomal characteristics were analyzed by Giemsa staining, Ag-NOR banding as well as fluorescence in situ hybridization (FISH) using microsatellites d(GC)15 probe. The results showed that the number of diploid chromosomes is 2n=34, while the fundamental number (NF) is 46 in both males and females. The types of chromosomes were 4 large metacentric, 6 large submetacentric, 2 medium telocentric, 2 small metacentric and 20 small telocentric chromosomes. The results of conventional Giemsa staining presented the diploid chromosome number differentiation even in the same genus. NORs are located at the secondary constriction to the telomere on the long arm of chromosome pair 5. There are no sex differences in karyotypes between males and females. FISH with d(GC)15 sequences were also displayed at the telomeres of most other chromosomes. We found that during metaphase I the homologous chromosomes showed synapsis, which can be defined as 19 ring bivalents and 17 haploid chromosomes (n=17) at metaphase II as a diploid species. The karyotype formula is as follows: 2n (34) = L4m+L6sm+M2t+S2m+S20t.

Genome-Wide Comparative Analysis of the R2R3-MYB Gene Family in Six Ipomoea Species and the Identification of Anthocyanin-Related Members in Sweet Potatoes.

Sweet potatoes (Ipomoea batatas) are one of the important tuberous root crops cultivated worldwide, and thier storage roots are rich in antioxidants, such as anthocyanins. R2R3-MYB is a large gene family involved in various biological processes, including anthocyanin biosynthesis. However, few reports about the R2R3-MYB gene family of sweet potatoes have been released to date. In the present study, a total of 695 typical R2R3-MYB genes were identified in six Ipomoea species, including 131 R2R3-MYB genes in sweet potatoes. A maximum likelihood phylogenetic analysis divided these genes into 36 clades, referring to the classification of 126 R2R3-MYB proteins of Arabidopsis. Clade C25(S12) has no members in six Ipomoea species, whereas four clades (i.e., clade C21, C26, C30, and C36), including 102 members, had no members in Arabidopsis, and they were identified as Ipomoea-specific clades. The identified R2R3-MYB genes were unevenly distributed on all chromosomes in six Ipomoea species genomes, and the collinearity analysis among hexaploid I. batatas and another five diploid Ipomoea species suggested that the sweet potato genome might have undergone a larger chromosome rearrangement during the evolution process. Further analyses of gene duplication events showed that whole-genome duplication, transposed duplication, and dispersed duplication events were the primary forces driving the R2R3-MYB gene family expansion of Ipomoea plants, and these duplicated genes experienced strong purifying selection because of their Ka/Ks ratio, which is less than 1. Additionally, the genomic sequence length of 131 IbR2R3-MYBs varied from 923 bp to ~12.9 kb with a mean of ~2.6 kb, and most of them had more than three exons. The Motif 1, 2, 3, and 4 formed typical R2 and R3 domains and were identified in all IbR2R3-MYB proteins. Finally, based on multiple RNA-seq datasets, two IbR2R3-MYB genes (IbMYB1/g17138.t1 and IbMYB113/g17108.t1) were relatively highly expressed in pigmented leaves and tuberous root flesh and skin, respectively; thus, they were identified to regulate tissue-specific anthocyanin accumulation in sweet potato. This study provides a basis for the evolution and function of the R2R3-MYB gene family in sweet potatoes and five other Ipomoea species.

Quantifying the fraction of new mutations that are recessive lethal.

The presence and impact of recessive lethal mutations have been widely documented in diploid outcrossing species. However, precise estimates of the proportion of new mutations that are recessive lethal remain limited. Here, we evaluate the performance of Fit∂a∂i, a commonly used method for inferring the distribution of fitness effects (DFE), in the presence of lethal mutations. Using simulations, we demonstrate that in both additive and recessive cases, inference of the deleterious nonlethal portion of the DFE is minimally affected by a small proportion (<10%) of lethal mutations. Additionally, we demonstrate that while Fit∂a∂i cannot estimate the fraction of recessive lethal mutations, Fit∂a∂i can accurately infer the fraction of additive lethal mutations. Finally, as an alternative approach to estimate the proportion of mutations that are recessive lethal, we employ models of mutation-selection-drift balance using existing genomic parameters and estimates of segregating recessive lethals for humans and Drosophila melanogaster. In both species, the segregating recessive lethal load can be explained by a very small fraction (<1%) of new nonsynonymous mutations being recessive lethal. Our results refute recent assertions of a much higher proportion of mutations being recessive lethal (4%-5%), while highlighting the need for additional information on the joint distribution of selection and dominance coefficients.

Genome-wide characterization of major latex protein gene family in peanut and expression analyses under drought and waterlogging stress.

Peanut is an important oilseed crop around the world which provides vegetable oil, protein and vitamins for humans. Major latex-like proteins (MLPs) play important roles in plant growth and development, as well as responses to biotic and abiotic stresses. However, their biological function in peanut is still unclear. In this study, a genome-wide identification of MLP genes in cultivated peanut and two diploid ancestor species was analyzed to determine their molecular evolutionary characteristics and the expression profile under drought and waterlogging stress conditions. Firstly, a total of 135 MLP genes were identified from the genome of tetraploid peanut (Arachis hypogaea) and two diploid species Arachis. duranensis and Arachis. ipaensis. Then, phylogenetic analysis revealed that MLP proteins were divided into five different evolutionary groups. These genes were distributed unevenly at the ends of chromosomes 3, 5, 7, 8, 9 and 10 in three Arachis species. The evolution of MLP gene family in peanut was conserved and led by tandem and segmental duplication. The prediction analysis of cis-acting elements showed that the promoter region of peanut MLP genes contained different proportions of transcription factors, plant hormones-responsive elements and so on. The expression pattern analysis showed that they were differentially expressed under waterlogging and drought stress. These results of this study provide a foundation for further research on the function of the important MLP genes in peanut.

Investigation of NLR Genes Reveals Divergent Evolution on NLRome in Diploid and Polyploid Species in Genus Trifolium.

Crop wild relatives contain a greater variety of phenotypic and genotypic diversity compared to their domesticated counterparts. Trifolium crop species have limited genetic diversity to cope with biotic and abiotic stresses due to artificial selection for consumer preferences. Here, we investigated the distribution and evolution of nucleotide-binding site leucine-rich repeat receptor (NLR) genes in the genus of Trifolium with the objective to identify reference NLR genes. We identified 412, 350, 306, 389 and 241 NLR genes were identified from Trifolium. subterraneum, T. pratense, T. occidentale, subgenome-A of T. repens and subgenome-B of T. repens, respectively. Phylogenetic and clustering analysis reveals seven sub-groups in genus Trifolium. Specific subgroups such as G4-CNL, CCG10-CNL and TIR-CNL show distinct duplication patterns in specific species, which suggests subgroup duplications that are the hallmarks of their divergent evolution. Furthermore, our results strongly suggest the overall expansion of NLR repertoire in T. subterraneum is due to gene duplication events and birth of gene families after speciation. Moreover, the NLRome of the allopolyploid species T. repens has evolved asymmetrically, with the subgenome -A showing expansion, while the subgenome-B underwent contraction. These findings provide crucial background data for comprehending NLR evolution in the Fabaceae family and offer a more comprehensive analysis of NLR genes as disease resistance genes.

The morphological diversity of pollen in the genus Gossypium

BackgroundPlant pollen has diverse morphological characteristics that can be consistently passed down from generation to generation. Information on pollen morphology is thus immensely important for plant classification and identification. In the genus Gossypium, however, in-depth research on pollen morphology is lacking, with only few reports on limited cotton species. To evaluate the diversity of pollen in Gossypium, we therefore conducted a comprehensive analysis of the pollen morphology of 33 cotton species and varieties using scanning electron microscopy.ResultsThe 33 analyzed cotton samples exhibited common pollen morphological features, including spherical shapes, radial symmetry, echination, panporation, and operculation, while the pollen size, spine shape, spine density and length showed distinctive features. Pollen size varied significantly among species, with diameters ranging from 62.43 μm in G. harknessii to 103.41 μm in G. barbadense. The exine had an echinate sculptural texture, and spines were mostly conical or sharply conical but occasionally rod-like. Spine density varied from 173 in G. incanum to 54 in G. gossypioides, while spine length ranged from 3.53 μm in G. herbaceum to 9.47 μm in G. barbadense. In addition, the 33 cotton species and varieties were grouped at a genetic distance of 3.83 into three clusters. Cluster I comprised five allotetraploid AD-genome cotton species, four D-genome species, and one K-genome species. Cluster II included 13 diploid species from A, B, D, E, and G genomes, whereas Cluster III only consisted one E-genome species G. incanum.ConclusionsAlthough pollen characteristics alone are not enough to resolve taxonomic and systematic relationships within the genus Gossypium, our results add to knowledge on palynomorphology and contribute to phenological information on these taxa. Our findings should aid future systematic and phylogenetic studies of the Gossypium genus.

A complete gap-free diploid genome in Saccharum complex and the genomic footprints of evolution in the highly polyploid Saccharum genus.

A diploid genome in the Saccharum complex facilitates our understanding of evolution in the highly polyploid Saccharum genus. Here we have generated a complete, gap-free genome assembly of Erianthus rufipilus, a diploid species within the Saccharum complex. The complete assembly revealed that centromere satellite homogenization was accompanied by the insertions of Gypsy retrotransposons, which drove centromere diversification. An overall low rate of gene transcription was observed in the palaeo-duplicated chromosome EruChr05 similar to other grasses, which might be regulated by methylation patterns mediated by homologous 24 nt small RNAs, and potentially mediating the functions of many nucleotide-binding site genes. Sequencing data for 211 accessions in the Saccharum complex indicated that Saccharum probably originated in the trans-Himalayan region from a diploid ancestor (x = 10) around 1.9-2.5 million years ago. Our study provides new insights into the origin and evolution of Saccharum and accelerates translational research in cereal genetics and genomics.

Characterization of centromeric DNA of Gossypium anomalum reveals sequence-independent enrichment dynamics of centromeric repeats.

Centromeres in eukaryotes are composed of highly repetitive DNAs, which evolve rapidly and are thought to achieve a favorable structure in mature centromeres. However, how the centromeric repeat evolves into an adaptive structure is largely unknown. We characterized the centromeric sequences of Gossypium anomalum through chromatin immunoprecipitation against CENH3 antibodies. We revealed that the G. anomalum centromeres contained only retrotransposon-like repeats but were depleted in long arrays of satellites. These retrotransposon-like centromeric repeats were present in the African-Asian and Australian lineage species, suggesting that they might have arisen in the common ancestor of these diploid species. Intriguingly, we observed a substantial increase and decrease in copy numbers among African-Asian and Australian lineages, respectively, for the retrotransposon-derived centromeric repeats without apparent structure or sequence variation in cotton. This result indicates that the sequence content is not a decisive aspect of the adaptive evolution of centromeric repeats or at least retrotransposon-like centromeric repeats. In addition, two active genes with potential roles in gametogenesis or flowering were identified in CENH3 nucleosome-binding regions. Our results provide new insights into the constitution of centromeric repetitive DNA and the adaptive evolution of centromeric repeats in plants.

Genome-wide identification of CBL family genes in Nicotiana tabacum and the functional analysis of NtCBL4A-1 under salt stress

The CALCINEURIN B-LIKE PROTEIN family genes (CBLs) encode a group of plant-specific calcium sensor proteins that play critical roles in plant development and stress response. Although genome-wide analysis of CBL family genes has been conducted in several diploid plant species, less is known about their functions in allopolyploid plant species. This study conducted a detailed analysis of tobacco CBL family genes. Of the 24 NtCBL genes identified, 14 were derived from the maternal genome donor N. sylvestris, and 10 were from the paternal genome donor N. tomentosiformis. Phylogenetic analysis revealed that NtCBLs could be clustered into three different groups with distinct N-terminal characteristics. All NtCBLs contain EF-hand motifs, but amino acid substitutions and configuration changes may be responsible for variations in their Ca2+-binding abilities. Alternative splicing is common for NtCBL gene transcripts, which may contribute to variations in translation efficiency or protein interaction preferences. Several NtCBL genes exhibited shoot- or root-predominant expression patterns, and some were responsive to salt and/or drought stress. NtCBL4A-1, an ortholog of Arabidopsis thaliana AtCBL4 (SOS3), was predominantly expressed in roots and exhibited a lower gene expression level under salt stress. Further functional analysis showed that overexpression of NtCBL4A-1 increased the sensitivity of transgenic tobacco to Na+-induced salt stress, which is inconsistent with the role of AtCBL4 in contributing to salt tolerance. This study provides data for identifying CBL gene functions and selecting candidate stress-tolerant genes in N. tabacum. Additionally, the results may aid gene function studies in other allopolyploid plant species.

Genome-wide identification of the MIKCc-type MADS-box gene family in Chrysanthemum lavandulifolium reveals their roles in the capitulum development.

Chrysanthemum ×morifolium is well known throughout the world for its diverse and exquisite flower types. However, due to the complicated genetic background of C. ×morifolium, it is difficult to understand the molecular mechanism of its flower development. And it limits the molecular breeding of improving chrysanthemum flower types. C. ×morifolium has the typical radial capitulum, and many researches showed that the members of the MIKCc-type MADS box gene family play a key role in the formation and development of the capitulum. However, it has been difficult to isolate the important MIKCc and investigate their roles in this process due to the lack of genomic information in chrysanthemum. Here, we identified MIKCc-type MADS box genes at whole genome-wide level in C. lavandulifolium, a diploid species closely related to C. ×morifolium, and investigated their roles in capitulum development by gene expression pattern analysis and protein interaction analysis. A total of 40 ClMIKCc were identified and were phylogenetically grouped into 12 clades. Members of all clades showed different enriched expression patterns during capitulum formation. We speculate that the E-class genes in C. lavandulifolium underwent subfunctionalization because they have a significantly expanded, more diverse expression patterns, and specifically tissue expression than AtSEPs. Meanwhile, we detected the C-class expressed in disc floret corolla, which could be the clue to explore the morphological differences between disc and ray floret corolla. In addition, the potential roles of some MIKCcs in complex inflorescence formation were explored by comparing the number and phylogenetic relationship of MIKCc subfamily members in Asteraceae with different capitulum types. Members of the FLC branch in Asteraceae were found to be possibly related to the differentiation and development of the ray floret.

Phylogenetic Analysis of Allotetraploid Species Using Polarized Genomic Sequences.

Phylogenetic analysis of polyploid hybrid species has long posed a formidable challenge as it requires the ability to distinguish between alleles of different ancestral origins in order to disentangle their individual evolutionary history. This problem has been previously addressed by conceiving phylogenies as reticulate networks, using a two-step phasing strategy that first identifies and segregates homoeologous loci and then, during a second phasing step, assigns each gene copy to one of the subgenomes of an allopolyploid species. Here, we propose an alternative approach, one that preserves the core idea behind phasing-to produce separate nucleotide sequences that capture the reticulate evolutionary history of a polyploid-while vastly simplifying its implementation by reducing a complex multistage procedure to a single phasing step. While most current methods used for phylogenetic reconstruction of polyploid species require sequencing reads to be pre-phased using experimental or computational methods-usually an expensive, complex, and/or time-consuming endeavor-phasing executed using our algorithm is performed directly on the multiple-sequence alignment (MSA), a key change that allows for the simultaneous segregation and sorting of gene copies. We introduce the concept of genomic polarization that, when applied to an allopolyploid species, produces nucleotide sequences that capture the fraction of a polyploid genome that deviates from that of a reference sequence, usually one of the other species present in the MSA. We show that if the reference sequence is one of the parental species, the polarized polyploid sequence has a close resemblance (high pairwise sequence identity) to the second parental species. This knowledge is harnessed to build a new heuristic algorithm where, by replacing the allopolyploid genomic sequence in the MSA by its polarized version, it is possible to identify the phylogenetic position of the polyploid's ancestral parents in an iterative process. The proposed methodology can be used with long-read and short-read high-throughput sequencing data and requires only one representative individual for each species to be included in the phylogenetic analysis. In its current form, it can be used in the analysis of phylogenies containing tetraploid and diploid species. We test the newly developed method extensively using simulated data in order to evaluate its accuracy. We show empirically that the use of polarized genomic sequences allows for the correct identification of both parental species of an allotetraploid with up to 97% certainty in phylogenies with moderate levels of incomplete lineage sorting (ILS) and 87% in phylogenies containing high levels of ILS. We then apply the polarization protocol to reconstruct the reticulate histories of Arabidopsis kamchatica and Arabidopsis suecica, two allopolyploids whose ancestry has been well documented. [Allopolyploidy; Arabidopsis; genomic polarization; homoeologs; incomplete lineage sorting; phasing; polyploid phylogenetics; reticulate evolution.].

Endangered spined loach Cobitis species in Southwest Germany: A complex network of native, hybrid and newly emerged populations

AbstractSuccessful conservation of cryptic species such as representatives of the genus Cobitis presents a challenge for fishery managers. Only Cobitis taenia L. is currently assumed to occur in southwest Germany, where it is classified as critically endangered. Established genetic markers and eDNA analysis identified three Cobitis species in the study area. Cobitis taenia and Cobitis elongatoides were widely distributed in the Rhine River system, primarily forming hybrid species complexes in which polyploid individuals dominate but coexist with diploid parent species. A third non‐native species, Cobitis bilineata, has recently become established in the southern part of the Rhine to which it probably migrated naturally via the connected Aare River system. In the Danube River catchment, two populations, including one previously thought to be extinct, comprise both diploid and polyploid representatives of C. elongatoides. Our results provide an important basis for more targeted conservation strategies for Cobitis species currently found in the region.

Intrinsic Allergenicity Potential of Salt-Soluble Protein Extracts from the Diploid, Tetraploid and Hexaploid Wheats: Validation Using an Adjuvant-Free Mouse Model

Wheat allergies are potentially life-threatening and, therefore, have become a major health concern at the global level. It is largely unknown at present whether genetic variation in allergenicity potential exists among hexaploid, tetraploid and diploid wheat species. Such information is critical in establishing a baseline allergenicity map to inform breeding efforts to identify hyper-, hypo- and non-allergenic varieties. We recently reported a novel mouse model of intrinsic allergenicity using the salt-soluble protein extract (SSPE) from durum, a tetraploid wheat (Triticum durum). Here, we validated the model for three other wheat species [hexaploid common wheat (Triticum aestivum), diploid einkorn wheat (Triticum monococcum), and the ancient diploid wheat progenitor, Aegilops tauschii], and then tested the hypothesis that the SSPEs from wheat species will exhibit differences in relative allergenicities. Balb/c mice were repeatedly exposed to SSPEs via the skin. Allergic sensitization potential was assessed by specific (s) IgE antibody responses. Oral anaphylaxis was quantified by the hypothermic shock response (HSR). The mucosal mast cell response (MMCR) was determined by measuring mast cell protease in the blood. While T. monococcum elicited the least, but significant, sensitization, others were comparable. Whereas Ae. taushcii elicited the least HSR, the other three elicited much higher HSRs. Similarly, while Ae. tauschii elicited the least MMCR, the other wheats elicited much higher MMCR as well. In conclusion, this pre-clinical comparative mapping strategy may be used to identify potentially hyper-, hypo- and non-allergenic wheat varieties via crossbreeding and genetic engineering methods.

Evolutionary and Genomic Diversity of True Polyploidy in Tetrapods.

True polyploid organisms have more than two chromosome sets in their somatic and germline cells. Polyploidy is a major evolutionary force and has played a significant role in the early genomic evolution of plants, different invertebrate taxa, chordates, and teleosts. However, the contribution of polyploidy to the generation of new genomic, ecological, and species diversity in tetrapods has traditionally been underestimated. Indeed, polyploidy represents an important pathway of genomic evolution, occurring in most higher-taxa tetrapods and displaying a variety of different forms, genomic configurations, and biological implications. Herein, we report and discuss the available information on the different origins and evolutionary and ecological significance of true polyploidy in tetrapods. Among the main tetrapod lineages, modern amphibians have an unparalleled diversity of polyploids and, until recently, they were considered to be the only vertebrates with closely related diploid and polyploid bisexual species or populations. In reptiles, polyploidy was thought to be restricted to squamates and associated with parthenogenesis. In birds and mammals, true polyploidy has generally been considered absent (non-tolerated). These views are being changed due to an accumulation of new data, and the impact as well as the different evolutionary and ecological implications of polyploidy in tetrapods, deserve a broader evaluation.

Wheat grain quality changes with water stress, zinc, and iron applications predicted by the solvent retention capacity (SRC)

Due to a lack of assessment tools for large samples, the effects of foliar Zn and Fe applications on grain nutritional and quality traits for wheat genotypes from diverse species and poloidal levels under water-limited conditions are less addressed. With the availability of simple/rapid solvent retention capacity (SRC) tests to assess a large number of samples for grain quality attributes, thirty-five wheat genotypes from 18 different species were field evaluated under four treatments of control, zinc (Zn), iron (Fe), and Zn + Fe foliar applications in two moisture regimes. Grain yield (GY) and thousand kernel weight (TKW) decreased due to water stress, while the grain quality traits including SRC values, flour protein content (FPC), Zn, and Fe content increased. Although diploid species showed higher FPC, grain Zn, and Fe than tetraploid (4x) and hexaploid (6x) species across all treatments, they presented much lower SRC values attributed to protein quality. Zinc application increased SRC profiles for 4x and 6x wheat and alleviated GY and TKW reduction due to water stress. A significant correlation between various grain and flour parameters confirmed the validity of the results and SRC quality predictions.

Geometric Morphometric Versus Genomic Patterns in a Large Polyploid Plant Species Complex

Plant species complexes represent a particularly interesting example of taxonomically complex groups (TCGs), linking hybridization, apomixis, and polyploidy with complex morphological patterns. In such TCGs, mosaic-like character combinations and conflicts of morphological data with molecular phylogenies present a major problem for species classification. Here, we used the large polyploid apomictic European Ranunculus auricomus complex to study relationships among five diploid sexual progenitor species and 75 polyploid apomictic derivate taxa, based on geometric morphometrics using 11,690 landmarked objects (basal and stem leaves, receptacles), genomic data (97,312 RAD-Seq loci, 48 phased target enrichment genes, 71 plastid regions) from 220 populations. We showed that (1) observed genomic clusters correspond to morphological groupings based on basal leaves and concatenated traits, and morphological groups were best resolved with RAD-Seq data; (2) described apomictic taxa usually overlap within trait morphospace except for those taxa at the space edges; (3) apomictic phenotypes are highly influenced by parental subgenome composition and to a lesser extent by climatic factors; and (4) allopolyploid apomictic taxa, compared to their sexual progenitor, resemble a mosaic of ecological and morphological intermediate to transgressive biotypes. The joint evaluation of phylogenomic, phenotypic, reproductive, and ecological data supports a revision of purely descriptive, subjective traditional morphological classifications.

Smooth Descent: A ploidy-aware algorithm to improve linkage mapping in the presence of genotyping errors.

Linkage mapping is an approach to order markers based on recombination events. Mapping algorithms cannot easily handle genotyping errors, which are common in high-throughput genotyping data. To solve this issue, strategies have been developed, aimed mostly at identifying and eliminating these errors. One such strategy is SMOOTH, an iterative algorithm to detect genotyping errors. Unlike other approaches, SMOOTH can also be used to impute the most probable alternative genotypes, but its application is limited to diploid species and to markers heterozygous in only one of the parents. In this study we adapted SMOOTH to expand its use to any marker type and to autopolyploids with the use of identity-by-descent probabilities, naming the updated algorithm Smooth Descent (SD). We applied SD to real and simulated data, showing that in the presence of genotyping errors this method produces better genetic maps in terms of marker order and map length. SD is particularly useful for error rates between 5% and 20% and when error rates are not homogeneous among markers or individuals. With a starting error rate of 10%, SD reduced it to ∼5% in diploids, ∼7% in tetraploids and ∼8.5% in hexaploids. Conversely, the correlation between true and estimated genetic maps increased by 0.03 in tetraploids and by 0.2 in hexaploids, while worsening slightly in diploids (∼0.0011). We also show that the combination of genotype curation and map re-estimation allowed us to obtain better genetic maps while correcting wrong genotypes. We have implemented this algorithm in the R package Smooth Descent.

Microscopical palynology: Birch woodland expansion and species hybridisation coincide with periods of climate warming during the Holocene epoch in Iceland.

Introgressive hybridisation between arctic dwarf birch (Betula nana) and European downy birch (B.pubescens) is relatively common in natural woodlands in Iceland. As dwarf birch is a diploid species and downy tree birch a tetraploid species, their hybrids are triploid. In the introgression process, triploid hybrids, which are partially fertile, can backcross with the parental species, producing progenies comprising introgressed diploid, triploid and tetraploid plants. Triploid plants produce both normal triporate pollen (with three pores) and abnormal, aborted pollen, due to dysfunctional meiosis. The type of pollen abnormality that can be detected and quantified is non-triporate pollen (with four or more pores in the pollen wall). We therefore used the presence of non-triporate pollen as a marker to trace birch hybridisation in the past. In the current study we examined fossil pollen in samples from Holocene sediments from three locations: Grímsnes (SW), Eyjafjördur (N) and Thistilsfjördur (NE Iceland). All three peat monoliths had the starting age of 10.3cal. ka BP. Ages were calibrated based on known tephra layers and by radiocarbon dating. The size of Betula pollen grains was measured, and the species proportions calculated from size. Non-triporate grains were detected in samples from all three locations and throughout the Holocene, but with different frequencies. The peaks of intense hybridisation followed birch woodland expansion in two major periods of the Holocene, both coinciding with a warming of climate. The first period occurred in the Early Holocene, around 9.5-7cal. ka BP, when the climate warmed rapidly after deglaciation. The second period occurred around 5-3.5cal. ka BP, well within the mid-Holocene Northern Hemisphere warming. A new wave of birch hybridisation appears to have started in the last few decades as the climate has warmed. Birch woodlands are likely to become more widespread. Introgressed shrub birch is likely to be more competitive over dwarf birch.