Recent Duplication Research Articles

Identification of a candidate sex determination gene in Culaea inconstans suggests convergent recruitment of an Amh duplicate in two lineages of stickleback.

Sex chromosomes vary greatly in their age and levels of differentiation across the tree of life. This variation is largely due to the rates of sex chromosome turnover in different lineages; however, we still lack an explanation for why sex chromosomes are so conserved in some lineages (e.g. mammals, birds) but so labile in others (e.g. teleosts, amphibians). To identify general mechanisms driving transitions in sex determination systems or forces which favour their conservation, we first require empirical data on sex chromosome systems from multiple lineages. Stickleback fishes are a valuable model lineage for the study of sex chromosome evolution due to variation in sex chromosome systems between closely-related species. Here, we identify the sex chromosome and a strong candidate for the master sex determination gene in the brook stickleback, Culaea inconstans. Using whole-genome sequencing of wild-caught samples and a lab cross, we identify AmhY, a male specific duplication of the gene Amh, as the candidate master sex determination gene. AmhY resides on Chromosome 20 in C. inconstans and is likely a recent duplication, as both AmhY and the sex-linked region of Chromosome 20 show little sequence divergence. Importantly, this duplicate AmhY represents the second independent duplication and recruitment of Amh as the sex determination gene in stickleback and the eighth example known across teleosts. We discuss this convergence in the context of sex chromosome turnovers and the role that the Amh/AmhrII pathway, which is crucial for sex determination, may play in the evolution of sex chromosomes in teleosts.

Genome-Wide Identification and Expression Analysis of CsCaM/CML Gene Family in Response to Low-Temperature and Salt Stresses in Chrysanthemum seticuspe.

Calmodulin (CaM) and calmodulin-like proteins (CML) act as significant Ca2+ sensors binding Ca2+ with EF-hand motifs and have been reported to be involved in various environmental stresses in plants. In this study, calmodulin CsCaM/CML gene family members were identified based on the genome of Chrysanthemum seticuspe published recently; a phylogenetic tree was constructed; gene structures and chromosomal locations of CsCaM/CML were depicted; cis-acting regulatory elements were predicted; collinearity and duplicate events of CaM/CML were analyzed using MCScanX software; and the expression levels of CsCaM/CML in response to abiotic stress were analyzed, based on the published RNA-seq data. We identified 86 CsCaM/CML (4 CsCaMs and 82 CsCMLs) genes in total. Promoter sequences of CsCaM/CML contained elements related to abiotic stresses (including low-temperature and anaerobic stresses) and plant hormones (including abscisic acid (ABA), MeJA, and salicylic acid). CsCaM/CML genes were distributed on nine chromosomes unevenly. Collinearity analysis indicated that recent segmental duplications significantly enlarged the scale of the CML family in C. seticuspe. Four CsCMLs (CsCML14, CsCML50, CsCML65, and CsCML79) were statistically differentially regulated under low-temperature and salt stress compared with those in the normal condition. These results indicate diverse roles of CsCaM/CML in plant development and in response to environmental stimuli in C. seticuspe.

The Antarctic Moss Pohlia nutans Genome Provides Insights Into the Evolution of Bryophytes and the Adaptation to Extreme Terrestrial Habitats.

The Antarctic continent has extreme natural environment and fragile ecosystem. Mosses are one of the dominant floras in the Antarctic continent. However, their genomic features and adaptation processes to extreme environments remain poorly understood. Here, we assembled the high-quality genome sequence of the Antarctic moss (Pohlia nutans) with 698.20 Mb and 22 chromosomes. We found that the high proportion of repeat sequences and a recent whole-genome duplication (WGD) contribute to the large size genome of P. nutans when compared to other bryophytes. The genome of P. nutans harbors the signatures of massive segmental gene duplications and large expansions of gene families, likely facilitating neofunctionalization. Genomic characteristics that may support the Antarctic lifestyle of this moss comprise expanded gene families involved in phenylpropanoid biosynthesis, unsaturated fatty acid biosynthesis, and plant hormone signal transduction. Additional contributions include the significant expansion and upregulation of several genes encoding DNA photolyase, antioxidant enzymes, flavonoid biosynthesis enzymes, possibly reflecting diverse adaptive strategies. Notably, integrated multi-omic analyses elucidate flavonoid biosynthesis may function as the reactive oxygen species detoxification under UV-B radiation. Our studies provide insight into the unique features of the Antarctic moss genome and their molecular responses to extreme terrestrial environments.

