Whole Genome Duplication Events Research Articles

The genome of okra (Abelmoschus esculentus) provides insights into its genome evolution and high nutrient content.

Okra (Abelmoschus esculentus) is an important vegetable crop with high nutritional value. However, the mechanism underlying its high nutrient content remains poorly understood. Here, we present a chromosome-scale genome of okra with a size of 1.19Gb. Comparative genomics analysis revealed the phylogenetic status of A. esculentus, as well as whole-genome duplication (WGD) events that have occurred widely across the Malvaceae species. We found that okra has experienced three additional WGDs compared with the diploid cotton Gossypium raimondii, resulting in a large chromosome number (2n = 130). After three WGDs, okra has undergone extensive genomic deletions and retained substantial numbers of genes related to secondary metabolite biosynthesis and environmental adaptation, resulting in significant differences between okra and G. raimondii in the gene families related to cellulose synthesis. Combining transcriptomic and metabolomic analysis, we revealed the relationship between gene expression and metabolite content change across different okra developmental stages. Furthermore, the sinapic acid/S-lignin biosynthesis-related gene families have experienced remarkable expansion in okra, and the expression of key enzymes involved in the sinapic acid/S-lignin biosynthesis pathway vary greatly across developmental periods, which partially explains the differences in metabolite content across the different stages. Our study gains insights into the comprehensive evolutionary history of Malvaceae species and the genetic basis that underlies the nutrient content changes in okra, which will facilitate the functional study and genetic improvement of okra varieties.

Comprehensive Analysis of YTH Domain-Containing Genes, Encoding m6A Reader and Their Response to Temperature Stresses and Yersinia ruckeri Infection in Rainbow Trout (Oncorhynchus mykiss)

YTH domain-containing genes are important readers of N6-methyladenosine (m6A) modifications with ability to directly affect the fates of distinct RNAs in organisms. Despite their importance, little is known about YTH domain-containing genes in teleosts until now. In the present study, a total of 10 YTH domain-containing genes have been systematically identified and functionally characterized in rainbow trout (Oncorhynchus mykiss). According to the phylogenetic tree, gene structure and syntenic analysis, these YTH domain-containing genes could be classified into three evolutionary subclades, including YTHDF, YTHDC1 and YTHDC2. Of them, the copy number of OmDF1, OmDF2, OmDF3, and OmDC1 were duplicated or even triplicated in rainbow trout due to the salmonid-specific whole-genome duplication event. The three-dimensional protein structure analysis revealed that there were similar structures and the same amino acid residues that were associated with cage formation between humans and rainbow trout, implying their similar manners in binding to m6A modification. Additionally, the results of qPCR experiment indicated that the expression patterns of a few YTH domain-containing genes, especially OmDF1b, OmDF3a and OmDF3b, were significantly different in liver tissue of rainbow trout under four different temperatures (7 °C, 11 °C, 15 °C, and 19 °C). The expression levels of OmDF1a, OmDF1b and OmDC1a were obviously repressed in spleen tissue of rainbow trout at 24 h after Yersinia ruckeri infection, while increased expression was detected in OmDF3b. This study provides a systemic overview of YTH domain-containing genes in rainbow trout and reveals their biological roles in responses to temperature stress and bacterial infection.

Genome-wide identification and characterization of the PbrATG family in Pyrus bretschneideri and functional analysis of PbrATG1a in response to Botryosphaeria dothidea

