Species-specific Duplications Research Articles

Chromosome-level genome assembly provides insights into the genome evolution and functional importance of the phenylpropanoid–flavonoid pathway in Thymus mongolicus

BackgroundThymus mongolicus (family Lamiaceae) is a Thyme subshrub with strong aroma and remarkable environmental adaptability. Limited genomic information limits the use of this plant.ResultsChromosome-level 605.2 Mb genome of T. mongolicus was generated, with 96.28% anchored to 12 pseudochromosomes. The repetitive sequences were dominant, accounting for 70.98%, and 32,593 protein-coding genes were predicted. Synteny analysis revealed that Lamiaceae species generally underwent two rounds of whole genome duplication; moreover, species-specific genome duplication was identified. A recent LTR retrotransposon burst and tandem duplication might play important roles in the formation of the Thymus genome. Using comparative genomic analysis, phylogenetic tree of seven Lamiaceae species was constructed, which revealed that Thyme plants evolved recently in the family. Under the phylogenetic framework, we performed functional enrichment analysis of the genes on nodes that contained the most gene duplication events (> 50% support) and of relevant significant expanded gene families. These genes were highly associated with environmental adaptation and biosynthesis of secondary metabolites. Combined transcriptome and metabolome analyses revealed that Peroxidases, Hydroxycinnamoyl-CoA shikimate/quinate hydroxycinnamoyl transferases, and 4-coumarate-CoA ligases genes were the essential regulators of the phenylpropanoid–flavonoid pathway. Their catalytic products (e.g., apigenin, naringenin chalcone, and several apigenin-related compounds) might be responsible for the environmental tolerance and aromatic properties of T. mongolicus.ConclusionThis study enhanced the understanding of the genomic evolution of T. mongolicus, enabling further exploration of its unique traits and applications, and contributed to the understanding of Lamiaceae genomics and evolutionary biology.

From genomics to metabolomics: Deciphering sanguinarine biosynthesis in Dicranostigma leptopodum

Dicranostigma leptopodum (Maxim) Fedde (DLF) is a renowned medicinal plant in China, known to be rich in alkaloids. However, the unavailability of a reference genome has impeded investigation into its plant metabolism and genetic breeding potential. Here we present a high-quality chromosomal-level genome assembly for DLF, derived using a combination of Nanopore long-read sequencing, Illumina short-read sequencing and Hi-C technologies. Our assembly genome spans a size of 621.81 Mb with an impressive contig N50 of 93.04 Mb. We show that the species-specific whole-genome duplication (WGD) of DLF and Papaver somniferum corresponded to two rounds of WGDs of Papaver setigerum. Furthermore, we integrated comprehensive homology searching, gene family analyses and construction of a gene-to-metabolite network. These efforts led to the discovery of co-expressed transcription factors, including NAC and bZIP, alongside sanguinarine (SAN) pathway genes CYP719 (CFS and SPS). Notably, we identified P6H as a promising gene for enhancing SAN production. By providing the first reference genome for Dicranostigma, our study confirms the genomic underpinning of SAN biosynthesis and establishes a foundation for advancing functional genomic research on Papaveraceae species. Our findings underscore the pivotal role of high-quality genome assemblies in elucidating genetic variations underlying the evolutionary origin of secondary metabolites.

Integrated global analysis in spider flowers illuminates features underlying the evolution and maintenance of C4 photosynthesis.

