Genome sequencing has become a routine task for biologists, but the challenge of gene structure annotation persists, impeding accurate genomic and genetic research. Here, we present a bioinformatics toolkit, SynGAP (Synteny-based Gene structure Annotation Polisher), which uses gene synteny information to accomplish precise and automated polishing of gene structure annotation of genomes. SynGAP offers exceptional capabilities in the improvement of gene structure annotation quality and the profiling of integrative gene synteny between species. Furthermore, an expression variation index is designed for comparative transcriptomics analysis to explore candidate genes responsible for the development of distinct traits observed in phylogenetically related species.

BACKGROUND: Obtaining de novo chromosome-level genome assemblies greatly enhances conservation and evolutionary biology studies. For many research teams, long-read sequencing technologies (that produce highly contiguous assemblies) remain unaffordable or unpractical. For the groups that display high synteny conservation, these limitations can be overcome by a reference-guided assembly using a close relative genome. Among chelonians, tortoises (Testudinidae) are considered one of the most endangered taxa, which calls for more genomic resources. Here we make the most of high synteny conservation in chelonians to produce the first chromosome-level genome assembly of the genus Testudo with one of the most iconic tortoise species in the Mediterranean basin: Testudo graeca.
RESULTS: We used high-quality, paired-end Illumina sequences to build a reference-guided assembly with the chromosome-level reference of Gopherus evgoodei. We reconstructed a 2.29 Gb haploid genome with a scaffold N50 of 107.598 Mb and 5.37% gaps. We sequenced 25,998 protein-coding genes, and identified 41.2% of the assembly as repeats. Demographic history reconstruction based on the genome revealed two events (population decline and recovery) that were consistent with previously suggested phylogeographic patterns for the species. This outlines the value of such reference-guided assemblies for phylogeographic studies.
CONCLUSIONS: Our results highlight the value of using close relatives to produce de novo draft assemblies in species where such resources are unavailable. Our annotated genome of T. graeca paves the way to delve deeper into the species\' evolutionary history and provides a valuable resource to enhance direct conservation efforts on their threatened populations.

Adenosine kinase (ADK) is a key enzyme widely distributed in plants, playing an important role in maintaining cellular energy homeostasis and regulating plant growth, development, and responses to environmental stresses. However, research on ADK genes in cotton (Gossypium hirsutum), an economically significant crop, has been limited. This study identified 92 ADK genes from four cotton species (G. arboreum, G. raimondii, G. hirsutum, and G. barbadense) using HMMER and Local BLASTP methods and classified them into six groups. Chromosomal localization revealed a random distribution of ADK genes in G. hirsutum, with 13 genes located on the At subgenome and 14 genes on the Dt subgenome. Gene structure analysis showed consistency in exon-intron organization within subgroups, while conserved motif analysis identified subgroup-specific motifs, indicating functional diversity. Synteny and collinearity mapping analysis revealed that the primary expansion mechanisms of the ADK gene family in cotton are polyploidy and segmental duplication. Cis-regulatory elements in GhADK promoters were classified into light response, hormone response, developmental regulation, and stress response. We also analyzed the expression patterns of GhADK genes under a low temperature (4 °C) and drought conditions. Most GhADK genes responded to cold stress with different expression patterns, indicating their roles in rapid response and long-term cold adaptation. Under drought stress, expression patterns varied, with some genes showing sustained high expression levels. The qRT-PCR validation of transcriptomic data confirmed the stress-induced expression patterns of selected GhADK genes. Functional analysis through the VIGS silencing of GhADK25 demonstrated its importance in cold and drought stress responses, with silencing resulting in poor growth under stress, highlighting its significance in stress tolerance. This study provides a basis for further understanding the evolutionary relationships and functions of the cotton ADK gene family.

The fructose-1,6-bisphosphate aldolase (FBA) gene family exists in higher plants, with the genes of this family playing significant roles in plant growth and development, as well as response to abiotic stresses. However, systematic reports on the FBA gene family and its functions in cucumber are lacking. In this study, we identified five cucumber FBA genes, named CsFBA1-5, that are distributed randomly across chromosomes. Phylogenetic analyses involving these cucumber FBAs, alongside eight Arabidopsis FBA proteins and eight tomato FBA proteins, were conducted to assess their homology. The CsFBAs were grouped into two clades. We also analyzed the physicochemical properties, motif composition, and gene structure of the cucumber FBAs. This analysis highlighted differences in the physicochemical properties and revealed highly conserved domains within the CsFBA family. Additionally, to explore the evolutionary relationships of the CsFBA family further, we constructed comparative syntenic maps with Arabidopsis and tomato, which showed high homology but only one segmental duplication event within the cucumber genome. Expression profiles indicated that the CsFBA gene family is responsive to various abiotic stresses, including low temperature, heat, and salt. Taken together, the results of this study provide a theoretical foundation for understanding the evolution of and future research into the functional characterization of cucumber FBA genes during plant growth and development.