Duplication of horizontally acquired GH5_2 enzymes played a central role in the evolution of longhorned beetles.

The rise of functional diversity through gene duplication contributed to the adaption of organisms to various environments. Here we investigate the evolution of putative cellulases of the subfamily 2 of glycoside hydrolase family 5 (GH5_2) in the Cerambycidae (longhorned beetles), a megadiverse assemblage of mostly xylophagous beetles. Cerambycidae originally acquired GH5_2 from a bacterial donor through horizontal gene transfer (HGT), and extant species harbor multiple copies that arose from gene duplication. We ask how these digestive enzymes contributed to the ability of these beetles to feed on wood. We analyzed 113 GH5_2, including the functional characterization of 52 of them, derived from 25 species covering most subfamilies of Cerambycidae. Ancestral gene duplications led to five well-defined groups with distinct substrate specificity, allowing these beetles to break down, in addition to cellulose, polysaccharides that are abundant in plant cell walls (PCWs), namely, xyloglucan, xylan, and mannans. Resurrecting the ancestral enzyme originally acquired by HGT, we show it was a cellulase that was able to break down glucomannan and xylan. Finally, recent gene duplications further expanded the catalytic repertoire of cerambycid GH5_2, giving rise to enzymes that favor transglycosylation over hydrolysis. We suggest that HGT and gene duplication, which shaped the evolution of GH5_2, played a central role in the ability of cerambycid beetles to use a PCW-rich diet and may have contributed to their successful radiation.

Comparative transcriptomic analysis unveils the deep phylogeny and secondary metabolite evolution of 116 Camellia plants.

Camellia plants include more than 200 species of great diversity and immense economic, ornamental, and cultural values. We sequenced the transcriptomes of 116 Camellia plants from almost all sections of the genus Camellia. We constructed a pan-transcriptome of Camellia plants with 89 394 gene families and then resolved the phylogeny of genus Camellia based on 405 high-quality low-copy core genes. Most of the inferred relationships are well supported by multiple nuclear gene trees and morphological traits. We provide strong evidence that Camellia plants shared a recent whole genome duplication event, followed by large expansions of transcription factor families associated with stress resistance and secondary metabolism. Secondary metabolites, particularly those associated with tea quality such as catechins and caffeine, were preferentially heavily accumulated in the Camellia plants from section Thea. We thoroughly examined the expression patterns of hundreds of genes associated with tea quality, and found that some of them exhibited significantly high expression and correlations with secondary metabolite accumulations in Thea species. We also released a web-accessible database for efficient retrieval of Camellia transcriptomes. The reported transcriptome sequences and obtained novel findings will facilitate the efficient conservation and utilization of Camellia germplasm towards a breeding program for cultivated tea, camellia, and oil-tea plants.

Repetitive Elements, Sequence Turnover and Cyto-Nuclear Gene Transfer in Gymnosperm Mitogenomes.

Among the three genomes in plant cells, the mitochondrial genome (mitogenome) is the least studied due to complex recombination and intergenomic transfer. In gymnosperms only ∼20 mitogenomes have been released thus far, which hinders a systematic investigation into the tempo and mode of mitochondrial DNA evolution in seed plants. Here, we report the complete mitogenome sequence of Platycladus orientalis (Cupressaceae). This mitogenome is assembled as two circular-mapping chromosomes with a size of ∼2.6 Mb and which contains 32 protein-coding genes, three rRNA and seven tRNA genes, and 1,068 RNA editing sites. Repetitive sequences, including dispersed repeats, transposable elements (TEs), and tandem repeats, made up 23% of the genome. Comparative analyses with 17 other mitogenomes representing the five gymnosperm lineages revealed a 30-fold difference in genome size, 80-fold in repetitive content, and 230-fold in substitution rate. We found dispersed repeats are highly associated with mitogenome expansion (r = 0.99), and most of them were accumulated during recent duplication events. Syntenic blocks and shared sequences between mitogenomes decay rapidly with divergence time (r = 0.53), with the exceptions of Ginkgo and Cycads which retained conserved genome structure over long evolutionary time. Our phylogenetic analysis supports a sister group relationship of Cupressophytes and Gnetophytes; both groups are unique in that they lost 8–12 protein-coding genes, of which 4–7 intact genes are likely transferred to nucleus. These two clades also show accelerated and highly variable substitution rates relative to other gymnosperms. Our study highlights the dynamic and enigmatic evolution of gymnosperm mitogenomes.