The carbon concentrating mechanism-C4 photosynthesis-represents a classic example of convergent evolution, but how this important trait originated and evolved remains largely enigmatic. The spider flower Gynandropsis gynandra is a valuable leafy vegetable crop and medicinal plant that has also been recognized as a C4 model species. Here we present a high-quality chromosome-scale annotated genome assembly of G. gynandra through a combination of Oxford Nanopore Technology (ONT), HiFi and Hi-C technology. The 17 super-scaffolds cover 98.66% of the estimated genome (997.61Mb), with a contig N50 of 11.43Mb and a scaffold N50 of 51.02Mb. Repetitive elements occupy up to 71.91% of its genome, and over half are long terminal repeat retrotransposons (LTR-RTs) derived from recent bursts, contributing to genome size expansion. Strikingly, LTR-RT explosion also played a critical role in C4 evolution by altering expression features of photosynthesis-associated genes via preferential insertion in promoters. Integrated multiomics analyses of G. gynandra and the ornamental horticulture C3 relative Tarenaya hassleriana reveal that species-specific whole-genome duplication, gene family expansion, recent LTR-RT amplification, and more recent tandem duplication events have all facilitated the evolution of C4 photosynthesis, revealing uniqueness of C4 evolution in the Cleome genus. Moreover, high leaf vein density and heat stress resilience are associated with shifted gene expression patterns. The mode of C3-to-C4 transition found here yields new insights into evolutionary convergence of a complex plant trait. The availability of this reference-grade genomic resource makes G. gynandra an ideal model system facilitating efforts toward C4-aimed crop engineering.

Chromosome-level genome assembly of Murraya paniculata sheds light on biosynthesis of floral volatiles

BackgroundMurraya paniculata (L.) Jack, commonly called orange jessamine in the family Rutaceae, is an important ornamental plant in tropical and subtropical regions which is famous for its strong fragrance. Although genome assemblies have been reported for many Rutaceae species, mainly in the genus Citrus, full genomic information has not been reported for M. paniculata, which is a prerequisite for in-depth genetic studies on Murraya and manipulation using genetic engineering techniques. Here, we report a high-quality chromosome-level genome assembly of M. paniculata and aim to provide insights on the molecular mechanisms of flower volatile biosynthesis.ResultsThe genome assembly with a contig N50 of 18.25 Mb consists of 9 pseudomolecules and has a total length of 216.86 Mb. Phylogenetic analysis revealed that M. paniculata diverged from the common ancestor approximately 25 million years ago and has not undergone any species-specific whole genome duplication events. Genome structural annotation and comparative genomics analysis revealed that there are obvious differences in transposon contents among the genomes of M. paniculata and Citrus species, especially in the upstream regions of genes. Research on the flower volatiles of M. paniculata and C. maxima at three flowering stages revealed significant differences in volatile composition with the flowers of C. maxima lacking benzaldehyde and phenylacetaldehyde. Notably, there are transposons inserted in the upstream region of the phenylacetaldehyde synthase (PAAS) genes Cg1g029630 and Cg1g029640 in C. maxima, but not in the upstream region of three PAAS genes Me2G_2379, Me2G_2381, and Me2G_2382 in M. paniculata. Our results indicated that compared to the low expression levels of PAAS genes in C. maxima, the higher expression levels of the three PAAS genes in M. paniculata are the main factor affecting the phenylacetaldehyde biosynthesis and causing the content difference of phenylacetaldehyde. The phenylacetaldehyde synthetic activities of the enzymes encoded by M. paniculata PAAS genes were validated by in vitro analyses.ConclusionsOur study provides useful genomic resources of M. paniculata for further research on Rutaceae plants, identifies new PAAS genes, and provides insights into how transposons contribute to variations in flower volatiles among Murraya and Citrus plants.

FTR33, a member of fish-specific TRIM (finTRIM) subfamily, regulates negatively type I IFN antiviral immunity in zebrafish

In mammals, the tripartite motif (TRIM) proteins have been identified as critical factors involved in various cellular processes, including antiviral immunity. In teleost fish, a subfamily of fish-specific TRIM (finTRIM, FTR) has emerged in genus- or species-specific duplication. In this study, a finTRIM gene, called ftr33, was identified in zebrafish (Danio rerio), and phylogenic analysis revealed that FTR33 is closely related with zebrafish FTR14. The FTR33 protein contains all conservative domains reported in other finTRIMs. The ftr33 has a constitutive expression in embryos and in tissues/organs of adult fish, and its expression can be induced following spring viremia of carp virus (SVCV) infection and interferon (IFN) stimulation. The overexpression of FTR33 significantly downregulated the expression of type I IFNs and IFN-stimulated genes (ISGs) both in vitro and in vivo, respectively, leading to the increased replication of SVCV. It was also found that FTR33 interacted with melanoma differentiation associated gene 5 (MDA5) or mitochondrial anti-viral signaling protein (MAVS) to weaken the promoter activity of type I IFN. It is thus concluded that the FTR33, as an ISG, in zebrafish can negatively regulate IFN-mediated antiviral response.

