Hybridization and a loss of sex shape genome-wide diversity and the origin of species in the evening primroses (Oenothera, Onagraceae)

Haplotype-phased genome unveils the butylphthalide biosynthesis and homoploid hybrid origin of Ligusticum chuanxiong

Butylphthalide is one of the first-line drugs for ischemic stroke therapy, while no biosynthetic enzyme for butylphthalide has been reported. Here, we present a haplotype-resolved genome of Ligusticum chuanxiong, a long-cultivated and phthalide-rich medicinal plant in Apiaceae. On the basis of comprehensive screening, four Fe(II)- and 2-oxoglutarate-dependent dioxygenases and two CYPs were mined and further biochemically verified as phthalide C-4/C-5 desaturases (P4,5Ds) that effectively promoted the forming of (S)-3-n-butylphthalide and butylidenephthalide. The substrate promiscuity and functional redundancy featured for P4,5Ds may contribute to the high phthalide diversity in L. chuanxiong. Notably, comparative genomic evidence supported L. chuanxiong as a homoploid hybrid with Ligusticum sinense as a potential parent. The two haplotypes demonstrated exceptional structure variance and diverged around 3.42 million years ago. Our study is an icebreaker for the dissection of phthalide biosynthetic pathway and reveals the hybrid origin of L. chuanxiong, which will facilitate the metabolic engineering for (S)-3-n-butylphthalide production and breeding for L. chuanxiong.

Migration of repetitive DNAs during evolution of the permanent translocation heterozygosity in the oyster plant (Tradescantia section Rhoeo)

Due to translocation heterozygosity for all chromosomes in the cell complement, the oyster plant (Tradescantia spathacea) forms a complete meiotic ring. It also shows Rabl-arrangement at interphase, featured by polar centromere clustering. We demonstrate that the pericentromeric regions of the oyster plant are homogenized in concert by three subtelomeric sequences: 45S rDNA, (TTTAGGG)n motif, and TSrepI repeat. The Rabl-based clustering of pericentromeric regions may have been an excellent device to combine the subtelomere-pericentromere sequence migration (via inversions) with the pericentromere-pericentromere DNA movement (via whole arm translocations) that altogether led to the concerted homogenization of all the pericentromeric domains by the subtelomeric sequences. We also show that the repetitive sequence landscape of interstitial chromosome regions contains many loci consisting of Arabidopsis-type telomeric sequence or of TSrepI repeat, and it is extensively heterozygous. However, the sequence arrangement on some chromosomal arms suggest segmental inversions that are fully or partially homozygous, a fact that could be explained if the inversions started to create linkages already in a bivalent-forming ancestor. Remarkably, the subterminal TSrepI loci reside exclusively on the longer arms that could be due to sharing sequences between similarly-sized chromosomal arms in the interphase nucleus. Altogether, our study spotlights the supergene system of the oyster plant as an excellent model to link complex chromosome rearrangements, evolution of repetitive sequences, and nuclear architecture.

The evolution of multi-gene families and metabolic pathways in the evening primroses (Oenothera: Onagraceae): A comparative transcriptomics approach

The plant genus Oenothera has played an important role in the study of plant evolution of genomes and plant defense and reproduction. Here, we build on the 1kp transcriptomic dataset by creating 44 new transcriptomes and analyzing a total of 63 transcriptomes to present a large-scale comparative study across 29 Oenothera species. Our dataset included 30.4 million reads per individual and 2.3 million transcripts on average. We used this transcriptome resource to examine genome-wide evolutionary patterns and functional diversification by searching for orthologous genes and performed gene family evolution analysis. We found wide heterogeneity in gene family evolution across the genus, with section Oenothera exhibiting the most pronounced evolutionary changes. Overall, more significant gene family expansions occurred than contractions. We also analyzed the molecular evolution of phenolic metabolism by retrieving proteins annotated for phenolic enzymatic complexes. We identified 1,568 phenolic genes arranged into 83 multigene families that varied widely across the genus. All taxa experienced rapid phenolic evolution (fast rate of genomic turnover) involving 33 gene families, which exhibited large expansions, gaining about 2-fold more genes than they lost. Upstream enzymes phenylalanine ammonia-lyase (PAL) and 4-coumaroyl: CoA ligase (4CL) accounted for most of the significant expansions and contractions. Our results suggest that adaptive and neutral evolutionary processes have contributed to Oenothera diversification and rapid gene family evolution.