The pear is an economic fruit that is widely planted around the world and is loved by people for its rich nutritional value. Autophagy is a self-protection mechanism in eukaryotes, and its occurrence often accompanied by the degradation of damaged substances in cells and the recycling of nutrients. Autophagy is one of the mechanisms through which plants respond to environmental stress and plays an important role in plant development and stress resistance. Functional studies of autophagy-related genes (ATGs) have been performed on a variety of plant species, but little information is available on the ATG family in pear (Pyrus bretschneideri Rehd). Therefore, we analyzed the evolutionary dynamics and performed a genome-wide characterization of the PbrATG gene family. A total of 28 PbrATG members were identified. Phylogenetic analysis showed that PbrATGs were more closely related to ATGs of European pear and apple. Evolutionary analysis revealed that whole-genome duplication (WGD) and dispersed duplication events were the main driving forces of PbrATG family expansion. Expression analysis of different pear tissues showed that all the genes were expressed in different pear tissues, and different PbrATGs are expressed at different times and in different locations. Moreover, all PbrATGs also responded to different abiotic stresses, especially salt and drought stress, which elicited the highest expression levels. Pear seedlings were subsequently infected with Botryosphaeria dothidea (B. dothidea). The results showed that different PbrATGs had different expression patterns at different infection stages. According to the gene expression data, PbrATG1a was selected as a key autophagy gene for further analysis. Silencing of PbrATG1a reduced the resistance of pear to B. dothidea, which resulted in increased lesions, reactive oxygen species (ROS) contents, antioxidant enzyme activity, and gene expression levels in the silenced pear seedlings after B. dothidea inoculation. In this study, a comprehensive bioinformatic analysis of ATGs was conducted, and the functions of PbrATGs in pear development and in response to stress were elucidated, which laid a foundation for further study of the molecular mechanism of autophagy and a new strategy for pear resistance breeding.

Independent rediploidization masks shared whole genome duplication in the sturgeon-paddlefish ancestor

Whole genome duplication (WGD) is a dramatic evolutionary event generating many new genes and which may play a role in survival through mass extinctions. Paddlefish and sturgeon are sister lineages that both show genomic evidence for ancient WGD. Until now this has been interpreted as two independent WGD events due to a preponderance of duplicate genes with independent histories. Here we show that although there is indeed a plurality of apparently independent gene duplications, these derive from a shared genome duplication event occurring well over 200 million years ago, likely close to the Permian-Triassic mass extinction period. This was followed by a prolonged process of reversion to stable diploid inheritance (rediploidization), that may have promoted survival during the Triassic-Jurassic mass extinction. We show that the sharing of this WGD is masked by the fact that paddlefish and sturgeon lineage divergence occurred before rediploidization had proceeded even half-way. Thus, for most genes the resolution to diploidy was lineage-specific. Because genes are only truly duplicated once diploid inheritance is established, the paddlefish and sturgeon genomes are thus a mosaic of shared and non-shared gene duplications resulting from a shared genome duplication event.

Dynamics of Gene Loss following Ancient Whole-Genome Duplication in the Cryptic Paramecium Complex.

Whole-genome duplications (WGDs) have shaped the gene repertoire of many eukaryotic lineages. The redundancy created by WGDs typically results in a phase of massive gene loss. However, some WGD-derived paralogs are maintained over long evolutionary periods, and the relative contributions of different selective pressures to their maintenance are still debated. Previous studies have revealed a history of three successive WGDs in the lineage of the ciliate Paramecium tetraurelia and two of its sister species from the Paramecium aurelia complex. Here, we report the genome sequence and analysis of 10 additional P. aurelia species and 1 additional out group, revealing aspects of post-WGD evolution in 13 species sharing a common ancestral WGD. Contrary to the morphological radiation of vertebrates that putatively followed two WGD events, members of the cryptic P. aurelia complex have remained morphologically indistinguishable after hundreds of millions of years. Biases in gene retention compatible with dosage constraints appear to play a major role opposing post-WGD gene loss across all 13 species. In addition, post-WGD gene loss has been slower in Paramecium than in other species having experienced genome duplication, suggesting that the selective pressures against post-WGD gene loss are especially strong in Paramecium. A near complete lack of recent single-gene duplications in Paramecium provides additional evidence for strong selective pressures against gene dosage changes. This exceptional data set of 13 species sharing an ancestral WGD and 2 closely related out group species will be a useful resource for future studies on Paramecium as a major model organism in the evolutionary cell biology.

Expansion and adaptive evolution of the WRKY transcription factor family in Avicennia mangrove trees