Species-specific gene duplication in Clostridia produces variations of cholesterol-dependent cytolysin with different cytotoxicity

Cholesterol-dependent cytolysin (CDC) is a bacterial toxin that binds to eukaryotic cholesterol-containing membranes, forms oligomeric complexes, and is inserted into the bilayer to create large aqueous pores. Recently, we reported a species-specific duplication of the hemolysin gene in group III Clostridium botulinum. The duplicated genes (bly1 and bly2) encoded two separate CDC proteins (botulinolysins; BLY1 and BLY2). Here, we aimed to investigate whether BLY1 and BLY2 exert differential cytotoxicity. We isolated two bly genes from C. botulinum and evaluated the cytotoxicity of two recombinant BLY proteins (rBLY1 and rBLY2) in HeLa, IEC-6, and NRK cells. rBLYs were cytotoxic to equine erythrocytes. rBLY1 showed higher hemolytic activity than rBLY2. rBLY2 showed no or very weak cytotoxicity to the HeLa, IEC-6, and NRK cells, whereas rBLY1 showed high cytotoxicity to these cells. The comparison of the amino acid sequence of BLYs with those of other CDCs revealed that the already-known amino acid residues involved in cholesterol-containing membrane recognition, oligomerization, and insertion into membranes are well conserved in both BLYs. However, several amino acid substitutions were observed in the conserved regions, particularly in L2 and L3 regions involved in cell binding. These findings suggest that gene duplication in group III C. botulinum evolved distinct functional specializations, and differential cytotoxicity of BLY1 and BLY2 could be due to the amino acid substitution in the conserved regions. However, the structural and functional comparisons of the two BLYs are essential to gain insights into the function of the CDCs.

A chromosome-level genome assembly of the Chinese cork oak (Quercus variabilis).

Quercus variabilis (Fagaceae) is an ecologically and economically important deciduous broadleaved tree species native to and widespread in East Asia. It is a valuable woody species and an indicator of local forest health, and occupies a dominant position in forest ecosystems in East Asia. However, genomic resources from Q. variabilis are still lacking. Here, we present a high-quality Q. variabilis genome generated by PacBio HiFi and Hi-C sequencing. The assembled genome size is 787 Mb, with a contig N50 of 26.04 Mb and scaffold N50 of 64.86 Mb, comprising 12 pseudo-chromosomes. The repetitive sequences constitute 67.6% of the genome, of which the majority are long terminal repeats, accounting for 46.62% of the genome. We used ab initio, RNA sequence-based and homology-based predictions to identify protein-coding genes. A total of 32,466 protein-coding genes were identified, of which 95.11% could be functionally annotated. Evolutionary analysis showed that Q. variabilis was more closely related to Q. suber than to Q. lobata or Q. robur. We found no evidence for species-specific whole genome duplications in Quercus after the species had diverged. This study provides the first genome assembly and the first gene annotation data for Q. variabilis. These resources will inform the design of further breeding strategies, and will be valuable in the study of genome editing and comparative genomics in oak species.

A chromosome-scale genome assembly of Quercus gilva: Insights into the evolution of Quercus section Cyclobalanopsis (Fagaceae).