Chromosome-level and haplotype-resolved genome provides insight into the tetraploid hybrid origin of patchouli

Patchouli (Pogostemon cablin (Blanco) Benth.), a member of the Lamiaceae family, is an important aromatic plant that has been widely used in medicine and perfumery. Here, we report a 1.94 Gb chromosome-scale assembly of the patchouli genome (contig N50 = 7.97 Mb). The gene annotation reveals that tandem duplication of sesquiterpene biosynthetic genes may be a major contributor to the biosynthesis of patchouli bioactivity components. We further phase the genome into two distinct subgenomes (A and B), and identify a chromosome substitution event that have occurred between them. Further investigations show that a burst of universal LTR-RTs in the A subgenome lead to the divergence between two subgenomes. However, no significant subgenome dominance is detected. Finally, we track the evolutionary scenario of patchouli including whole genome tetraploidization, subgenome divergency, hybridization, and chromosome substitution, which are the key forces to determine the complexity of patchouli genome. Our work sheds light on the evolutionary history of patchouli and offers unprecedented genomic resources for fundamental patchouli research and elite germplasm development.

The ploidy level of patchouli, an aromatic plant in the Lamiaceae family, remain unclear. Here, the authors assemble a chromosome-level and haplotype-resolved genome for patchouli and reveal that it is tetraploid hybrid as well as compensated aneuploidy.

Target enrichment and extensive population sampling help untangle the recent, rapid radiation of Oenothera sect. Calylophus

Oenothera sect. Calylophus is a North American group of 13 recognized taxa in the evening primrose family (Onagraceae) with an evolutionary history that may include independent origins of bee pollination, edaphic endemism, and permanent translocation heterozygosity. Like other groups that radiated relatively recently and rapidly, taxon boundaries within Oenothera sect. Calylophus have remained challenging to circumscribe. In this study, we used target enrichment, flanking non-coding regions, gene tree/species tree methods, tests for gene flow modified for target-enrichment data, and morphometric analysis to reconstruct phylogenetic hypotheses, evaluate current taxon circumscriptions, and examine character evolution in Oenothera sect. Calylophus. Because sect. Calylophus comprises a clade with a relatively restricted geographic range, we were able to extensively sample across the range of geographic, edaphic and morphological diversity in the group. We found that the combination of exons and flanking non-coding regions led to improved support for species relationships. We reconstructed potential hybrid origins of some accessions and note that if processes such as hybridization are not taken into account, the number of inferred evolutionary transitions may be artificially inflated. We recovered strong evidence for multiple evolutionary origins of bee pollination from ancestral hawkmoth pollination, edaphic specialization on gypsum, and permanent translocation heterozygosity. This study applies newly emerging techniques alongside dense infraspecific sampling and morphological analyses to effectively reconstruct the recalcitrant history of a rapid radiation.

A photosynthesis operon in the chloroplast genome drives speciation in evening primroses