Mangroves are adapted to intertidal zones, which present extreme environmental conditions. WRKYs are among the most prominent transcription factors (TFs) in higher plants and act through various interconnected networks to regulate responses to multiple abiotic stressors. Here, based on omic data, we investigated the landscape and evolutionary patterns of WRKYs in the main mangrove genus Avicennia. We found that both the number and the proportion of TFs and WRKYs in Avicennia species exceeded their inland relatives, indicating a significant expansion of WRKYs in Avicennia. We identified 109 WRKY genes in the representative species Avicennia marina. Comparative genomic analysis showed that two recent whole-genome duplication (WGD) events played a critical role in the expansion of WRKYs, and 88% of Avicennia marina WRKYs (AmWRKYs) have been retained following these WGDs. Applying comparative transcriptomics on roots under experimental salt gradients, we inferred that there is high divergence in the expression of WGD-retained AmWRKYs. Moreover, we found that the expression of 16 AmWRKYs was stable between freshwater and moderately saline water but increased when the trees were exposed to high salinity. In particular, 14 duplicates were retained following the two recent WGD events, indicating potential neo- and sub-functionalization. We also found that WRKYs could interact with other upregulated genes involved in signalling pathways and natural antioxidant biosynthesis to enhance salt tolerance, contributing to the adaptation to intertidal zones. Our omic data of the WRKY family in A. marina broadens the understanding of how a TF family relates to the adaptive evolution of mangroves.

E. urophylla × E. grandis high-quality genome and comparative genomics provide insights on evolution and diversification of eucalyptus

BackgroundEucalyptus urophylla × Eucalyptus grandis, an economically important forest tree, provides important raw material for energy and reduces damage to native forests. However, the absence of a high-quality E. urophylla × E. grandis reference genome has significantly hindered its evolution and genetic analysis.ResultsWe successfully presented a high-quality reference genome of E. urophylla × E. grandis (545.75 Mb; scaffold N50, 51.62 Mb) using a combination of the Illumina, PacBio HiFi, and Hi-C sequencing platforms. A total of 34,502 genes and 58.56% of the repetitive sequences in this genome were annotated. Using genome evolution analyses, we identified a recent whole-genome duplication (WGD) event in E. urophylla × E. grandis. We further found that gene families associated with starch and sucrose metabolism, flavonoid biosynthesis, and plant–pathogen interaction were significantly expanded in E. urophylla × E. grandis. Moreover, comparative genomic and evolutionary analyses showed large structural variations among the different chromosomes of the 34 Eucalyptus accessions, which were divided into six clades.ConclusionsOverall, our findings provide a valuable resource for expanding our understanding of the E. urophylla × E. grandis genome evolution, genetic improvement, and its comparative biology.

Tracing 100 million years of grass genome evolutionary plasticity.

Grasses derive from a family of monocotyledonous plants that includes crops of major economic importance such as wheat, rice, sorghum and barley, sharing a common ancestor some 100 million years ago. The genomic attributes of plant adaptation remain obscure and the consequences of recurrent whole genome duplications (WGD) or polyploidization events, a major force in plant evolution, remain largely speculative. We conducted a comparative analysis of omics data from ten grass species to unveil structural (inversions, fusions, fissions, duplications, substitutions) and regulatory (expression and methylation) basis of genome plasticity, as possible attributes of plant long lasting evolution and adaptation. The present study demonstrates that diverged polyploid lineages sharing a common WGD event often present the same patterns of structural changes and evolutionary dynamics, but these patterns are difficult to generalize across independent WGD events as a result of non-WGD factors such as selection and domestication of crops. Polyploidy is unequivocally linked to the evolutionary success of grasses during the past 100 million years, although it remains difficult to attribute this success to particular genomic consequences of polyploidization, suggesting that polyploids harness the potential of genome duplication, at least partially, in lineage-specific ways. Overall, the present study clearly demonstrates that post-polyploidization reprogramming is more complex than traditionally reported in investigating single species and calls for a critical and comprehensive comparison across independently polyploidized lineages.

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.

Large-scale long terminal repeat insertions produced a significant set of novel transcripts in cotton.

