Synteny Research Articles

Genome-wide exploration of the CONSTANS-like (COL) gene family and its potential role in regulating plant flowering time in foxtail millet (Setaria italica)

In photoperiod-sensitive plants, the CO-like gene (CONSTANS-like, COL) has a crucial function in regulating the timing of flowering. The blooming period in foxtail millet is greatly influenced by the duration of daylight; however, there is a scarcity of data regarding the molecular properties of the COL genes in the foxtail millet. In this study, we conducted a comprehensive analysis of the COL gene family in foxtail millet at the genome-wide level. We identified 11 SiCOL genes and performed gene structure analysis, which showed pronounced variation in gene length and intron number among the genes. The examination of COL proteins in foxtail millet and other plant species using phylogenetic analysis revealed that they could be clustered into three distinct groups. Cis-element analysis identified elements related to light-responsiveness, hormones, and abiotic stress in the promoter region of the SiCOL gene. Furthermore, tissue-specific expression analysis showed widespread expression of all 11 SiCOL genes in various foxtail millet tissues and organs, particularly in leaves and panicles. Collinearity analysis identified 14 syntenic gene pairs in both foxtail millet and rice. The results also revealed diurnal oscillations in the transcription levels of SiCOL genes under different light conditions. Moreover, among the 11 genes, SiCO, SiCOL1, and SiCOL6 expression levels were negatively correlated with flowering time variation in two foxtail cultivars. Additionally, upon constructing a network of predicted molecular interactions, FLOWER LOCUS-like (FTL) and Phytochromea A (PHY A) were suggested to potentially interact with SiCO, SiCOL1, and SiCOL6. SiCO, SiCOL1, and SiCOL6 have the potential for flowering and heading in foxtail millet. This research enhances our comprehension of the role and control of the SiCOL gene family constituents in foxtail millet, establishing a basis for future investigations.

Genome-wide analysis of the MYB gene family and functional analysis of BhMYB79 in wax gourd

Wax gourd (Benincasa hispida) is an important cucurbit crop with economic and medicinal value. The myeloblastosis (MYB) gene family is one of the largest gene families in plants and regulates various biological processes, whereas the MYB gene family has not been systematically studied in wax gourd. In this study, we performed genome-wide identification of the MYB gene family in wax gourd and analyzed their phylogenetic relationship, MYB DNA-binding domain (MYB DBD), gene structure, protein motif, synteny, duplication mode and expression pattern. As a result, a total of 215 BhMYB genes (BhMYBs) were identified, belonging to four subfamilies: 1R-, 2R-, 3R- and 4R-MYB subfamilies. Genes of 1R-MYB subfamily and 2R-MYB subfamily were subdivided into different subgroups respectively. The analysis of MYB DBD, gene structure and protein motif showed that the most genes in the same subgroup have similar characteristics and the 2R-MYB genes were more conserved than the 1R-MYB genes. Interestingly, the long terminal retrotransposons (LTR-RTs) were found in the long introns of several BhMYBs. The results of synteny analysis showed that there were more syntenic gene pairs between wax gourd and other cucurbit crops, while the least number of syntenic gene pairs existed between wax gourd and rice. Gene duplication was the main reason for the expansion of the MYB gene family in wax gourd, with the transposed duplication (TRD) mode contributing more. All duplication BhMYB genes were under purifying selection pressure. Further expression analysis showed that many BhMYBs exhibited obvious tissue-specific expression and several BhMYBs were significantly induced by one or more abiotic stresses. BhMYB79 was particularly expressed in roots and significantly induced by salt, drought, cold and heat stresses, overexpression of which led to reduced tolerance to salt stress in Arabidopsis. In conclusion, our results provide a systematic analysis of wax gourd MYB gene family and facilitate the biological role study of BhMYB79 during wax gourd salt stress response process.

Aconitate decarboxylase 1 (Acod1) suppresses inflammatory stimuli-induced immune responses in Japanese flounder (Paralichthys olivaceus)