Chromosome-level genome assembly of Indian mangrove (Ceriops tagal) revealed a genome-wide duplication event predating the divergence of Rhizophoraceae mangrove species.

Mangrove ecosystems are unique, highly diverse, provide benefits to humans, and aid in coastal protection. The Indian mangrove, or spurred mangrove, [Ceriops tagal (Perr.) C. B. Rob.] is a member of the Rhizophoraceae family and is commonly found along the intertidal zones in tropical regions in Southeast Asia, southern Asia, and Africa. Here, we present the first high-quality reference genome assembly of the Ceriops species. A preliminary draft assembly, generated from the 10× Genomics linked-read library, was scaffolded using the proximity ligation chromatin contact mapping technique (Hi-C) to obtain a chromosome-scale assembly of 231,919,005 bases with an N50 length of 11,408,429 bases. The benchmarking universal single-copy orthologs (BUSCO) analysis revealed that C. tagal gene predictions recovered 95.8% of the highly conserved orthologs. Phylogenetic analyses suggested that C. tagal diverged from the last common ancestor of flat-leaf spurred mangrove [C. decandra (Griff.) Ding Hou] and C. zippeliana Blume ∼10.4 million yr ago (MYA), and the last common ancestor of genera Ceriops, Kandelia, and Rhizophora diverged from that of genus Bruguiera ∼49.4 MYA. In addition, our analysis of the transversion rate at fourfold-degenerate sites from orthologous gene pairs provided evidence supporting a recent whole-genome duplication in C. tagal. The STRUCTURE and principal component analyses illustrated that C. tagal individuals investigated in this study were the admixture of two subpopulations, the genetic background of which was influenced primarily by location. The availability of genomic and transcriptomic resources and biodiversity data reported in this work will be useful for future studies that may shed light on adaptive evolutions of mangrove species.

Genome-Wide Identification and Evolution Analysis of the Gibberellin Oxidase Gene Family in Six Gramineae Crops.

The plant hormones gibberellins (GAs) regulate plant growth and development and are closely related to the yield of cash crops. The GA oxidases (GAoxs), including the GA2ox, GA3ox, and GA20ox subfamilies, play pivotal roles in GAs’ biosynthesis and metabolism, but their classification and evolutionary pattern in Gramineae crops remain unclear. We thus conducted a comparative genomic study of GAox genes in six Gramineae representative crops, namely, Setaria italica (Si), Zea mays (Zm), Sorghum bicolor (Sb), Hordeum vulgare (Hv), Brachypodium distachyon (Bd), and Oryza sativa (Os). A total of 105 GAox genes were identified in these six crop genomes, belonging to the C19-GA2ox, C20-GA2ox, GA3ox, and GA20ox subfamilies. Based on orthogroup (OG) analysis, GAox genes were divided into nine OGs and the number of GAox genes in each of the OGs was similar among all tested crops, which indicated that GAox genes may have completed their family differentiations before the species differentiations of the tested species. The motif composition of GAox proteins showed that motifs 1, 2, 4, and 5, forming the 2OG-FeII_Oxy domain, were conserved in all identified GAox protein sequences, while motifs 11, 14, and 15 existed specifically in the GA20ox, C19-GA2ox, and C20-GA2ox protein sequences. Subsequently, the results of gene duplication events suggested that GAox genes mainly expanded in the form of WGD/SD and underwent purification selection and that maize had more GAox genes than other species due to its recent duplication events. The cis-acting elements analysis indicated that GAox genes may respond to growth and development, stress, hormones, and light signals. Moreover, the expression profiles of rice and maize showed that GAox genes were predominantly expressed in the panicles of the above two plants and the expression of several GAox genes was significantly induced by salt or cold stresses. In conclusion, our results provided further insight into GAox genes’ evolutionary differences among six representative Gramineae and highlighted GAox genes that may play a role in abiotic stress.