Quercus gilva is an ecologically and economically important species of Quercus section Cyclobalanopsis and is a dominant species in evergreen broad-leaved forests in subtropical regions of East Asia. In the present study, we reported a high-quality chromosome-scale genome assembly of Q. gilva, the first reference genome for section Cyclobalanopsis, using the combination of Illumina and PacBio sequencing with Hi-C technologies. The assembled genome size of Q. gilva was 889.71 Mb, with a contig number of 773 and a contig N50 of 28.32 Mb. Hi-C scaffolding anchored 859.07 Mb contigs (96.54% of the assembled genome) onto 12 pseudochromosomes, with a scaffold N50 of 70.35 Mb. A combination of de novo, homology-based, and transcript-based predictions predicted a final set of 36,442 protein-coding genes distributed on 12 pseudochromosomes, and 97.73% of them were functionally annotated. A total of 535.64 Mb (60.20%) of repetitive sequences were identified. Genome evolution analysis revealed that Q. gilva was most closely related to Q. suber and they diverged at 40.35 Ma, and Q. gilva did not experience species-specific whole-genome duplication in addition to the ancient gamma (γ) whole-genome triplication event shared by core eudicot plants. Q. gilva underwent considerable gene family expansion and contraction, with 598 expanded and 6,509 contracted gene families detected. The first chromosome-scale genome of Q. gilva will promote its germplasm conservation and genetic improvement and provide essential resources for better studying the evolution of Quercus section Cyclobalanopsis.

Chromosome-level genome assembly and characterization of Sophora Japonica.

Sophora japonica is a medium-size deciduous tree belonging to Leguminosae family and famous for its high ecological, economic and medicinal value. Here, we reveal a draft genome of S. japonica, which was ∼511.49 Mb long (contig N50 size of 17.34 Mb) based on Illumina, Nanopore and Hi-C data. We reliably assembled 110 contigs into 14 chromosomes, representing 91.62% of the total genome, with an improved N50 size of 31.32 Mb based on Hi-C data. Further investigation identified 271.76 Mb (53.13%) of repetitive sequences and 31,000 protein-coding genes, of which 30,721 (99.1%) were functionally annotated. Phylogenetic analysis indicates that S. japonica separated from Arabidopsis thaliana and Glycine max ∼107.53 and 61.24 million years ago, respectively. We detected evidence of species-specific and common-legume whole-genome duplication events in S. japonica. We further found that multiple TF families (e.g. BBX and PAL) have expanded in S. japonica, which might have led to its enhanced tolerance to abiotic stress. In addition, S. japonica harbours more genes involved in the lignin and cellulose biosynthesis pathways than the other two species. Finally, population genomic analyses revealed no obvious differentiation among geographical groups and the effective population size continuously declined since 2 Ma. Our genomic data provide a powerful comparative framework to study the adaptation, evolution and active ingredients biosynthesis in S. japonica. More importantly, our high-quality S. japonica genome is important for elucidating the biosynthesis of its main bioactive components, and improving its production and/or processing.

Recurrent Duplication and Diversification of Acrosomal Fertilization Proteins in Abalone.

Reproductive proteins mediating fertilization commonly exhibit rapid sequence diversification driven by positive selection. This pattern has been observed among nearly all taxonomic groups, including mammals, invertebrates, and plants, and is remarkable given the essential nature of the molecular interactions mediating fertilization. Gene duplication is another important mechanism that facilitates the generation of molecular novelty through functional divergence. Following duplication, paralogs may partition ancestral gene function (subfunctionalization) or acquire new roles (neofunctionalization). However, the contributions of duplication followed by sequence diversification to the molecular diversity of gamete recognition genes has been understudied in many models of fertilization. The marine gastropod mollusk abalone is a classic model for fertilization. Its two acrosomal proteins (lysin and sp18) are ancient gene duplicates with unique gamete recognition functions. Through detailed genomic and bioinformatic analyses we show how duplication events followed by sequence diversification has played an ongoing role in the evolution of abalone acrosomal proteins. The common ancestor of abalone had four members of its acrosomal protein family in a tandem gene array that repeatedly experienced positive selection. We find that both sp18 paralogs contain positively selected sites located in different regions of the paralogs, suggestive of functional divergence where selection acted upon distinct binding interfaces in each paralog. Further, a more recent species-specific duplication of both lysin and sp18 in the European abalone H. tuberculata is described. Despite clade-specific acrosomal protein paralogs, there are no concomitant duplications of egg coat proteins in H. tuberculata, indicating that duplication of egg proteins per se is not responsible for retention of duplicated acrosomal proteins. We hypothesize that, in a manner analogous to host/pathogen evolution, sperm proteins are selected for increased diversity through extensive sequence divergence and recurrent duplication driven by conflict mechanisms.

