First chromosome scale genomes of ithomiine butterflies (Nymphalidae: Ithomiini): Comparative models for mimicry genetic studies

The ithomiine butterflies (Nymphalidae: Danainae) represent the largest known radiation of Müllerian mimetic butterflies. They dominate by number the mimetic butterfly communities, which include species such as the iconic neotropical Heliconius genus. Recent studies on the ecology and genetics of speciation in Ithomiini have suggested that sexual pheromones, colour pattern and perhaps hostplant could drive reproductive isolation. However, no reference genome was available for Ithomiini, which has hindered further exploration on the genetic architecture of these candidate traits, and more generally on the genomic patterns of divergence. Here, we generated high-quality, chromosome-scale genome assemblies for two Melinaea species, M. marsaeus and M. menophilus, and a draft genome of the species Ithomia salapia. We obtained genomes with a size ranging from 396 to 503 Mb across the three species and scaffold N50 of 40.5 and 23.2 Mb for the two chromosome-scale assemblies. Using collinearity analyses we identified massive rearrangements between the two closely related Melinaea species. An annotation of transposable elements and gene content was performed, as well as a specialist annotation to target chemosensory genes, which is crucial for host plant detection and mate recognition in mimetic species. A comparative genomic approach revealed independent gene expansions in ithomiines and particularly in gustatory receptor genes. These first three genomes of ithomiine mimetic butterflies constitute a valuable addition and a welcome comparison to existing biological models such as Heliconius, and will enable further understanding of the mechanisms of adaptation in butterflies.     

DNA barcoding of Cymbidium by genome skimming: Call for next-generation nuclear barcodes

Cymbidium is an orchid genus that has undergone rapid radiation and has high ornamental, economic, ecological and cultural importance, but its classification based on morphology is controversial. The plastid genome (plastome), as an extension of plant standard DNA barcodes, has been widely used as a potential molecular marker for identifying recently diverged species or complicated plant groups. In this study, we newly generated 237 plastomes of 50 species (at least two individuals per species) by genome skimming, covering 71.4% of members of the genus Cymbidium. Sequence-based analyses (barcoding gaps and automatic barcode gap discovery) and tree-based analyses (maximum likelihood, Bayesian inference and multirate Poisson tree processes model) were conducted for species identification of Cymbidium. Our work provides a comprehensive DNA barcode reference library for Cymbidium species identification. The results show that compared with standard DNA barcodes (rbcL + matK) as well as the plastid trnH-psbA, the species identification rate of the plastome increased moderately from 58% to 68%. At the same time, we propose an optimized identification strategy for Cymbidium species. The plastome cannot completely resolve the species identification of Cymbidium, the main reasons being incomplete lineage sorting, artificial cultivation, natural hybridization and chloroplast capture. To further explore the potential use of nuclear data in identifying species, the Skmer method was adopted and the identification rate increased to 72%. It appears that nuclear genome data have a vital role in species identification and are expected to be used as next-generation nuclear barcodes.     

ploidyfrost: Reference-free estimation of ploidy level from whole genome sequencing data based on de Bruijn graphs

Polyploidy is ubiquitous and its consequences are complex and variable. A change of ploidy level generally influences genetic diversity and results in morphological, physiological and ecological differences between cells or organisms with different ploidy levels. To avoid cumbersome experiments and take advantage of the less biased information provided by the vast amounts of genome sequencing data, computational tools for ploidy estimation are urgently needed. Until now, although a few such tools have been developed, many aspects of this estimation, such as the requirement of a reference genome, the lack of informative results and objective inferences, and the influence of false positives from errors and repeats, need further improvement. We have developed ploidyfrost , a de Bruijn graph-based method, to estimate ploidy levels from whole genome sequencing data sets without a reference genome. ploidyfrost provides a visual representation of allele frequency distribution generated using the ggplot2 package as well as quantitative results using the Gaussian mixture model. In addition, it takes advantage of colouring information encoded in coloured de Bruijn graphs to analyse multiple samples simultaneously and to flexibly filter putative false positives. We evaluated the performance of ploidyfrost by analysing highly heterozygous or repetitive samples of Cyclocarya paliurus and a complex allooctoploid sample of Fragaria × ananassa. Moreover, we demonstrated that the accuracy of analysis results can be improved by constraining a threshold such as Cramér's V coefficient on variant features, which may significantly reduce the side effects of sequencing errors and annoying repeats on the graphical structure constructed.     