Genomic analysis has revealed that the 1,637-Mb Gossypium arboreum genome contains approximately 81% transposable elements (TEs), while only 57% of the 735-Mb G. raimondii genome is occupied by TEs. In this study, we investigated whether there were unknown transcripts associated with TE or TE fragments and, if so, how these new transcripts were evolved and regulated. As sequence depths increased from 4 to 100 G, a total of 10,284 novel intergenic transcripts (intergenic genes) were discovered. On average, approximately 84% of these intergenic transcripts possibly overlapped with the long terminal repeat (LTR) insertions in the otherwise untranscribed intergenic regions and were expressed at relatively low levels. Most of these intergenic transcripts possessed no transcription activation markers, while the majority of the regular genic genes possessed at least one such marker. Genes without transcription activation markers formed their+1 and -1 nucleosomes more closely (only (117±1.4)bp apart), while twice as big spaces (approximately (403.5±46.0) bp apart) were detected for genes with the activation markers. The analysis of 183 previously assembled genomes across three different kingdoms demonstrated systematically that intergenic transcript numbers in a given genome correlated positively with its LTR content. Evolutionary analysis revealed that genic genes originated during one of the whole-genome duplication events around 137.7 million years ago (MYA) for all eudicot genomes or 13.7 MYA for the Gossypium family, respectively, while the intergenic transcripts evolved around 1.6 MYA, resultant of the last LTR insertion. The characterization of these low-transcribed intergenic transcripts can facilitate our understanding of the potential biological roles played by LTRs during speciation and diversifications.

Phylogenetic placement of whole genome duplications in yeasts through quantitative analysis of hierarchical orthologous groups

Background: Whole genome duplications (WGD) are genomic events leading to formation of polyploid organisms. Resulting duplicated genes play important roles in driving species evolution and diversification. After such events, the initial ploidy is usually restored, complicating their detection across evolution. With the advance of bioinformatics and the rising number of new well-assembled genomes, new detection methods are ongoingly being developed to overcome the weaknesses of different approaches. Results: Here we propose a novel method for detecting WGD in yeast lineages based on the quantitative and comparative analysis of hierarchical orthologous groups (HOGs) of duplicated genes for a given set of organisms. We reconstruct ancestral genomes to obtain evolutionary information for each phylogenetic branch. This reconstruction relies on the inference of HOGs from the selected species’ proteomes. To estimate WGD events, the number of HOGs of duplicated genes across all taxonomic ranges are adjusted according to the molecular clock hypothesis and by the average genome size. Branches with a significant increase in the adjusted number of duplicated gene families are kept as candidates for WGD placement. The developed method was tested on two real datasets and showed promising results in phylogenetic WGD placements on the yeast lineage.

Genome-Wide Identification and Expression Analysis of CAMTA Gene Family Implies PbrCAMTA2 Involved in Fruit Softening in Pear

CAMTA are calcium-modulating binding transcription factors that contribute to plant development. We identified 46 CAMTA genes from eight Rosaceae species and divided them into five subgroups based on a phylogenetic tree. Our analysis indicated that CAMTA is a highly conserved family among Rosaceae species, with a conserved DNA-binding domain (CG-1) and a conserved transcription factor immunoglobulin domain (TIG). Following a recent whole-genome duplication event, the genomes of Chinese white pear, European pear, and apple experienced significant expansion, resulting in the number of CAMTA genes being twice that of the other species. Cis-element identification showed that the distribution of the zein metabolism regulation-responsive element was different in the promoters of Chinese white pear (55.56%) and European pear (11.11%) CAMTA gene families. The gene expression results showed that PbrCAMTA1, 2, 6, 7 was highly expressed in pear fruit. Among them, PbrCAMTA2 may have a key influence on fruit softening, as observed in transient transformation experiments. In conclusion, our results provide crucial insights into the evolution of the CAMTA gene family in pear and other Rosaceae species and identify a candidate PbrCAMTA gene, which is involved in the dynamic development of pear fruits.

Chromosome-scale genome assembly of marigold (Tagetes erecta L.): An ornamental plant and feedstock for industrial lutein production