Domain-Dependent Evolution Explains Functional Homology of Protostome and Deuterostome Complement C3-Like Proteins.

Complement proteins emerged early in evolution but outside the vertebrate clade they are poorly characterized. An evolutionary model of C3 family members revealed that in contrast to vertebrates the evolutionary trajectory of C3-like genes in cnidarian, protostomes and invertebrate deuterostomes was highly divergent due to independent lineage and species-specific duplications. The deduced C3-like and vertebrate C3, C4 and C5 proteins had low sequence conservation, but extraordinarily high structural conservation and 2-chain and 3-chain protein isoforms repeatedly emerged. Functional characterization of three C3-like isoforms in a bivalve representative revealed that in common with vertebrates complement proteins they were cleaved into two subunits, b and a, and the latter regulated inflammation-related genes, chemotaxis and phagocytosis. Changes within the thioester bond cleavage sites and the a-subunit protein (ANATO domain) explained the functional differentiation of bivalve C3-like. The emergence of domain-related functions early during evolution explains the overlapping functions of bivalve C3-like and vertebrate C3, C4 and C5, despite low sequence conservation and indicates that evolutionary pressure acted to conserve protein domain organization rather than the primary sequence.

Evolution and expression of LEAFY genes in ferns and lycophytes

BackgroundThe LEAFY (LFY) transcription factors are present in algae and across land plants. The available expression and functional data of these genes in embryophytes suggest that LFY genes control a plethora of processes including the first zygotic cell division in bryophytes, shoot cell divisions of the gametophyte and sporophyte in ferns, cone differentiation in gymnosperms and floral meristem identity in flowering plants. However, their putative plesiomorphic role in plant reproductive transition in vascular plants remains untested.ResultsWe perform Maximum Likelihood (ML) phylogenetic analyses for the LFY gene lineage in embryophytes with expanded sampling in lycophytes and ferns. We recover the previously identified seed plant duplication that results in LEAFY and NEEDLY paralogs. In addition, we recover multiple species-specific duplications in ferns and lycophytes and large-scale duplications possibly correlated with the occurrence of whole genome duplication (WGD) events in Equisetales and Salviniales. To test putative roles in diverse ferns and lycophytes we perform LFY expression analyses in Adiantum raddianum, Equisetum giganteum and Selaginella moellendorffii. Our results show that LFY genes are active in vegetative and reproductive tissues, with higher expression in early fertile developmental stages and during sporangia differentiation.ConclusionsOur data point to previously unrecognized roles of LFY genes in sporangia differentiation in lycophytes and ferns and suggests that functions linked to reproductive structure development are not exclusive to seed plant LFY homologs.

Different scales of gene duplications occurring at different times have jointly shaped the NBS-LRR genes in Prunus species

In this study, genome-wide identification, phylogenetic relationships, duplication time and selective pressure of the NBS-LRR genes, an important group of plant disease-resistance genes (R genes), were performed to uncover their genetic evolutionary patterns in the six Prunus species. A total of 1946 NBS-LRR genes were identified; specifically, 589, 361, 284, 281, 318, and 113 were identified in Prunus yedoensis, P. domestica, P. avium, P. dulcis, P. persica and P. yedoensis var. nudiflora, respectively. Two NBS-LRR gene subclasses, TIR-NBS-LRR (TNL) and non-TIR-NBS-LRR (non-TNL), were also discovered. In total, 435 TNL and 1511 non-TNL genes were identified and could be classified into 30/55/75 and 103/158/191 multi-gene families, respectively, according to three different criteria. Higher Ks and Ka/Ks values were detected in TNL gene families than in non-TNL gene families. These results indicated that the TNL genes had more members involved in relatively ancient duplications and were affected by stronger selection pressure than the non-TNL genes. In general, the NBS-LRR genes were shaped by species-specific duplications, and lineage-specific duplications occurred at recent and relatively ancient periods among the six Prunus species. Therefore, different duplicated copies of NBS-LRRs can resist specific pathogens and will provide an R-gene library for resistance breeding in Prunus species.