Aconitate decarboxylase 1 (ACOD1, formerly termed as immunoresponsive gene 1) is an enzyme responsible for generation of anti-inflammatory metabolite itaconate through decarboxylation of cis-aconitate and plays a key role in inflammation regulation and anti-microbial host defense in mammals. However, the immune functions of Acod1 in fish remain largely undetermined. Here, we characterize a Acod1 homolog that is highly induced by LPS, poly(I:C) and extracellular ATP (eATP) challenges in the Japanese flounder Paralichthys olivaceus head kidney macrophages (HKMs). We find that Acod1 syntenic genes in P. olivaceus are similar with those in mammalian species and the deduced Japanese flounder ACOD1 protein shares similar protein domains and 3-D structures with its counterparts in mammals. We next show that itaconate derivatives dimethyl itaconate (DMI) and citraconate (Cit) can markedly down-regulate inflammatory stimuli, including LPS-, poly(I:C)- and eATP-induced inflammatory mediator nitric oxide (NO) production and gene expression of proinflammatory cytokines (IL-1β and TNF-α), chemokines (CXCL-8, CXCL-9 and CXCL-13), HMGB1 and inflammasome components (NLRP3, NLRC3, caspase1, ASC and GSDME), but up-regulate anti-inflammatory gene (IL-10 and SOCS1/3) expression in P. olivaceus HKMs. In addition, DMI and Cit treatment significantly inhibits Edwardsiella piscicida growth in P. olivaceus HKMs. Furthermore, knockdown of the endogenous Acod1 expression by siRNA silencing not only significantly increases LPS-, poly(I:C)- and eATP-induced proinflammatory gene expression but also decreases anti-inflammatory gene expression in the HKMs. In contrast, the LPS-, poly(I:C)- and eATP-indued proinflammatory gene expression and NO production are greatly attenuated in ACOD1 overexpressed Japanese flounder FG-9307 cells. Collectively, our studies demonstrate that the inducible fish Acod1 gene is a novel immune regulator in dampening inflammatory stimuli-induced immune responses and ACOD1-produced itaconate may play an essential role in anti-bacterial immunity in fish macrophages.

Convergent evolution of desiccation tolerance in grasses.

Desiccation tolerance has evolved repeatedly in plants as an adaptation to survive extreme environments. Plants use similar biophysical and cellular mechanisms to survive life without water, but convergence at the molecular, gene and regulatory levels remains to be tested. Here we explore the evolutionary mechanisms underlying the recurrent evolution of desiccation tolerance across grasses. We observed substantial convergence in gene duplication and expression patterns associated with desiccation. Syntenic genes of shared origin are activated across species, indicative of parallel evolution. In other cases, similar metabolic pathways are induced but using different gene sets, pointing towards phenotypic convergence. Species-specific mechanisms supplement these shared core mechanisms, underlining the complexity and diversity of evolutionary adaptations to drought. Our findings provide insight into the evolutionary processes driving desiccation tolerance and highlight the roles of parallel and convergent evolution in response to environmental challenges.

The biosynthesis of asiaticoside and madecassoside reveals tandem duplication–directed evolution of glycoside glycosyltransferases in the Apiales