Marigold (Tagetes erecta L., African marigold) is a widely grown ornamental plant and a main source of the carotenoid lutein for the industrial production of pharmaceuticals, food coloring, and feed additives. To gain a deeper understanding of the genetic mechanism of lutein in marigold, a chromosome-scale assembly of the marigold (T. erecta V-01) genome was completed based on Illumina, PacBio, and Hi-C reads. The 707.21-Mb assembled genome consisted of 35 834 annotated protein-coding genes, with 97.7% genomic integrity. We anchored 87.8% of the contigs (covering 621.20 Mb) to 12 pseudochromosomes, bringing the scaffold N50 length to 54.15 Mb. Phylogenetic analysis showed that marigold was closely related to the Asteraceae species bitter vine (Mikania micrantha) and sunflower (Helianthus annuus), all three of which originated in the Americas. Marigold diverged from the sunflower clade 23.57 million years ago (MYA) and from M. micrantha 19.59 MYA. Marigold has undergone three whole-genome duplication events, as well as a recent whole-genome duplication event (WGD-2) common to H. annuus and M. micrantha. Marigold gene families were significantly less expanded than those of M. micrantha or H. annuus, and the marigold genome contained significantly fewer interspersed repeats, which might account for its smaller genome. In addition, a range of candidate genes involved in the lutein biosynthetic pathway were identified. The high-quality reference genome obtained in this study provided a valuable genomic resource for studying the evolution of the Asteraceae family and for improving marigold breeding strategies.

Abstract 1690: Myc suppression permits entry into the cancer endocycle to evade toxic effects of chemotherapy

Abstract Recent advances in precision medicine and immunotherapy have made great strides in cancer treatment and care, but chemoresistance remains a major challenge. One method of chemoresistance employed by cancer cells is accession of the protective polyaneuploid cancer cell (PACC) state, in which cells enter a functionally dormant state upon environmental stress. Entry into the PACC state begins with a whole-genome duplication event followed by endocycling, where the cells are metabolically active and cycling through repeated G1, S and G2 phases of the cell cycle without dividing. While in this state, cells avoid the mitotic stress exerted by common chemotherapeutics and survive therapy, eventually undergoing a depolyploidization process to repopulate the tumor, driving cancer lethality. The mechanisms by which a cancer cell can switch between a mitotic cell cycle and an endocycle are unknown. Cancer is characterized by deregulated cell proliferation, often driven by mutation or amplification of genes that govern the cell cycle. C-myc, an oncogenic transcription factor, is dysregulated in roughly 70% of cancers and is clinically used as a negative prognostic factor for disease-free survival and metastasis. C-myc is responsible for the transcriptional activation of many cell cycle-related genes, and repression of cell cycle inhibitors. Our data indicate that endocycling cells in the PACC state downregulate c-myc expression. Transcriptomic analyses show repression of myc targets, along with supporting protein data that total c-myc protein levels are decreased. While myc downregulation is expected in non-proliferative cells, the compensatory mechanisms mediating the G1/S transition in the cancer endocycle remain unknown. These observations suggest the presence of underlying mechanisms to reprogram the mitotic cell cycle to access the PACC state. We hypothesize that myc downregulation is essential for cells to activate the cancer endocyle and enter the protective PACC state. Understanding the mechanisms that govern this transition is crucial to generate next-generation cancer therapeutics to block accession or disrupt this state. Citation Format: Michael Loycano, Kenneth J. Pienta, Sarah R. Amend. Myc suppression permits entry into the cancer endocycle to evade toxic effects of chemotherapy [abstract]. In: Proceedings of the American Association for Cancer Research Annual Meeting 2023; Part 1 (Regular and Invited Abstracts); 2023 Apr 14-19; Orlando, FL. Philadelphia (PA): AACR; Cancer Res 2023;83(7_Suppl):Abstract nr 1690.

A homology guide for Pacific salmon genus Oncorhynchus resolves patterns of ohnolog retention, resolution and local adaptation following the salmonid-specific whole-genome duplication event.