Genome-wide characterization, evolution, structure, and expression analysis of the F-box genes in Caenorhabditis

BackgroundF-box proteins represent a diverse class of adaptor proteins of the ubiquitin-proteasome system (UPS) that play critical roles in the cell cycle, signal transduction, and immune response by removing or modifying cellular regulators. Among closely related organisms of the Caenorhabditis genus, remarkable divergence in F-box gene copy numbers was caused by sizeable species-specific expansion and contraction. Although F-box gene number expansion plays a vital role in shaping genomic diversity, little is known about molecular evolutionary mechanisms responsible for substantial differences in gene number of F-box genes and their functional diversification in Caenorhabditis. Here, we performed a comprehensive evolution and underlying mechanism analysis of F-box genes in five species of Caenorhabditis genus, including C. brenneri, C. briggsae, C. elegans, C. japonica, and C. remanei.ResultsHerein, we identified and characterized 594, 192, 377, 39, 1426 F-box homologs encoding putative F-box proteins in the genome of C. brenneri, C. briggsae, C. elegans, C. japonica, and C. remanei, respectively. Our work suggested that extensive species-specific tandem duplication followed by a small amount of gene loss was the primary mechanism responsible for F-box gene number divergence in Caenorhabditis genus. After F-box gene duplication events occurred, multiple mechanisms have contributed to gene structure divergence, including exon/intron gain/loss, exonization/pseudoexonization, exon/intron boundaries alteration, exon splits, and intron elongation by tandem repeats. Based on high-throughput RNA sequencing data analysis, we proposed that F-box gene functions have diversified by sub-functionalization through highly divergent stage-specific expression patterns in Caenorhabditis species.ConclusionsMassive species-specific tandem duplications and occasional gene loss drove the rapid evolution of the F-box gene family in Caenorhabditis, leading to complex gene structural variation and diversified functions affecting growth and development within and among Caenorhabditis species. In summary, our findings outline the evolution of F-box genes in the Caenorhabditis genome and lay the foundation for future functional studies.

Evolution and Potential Function in Molluscs of Neuropeptide and Receptor Homologues of the Insect Allatostatins.

The allatostatins (ASTs), AST-A, AST-B and AST-C, have mainly been investigated in insects. They are a large group of small pleotropic alloregulatory neuropeptides that are unrelated in sequence and activate receptors of the rhodopsin G-protein coupled receptor family (GPCRs). The characteristics and functions of the homologue systems in the molluscs (Buccalin, MIP and AST-C-like), the second most diverse group of protostomes after the arthropods, and of high interest for evolutionary studies due to their less rearranged genomes remains to be explored. In the present study their evolution is deciphered in molluscs and putative functions assigned in bivalves through meta-analysis of transcriptomes and experiments. Homologues of the three arthropod AST-type peptide precursors were identified in molluscs and produce a larger number of mature peptides than in insects. The number of putative receptors were also distinct across mollusc species due to lineage and species-specific duplications. Our evolutionary analysis of the receptors identified for the first time in a mollusc, the cephalopod, GALR-like genes, which challenges the accepted paradigm that AST-AR/buccalin-Rs are the orthologues of vertebrate GALRs in protostomes. Tissue transcriptomes revealed the peptides, and their putative receptors have a widespread distribution in bivalves and in the bivalve Mytilus galloprovincialis, elements of the three peptide-receptor systems are highly abundant in the mantle an innate immune barrier tissue. Exposure of M. galloprovincialis to lipopolysaccharide or a marine pathogenic bacterium, Vibrio harveyi, provoked significant modifications in the expression of genes of the peptide precursor and receptors of the AST-C-like system in the mantle suggesting involvement in the immune response. Overall, our study reveals that homologues of the arthropod AST-systems in molluscs are potentially more complex due to the greater number of putative mature peptides and receptor genes. In bivalves they have a broad and varying tissue distribution and abundance, and the elements of the AST-C-like family may have a putative function in the immune response.