Evolution of enzymatic function through gene duplication in the metallocarboxypeptidase family

Changes in the genetic code occur frequently through, for example, DNA breaks and recombination events, often leading to the duplication of genes. The M14 family of metallocarboxypeptidases suggests such a history of gene duplication. In order to understand the functional outcome of such events, we investigated more recent duplication events in the M14 family of enzymes through analysis of genomic data found in the Ensembl genome database. Of the 23 M14 genes found in most vertebrates, four genes were present as two or more copies in many fish species. The AEBP1, CPXM1, and CPZ genes were duplicated in most fish with an arrangement suggesting a large‐scale event such as whole‐genome duplication. The CPO gene was also frequently duplicated in many fish species, but with copies tandemly arranged, suggestive of crossing‐over errors. A close analysis of the synteny of four CPO paralogs found in Xenopus tropicalis suggested potential transposon‐mediated duplication. Prediction of enzyme functions for the products of these genes showed that a greater amount of neofunctionalization was present in CPO paralogs as compared with other carboxypeptidase enzymes; purifying selection to maintain function remained detectable. To further examine function, the four CPO paralogs found in Xenopus tropicalis were examined biochemically. All were expressed in HEK293T cells, but no enzymatic activity was detected. Subcellular distribution suggested a possible function on lipid droplets similar to that previously shown for human CPO. Experiments to further examine enzyme activity using the Sf9 expression system are ongoing.

The Two Domains of the Avian Double-β-Defensin AvBD11 Have Different Ancestors, Common with Potential Monodomain Crocodile and Turtle Defensins.

Simple SummaryVertebrate defensins are a multigene family of antimicrobial peptides that evolved following a series of gene duplication and divergence events during the expansion of vertebrates. In birds, the repertoire of avian defensins contains an atypical defensin, namely AvBD11 (avian beta-defensin 11), which consists of two repeated but divergent defensin units (or domains) while most vertebrate defensins only possess one unit. In this study, we investigated the evolutionary scenario leading to the formation of this double defensin in birds by comparing each defensin unit of AvBD11 with other defensins from birds and closely related reptiles (crocodile, turtles) predicted to have a single defensin unit. Our most outstanding results suggest that the double defensin AvBD11 probably appeared following a fusion of two ancestral genes or from an ancestral double defensin, but not from a recent internal duplication as it can be observed in other types of proteins with domain repeats.Beta-defensins are an essential group of cysteine-rich host-defence peptides involved in vertebrate innate immunity and are generally monodomain. Among bird defensins, the avian β-defensin 11 (AvBD11) is unique because of its peculiar structure composed of two β-defensin domains. The reasons for the appearance of such ‘polydefensins’ during the evolution of several, but not all branches of vertebrates, still remain an open question. In this study, we aimed at exploring the origin and evolution of the bird AvBD11 using a phylogenetic approach. Although they are homologous, the N- and C-terminal domains of AvBD11 share low protein sequence similarity and possess different cysteine spacing patterns. Interestingly, strong variations in charge properties can be observed on the C-terminal domain depending on bird species but, despite this feature, no positive selection was detected on the AvBD11 gene (neither on site nor on branches). The comparison of AvBD11 protein sequences in different bird species, however, suggests that some amino acid residues may have undergone convergent evolution. The phylogenetic tree of avian defensins revealed that each domain of AvBD11 is distant from ovodefensins (OvoDs) and may have arisen from different ancestral defensins. Strikingly, our phylogenetic analysis demonstrated that each domain of AvBD11 has common ancestors with different putative monodomain β-defensins from crocodiles and turtles and are even more closely related with these reptilian defensins than with their avian paralogs. Our findings support that AvBD11′s domains, which differ in their cysteine spacing and charge distribution, do not result from a recent internal duplication but most likely originate from a fusion of two different ancestral genes or from an ancestral double-defensin arisen before the Testudines-Archosauria split.