Salmonid fishes have emerged as a tractable model to study whole-genome duplications (WGDs) as this group has undergone four rounds of WGDs. While most of the salmonid genome has returned to a diploid state, a significant proportion of genes are maintained as duplicates and are referred to as ohnologs. The fact that much of the modern salmonid gene repertoire is comprised of ohnologs, while other genes have returned to their singleton state creates complications for genetic studies by obscuring homology relationships. The difficulty this creates is particularly prominent in Pacific salmonids belonging to genus Oncorhynchus who are the focus of intense genetics-based conservation and management efforts owing to the important ecological and cultural roles these fish play. To address this gap, we generated a homology guide for six species of Oncorhynchus with available genomes and used this guide to describe patterns of ohnolog retention and resolution. Overall, we find that ohnologs comprise approximately half of each species modern gene repertoires, which are functionally enriched for genes involved in DNA binding, while the less numerous singleton genes are heavily enriched in dosage-sensitive processes such as mitochondrial metabolism. Additionally, by reanalyzing published expression data from locally adapted strains of O. mykiss, we show that numerous ohnologs exhibit adaptive expression profiles; however, ohnologs are not more likely to display adaptive signatures than either paralogs or singletons. Finally, we demonstrate the utility of our homology guide by investigating the evolutionary relationship among genes highlighted as playing a role in salmonid life-history traits or gene editing targets.

The chromosome-level genome of Eucommia ulmoides provides insights into sex differentiation and α-linolenic acid biosynthesis.

Eucommia ulmoides Oliver is a typical dioecious plant endemic to China that has great medicinal and economic value. Here, we report a high-quality chromosome-level female genome of E. ulmoides obtained by PacBio and Hi-C technologies. The size of the female genome assembly was 1.01 Gb with 17 pseudochromosomes and 31,665 protein coding genes. In addition, Hi-C technology was used to reassemble the male genome released in 2018. The reassembled male genome was 1.24 Gb with the superscaffold N50 (48.30 Mb), which was increased 25.69 times, and the number of predicted genes increased by 11,266. Genome evolution analysis indicated that E. ulmoides has undergone two whole-genome duplication events before the divergence of female and male, including core eudicot γ whole-genome triplication event (γ-WGT) and a recent whole genome duplication (WGD) at approximately 27.3 million years ago (Mya). Based on transcriptome analysis, EuAP3 and EuAG may be the key genes involved in regulating the sex differentiation of E. ulmoides. Pathway analysis showed that the high expression of ω-3 fatty acid desaturase coding gene EU0103017 was an important reason for the high α-linolenic acid content in E. ulmoides. The genome of female and male E. ulmoides presented here is a valuable resource for the molecular biological study of sex differentiation of E. ulmoides and also will provide assistance for the breeding of superior varieties.

The evolution and immunomodulatory role of Zc3h12 proteins in zebrafish (Danio rerio)

Zc3h12 family is an important RNA-binding protein family regulating mRNA of inflammatory cytokines in mammals. However, there are few studies on their post-transcriptional level regulation of inflammatory cytokines in fish. Here, we investigated the evolution of zebrafish Zc3h12 family and explored their immunomodulatory role. Phylogenetic and syntenic analysis indicated the number of zc3h12 family members had increased ranging from a single member in invertebrates to a single copy of four members in mammals. As the most evolutionarily diverse group of vertebrates, the number of zc3h12 family members was more complex and diverse in the teleost, each member experienced different fates and followed different rules in multiple rounds of whole-genome duplication events. Thereinto, zebrafish contained three zc3h12 genes, among which zc3h12aa and zc3h12ab were duplicated from the same gene. Zebrafish Zc3h12 family could recognize the 3′-UTR regions of inflammatory cytokines through binding to the specific RNA secondary structure and negatively regulate their expression. Deletion of either Zc3h12 domains or mutation of the key amino acid in RNAase domain attenuated their modulatory effect, suggesting both domain and RNAase activity are important to the immunomodulatory role. These results elucidated the evolution of Zc3h12 family and uncovered Zc3h12-mediated post-transcriptional regulation of cytokines in zebrafish.

Nuclear phylogeny and insights into whole-genome duplications and reproductive development of Solanaceae plants