Certain plant species within the Apiales order accumulate triterpenoid saponins that feature a distinctive glucose-glucose-rhamnose (G-G-R) sugar chain attached at the C-28 position of the pentacyclic triterpene skeleton. Until recently, the genomic basis underlying the biosynthesis and evolution of this sugar chain has remained elusive. In this study, we identified two novel glycoside glycosyltransferases (GGTs) that can sequentially install the sugar chain’s second D-glucose and third L-rhamnose during the biosynthesis of asiaticoside and madecassoside, two representative G-G-R sugar chain–containing triterpenoid saponins produced by Centella asiatica. Enzymatic assays revealed the remarkable substrate promiscuity of the two GGTs and the key residues crucial for sugar-donor selectivity of the glucosyltransferase and rhamnosyltransferase. We further identified syntenic tandem gene duplicates of the two GGTs in the Apiaceae and Araliaceae families, suggesting a well-conserved genomic basis underlying sugar chain assembly that likely has evolved in the early ancestors of the Apiales order. Moreover, expression patterns of the two GGTs in pierced leaves of C. asiatica were found to be correlated with the production of asiaticoside and madecassoside, implying their involvement in host defense against herbivores and pathogens. Our work sheds light on the biosynthesis and evolution of complex saponin sugars, paving the way for future engineering of diverse bioactive triterpenoids with unique glycoforms.

Transposable elements in Rosaceae: insights into genome evolution, expression dynamics, and syntenic gene regulation.

Transposable elements (TEs) exert significant influence on plant genomic structure and gene expression. Here, we explored TE-related aspects across 14 Rosaceae genomes, investigating genomic distribution, transposition activity, expression patterns, and nearby differentially expressed genes (DEGs). Analyses unveiled distinct long terminal repeat retrotransposon (LTR-RT) evolutionary patterns, reflecting varied genome size changes among nine species over the past million years. In the past 2.5 million years, Rubus idaeus showed a transposition rate twice as fast as Fragaria vesca, while Pyrus bretschneideri displayed significantly faster transposition compared with Crataegus pinnatifida. Genes adjacent to recent TE insertions were linked to adversity resistance, while those near previous insertions were functionally enriched in morphogenesis, enzyme activity, and metabolic processes. Expression analysis revealed diverse responses of LTR-RTs to internal or external conditions. Furthermore, we identified 3695 pairs of syntenic DEGs proximal to TEs in Malus domestica cv. 'Gala' and M. domestica (GDDH13), suggesting TE insertions may contribute to varietal trait differences in these apple varieties. Our study across representative Rosaceae species underscores the pivotal role of TEs in plant genome evolution within this diverse family. It elucidates how these elements regulate syntenic DEGs on a genome-wide scale, offering insights into Rosaceae-specific genomic evolution.

Dynamic changes in the plastid and mitochondrial genomes of the angiosperm Corydalis pauciovulata (Papaveraceae)

BackgroundCorydalis DC., the largest genus in the family Papaveraceae, comprises > 465 species. Complete plastid genomes (plastomes) of Corydalis show evolutionary changes, including syntenic arrangements, gene losses and duplications, and IR boundary shifts. However, little is known about the evolution of the mitochondrial genome (mitogenome) in Corydalis. Both the organelle genomes and transcriptomes are needed to better understand the relationships between the patterns of evolution in mitochondrial and plastid genomes.ResultsWe obtained complete plastid and mitochondrial genomes from Corydalis pauciovulata using a hybrid assembly of Illumina and Oxford Nanopore Technologies reads to assess the evolutionary parallels between the organelle genomes. The mitogenome and plastome of C. pauciovulata had sizes of 675,483 bp and 185,814 bp, respectively. Three ancestral gene clusters were missing from the mitogenome, and expanded IR (46,060 bp) and miniaturized SSC (202 bp) regions were identified in the plastome. The mitogenome and plastome of C. pauciovulata contained 41 and 67 protein-coding genes, respectively; the loss of genes was a plastid-specific event. We also generated a draft genome and transcriptome for C. pauciovulata. A combination of genomic and transcriptomic data supported the functional replacement of acetyl-CoA carboxylase subunit β (accD) by intracellular transfer to the nucleus in C. pauciovulata. In contrast, our analyses suggested a concurrent loss of the NADH-plastoquinone oxidoreductase (ndh) complex in both the nuclear and plastid genomes. Finally, we performed genomic and transcriptomic analyses to characterize DNA replication, recombination, and repair (DNA-RRR) genes in C. pauciovulata as well as the transcriptomes of Liriodendron tulipifera and Nelumbo nuicifera. We obtained 25 DNA-RRR genes and identified their structure in C. pauciovulata. Pairwise comparisons of nonsynonymous (dN) and synonymous (dS) substitution rates revealed that several DNA-RRR genes in C. pauciovulata have higher dN and dS values than those in N. nuicifera.ConclusionsThe C. pauciovulata genomic data generated here provide a valuable resource for understanding the evolution of Corydalis organelle genomes. The first mitogenome of Papaveraceae provides an example that can be explored by other researchers sequencing the mitogenomes of related plants. Our results also provide fundamental information about DNA-RRR genes in Corydalis and their related rate variation, which elucidates the relationships between DNA-RRR genes and organelle genome stability.