Genomic Analysis Based on Chromosome-Level Genome Assembly Reveals an Expansion of Terpene Biosynthesis of Azadirachta indica.

Azadirachta indica (neem), an evergreen tree of the Meliaceae family, is a source of the potent biopesticide azadirachtin. The lack of a chromosome-level assembly impedes an in-depth understanding of its genome architecture and the comparative genomic analysis of A. indica. Here, a high-quality genome assembly of A. indica was constructed using a combination of data from Illumina, PacBio, and Hi-C technology, which is the first chromosome-scale genome assembly of A. indica. Based on the length of our assembly, the genome size of A. indica is estimated to be 281 Mb anchored to 14 chromosomes (contig N50 = 6 Mb and scaffold N50 = 19 Mb). The genome assembly contained 115 Mb repetitive elements and 25,767 protein-coding genes. Evolutional analysis revealed that A. indica didn’t experience any whole-genome duplication (WGD) event after the core eudicot γ event, but some genes and genome segment might likely experienced recent duplications. The secondary metabolite clusters, TPS genes, and CYP genes were also identified. Comparative genomic analysis revealed that most of the A. indica-specific TPS genes and CYP genes were located on the terpene-related clusters on chromosome 13. It is suggested that chromosome 13 may play an important role in the specific terpene biosynthesis of A. indica. The gene duplication events may be responsible for the terpene biosynthesis expansion in A. indica. The genomic dataset and genomic analysis created for A. indica will shed light on terpene biosynthesis in A. indica and facilitate comparative genomic research of the family Meliaceae.

A chromosome-level genome assembly of chia provides insights into high omega-3 content and coat color variation of its seeds

Chia (Salvia hispanica) is a functional food crop for humans. Although its seeds contain high omega-3 fatty acids, the seed yield of chia is still low. Genomic resources available for this plant are limited. We report the first high-quality chromosome-level genome sequence of chia. The assembled genome size was 347.6 Mb and covered 98.1% of the estimated genome size. A total of 31 069 protein-coding genes were predicted. The absence of recent whole-genome duplication and the relatively low intensity of transposable element expansion in chia compared to its sister species contribute to its small genome size. Transcriptome sequencing and gene duplication analysis reveal that the expansion of the fab2 gene family is likely to be related to the high content of omega-3 in seeds. The white seed coat color is determined by a single locus on chromosome 4. This study provides novel insights into the evolution of Salvia species and high omega-3 content, as well as valuable genomic resources for genetic improvement of important commercial traits of chia and its related species.

Multi-omics reveal differentiation and maintenance of dimorphic flowers in an alpine plant on the Qinghai-Tibet Plateau.

Dimorphic flowers growing on a single individual plant play a critical role in extreme adaption and reproductive assurance in plants and have high ecological and evolutionary significance. However, the omics bases underlying such a differentiation and maintenance remain largely unknown. We aimed to investigate this through genomic, transcriptome and metabolomic analyses of dimorphic flowers in an alpine biennial, Sinoswertia tetraptera (Gentianaceae). A high-quality chromosome-level genome sequence (903Mb) was first assembled for S.tetraptera with 31,359 protein-coding genes annotated. Two rounds of recent independent whole-genome duplication (WGD) were revealed. Numerous genes from the recent species-specific WGD were found to be differentially expressed in the two types of flowers, and this may have helped contribute to the origin of this innovative trait. The genes with contrasting expressions between flowers were related to biosynthesis of hormones, floral pigments (carotenoids and flavonoids) and iridoid compounds, which are involved in both flower development and colour. Metabolomic analyses similarly suggested differential concentrations of these chemicals in the two types of flowers. The expression interactions between multiple genes may together lead to contrasting morphology and chemical concentration and open versus closed pollination of the dimorphic flowers in this species for reproductive assurance.

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.

The first genome for the Cape Primrose Streptocarpus rexii (Gesneriaceae), a model plant for studying meristem-driven shoot diversity.