A chromosome-level genome assembly enables the identification of the follicule stimulating hormone receptor as the master sex-determining gene in the flatfish Solea senegalensis

Sex determination (SD) shows huge variation among fish and a high evolutionary rate, as illustrated by the Pleuronectiformes (flatfishes). This order is characterized by its adaptation to demersal life, compact genomes and diversity of SD mechanisms. Here, we assembled the Solea senegalensis genome, a flatfish of great commercial value, into 82 contigs (614 Mb) combining long- and short-read sequencing, which were next scaffolded using a highly dense genetic map (28,838 markers, 21 linkage groups), representing 98.9% of the assembly. Further, we established the correspondence between the assembly and the 21 chromosomes by using BAC-FISH. Whole genome resequencing of six males and six females enabled the identification of 41 single nucleotide polymorphism variants in the follicle stimulating hormone receptor (fshr) consistent with an XX/XY SD system. The observed sex association was validated in a broader independent sample, providing a novel molecular sexing tool. The fshr gene displayed differential expression between male and female gonads from 86 days post-fertilization, when the gonad is still an undifferentiated primordium, concomitant with the activation of amh and cyp19a1a, testis and ovary marker genes, respectively, in males and females. The Y-linked fshr allele, which included 24 nonsynonymous variants and showed a highly divergent 3D protein structure, was overexpressed in males compared to the X-linked allele at all stages of gonadal differentiation. We hypothesize a mechanism hampering the action of the follicle stimulating hormone driving the undifferentiated gonad toward testis.     

Mitochondrial genome phylogeny reveals the deep-time origin of Gomphomastacinae (Orthoptera: Eumastacidae) and its alpine genera in China

Gomphomastacinae is a grasshopper subfamily in Eumastacidae, with a morphology and distribution distinct from other subfamilies. The alpine genera of Gomphomastacinae that inhabit the Qinghai–Tibet Plateau in China show unique characteristics adapted to high-altitude life. However, their phylogenetic position and biogeographic history remain controversial. Thus, to determine the diversification history of these alpine genera and the origin of the subfamily, we obtained mitochondrial genome sequences from all seven Gomphomastacinae genera distributed in China. The reconstructed phylogeny was well supported and confirmed the phylogenetic position of Gomphomastacinae within Eumastacidae. Time calibration revealed a deep-time origin of the subfamily dating back to the Cretaceous period, and the diversification among alpine genera was also an ancient pre-Miocene event (30–50 Ma). Based on phylogeny and time estimates, the most likely biogeographic scenario is that Gomphomastacinae originated from an ancestral lineage that lived in East Gondwana and dispersed to Central and Western Asia through India. Subsequently, the alpine genera likely diverged along with the uplift of the Qinghai–Tibet Plateau and survived drastic climate change by in situ adaptation to high-altitude dwellings.     

Comparative analysis of 343 plastid genomes of Solanum section Petota: Insights into potato diversity,phylogeny, and species discrimination

Common potato (Solanum tuberosum L.) and its wild relatives belong to Solanum section Petota. This section's phylogeny and species delimitation are complicated due to various ploidy levels, high heterozygosity, and frequent interspecific hybridization. Compared to the nuclear genome, the plastid genome is more conserved, has a haploid nature, and has a lower nucleotide substitution rate, providing informative alternative insights into the phylogenetic study of section Petota. Here, we analyzed 343 potato plastid genomes from 53 wild and four cultivated species. The diversity of sequences and genomes was comprehensively analyzed. A total of 24 species were placed in a phylogenetic tree based on genomic data for the first time. Overall, our results not only confirmed most existing clades and species boundaries inferred by nuclear evidence but also provided some distinctive species clade belonging and the maternally inherited evidence supporting the hybrid origin of some species. Furthermore, the divergence times between the major potato clades were estimated. In addition, the species discriminatory power of universal barcodes, nuclear ribosomal DNA, and whole and partial plastid genomes and their combinations were thoroughly evaluated; the plastid genome performed best but had limited discriminatory power for all survey species (40%). Overall, our study provided not only new insights into phylogeny and DNA barcoding of potato but also provided valuable genetic data resources for further systematical research of Petota.     