Identification of the ALMT gene family in the potato (Solanum tuberosum L.) and analysis of the function of StALMT6/10 in response to aluminum toxicity.

Aluminum (Al)-activated malate transporters (ALMTs) play an important role in the response to Al toxicity, maintenance of ion homeostasis balance, mineral nutrient distribution, and fruit quality development in plants. However, the function of the StALMT gene family in potato remains unknown. In this study, 14 StALMT genes were identified from the potato genome, unevenly distributed on seven different chromosomes. Collinearity and synteny analyses of ALMT genes showed that potatoes had 6 and 22 orthologous gene pairs with Arabidopsis and tomatoes, respectively, and more than one syntenic gene pair was identified for some StALMT gene family members. Real-time quantitative polymerase chain reaction (qPCR) results showed differential expression levels of StALMT gene family members in different tissues of the potato. Interestingly, StALMT1, StALMT6, StALMT8, StALMT10, and StALMT12 had higher expression in the root of the potato cultivar Qingshu No. 9. After being subjected to Al3+ stress for 24 h, the expression of StALMT6 and StALMT10 in roots was evidently increased, displaying their decisive role in Al3+ toxicity. In addition, overexpression of StALMT6 and StALMT10 in Arabidopsis enhanced its tolerance to Al toxicity. These results indicate that StALMT6 and StALMT10 impart Al3+ resistance in the potato, establishing the foundation for further studies of the biological functions of these genes.

Genome-wide identification of the Pinus tabuliformis CONSTANS-like gene family and their potential roles in reproductive cone development

The CONSTANS-like (COL) genes, as a core transcription factor in the photoperiod regulation pathway, play a key role in plant reproduction development. However, their molecular characterization has rarely been studied in Pinus tabuliformis. Here, 10 PtCOL genes were identified in the P. tabuliformis genome and multiple sequence alignments have indicated that the PtCOL proteins contained highly conserved B-BOX1 and CCT domains. Sequence similarity analysis showed that PtCOL1 and PtCOL3 had the higher similarity with Norway spruce COLs (PaCOL2 and PaCOL1) and Arabidopsis COLs (AtCOL3, 4 and 5), respectively. Phylogeny and gene structure analyses revealed that PtCOLs were divided into three subgroups, each with identical or similar distributions of exons, introns, and motifs. Moreover, 10 PtCOLs were distributed on 6 chromosomes and PtCOL9 has syntenic gene pairs in both Ginkgo biloba and Sequoiadendron giganteum. Interestingly, in transcriptome profiles, most PtCOLs exhibited a diurnal oscillation pattern under both long (LD) and short (SD) day conditions. Additionally, PtCOLs were highly expressed in needles and female cones, and showed different spatial expression patterns. Among the ten PtCOLs, PtCOL1/3 heterologous overexpression Arabidopsis displayed a delayed-flowering phenotype under SD, indicating that they are likely to play a crucial role in the reproductive development. Additionally, PtCOL1 and PtCOL3 were not only capable of interacting with each other, but they were each capable of interacting with themselves. Furthermore, PtCOL1 and PtCOL3 were also involved in the MADS-box protein-protein interaction (PPI) network in P. tabuliformis cone development. Direct interactions of PtDAL11 with PtCOL1/3 impeded PtCOL1/3 translocation into the nucleus. In summary, this study provided comprehensive understanding for the functions of the PtCOL gene family and revealed their biological roles in the photoperiod-dependent P. tabuliformis cone development.