Genetic incompatibility between the cytoplasm and the nucleus is thought to be a major factor in species formation, but mechanistic understanding of this process is poor. In evening primroses (Oenothera spp.), a model plant for organelle genetics and population biology, hybrid offspring regularly display chloroplast-nuclear incompatibility. This usually manifests in bleached plants, more rarely in hybrid sterility or embryonic lethality. Hence, most of these incompatibilities affect photosynthetic capability, a trait that is under selection in changing environments. Here we show that light-dependent misregulation of the plastid psbB operon, which encodes core subunits of photosystem II and the cytochrome b6f complex, can lead to hybrid incompatibility, and this ultimately drives speciation. This misregulation causes an impaired light acclimation response in incompatible plants. Moreover, as a result of their different chloroplast genotypes, the parental lines differ in photosynthesis performance upon exposure to different light conditions. Significantly, the incompatible chloroplast genome is naturally found in xeric habitats with high light intensities, whereas the compatible one is limited to mesic habitats. Consequently, our data raise the possibility that the hybridization barrier evolved as a result of adaptation to specific climatic conditions.

A nuclear phylogenomic study of the angiosperm order Myrtales, exploring the potential and limitations of the universal Angiosperms353 probe set

PREMISE

To further advance the understanding of the species-rich, economically and ecologically important angiosperm order Myrtales in the rosid clade, comprising nine families, approximately 400 genera and almost 14,000 species occurring on all continents (except Antarctica), we tested the Angiosperms353 probe kit.

METHODS

We combined high-throughput sequencing and target enrichment with the Angiosperms353 probe kit to evaluate a sample of 485 species across 305 genera (76% of all genera in the order).

RESULTS

Results provide the most comprehensive phylogenetic hypothesis for the order to date. Relationships at all ranks, such as the relationship of the early-diverging families, often reflect previous studies, but gene conflict is evident, and relationships previously found to be uncertain often remain so. Technical considerations for processing HTS data are also discussed.

CONCLUSIONS

High-throughput sequencing and the Angiosperms353 probe kit are powerful tools for phylogenomic analysis, but better understanding of the genetic data available is required to identify genes and gene trees that account for likely incomplete lineage sorting and/or hybridization events.

Identification of tandem repeat families from long-read sequences of Humulus lupulus

Hop (Humulus lupulus L.) is known for its use as a bittering agent in beer and has a rich history of cultivation, beginning in Europe and now spanning the globe. There are five wild varieties worldwide, which may have been introgressed with cultivated varieties. As a dioecious species, its obligate outcrossing, non-Mendelian inheritance, and genomic structural variability have confounded directed breeding efforts. Consequently, understanding the hop genome represents a considerable challenge, requiring additional resources. In order to facilitate investigations into the transmission genetics of hop, we report here a tandem repeat discovery pipeline developed using k-mer filtering and dot plot analysis of PacBio long-read sequences from the hop cultivar Apollo. From this we identified 17 new and distinct tandem repeat sequence families, which represent candidates for FISH probe development. For two of these candidates, HuluTR120 and HuluTR225, we produced oligonucleotide FISH probes from conserved regions of and demonstrated their utility by staining meiotic chromosomes from wild hop, var. neomexicanus to address, for example, questions about hop transmission genetics. Collectively, these tandem repeat sequence families represent new resources suitable for development of additional cytogenomic tools for hop research.

Genetic diversity, demographic history and neo-sex chromosomes in the Critically Endangered Raso lark

Small effective population sizes could expose island species to inbreeding and loss of genetic variation. Here, we investigate factors shaping genetic diversity in the Raso lark, which has been restricted to a single islet for approximately 500 years, with a population size of a few hundred. We assembled a reference genome for the related Eurasian skylark and then assessed diversity and demographic history using RAD-seq data (75 samples from Raso larks and two related mainland species). We first identify broad tracts of suppressed recombination in females, indicating enlarged neo-sex chromosomes. We then show that genetic diversity across autosomes in the Raso lark is lower than in its mainland relatives, but inconsistent with long-term persistence at its current population size. Finally, we find that genetic signatures of the recent population contraction are overshadowed by an ancient expansion and persistence of a very large population until the human settlement of Cape Verde. Our findings show how genome-wide approaches to study endangered species can help avoid confounding effects of genome architecture on diversity estimates, and how present-day diversity can be shaped by ancient demographic events.