SINA (Seven in absentia) E3 ubiquitin ligases are a family of RING (really interesting new gene) E3 ubiquitin ligases, and they play a crucial role in regulating plant growth and development, hormone response, and abiotic and biotic stress. However, there is little research on the SINA gene family in U. rhynchophylla. In this study, a total of 10 UrSINA genes were identified from the U. rhynchophylla genome. The results of multiple sequence alignments and chromosomal locations show that 10 UrSINA genes were unevenly located on 22 chromosomes, and each UrSINA protein contained a SINA domain at the N-terminal and RING domains at the C-terminal. Synteny analysis showed that there are no tandem duplication gene pairs and there are four segmental gene pairs in U. rhynchophylla, contributing to the expansion of the gene family. Furthermore, almost all UrSINA genes contained the same gene structure, with three exons and two introns, and there were many cis-acting elements relating to plant hormones, light responses, and biotic and abiotic stress. The results of qRT-PCR show that most UrSINA genes were expressed in stems, with the least expression in roots; meanwhile, most UrSINA genes and key enzyme genes were responsive to ABA and MeJA hormones with overlapping but different expression patterns. Co-expression analysis showed that UrSINA1 might participate in the TIA pathway under ABA treatment, and UrSINA5 and UrSINA6 might participate in the TIA pathway under MeJA treatment. The mining of UrSINA genes in the U. rhynchophylla provided novel information for understanding the SINA gene and its function in plant secondary metabolites, growth, and development.

Dendrobium loddigesii is a precious traditional Chinese medicine with high medicinal and ornamental value. However, the characterization of its mitochondrial genome is still pending. Here, we assembled the complete mitochondrial genome of D. loddigesii and discovered that its genome possessed a complex multi-chromosome structure. The mitogenome of D. loddigesii consisted of 17 circular subgenomes, ranging in size from 16,323 bp to 56,781 bp. The total length of the mitogenome was 513,356 bp, with a GC content of 43.41%. The mitogenome contained 70 genes, comprising 36 protein-coding genes (PCGs), 31 tRNA genes, and 3 rRNA genes. Furthermore, we detected 403 repeat sequences as well as identified 482 RNA-editing sites and 8154 codons across all PCGs. Following the sequence similarity analysis, 27 fragments exhibiting homology to both the mitogenome and chloroplast genome were discovered, accounting for 9.86% mitogenome of D. loddigesii. Synteny analysis revealed numerous sequence rearrangements in D. loddigesii and the mitogenomes of related species. Phylogenetic analysis strongly supported that D. loddigesii and D. Amplum formed a single clade with 100% bootstrap support. The outcomes will significantly augment the orchid mitochondrial genome database, offering profound insights into Dendrobium\'s intricate mitochondrial genome architecture.

OBJECTIVE: The history of angiosperms is marked by repeated rounds of ancient whole-genome duplications (WGDs). Here we used state-of-the-art methods to provide an up-to-date view of the distribution of WGDs in the history of angiosperms that considers both uncertainty introduced by different WGD inference methods and different underlying species-tree hypotheses.
METHODS: We used the distribution synonymous divergences (Ks) of paralogs and orthologs from transcriptomic and genomic data to infer and place WGDs across two hypothesized angiosperm phylogenies. We further tested these WGD hypotheses with syntenic inferences and Bayesian models of duplicate gene gain and loss.
RESULTS: The predicted number of WGDs in the history of angiosperms (~170) based on the current taxon sampling is largely similar across different inference methods, but varies in the precise placement of WGDs on the phylogeny. Ks-based methods often yield alternative hypothesized WGD placements due to variation in substitution rates among lineages. Phylogenetic models of duplicate gene gain and loss are more robust to topological variation. However, errors in species-tree inference can still produce spurious WGD hypotheses, regardless of method used.
CONCLUSIONS: Here we showed that different WGD inference methods largely agree on an average of 3.5 WGD in the history of individual angiosperm species. However, the precise placement of WGDs on the phylogeny is subject to the WGD inference method and tree topology. As researchers continue to test hypotheses regarding the impacts ancient WGDs have on angiosperm evolution, it is important to consider the uncertainty of the phylogeny as well as WGD inference methods.