Genome-Wide Identification and Characterization of WRKY Transcription Factors in Betula platyphylla Suk. and Their Responses to Abiotic Stresses.

The WRKY transcription factor (TF) family is one the largest plant-specific transcription factor families. It has been proven to play significant roles in multiple plant biological processes, especially stress response. Although many WRKY TFs have been identified in various plant species, WRKYs in white birch (Betula platyphylla Suk.) remain to be studied. Here, we identified a total of 68 BpWRKYs, which could be classified into four main groups. The basic physiochemical properties of these TFs were analyzed using bioinformatics tools, including molecular weight, isoelectric point, chromosome location, and predicted subcellular localization. Most BpWRKYs were predicted to be located in the nucleus. Synteny analysis found 17 syntenic gene pairs among BpWRKYs and 52 syntenic gene pairs between BpWRKYs and AtWRKYs. The cis-acting elements in the promoters of BpWRKYs could be enriched in multiple plant biological processes, including stress response, hormone response, growth and development, and binding sites. Tissue-specific expression analysis using qRT-PCR showed that most BpWRKYs exhibited highest expression levels in the root. After ABA, salt (NaCl), or cold treatment, different BpWRKYs showed different expression patterns at different treatment times. Furthermore, the results of the Y2H assay proved the interaction between BpWRKY17 and a cold-responsive TF, BpCBF7. By transient expression assay, BpWRKY17 and BpWRKY67 were localized in the nucleus, consistent with the previous prediction. Our study hopes to shed light for research on WRKY TFs and plant stress response.

GET_PANGENES: calling pangenes from plant genome alignments confirms presence-absence variation

Crop pangenomes made from individual cultivar assemblies promise easy access to conserved genes, but genome content variability and inconsistent identifiers hamper their exploration. To address this, we define pangenes, which summarize a species coding potential and link back to original annotations. The protocol get_pangenes performs whole genome alignments (WGA) to call syntenic gene models based on coordinate overlaps. A benchmark with small and large plant genomes shows that pangenes recapitulate phylogeny-based orthologies and produce complete soft-core gene sets. Moreover, WGAs support lift-over and help confirm gene presence-absence variation. Source code and documentation: https://github.com/Ensembl/plant-scripts.

Integrating Genomic and Chromosomal Data: A Cytogenetic Study of Transancistrus santarosensis (Loricariidae: Hypostominae) with Characterization of a ZZ/ZW Sex Chromosome System.

The plecos (Loricariidae) fish represent a great model for cytogenetic investigations due to their variety of karyotypes, including diploid and polyploid genomes, and different types of sex chromosomes. In this study we investigate Transancistrus santarosensis a rare loricariid endemic to Ecuador, integrating cytogenetic methods with specimens' molecular identification by mtDNA, to describe the the species karyotype. We aim to verify whether sex chromosomes are cytologically identifiable and if they are associated with the accumulation of repetitive sequences present in other species of the family. The analysis of the karyotype (2n = 54 chromosomes) excludes recent centric fusion and pericentromeric inversion and suggests the presence of a ZZ/ZW sex chromosome system at an early stage of differentiation: the W chromosome is degenerated but is not characterized by the presence of differential sex-specific repetitive DNAs. Data indicate that although T. santarosensis has retained the ancestral diploid number of Loricariidae, it accumulated heterochromatin and shows non-syntenic ribosomal genes localization, chromosomal traits considered apomorphic in the family.

Characterization of the male-specific region containing the candidate sex-determining gene in Amur catfish (Silurus asotus) using third-generation- and pool-sequencing data