Multi-tissue transcriptome analysis of two Begonia species reveals dynamic patterns of evolution in the chalcone synthase gene family

Begonia is an important horticultural plant group, as well as one of the most speciose Angiosperm genera, with over 2000 described species. Genus wide studies of genome size have shown that Begonia has a highly variable genome size, and analysis of paralog pairs has previously suggested that Begonia underwent a whole genome duplication. We address the contribution of gene duplication to the generation of diversity in Begonia using a multi-tissue RNA-seq approach. We chose to focus on chalcone synthase (CHS), a gene family having been shown to be involved in biotic and abiotic stress responses in other plant species, in particular its importance in maximising the use of variable light levels in tropical plants. We used RNA-seq to sample six tissues across two closely related but ecologically and morphologically divergent species, Begonia conchifolia and B. plebeja, yielding 17,012 and 19,969 annotated unigenes respectively. We identified the chalcone synthase gene family members in our Begonia study species, as well as in Hillebrandia sandwicensis, the monotypic sister genus to Begonia, Cucumis sativus, Arabidopsis thaliana, and Zea mays. Phylogenetic analysis suggested the CHS gene family has high duplicate turnover, all members of CHS identified in Begonia arising recently, after the divergence of Begonia and Cucumis. Expression profiles were similar within orthologous pairs, but we saw high inter-ortholog expression variation. Sequence analysis showed relaxed selective constraints on some ortholog pairs, with substitutions at conserved sites. Evidence of pseudogenisation and species specific duplication indicate that lineage specific differences are already beginning to accumulate since the divergence of our study species. We conclude that there is evidence for a role of gene duplication in generating diversity through sequence and expression divergence in Begonia.

The complex evolution of the metazoan HSP70 gene family

The metazoan 70-kDa heat shock protein (HSP70) family contains several members localized in different subcellular compartments. The cytosolic members have been classified into inducible HSP70s and constitutive heat shock cognates (HSC70s), but their distinction and evolutionary relationship remain unclear because of occasional reports of “constitutive HSP70s” and the lack of cross-phylum comparisons. Here we provide novel insights into the evolution of these important molecular chaperones. Phylogenetic analyses of 125 full-length HSP70s from a broad range of phyla revealed an ancient duplication that gave rise to two lineages from which all metazoan cytosolic HSP70s descend. One lineage (A) contains a relatively small number of genes from many invertebrate phyla, none of which have been shown to be constitutively expressed (i.e., either inducible or unknown). The other lineage (B) included both inducible and constitutive genes from diverse phyla. Species-specific duplications are present in both lineages, and Lineage B contains well-supported phylum-specific clades for Platyhelminthes, Rotifera, Nematoda, Porifera/Cnidaria, and Chordata. Some genes in Lineage B have likely independently acquired inducibility, which may explain the sporadic distribution of “HSP70” or “HSC70” in previous phylogenetic analyses. Consistent with the diversification history within each group, inducible members show lower purifying selection pressure compared to constitutive members. These results illustrate the evolutionary history of the HSP70 family, encouraging us to propose a new nomenclature: “HSP70 + subcellular localization + linage + copy number in the organism + inducible or constitutive, if known.” e.g., HSP70cA1i for cytosolic Lineage A, copy 1, inducible.

Evolution of hes gene family in vertebrates: the hes5 cluster genes have specifically increased in frogs