Genomic analyses reveal natural selection on reproduction related genes between two closely related Populus (Salicaceae) species

Identifying the factors that cause reproductive isolation and their relative importance in species divergence is crucial to our understanding of speciation processes. In most species, natural selection is commonly considered to play a large role in driving speciation. Based on whole genome re-sequencing data from 27 Populus alba and 28 Populus adenopoda individuals, we explored the factors related to reproductive isolation of these two closely related species. The results showed that the two species diverged ~5–10 million years ago (Ma), when the Qinghai–Tibet Plateau reached a certain height and the inland climate of the Asian continent became arid. In highly differentiated genomic regions, the relative divergence (F_ST) and absolute divergence (d_xy) were significantly higher than the genomic background, θπ and shared polymorphisms decreased whereas fixed differences increased, which indicated that natural selection played a key role in the reproductive isolation of the two species. In addition, we found several genes that were related to reproduction that may be involved in explaining the reproductive isolation. Using phylogenetic trees resolved from haplotype data of Populus tomentosa and P. adenopoda, the maternal origin of P. tomentosa from P. adenopoda was likely to be located in Hubei and Chongqing Provinces.     

Phylogenomic study of Amorphophallus (Alismatales; Araceae): When plastid DNA gene sequences help to resolve the backbone subgeneric delineation

Encompassing ca. 200 species distributed in paleotropical Africa and Asia, Amorphophallus is one of the largest genera of Araceae. In spite of the great economic interest in its glucomannan production, only a few studies have attempted to grasp the evolutionary history of this genus. In the current state of knowledge, four main clades, mostly linked to biogeographical delineation, have been identified from phylogenies based on a few genes. However, relationships among and within these clades still remain unclear, due to the rapid radiation that occurred during the early evolutionary history of the genus. Here, we generated genome skimming libraries for 43 specimens from 36 species distributed across the 4 clades, which allowed us to produce a phylogenetic matrix for a set of 71 plastid genes. Our phylogenies confirm the monophyly of these clades but show a new and well-resolved arrangement among these clades. Our analyses therefore provide a new scenario and timeline for the evolution of the main Amorphophallus clades, consistent with the morphological characteristics of the clades. The inferred scenario is also in agreement with climate dynamics and the onset of long-distance dispersal by the earliest migratory birds near the Oligocene/Miocene transition around 23 million years ago. Our study provides an up-to-date baseline to understand biogeographic and ecological processes that shaped the current diversity and distribution of Amorphophallus, paving the way for larger-scale phylogenomic studies based on plastid and nuclear genomes.     

EVOLUTIONARY RELATIONSHIPS BETWEEN FOUR SPECIES OF CLADOPHORA (CLADOPHORALES,CHLOROPHYTA) BASED ON DNA–DNA HYBRIDIZATION1

Analysis of the reassociation kinetics of the DNA from Cladophora pellucida (Huds.) Kütz. indicates that the genome of this benthic alga is comprised of approximately 75% repetitive sequences. Single-copy sequences reassociated with a rate constant of 1.8 × 10^?3 M^?1· s^?1, which corresponds to a haploid genome size of 4.7 × 10⁸ bp. Genotypic relationships between members of the form section Longiarticulatae were determined by the method of DNA–DNA hybridization. No significant divergence was observed between the single-copy sequences of C. pellucida isolates from the East Atlantic coast and Mediterranean Sea. Cladophora feredayi Harv. and C. att. ad pellucida from Australia and C. pellucidoidea van den Hoek from the West Atlantic coast were highly and about equally divergent from C. pellucida. The data support the hypothesis that the West Atlantic–West Pacific divergence reflects the middle Miocene closure of the Mediterranean–Indo-Pacific seaways, and the hypothesis that the Northwest Atlantic–Northeast Atlantic divergence reflects the middle Miocene thermal separation of these coasts.