Amur catfish (Silurus asotus) is an ecologically and economically important fish species in Asia. Here, we assembled the female and male Amur catfish genomes, with genome sizes of 757.15 and 755.44 Mb, respectively, at the chromosome level using nanopore and Hi-C technologies. Consistent with the known diploid chromosome count, both genomes contained 29 chromosome-size scaffolds covering 98.80 and 98.73 % of the complete haplotypic assembly with scaffold N50 of 28.87 and 27.29 Mb, respectively. The female (n = 40) and male (n = 40) pools were re-sequenced. Comparative analysis of sequencing and re-sequencing data from both sexes confirmed the presence of an XX/XY sex determination system in Amur catfish and revealed Chr5 as the sex chromosome containing an approximately 400 kb Y-specific region (MSY). Gene annotation revealed a male-specific duplicate of amhr2, namely amhr2y, in MSY, which is male-specific in different wild populations and expressed only in the testes. Amur catfish shared partially syntenic MSY and amhr2y genes with the southern catfish (S. meridionalis, Chr24), which were located on different chromosomes. High sequence divergence between amhr2y and amhr2 and high sequence similarity with amhr2y were observed in both species. These results indicate the common origin of the sex-determining (SD) gene and transition of amhr2y in the two Silurus species. Accumulation of repetitive elements in the MSY of both species may be the main driver of the transition of amhr2y. Overall, our study provides valuable catfish genomic resources. Moreover, determination of amhr2y as the candidate SD gene in Amur catfish provides another example of amhr2 as the SD gene in fish.

Genome-Wide Identification and Characterization of the Phytochrome Gene Family in Peanut.

To investigate the potential role of phytochrome (PHY) in peanut growth and its response to environmental fluctuations, eight candidate AhPHY genes were identified via genome-wide analysis of cultivated peanut. These AhPHY polypeptides were determined to possess acidic and hydrophilic physiochemical properties and exhibit subcellular localization patterns consistent with residence in the nucleus and cytoplasm. Phylogenetic analysis revealed that the AhPHY gene family members were classified into three subgroups homologous to the PHYA/B/E progenitors of Arabidopsis. AhPHY genes within the same clade largely displayed analogous gene structure, conserved motifs, and phosphorylation sites. AhPHY exhibited symmetrical distribution across peanut chromosomes, with 7 intraspecific syntenic gene pairs in peanut, as well as 4 and 20 interspecific PHY syntenic gene pairs in Arabidopsis and soybean, respectively. A total of 42 cis-elements were predicted in AhPHY promoters, including elements implicated in phytohormone regulation, stress induction, physiology, and photoresponse, suggesting putative fundamental roles across diverse biological processes. Moreover, spatiotemporal transcript profiling of AhPHY genes in various peanut tissues revealed distinct expression patterns for each member, alluding to putative functional specialization. This study contributes novel insights into the classification, structure, molecular evolution, and expression profiles of the peanut phytochrome gene family, and also provides phototransduction gene resources for further mechanistic characterization.

Deep genome skimming reveals the hybrid origin of Pseudosasa gracilis (Poaceae: Bambusoideae)

Pseudosasa gracilis (Poaceae: Bambusoideae) is a temperate woody bamboo species endemic to South-central China with a narrow distribution. Previous phylogenetic studies revealed an unexpected, isolated phylogenetic position of Ps. gracilis. Here we conducted phylogenomic analysis by sampling populations of Ps. gracilis and its sympatric species Ps. nanunica and Sinosasa polytricha reflecting different genomic signals, by deep genome skimming. Integrating molecular evidence from chloroplast genes and genome-wide SNPs, we deciphered the phylogenetic relationships of Ps. gracilis. Both plastid and nuclear data indicate that Ps. gracilis is more closely related to Sinosasa, which is discordant with the taxonomic treatment. To further explore this molecular-morphological conflict, we screened 411 “perfect-copy” syntenic genes to reconstruct phylogenies using both the concatenation and coalescent methods. We observed extensive discordance between gene trees and the putative species tree. A significant hybridization event was detected based on 411 genes from the D subgenome, showing Ps. gracilis was a hybrid descendant between Sinosasa longiligulata and Ps. nanunica, with 63.56% and 36.44% inheritance probabilities of each parent. Moreover, introgression events were detected in the C subgenome between Ps. gracilis and S. polytricha in the same distribution region. Our findings suggest that sympatric hybridization and introgression play a crucial role in the origin of Ps. gracilis. By providing an empirical example of bamboo of hybrid origin using comprehensive analyses based on genomic data from different inheritance systems and morphological characters, our study represents a step forward in understanding of reticulate evolution of bamboos.