Backgroundhes genes are chordate homologs of Drosophila genes, hairy and enhancer of split, which encode a basic helix-loop-helix (bHLH) transcriptional repressor with a WRPW motif. Various developmental functions of hes genes, including early embryogenesis and neurogenesis, have been elucidated in vertebrates. However, their orthologous relationships remain unclear partly because of less conservation of relatively short amino acid sequences, the fact that the genome was not analyzed as it is today, and species-specific genome duplication. This results in complicated gene names in vertebrates, which are not consistent in orthologs. We previously revealed that Xenopus frogs have two clusters of hes5, named “the hes5.1 cluster” and “the hes5.3 cluster”, but the origin and the conservation have not yet been revealed.ResultsHere, we elucidated the orthologous and paralogous relationships of all hes genes of human, mouse, chicken, gecko, zebrafish, medaka, coelacanth, spotted gar, elephant shark and three species of frogs, Xenopus tropicalis (X. tropicalis), X. laevis, Nanorana parkeri, by phylogenetic and synteny analyses. Any duplicated hes5 were not found in mammals, whereas hes5 clusters in teleost were conserved although not as many genes as the three frog species. In addition, hes5 cluster-like structure was found in the elephant shark genome, but not found in cyclostomata.ConclusionThese data suggest that the hes5 cluster existed in the gnathostome ancestor but became a single gene in mammals. The number of hes5 cluster genes were specifically large in frogs.

Evolution of the Subgroup 6 R2R3-MYB Genes and Their Contribution to Floral Color in the Perianth-Bearing Piperales.

Flavonoids, carotenoids, betalains, and chlorophylls are the plant pigments responsible for floral color. Anthocyanins, a class of flavonoids, are largely responsible for the red, purple, pink, and blue colors. R2R3-MYB genes belonging to subgroup 6 (SG6) are the upstream regulatory factors of the anthocyanin biosynthetic pathway. The canonical members of these genes in Arabidopsis include AtMYB75, AtMYB90, AtMYB113, and AtMYB114. The Aristolochiaceae is an angiosperm lineage with diverse floral groundplans and perianth colors. Saruma henryi exhibits a biseriate perianth with green sepals and yellow petals. All other genera have sepals only, with colors ranging from green (in Lactoris) to a plethora of yellow to red and purple mixtures. Here, we isolated and reconstructed the SG6 R2R3-MYB gene lineage evolution in angiosperms with sampling emphasis in Aristolochiaceae. We found numerous species-specific duplications of this gene lineage in core eudicots and local duplications in Aristolochiaceae for Saruma and Asarum. Expression of SG6 R2R3-MYB genes examined in different developmental stages and plant organs of four Aristolochiaceae species, largely overlaps with red and purple pigments, suggesting a role in anthocyanin and flavonoid synthesis and accumulation. A directed RNA-seq analysis corroborated our RT-PCR analyses, by showing that these structural enzymes activate during perianth development in Aristolochia fimbriata and that the regulatory genes are expressed in correlation with color phenotype. Finally, the reconstruction of the flavonoid and anthocyanin metabolic pathways using predicted peptides from transcriptomic data show that all pivotal enzymes are present in the analyzed species. We conclude that the regulatory genes as well as the biosynthetic pathway are largely conserved across angiosperms. In addition, the Aristolochiaceae emerges as a remarkable group to study the genetic regulatory network for floral color, as their members exhibit an outstanding floral diversity with elaborate color patterns and the genetic complement for SG6 R2R3-MYB genes is simpler than in core eudicot model species.

The American Paddlefish Genome Provides Novel Insights into Chromosomal Evolution and Bone Mineralization in Early Vertebrates.

Sturgeons and paddlefishes (Acipenseriformes) occupy the basal position of ray-finned fishes, although they have cartilaginous skeletons as in Chondrichthyes. This evolutionary status and their morphological specializations make them a research focus, but their complex genomes (polyploidy and the presence of microchromosomes) bring obstacles and challenges to molecular studies. Here, we generated the first high-quality genome assembly of the American paddlefish (Polyodon spathula) at a chromosome level. Comparative genomic analyses revealed a recent species-specific whole-genome duplication event, and extensive chromosomal changes, including head-to-head fusions of pairs of intact, large ancestral chromosomes within the paddlefish. We also provide an overview of the paddlefish SCPP (secretory calcium-binding phosphoprotein) repertoire that is responsible for tissue mineralization, demonstrating that the earliest flourishing of SCPP members occurred at least before the split between Acipenseriformes and teleosts. In summary, this genome assembly provides a genetic resource for understanding chromosomal evolution in polyploid nonteleost fishes and bone mineralization in early vertebrates.