Sex chromosomes are evolutionarily labile in many animals and sometimes fuse with autosomes, creating so-called neo-sex chromosomes. Fusions between sex chromosomes and autosomes have been proposed to reduce sexual conflict and to promote adaptation and reproductive isolation among species. Recently, advances in genomics have fuelled the discovery of such fusions across the tree of life. Here, we discovered multiple fusions leading to neo-sex chromosomes in the sapho subclade of the classical adaptive radiation of Heliconius butterflies. Heliconius butterflies generally have 21 chromosomes with very high synteny. However, the five Heliconius species in the sapho subclade show large variation in chromosome number ranging from 21 to 60. We find that the W chromosome is fused with chromosome 4 in all of them. Two sister species pairs show subsequent fusions between the W and chromosomes 9 or 14, respectively. These fusions between autosomes and sex chromosomes make Heliconius butterflies an ideal system for studying the role of neo-sex chromosomes in adaptive radiations and the degeneration of sex chromosomes over time. Our findings emphasize the capability of short-read resequencing to detect genomic signatures of fusion events between sex chromosomes and autosomes even when sex chromosomes are not explicitly assembled.

Islands are crucial evolutionary hotspots, providing unique opportunities for differentiation of novel biodiversity and long-term segregation of endemic species. Islands are also fragile ecosystems, where biodiversity is more exposed to environmental and anthropogenic pressures than on continents. The Ponza grayling, Hipparchia sbordonii, is an endemic butterfly species that is currently found only in two tiny islands of the Pontine archipelago, off the coast of Italy, occupying an area smaller than 10 km2. It has been classified as Endangered (IUCN) because of the extremely limited area of occurrence, population fragmentation, and the recent demographic decline. Thanks to a combination of different assemblers of long and short genomic reads, bulk transcriptome RNAseq, and synteny analysis with phylogenetically close butterflies, we produced a highly contiguous, chromosome-scale annotated reference genome for the Ponza grayling, including 28 autosomes and the Z sexual chromosomes. The final assembly spanned 388.61 Gb with a contig N50 of 14.5 Mb and a BUSCO completeness score of 98.5%. Synteny analysis using four other butterfly species revealed high collinearity with Hipparchia semele and highlighted 10 intrachromosomal inversions longer than 10 kb, of which two appeared on the lineage leading to H. sbordonii. Our results show that a chromosome-scale reference genome is attainable also when chromatin conformation data may be impractical or present specific technical challenges. The high-quality genomic resource for H. sbordonii opens up new opportunities for the accurate assessment of genetic diversity and genetic load and for the investigations of the genomic novelties characterizing the evolutionary path of this endemic island species.

GAI-RGA-and-SCR (GRAS) transcription factors can regulate many biological processes such as plant growth and development and stress defense, but there are few related studies in sugar beet. Salt stress can seriously affect the yield and quality of sugar beet (Beta vulgaris). Therefore, this study used bioinformatics methods to identify GRAS transcription factors in sugar beet and analyzed their structural characteristics, evolutionary relationships, regulatory networks and salt stress response patterns. A total of 28 BvGRAS genes were identified in the whole genome of sugar beet, and the sequence composition was relatively conservative. According to the topology of the phylogenetic tree, BvGRAS can be divided into nine subfamilies: LISCL, SHR, PAT1, SCR, SCL3, LAS, SCL4/7, HAM and DELLA. Synteny analysis showed that there were two pairs of fragment replication genes in the BvGRAS gene, indicating that gene replication was not the main source of BvGRAS family members. Regulatory network analysis showed that BvGRAS could participate in the regulation of protein interaction, material transport, redox balance, ion homeostasis, osmotic substance accumulation and plant morphological structure to affect the tolerance of sugar beet to salt stress. Under salt stress, BvGRAS and its target genes showed an up-regulated expression trend. Among them, BvGRAS-15, BvGRAS-19, BvGRAS-20, BvGRAS-21, LOC104892636 and LOC104893770 may be the key genes for sugar beet\'s salt stress response. In this study, the structural characteristics and biological functions of BvGRAS transcription factors were analyzed, which provided data for the further study of the molecular mechanisms of salt stress and molecular breeding of sugar beet.