Genome-Wide Identification and Functional Analysis of the AP2/ERF Transcription Factor Family in Citrus Rootstock under Waterlogging Stress.

Citrus plants are sensitive to waterlogging, and the roots are the first plant organ affected by hypoxic stress. The AP2/ERF (APETALA2/ethylene-responsive element binding factors) can modulate plant growth and development. However, the information on AP2/ERF genes in citrus rootstock and their involvement in waterlogging conditions is limited. Previously, a rootstock cultivar, Citrus junos cv. Pujiang Xiangcheng was found to be highly tolerant to waterlogging stress. In this study, a total of 119 AP2/ERF members were identified in the C. junos genome. Conserved motif and gene structure analyses indicated the evolutionary conservation of PjAP2/ERFs. Syntenic gene analysis revealed 22 collinearity pairs among the 119 PjAP2/ERFs. The expression profiles under waterlogging stress showed differential expression of PjAP2/ERFs, of which, PjERF13 was highly expressed in both root and leaf. Furthermore, the heterologous expression of PjERF13 significantly enhanced the tolerance of transgenic tobacco to waterlogging stress. The overexpression of PjERF13 decreased the oxidative damage in the transgenic plants by reducing the H2O2 and MDA contents and increasing the antioxidant enzyme activities in the root and leaf. Overall, the current study provided basic information on the AP2/ERF family in the citrus rootstock and uncovered their potential function in positively regulating the waterlogging stress response.

Systematic Analysis of NRAMP Family Genes in Areca catechu and Its Response to Zn/Fe Deficiency Stress.

Areca catechu is a commercially important medicinal plant widely cultivated in tropical regions. The natural resistance-associated macrophage protein (NRAMP) is widespread in plants and plays critical roles in transporting metal ions, plant growth, and development. However, the information on NRAMPs in A. catechu is quite limited. In this study, we identified 12 NRAMPs genes in the areca genome, which were classified into five groups by phylogenetic analysis. Subcellular localization analysis reveals that, except for NRAMP2, NRAMP3, and NRAMP11, which are localized in chloroplasts, all other NRAMPs are localized on the plasma membrane. Genomic distribution analysis shows that 12 NRAMPs genes are unevenly spread on seven chromosomes. Sequence analysis shows that motif 1 and motif 6 are highly conserved motifs in 12 NRAMPs. Synteny analysis provided deep insight into the evolutionary characteristics of AcNRAMP genes. Among the A. catechu and the other three representative species, we identified a total of 19 syntenic gene pairs. Analysis of Ka/Ks values indicates that AcNRAMP genes are subjected to purifying selection in the evolutionary process. Analysis of cis-acting elements reveals that AcNRAMP genes promoter sequences contain light-responsive elements, defense- and stress-responsive elements, and plant growth/development-responsive elements. Expression profiling confirms distinct expression patterns of AcNRAMP genes in different organs and responses to Zn/Fe deficiency stress in leaves and roots. Taken together, our results lay a foundation for further exploration of the AcNRAMPs regulatory function in areca response to Fe and Zn deficiency.

Synteny Identifies Reliable Orthologs for Phylogenomics and Comparative Genomics of the Brassicaceae.