Cape Primroses (Streptocarpus, Gesneriaceae) are an ideal study system for investigating the genetics underlying species diversity in angiosperms. Streptocarpus rexii has served as a model species for plant developmental research for over five decades due to its unusual extended meristem activity present in the leaves. In this study, we sequenced and assembled the complete nuclear, chloroplast, and mitochondrial genomes of S. rexii using Oxford Nanopore Technologies long read sequencing. Two flow cells of PromethION sequencing resulted in 32 billion reads and were sufficient to generate a draft assembly including the chloroplast, mitochondrial and nuclear genomes, spanning 776 Mbp. The final nuclear genome assembly contained 5,855 contigs, spanning 766 Mbp of the 929‐Mbp haploid genome with an N50 of 3.7 Mbp and an L50 of 57 contigs. Over 70% of the draft genome was identified as repeats. A genome repeat library of Gesneriaceae was generated and used for genome annotation, with a total of 45,045 genes annotated in the S. rexii genome. Ks plots of the paranomes suggested a recent whole genome duplication event, shared between S. rexii and Primulina huaijiensis. A new chloroplast and mitochondrial genome assembly method, based on contig coverage and identification, was developed, and successfully used to assemble both organellar genomes of S. rexii . This method was developed into a pipeline and proved widely applicable. The nuclear genome of S. rexii and other datasets generated and reported here will be invaluable resources for further research to aid in the identification of genes involved in morphological variation underpinning plant diversification.

The complete sequence of a human genome.

Since its initial release in 2000, the human reference genome has covered only the euchromatic fraction of the genome, leaving important heterochromatic regions unfinished. Addressing the remaining 8% of the genome, the Telomere-to-Telomere (T2T) Consortium presents a complete 3.055 billion base pair (bp) sequence of a human genome, T2T-CHM13, that includes gapless assemblies for all chromosomes except Y, corrects errors in the prior references, and introduces nearly 200 million bp of sequence containing 1,956 gene predictions, 99 of which are predicted to be protein coding. The completed regions include all centromeric satellite arrays, recent segmental duplications, and the short arms of all five acrocentric chromosomes, unlocking these complex regions of the genome to variational and functional studies.

Innovation and Emerging Roles of Populus trichocarpa TEOSINTE BRANCHED1/CYCLOIDEA/PROLIFERATING CELL FACTOR Transcription Factors in Abiotic Stresses by Whole-Genome Duplication.

The TEOSINTE BRANCHED1/CYCLOIDEA/PROLIFERATING CELL FACTOR (TCP) family proteins are plant-specific transcription factors that have been well-acknowledged for designing the architectures of plant branch, shoot, and inflorescence. However, evidence for their innovation and emerging role in abiotic stress has been lacking. In this study, we identified a total of 36 TCP genes in Populus trichocarpa, 50% more than that in Arabidopsis (i.e., 24). Comparative intra-genomes showed that such significant innovation was mainly due to the most recent whole genome duplication (rWGD) in Populus lineage around Cretaceous-Paleogene (K-Pg) boundary after the divergence from Arabidopsis. Transcriptome analysis showed that the expressions of PtrTCP genes varied among leaf, stem, and root, and they could also be elaborately regulated by abiotic stresses (e.g., cold and salt). Moreover, co-expression network identified a cold-associated regulatory module including PtrTCP31, PtrTCP10, and PtrTCP36. Of them, PtrTCP10 was rWGD-duplicated from PtrTCP31 and evolved a strong capability of cold induction, which might suggest a neofunctionalization of PtrTCP genes and contribute to the adaptation of Populus lineage during the Cenozoic global cooling. Evidentially, overexpression of PtrTCP10 into Arabidopsis increased freezing tolerance and salt susceptibility. Integrating co-expression network and cis-regulatory element analysis confirmed that PtrTCP10 can regulate the well-known cold- and salt-relevant genes (e.g., ZAT10, GolS2, and SOS1), proving that PtrTCP10 is an evolutionary innovation in P. trichocarpa response to environmental changes. Altogether, our results provide evidence of the rWGD in P. trichocarpa responsible for the innovation of PtrTCP genes and their emerging roles in environmental stresses.

Genomic insights into lineage-specific evolution of the oleosin family in Euphorbiaceae