Solanaceae, the nightshade family, have ∼2700 species, including the important crops potato and tomato, ornamentals, and medicinal plants. Several sequenced Solanaceae genomes show evidence for whole-genome duplication (WGD), providing an excellent opportunity to investigate WGD and its impacts. Here, we generated 93 transcriptomes/genomes and combined them with 87 public datasets, for a total of 180 Solanaceae species representing all four subfamilies and 14 of 15 tribes. Nearly 1700 nuclear genes from these transcriptomic/genomic datasets were used to reconstruct a highly resolved Solanaceae phylogenetic tree with six major clades. The Solanaceae tree supports four previously recognized subfamilies (Goetzeioideae, Cestroideae, Nicotianoideae, and Solanoideae) and the designation of three other subfamilies (Schizanthoideae, Schwenckioideae, and Petunioideae), with the placement of several previously unassigned genera. We placed a Solanaceae-specific whole-genome triplication (WGT1) at ∼81 million years ago (mya), before the divergence of Schizanthoideae from other Solanaceae subfamilies at ∼73 mya. In addition, we detected two gene duplication bursts (GDBs) supporting proposed WGD events and four other GDBs. An investigation of the evolutionary histories of homologs of carpel and fruit developmental genes in 14 gene (sub)families revealed that 21 gene clades have retained gene duplicates. These were likely generated by the Solanaceae WGT1 and may have promoted fleshy fruit development. This study presents a well-resolved Solanaceae phylogeny and a new perspective on retained gene duplicates and carpel/fruit development, providing an improved understanding of Solanaceae evolution.

Evolutionary history of two evergreen Rhododendron species as revealed by chromosome-level genome assembly.

The genus Rhododendron (Ericaceae), a species-rich and widely distributed genus of woody plants, is distinguished for the beautiful and diverse flowers. Rhododendron delavayi Franch. and Rhododendron irroratum Franch., are highly attractive species widely distributed in south-west China and abundant new varieties have been selected from their genetic resources. We constructed chromosome-scale genome assemblies for Rhododendron delavayi and Rhododendron irroratum. Phylogenetic and whole-genome duplication analyses were performed to elucidate the evolutionary history of Rhododendron. Further, different types of gene duplications were identified and their contributions to gene family expansion were investigated. Finally, comprehensive characterization and evolutionary analysis of R2R3-MYB and NBS-encoding genes were conducted to explore their evolutionary patterns. The phylogenetic analysis classified Rhododendron species into two sister clades, 'rhododendrons' and 'azaleas'. Whole-genome duplication (WGD) analysis unveiled only one WGD event that occurred in Rhododendron after the ancestral γ triplication. Gene duplication and gene family expansion analyses suggested that the younger tandem and proximal duplications contributed greatly to the expansion of gene families involved in secondary metabolite biosynthesis and stress response. The candidate R2R3-MYB genes likely regulating anthocyanin biosynthesis and stress tolerance in Rhododendron will facilitate the breeding for ornamental use. NBS-encoding genes had undergone significant expansion and experienced species-specific gain and loss events in Rhododendron plants. The reference genomes presented here will provide important genetic resources for molecular breeding and genetic improvement of plants in this economically important Rhododendron genus.

The chromosome-level genome of double-petal phenotype jasmine (Jasminum sambac Aiton) provides insights into the biosynthesis of floral scent

Jasmine (Jasminum sambac Aiton) is a well-known cultivated plant species for its fragrant flowers used in the perfume industry and cosmetics. However, the genetic basis of its floral scent is largely unknown. In this study, using PacBio, Illumina, 10× Genomics and high-throughput chromosome conformation capture (Hi-C) sequencing technologies, a high-quality chromosome-level reference genome for J. sambac was obtained, exploiting a double-petal phenotype cultivar ‘Shuangbanmoli’ (JSSB). The results showed that the final assembled genome of JSSB is 580.33 Mb in size (contig N50 = 1.05 Mb; scaffold N50 = 45.07 Mb) with a total of 39 618 predicted protein-coding genes. Our analyses revealed that the JSSB genome has undergone an ancient whole-genome duplication (WGD) event at 91.68 million years ago (Mya). It was estimated that J. sambac diverged from the lineage leading to Olea europaea and Osmanthus fragrans about 28.8 Mya. On the basis of a combination of genomic, transcriptomic and metabolomic analyses, a range of floral scent volatiles and genes were identified involved in the benzenoid/phenylpropanoid and terpenoid biosynthesis pathways. The results provide new insights into the molecular mechanism of its fragrance biosynthesis in jasmine.