Large genomic data sets are becoming the new normal in phylogenetic research, but the identification of true orthologous genes and the exclusion of problematic paralogs is still challenging when applying commonly used sequencing methods such as target enrichment. Here, we compared conventional ortholog detection using OrthoFinder with ortholog detection through genomic synteny in a data set of 11 representative diploid Brassicaceae whole-genome sequences spanning the entire phylogenetic space. Then, we evaluated the resulting gene sets regarding gene number, functional annotation, and gene and species tree resolution. Finally, we used the syntenic gene sets for comparative genomics and ancestral genome analysis. The use of synteny resulted in considerably more orthologs and also allowed us to reliably identify paralogs. Surprisingly, we did not detect notable differences between species trees reconstructed from syntenic orthologs when compared with other gene sets, including the Angiosperms353 set and a Brassicaceae-specific target enrichment gene set. However, the synteny data set comprised a multitude of gene functions, strongly suggesting that this method of marker selection for phylogenomics is suitable for studies that value downstream gene function analysis, gene interaction, and network studies. Finally, we present the first ancestral genome reconstruction for the Core Brassicaceae which predating the Brassicaceae lineage diversification ∼25 million years ago.

Chromosome genome assembly for the meagre, Argyrosomus regius, reveals species adaptations and sciaenid sex-related locus evolution.

The meagre, Argyrosomus regius, has recently become a species of increasing economic interest for the Mediterranean aquaculture and there is ongoing work to boost production efficiency through selective breeding. Access to the complete genomic sequence will provide an essential resource for studying quantitative trait-associated loci and exploring the genetic diversity of different wild populations and aquaculture stocks in more detail. Here, we present the first complete genome for A. regius, produced through a combination of long and short read technologies and an efficient in-house developed pipeline for assembly and polishing. Scaffolding using previous linkage map data allowed us to reconstruct a chromosome level assembly with high completeness, complemented with gene annotation and repeat masking. The 696Mb long assembly has an N50 = 27.87Mb and an L50 = 12, with 92.85% of its length placed in 24 chromosomes. We use this new resource to study the evolution of the meagre genome and other Sciaenids, via a comparative analysis of 25 high-quality teleost genomes. Combining a rigorous investigation of gene duplications with base-wise conservation analysis, we identify candidate loci related to immune, fat metabolism and growth adaptations in the meagre. Following phylogenomic reconstruction, we show highly conserved synteny within Sciaenidae. In contrast, we report rapidly evolving syntenic rearrangements and gene copy changes in the sex-related dmrt1 neighbourhood in meagre and other members of the family. These novel genomic datasets and findings will add important new tools for aquaculture studies and greatly facilitate husbandry and breeding work in the species.

Navigating the CoGe Online Software Suite for Polyploidy Research.

The CoGe software suite at genomevolution.org hosts a number of tools that facilitate genomic research on plant and animal whole-genome multiplication-polyploidy. SynMap permits analysis and visualization of two-way syntenic dotplot alignments of genomes, includes many options and data/graphics download possibilities, and even permits three-genome synteny maps and interactive views. FractBias is a tool that operates within SynMap that permits calculation and graphic display of genome fragments (such as chromosomes) of one species mapped to another, displaying both blockwise homology depths and the extent of syntenic gene (syntelog) loss following polyploidy events. SynMap macrosynteny results can segue into the microsynteny tool GEvo, which provides genome-browser-like views of homologous genome blocks. CoGe FeatView allows call-up of given gene features already stored in the CoGe resource, and CoGeBlast permits searches for additional features that can be analyzed or downloadedfurther. Links from these tools can be fed into SynFind, which can find syntenic blocks surrounding a feature across multiple specified genomes while also simultaneously providing overall genome-wide syntenic depth calculations that can be interpreted to reflect polyploidy levels. Here, we describe basic use of these tools on the CoGe software suite.