BackgroundLipid droplets (LDs) present in land plants serve as an essential energy and carbon reserve for seed germination and seedling development. Oleosins, the most abundant structural proteins of LDs, comprise a small family involved in LD formation, stabilization and degradation. Despite their importance, our knowledge on oleosins is still poor in Euphorbiaceae, a large plant family that contains several important oil-bearing species.ResultsTo uncover lineage-specific evolution of oleosin genes in Euphorbiaceae, in this study, we performed a genome-wide identification and comprehensive comparison of the oleosin family in Euphorbiaceae species with available genome sequences, i.e. castor bean (Ricinus communis), physic nut (Jatropha curcas), tung tree (Vernicia fordii), Mercurialis annua, cassava (Manihot esculenta) and rubber tree (Hevea brasiliensis), and a number of five, five, five, five, eight and eight members were found, respectively. Synteny analysis revealed one-to-one collinear relationship of oleosin genes between the former four (i.e. castor bean, physic nut, tung tree and M. annua) as well as latter two species (i.e. cassava and rubber tree), whereas one-to-one and one-to-two collinear relationships were observed between physic nut and cassava, reflecting the occurrence of one recent whole-genome duplication (WGD) in the last common ancestor of cassava and rubber tree. The presence of five ortholog groups representing three previously defined clades (i.e. U, SL and SH) dates back at least to the Malpighiales ancestor, because they are also conserved in poplar (Populus trichocarpa), a tree having experienced one Salicaceae-specific recent WGD. As observed in poplar, WGD was shown to be the main driver for the family expansion in both cassava and rubber tree. Nevertheless, same retention patterns of WGD-derived duplicates observed in cassava and rubber tree are somewhat different from that of poplar, though certain homologous fragments are still present in rubber tree. Further transcriptional profiling revealed an apparent seed-predominant expression pattern of oleosin genes in physic nut, castor bean and rubber tree. Moreover, structure and expression divergence of paralogous pairs were also observed in both cassava and rubber tree.ConclusionComparative genomics analysis of oleosin genes reported in this study improved our knowledge on lineage-specific family evolution in Euphorbiaceae, which also provides valuable information for further functional analysis and utilization of key members and their promoters.

Strategies of tolerance reflected in two North American maple genomes.

Significance StatementThe first chromosome‐scale assemblies for North American members of the Acer genus, sugar maple (Acer saccharum) and boxelder (Acer negundo), as well as transcriptomic evaluation of the abiotic stress response in A. saccharum are reported. This integrated study describes in‐depth aspects contributing to each species' approach to tolerance and applies current knowledge in many areas of plant genome biology with Acer physiology to help convey the genomic complexities underlying tolerance in broadleaf tree species.

Reference genome of lumpfish Cyclopterus lumpus Linnaeus provides evidence of male heterogametic sex determination through the AMH pathway.

Teleosts exhibit extensive diversity of sex determination (SD) systems and mechanisms, providing the opportunity to study the evolution of SD and sex chromosomes. Here we sequenced the genome of the common lumpfish (Cyclopterus lumpus Linnaeus), a species of increasing importance to aquaculture, and identified the SD region and master SD locus using a 70K single nucleotide polymorphism array and tissue-specific expression data. The chromosome-level assembly identified 25 diploid chromosomes with a total size of 572.89Mb, a scaffold N50 of 23.86Mb and genome annotation-predicted 21,480 protein-coding genes. Genome-wide association analysis located a highly sex-associated region on chromosome 13, suggesting that anti-Müllerian hormone (AMH) is the putative SD factor. Linkage disequilibrium and heterozygosity across chromosome 13support a proto-XX/XY system, with an absence of widespread chromosome divergence between sexes. We identified three copies of AMH in the lumpfish primary and alternate haplotype assemblies localized in the SD region. Comparison to sequences from other teleosts suggested a monophyletic relationship and conservation within the Cottioidei. One AMH copy showed similarity to AMH/AMHY in a related species and was also the only copy with expression in testis tissue, suggesting this copy may be the functional copy of AMH in lumpfish. The two other copies arranged in tandem inverted duplication were highly similar, suggesting a recent duplication event. This study provides a resource for the study of early sex chromosome evolution and novel genomic resources that benefits lumpfish conservation management and aquaculture.