Dfam: a database of repetitive DNA based on profile hidden Markov models

VIBES: a workflow for annotating and visualizing viral sequences integrated into bacterial genomes

Bacteriophages are viruses that infect bacteria. Many bacteriophages integrate their genomes into the bacterial chromosome and become prophages. Prophages may substantially burden or benefit host bacteria fitness, acting in some cases as parasites and in others as mutualists. Some prophages have been demonstrated to increase host virulence. The increasing ease of bacterial genome sequencing provides an opportunity to deeply explore prophage prevalence and insertion sites. Here we present VIBES (Viral Integrations in Bacterial genomES), a workflow intended to automate prophage annotation in complete bacterial genome sequences. VIBES provides additional context to prophage annotations by annotating bacterial genes and viral proteins in user-provided bacterial and viral genomes. The VIBES pipeline is implemented as a Nextflow-driven workflow, providing a simple, unified interface for execution on local, cluster and cloud computing environments. For each step of the pipeline, a container including all necessary software dependencies is provided. VIBES produces results in simple tab-separated format and generates intuitive and interactive visualizations for data exploration. Despite VIBES's primary emphasis on prophage annotation, its generic alignment-based design allows it to be deployed as a general-purpose sequence similarity search manager. We demonstrate the utility of the VIBES prophage annotation workflow by searching for 178 Pf phage genomes across 1072 Pseudomonas spp. genomes.

A high-quality chromosome-scale genome assembly of blood orange, an important pigmented sweet orange variety

Blood orange (BO) is a rare red-fleshed sweet orange (SWO) with a high anthocyanin content and is associated with numerous health-related benefits. Here, we reported a high-quality chromosome-scale genome assembly for Neixiu (NX) BO, reaching 336.63 Mb in length with contig and scaffold N50 values of 30.6 Mb. Furthermore, 96% of the assembled sequences were successfully anchored to 9 pseudo-chromosomes. The genome assembly also revealed the presence of 37.87% transposon elements and 7.64% tandem repeats, and the annotation of 30,395 protein-coding genes. A high level of genome synteny was observed between BO and SWO, further supporting their genetic similarity. The speciation event that gave rise to the Citrus species predated the duplication event found within them. The genome-wide variation between NX and SWO was also compared. This first high-quality BO genome will serve as a fundamental basis for future studies on functional genomics and genome evolution.

Chromosome-level genome assembly of Prunella vulgaris L. provides insights into pentacyclic triterpenoid biosynthesis

Prunella vulgaris is one of the bestselling and widely used medicinal herbs. It is recorded as an ace medicine for cleansing and protecting the liver in Chinese Pharmacopoeia and has been used as the main constitutions of many herbal tea formulas in China for centuries. It is also a traditional folk medicine in Europe and other countries of Asia. Pentacyclic triterpenoids are a major class of bioactive compounds produced in P. vulgaris. However, their biosynthetic mechanism remains to be elucidated. Here, we report a chromosome-level reference genome of P. vulgaris using an approach combining Illumina, ONT, and Hi-C technologies. It is 671.95 Mb in size with a scaffold N50 of 49.10 Mb and a complete BUSCO of 98.45%. About 98.31% of the sequence was anchored into 14 pseudochromosomes. Comparative genome analysis revealed a recent WGD in P. vulgaris. Genome-wide analysis identified 35 932 protein-coding genes (PCGs), of which 59 encode enzymes involved in 2,3-oxidosqualene biosynthesis. In addition, 10 PvOSC, 358 PvCYP, and 177 PvUGT genes were identified, of which five PvOSCs, 25 PvCYPs, and 9 PvUGTs were predicted to be involved in the biosynthesis of pentacyclic triterpenoids. Biochemical activity assay of PvOSC2, PvOSC4, and PvOSC6 recombinant proteins showed that they were mixed amyrin synthase (MAS), lupeol synthase (LUS), and β-amyrin synthase (BAS), respectively. The results provide a solid foundation for further elucidating the biosynthetic mechanism of pentacyclic triterpenoids in P. vulgaris.

Differential selection of yield and quality traits has shaped genomic signatures of cowpea domestication and improvement

Cowpeas (tropical legumes) are important in ensuring food and nutritional security in developing countries, especially in sub-Saharan Africa. Herein, we report two high-quality genome assemblies of grain and vegetable cowpeas and we re-sequenced 344 accessions to characterize the genomic variations landscape. We identified 39 loci for ten important agronomic traits and more than 541 potential loci that underwent selection during cowpea domestication and improvement. In particular, the synchronous selections of the pod-shattering loci and their neighboring stress-relevant loci probably led to the enhancement of pod-shattering resistance and the compromise of stress resistance during the domestication from grain to vegetable cowpeas. Moreover, differential selections on multiple loci associated with pod length, grain number per pod, seed weight, pod and seed soluble sugars, and seed crude proteins shaped the yield and quality diversity in cowpeas. Our findings provide genomic insights into cowpea domestication and improvement footprints, enabling further genome-informed cultivar improvement of cowpeas.

De novo genome assembly and population genomics of a shrub tree Barthea barthei (Hance) krass provide insights into the adaptive color variations

Flower color is a classic example of an ecologically important trait under selection in plants. Understanding the genetic mechanisms underlying shifts in flower color can provide key insights into ecological speciation. In this study, we investigated the genetic basis of flower color divergence in Barthea barthei, a shrub tree species exhibiting natural variation in flower color. We assembled a high-quality genome assembly for B. barthei with a contig N50 of 2.39 Mb and a scaffold N50 of 16.21 Mb. The assembly was annotated with 46,430 protein-coding genes and 1,560 non-coding RNAs. Genome synteny analysis revealed two recent tetraploidization events in B. barthei, estimated to have occurred at approximately 17 and 63 million years ago. These tetraploidization events resulted in massive duplicated gene content, with over 70% of genes retained in collinear blocks. Gene family members of the core regulators of the MBW complex were significantly expanded in B. barthei compared to Arabidopsis, suggesting that these duplications may have provided raw genetic material for the evolution of novel regulatory interactions and the diversification of anthocyanin pigmentation. Transcriptome profiling of B. barthei flowers revealed differential expression of 9 transcription factors related to anthocyanin biosynthesis between the two ecotypes. Six of these differentially expressed transcription factors were identified as high-confidence candidates for adaptive evolution based on positive selection signals. This study provides insights into the genetic basis of flower color divergence and the evolutionary mechanisms underlying ecological adaptation in plants.

Gene gain and loss from the Asian corn borer W chromosome

BACKGROUND

Sex-limited chromosomes Y and W share some characteristics, including the degeneration of protein-coding genes, enrichment of repetitive elements, and heterochromatin. However, although many studies have suggested that Y chromosomes retain genes related to male function, far less is known about W chromosomes and whether they retain genes related to female-specific function.

RESULTS

Here, we built a chromosome-level genome assembly of the Asian corn borer, Ostrinia furnacalis Guenée (Lepidoptera: Crambidae, Pyraloidea), an economically important pest in corn, from a female, including both the Z and W chromosome. Despite deep conservation of the Z chromosome across Lepidoptera, our chromosome-level W assembly reveals little conservation with available W chromosome sequence in related species or with the Z chromosome, consistent with a non-canonical origin of the W chromosome. The W chromosome has accumulated significant repetitive elements and experienced rapid gene gain from the remainder of the genome, with most genes exhibiting pseudogenization after duplication to the W. The genes that retain significant expression are largely enriched for functions in DNA recombination, the nucleosome, chromatin, and DNA binding, likely related to meiotic and mitotic processes within the female gonad.

CONCLUSIONS

Overall, our chromosome-level genome assembly supports the non-canonical origin of the W chromosome in O. furnacalis, which experienced rapid gene gain and loss, with the retention of genes related to female-specific function.

Study of Dispersed Repeats in the Cyanidioschyzon merolae Genome

In this study, we applied the iterative procedure (IP) method to search for families of highly diverged dispersed repeats in the genome of Cyanidioschyzon merolae, which contains over 16 million bases. The algorithm included the construction of position weight matrices (PWMs) for repeat families and the identification of more dispersed repeats based on the PWMs using dynamic programming. The results showed that the C. merolae genome contained 20 repeat families comprising a total of 33,938 dispersed repeats, which is significantly more than has been previously found using other methods. The repeats varied in length from 108 to 600 bp (522.54 bp in average) and occupied more than 72% of the C. merolae genome, whereas previously identified repeats, including tandem repeats, have been shown to constitute only about 28%. The high genomic content of dispersed repeats and their location in the coding regions suggest a significant role in the regulation of the functional activity of the genome.

Pangenome reconstruction in rats enhances genotype-phenotype mapping and novel variant discovery

The HXB/BXH family of recombinant inbred rat strains is a unique genetic resource that has been extensively phenotyped over 25 years, resulting in a vast dataset of quantitative molecular and physiological phenotypes. We built a pangenome graph from 10x Genomics Linked-Read data for 31 recombinant inbred rats to study genetic variation and association mapping. The pangenome includes 0.2Gb of sequence that is not present the reference mRatBN7.2, confirming the capture of substantial additional variation. We validated variants in challenging regions, including complex structural variants resolving into multiple haplotypes. Phenome-wide association analysis of validated SNPs uncovered variants associated with glucose/insulin levels and hippocampal gene expression. We propose an interaction between Pirl1l1, chromogranin expression, TNF-α levels, and insulin regulation. This study demonstrates the utility of linked-read pangenomes for comprehensive variant detection and mapping phenotypic diversity in a widely used rat genetic reference panel.

LINE-1 mRNA 3' end dynamics shape its biology and retrotransposition potential

LINE-1 (L1) retrotransposons are mobile genetic elements that create new genomic insertions by a copy-paste mechanism involving L1 RNA/RNP intermediates. L1 encodes two ORFs, of which L1-ORF2p nicks genomic DNA and reverse transcribes L1 mRNA using the nicked DNA as a primer which base-pairs with poly(A) tail of L1 mRNA. To better understand the importance of non-templated L1 3' ends' dynamics and the interplay between L1 3' and 5' ends, we investigated the effects of genomic knock-outs and temporal knock-downs of XRN1, DCP2, and other factors. We hypothesized that in the absence of XRN1, the major 5'→3' exoribonuclease, there would be more L1 mRNA and retrotransposition. Conversely, we observed that loss of XRN1 decreased L1 retrotransposition. This occurred despite slight stabilization of L1 mRNA, but with decreased L1 RNP formation. Similarly, loss of DCP2, the catalytic subunit of the decapping complex, lowered retrotransposition despite increased steady-state levels of L1 proteins. In both XRN1 and DCP2 depletions we observed shortening of L1 3' poly(A) tails and their increased uridylation by TUT4/7. We explain the observed reduction of L1 retrotransposition by the changed qualities of non-templated L1 mRNA 3' ends demonstrating the important role of L1 3' end dynamics in L1 biology.

WAS IT A MATch I SAW? Approximate palindromes lead to overstated false match rates in benchmarks using reversed sequences

Background

Software for labeling biological sequences typically produces a theory-based statistic for each match (the E-value) that indicates the likelihood of seeing that match's score by chance. E-values accurately predict false match rate for comparisons of random (shuffled) sequences, and thus provide a reasoned mechanism for setting score thresholds that enable high sensitivity with low expected false match rate. This threshold-setting strategy is challenged by real biological sequences, which contain regions of local repetition and low sequence complexity that cause excess matches between non-homologous sequences. Knowing this, tool developers often develop benchmarks that use realistic-seeming decoy sequences to explore empirical tradeoffs between sensitivity and false match rate. A recent trend has been to employ reversed biological sequences as realistic decoys, because these preserve the distribution of letters and the existence of local repeats, while disrupting the original sequence's functional properties. However, we and others have observed that sequences appear to produce high scoring alignments to their reversals with surprising frequency, leading to overstatement of false match risk that may negatively affect downstream analysis.

Results

We demonstrate that an alignment between a sequence S and its (possibly mutated) reversal tends to produce higher scores than alignment between truly unrelated sequences, even when S is a shuffled string with no notable repetitive or low-complexity regions. This phenomenon is due to the unintuitive fact that (even randomly shuffled) sequences contain palindromes that are on average longer than the longest common substrings (LCS) shared between permuted variants of the same sequence. Though the expected palindrome length is only slightly larger than the expected LCS, the distribution of alignment scores involving reversed sequences is strongly right-shifted, leading to greatly increased frequency of high-scoring alignments to reversed sequences.

Impact

Overestimates of false match risk can motivate unnecessarily high score thresholds, leading to potentially reduced true match sensitivity. Also, when tool sensitivity is only reported up to the score of the first matched decoy sequence, a large decoy set consisting of reversed sequences can obscure sensitivity differences between tools. As a result of these observations, we advise that reversed biological sequences be used as decoys only when care is taken to remove positive matches in the original (un-reversed) sequences, or when overstatement of false labeling is not a concern. Though the primary focus of the analysis is on sequence annotation, we also demonstrate that the prevalence of internal palindromes may lead to an overstatement of the rate of false labels in protein identification with mass spectrometry.

Chromosome-level genome of putative autohexaploid Actinidia deliciosa provides insights into polyploidisation and evolution

Actinidia ('Mihoutao' in Chinese) includes species with complex ploidy, among which diploid Actinidia chinensis and hexaploid Actinidia deliciosa are economically and nutritionally important fruit crops. Actinidia deliciosa has been proposed to be an autohexaploid (2n = 174) with diploid A. chinensis (2n = 58) as the putative parent. A CCS-based assembly anchored to a high-resolution linkage map provided a chromosome-resolved genome for hexaploid A. deliciosa yielded a 3.91-Gb assembly of 174 pseudochromosomes comprising 29 homologous groups with 6 members each, which contain 39 854 genes with an average of 4.57 alleles per gene. Here we provide evidence that much of the hexaploid genome matches diploid A. chinensis; 95.5% of homologous gene pairs exhibited >90% similarity. However, intragenome and intergenome comparisons of synteny indicate chromosomal changes. Our data, therefore, indicate that if A. deliciosa is an autoploid, chromosomal rearrangement occurred following autohexaploidy. A highly diversified pattern of gene expression and a history of rapid population expansion after polyploidisation likely facilitated the adaptation and niche differentiation of A. deliciosa in nature. The allele-defined hexaploid genome of A. deliciosa provides new genomic resources to accelerate crop improvement and to understand polyploid genome evolution.

The Idesia polycarpa genome provides insights into its evolution and oil biosynthesis

The deciduous tree Idesia polycarpa can provide premium edible oil with high polyunsaturated fatty acid contents. Here, we generate its high-quality reference genome, which is ∼1.21 Gb, comprising 21 pseudochromosomes and 42,086 protein-coding genes. Phylogenetic and genomic synteny analyses show that it diverged with Populus trichocarpa about 16.28 million years ago. Notably, most fatty acid biosynthesis genes are not only increased in number in its genome but are also highly expressed in the fruits. Moreover, we identify, through genome-wide association analysis and RNA sequencing, the I. polycarpa SUGAR TRANSPORTER 5 (IpSTP5) gene as a positive regulator of high oil accumulation in the fruits. Silencing of IpSTP5 by virus-induced gene silencing causes a significant reduction of oil content in the fruits, suggesting it has the potential to be used as a molecular marker to breed the high-oil-content cultivars. Our results collectively lay the foundation for breeding the elite cultivars of I. polycarpa.

Genome sequencing and analysis of black flounder (Paralichthys orbignyanus) reveals new insights into Pleuronectiformes genomic size and structure

Black flounder (Paralichthys orbignyanus, Pleuronectiformes) is a commercially significant marine fish with promising aquaculture potential in Argentina. Despite extensive studies on Black flounder aquaculture, its limited genetic information available hampers the crucial role genetics plays in the development of this activity. In this study, we first employed Illumina sequencing technology to sequence the entire genome of Black flounder. Utilizing two independent libraries-one from a female and another from a male-with 150 bp paired-end reads, a mean insert length of 350 bp, and over 35 X-fold coverage, we achieved assemblies resulting in a genome size of ~ 538 Mbp. Analysis of the assemblies revealed that more than 98% of the core genes were present, with more than 78% of them having more than 50% coverage. This indicates a somehow complete and accurate genome at the coding sequence level. This genome contains 25,231 protein-coding genes, 445 tRNAs, 3 rRNAs, and more than 1,500 non-coding RNAs of other types. Black flounder, along with pufferfishes, seahorses, pipefishes, and anabantid fish, displays a smaller genome compared to most other teleost groups. In vertebrates, the number of transposable elements (TEs) is often correlated with genome size. However, it remains unclear whether the sizes of introns and exons also play a role in determining genome size. Hence, to elucidate the potential factors contributing to this reduced genome size, we conducted a comparative genomic analysis between Black flounder and other teleost orders to determine if the small genomic size could be explained by repetitive elements or gene features, including the whole genome genes and introns sizes. We show that the smaller genome size of flounders can be attributed to several factors, including changes in the number of repetitive elements, and decreased gene size, particularly due to lower amount of very large and small introns. Thus, these components appear to be involved in the genome reduction in Black flounder. Despite these insights, the full implications and potential benefits of genome reduction in Black flounder for reproduction and aquaculture remain incompletely understood, necessitating further research.

A chromosome-level genome assembly of yellow stem borer (Scirpophaga incertulas)

The yellow stem borer Scirpophaga incertulas is the dominant pest of rice in tropical Asia. However, the lack of genomic resources makes it difficult to understand their invasiveness and ecological adaptation. A high-quality chromosome-level genome of S. incertulas, a monophagous rice pest, was assembled by combining Illumina short reads, PacBio HiFi long sequencing, and Hi-C scaffolding technology. The final genome size was 695.65 Mb, with a scaffold N50 of 28.02 Mb, and 93.50% of the assembled sequences were anchored to 22 chromosomes. BUSCO analysis demonstrated that this genome assembly had a high level of completeness, with 97.65% gene coverage. A total of 14,850 protein-coding genes and 366.98 Mb of transposable elements were identified. In addition, comparative genomic analyses indicated that chemosensory processes and detoxification capacity may play critical roles in the specialized host preference of S. incertulas. In summary, the chromosome-level genome assembly of S. incertulas provides a valuable genetic resource for understanding the biological characteristics of its invasiveness and developing an efficient management strategy.

Ancestral genome reconstruction enhances transposable element annotation by identifying degenerate integrants

Growing evidence indicates that transposable elements (TEs) play important roles in evolution by providing genomes with coding and non-coding sequences. Identification of TE-derived functional elements, however, has relied on TE annotations in individual species, which limits its scope to relatively intact TE sequences. Here, we report a novel approach to uncover previously unannotated degenerate TEs (degTEs) by probing multiple ancestral genomes reconstructed from hundreds of species. We applied this method to the human genome and achieved a 10.8% increase in coverage over the most recent annotation. Further, we discovered that degTEs contribute to various cis-regulatory elements and transcription factor binding sites, including those of a known TE-controlling family, the KRAB zinc-finger proteins. We also report unannotated chimeric transcripts between degTEs and human genes expressed in embryos. This study provides a novel methodology and a freely available resource that will facilitate the investigation of TE co-option events on a full scale.

A chromosome-level genome assembly of the Rhus gall aphid Schlechtendalia chinensis provides insight into the endogenization of Parvovirus-like DNA sequences

The Rhus gall aphid, Schlechtendalia chinensis, feeds on its primary host plant Rhus chinensis to induce galls, which have economic importance in medicines and the food industry. Rhus gall aphids have a unique life cycle and are economically beneficial but there is huge gap in genomic information about this group of aphids. Schlechtendalia chinensis induces rich-tannin galls on its host plant and is emerging as a model organism for both commercial applications and applied research in the context of gall production by insects. Here, we generated a high-quality chromosome-level assembly for the S. chinensis genome, enabling the comparison between S. chinensis and non-galling aphids. The final genome assembly is 344.59 Mb with 91.71% of the assembled sequences anchored into 13 chromosomes. We predicted 15,013 genes, of which 14,582 (97.13%) coding genes were annotated, and 99% of the predicted genes were anchored to the 13 chromosomes. This assembly reveals the endogenization of parvovirus-related DNA sequences (PRDs) in the S. chinensis genome, which could play a role in environmental adaptations. We demonstrated the characterization and classification of cytochrome P450s in the genome assembly, which are functionally crucial for sap-feeding insects and have roles in detoxification and insecticide resistance. This genome assembly also revealed the whole genome duplication events in S. chinensis, which can be considered in comparative evolutionary analysis. Our work represents a reference genome for gall-forming aphids that could be used for comparative genomic studies between galling and non-galling aphids and provides the first insight into the endogenization of PRDs in the genome of galling aphids. It also provides novel genetic information for future research on gall-formation and insect-plant interactions.

Populus cathayana genome and population resequencing provide insights into its evolution and adaptation

Populus cathayana Rehder, an indigenous poplar species of ecological and economic importance, is widely distributed in a high-elevation range from southwest to northeast China. Further development of this species as a sustainable poplar resource has been hindered by a lack of genome information the at the population level. Here, we produced a chromosome-level genome assembly of P. cathayana, covering 406.55 Mb (scaffold N50 = 20.86 Mb) and consisting of 19 chromosomes, with 35 977 protein-coding genes. Subsequently, we made a genomic variation atlas of 438 wild individuals covering 36 representative geographic areas of P. cathayana, which were divided into four geographic groups. It was inferred that the Northwest China regions served as the genetic diversity centers and a population bottleneck happened during the history of P. cathayana. By genotype-environment association analysis, 947 environment-association loci were significantly associated with temperature, solar radiation, precipitation, and altitude variables. We identified local adaptation genes involved in DNA repair and UV radiation response, among which UVR8, HY5, and CUL4 had key roles in high-altitude adaptation of P. cathayana. Predictions of adaptive potential under future climate conditions showed that P. cathayana populations in areas with drastic climate change were anticipated to have greater maladaptation risk. These results provide comprehensive insights for understanding wild poplar evolution and optimizing adaptive potential in molecular breeding.

Sensitive and error-tolerant annotation of protein-coding DNA with BATH

We present BATH, a tool for highly sensitive annotation of protein-coding DNA based on direct alignment of that DNA to a database of protein sequences or profile hidden Markov models (pHMMs). BATH is built on top of the HMMER3 code base, and simplifies the annotation workflow for pHMM-based annotation by providing a straightforward input interface and easy-to-interpret output. BATH also introduces novel frameshift-aware algorithms to detect frameshift-inducing nucleotide insertions and deletions (indels). BATH matches the accuracy of HMMER3 for annotation of sequences containing no errors, and produces superior accuracy to all tested tools for annotation of sequences containing nucleotide indels. These results suggest that BATH should be used when high annotation sensitivity is required, particularly when frameshift errors are expected to interrupt protein-coding regions, as is true with long read sequencing data and in the context of pseudogenes.

Polyploid Genome Assembly Provides Insights into Morphological Development and Ascorbic Acid Accumulation of Sauropus androgynus

Sauropus androgynus (S. androgynus) (2n = 4x = 52) is one of the most popular functional leafy vegetables in South and Southeast Asia. With its rich nutritional and pharmaceutical values, it has traditionally had widespread use for dietary and herbal purposes. Here, the genome of S. androgynus was sequenced and assembled, revealing a genome size of 1.55 Gb with 26 pseudo-chromosomes. Phylogenetic analysis traced back the divergence of Sauropus from Phyllanthus to approximately 29.67 million years ago (Mya). Genome analysis revealed that S. androgynus polyploidized around 20.51 Mya and shared a γ event about 132.95 Mya. Gene function analysis suggested that the expansion of pathways related to phloem development, lignin biosynthesis, and photosynthesis tended to result in the morphological differences among species within the Phyllanthaceae family, characterized by varying ploidy levels. The high accumulation of ascorbic acid in S. androgynus was attributed to the high expression of genes associated with the L-galactose pathway and recycling pathway. Moreover, the expanded gene families of S. androgynus exhibited multiple biochemical pathways associated with its comprehensive pharmacological activity, geographic adaptation and distinctive pleasurable flavor. Altogether, our findings represent a crucial genomic asset for S. androgynus, casting light on the intricate ploidy within the Phyllanthaceae family.

Conserved and derived expression patterns and positive selection on dental genes reveal complex evolutionary context of ever-growing rodent molars

Background

Continuously growing teeth are an important innovation in mammalian evolution, yet genetic regulation of continuous growth by stem cells remains incompletely understood. Dental stem cells are lost at the onset of tooth root formation, but this loss of continuous crown growth is difficult to study in the mouse because regulatory signaling overlaps with signals that pattern tooth size and shape. Within the voles (Cricetidae, Rodentia, Glires), species have evolved both rooted and unrooted molars that have similar size and shape. We assembled a de novo genome of Myodes glareolus, a vole with high-crowned, rooted molars, and performed genomic and transcriptomic analyses in a broad phylogenetic context of Glires (rodents and lagomorphs) to assess differential selection and evolution in tooth forming genes.

Results

Our de novo genome recovered 91% of single-copy orthologs for Euarchontoglires and had a total length of 2.44 Gigabases, enabling genomic and transcriptomic analyses. We identified six dental genes undergoing positive selection across Glires and two genes undergoing positive selection in species with unrooted molars, Dspp and Aqp1. Transcriptomics analyses demonstrated conserved patterns of dental gene expression with species-specific variation likely related to developmental timing and morphological differences between mouse and vole molars.

Conclusions

Our results support ongoing dental gene evolution in rodents with unrooted molars. We identify candidate genes for further functional analyses, particularly Dspp, which plays an important role in mineralizing tissues. Our expression results support conservation of dental genes between voles and model species like mice, while revealing significant effects of overall tooth morphology on gene expression.

Chromosomal-level genome and metabolome analyses of highly heterozygous allohexaploid Dendrocalamus brandisii elucidate shoot quality and developmental characteristics

Dendrocalamus brandisii (Munro) Kurz is a sympodial bamboo species with inimitable taste and flavorful shoots. Its rapid growth and use as high-quality material make this bamboo species highly valued for both food processing and wood applications. However, genome information for D. brandisii is lacking, primarily due to its polyploidy and large genome size. Here, we assembled a high-quality genome for hexaploid D. brandisii, which comprises 70 chromosomes with a total size of 2,756 Mb, using long-read HiFi sequencing. Furthermore, we accurately separated the genome into its three constituent subgenomes. We used Oxford Nanopore Technologies long reads to construct a transcriptomic dataset covering 15 tissues for gene annotation to complement our genome assembly, revealing differential gene expression and post-transcriptional regulation. By integrating metabolome analysis, we unveiled that well-balanced lignin formation, as well as abundant flavonoid and fructose contents, contribute to the superior quality of D. brandisii shoots. Integrating genomic, transcriptomic, and metabolomic datasets provided a solid foundation for enhancing bamboo shoot quality and developing efficient gene-editing techniques. This study should facilitate research on D. brandisii and enhance its use as a food source and wood material by providing crucial genomic resources.

Viral Diversity in Samples of Freshwater Gastropods Benedictia baicalensis (Caenogastropoda: Benedictiidae) Revealed by Total RNA-Sequencing

Previously, the main studies were focused on viruses that cause disease in commercial and farmed shellfish and cause damage to food enterprises (for example, Ostreavirusostreidmalaco1, Aurivirus haliotidmalaco1 and Aquabirnavirus tellinae). Advances in high-throughput sequencing technologies have extended the studies to natural populations of mollusks (and other invertebrates) as unexplored niches of viral diversity and possible sources of emerging diseases. These studies have revealed a huge diversity of mostly previously unknown viruses and filled gaps in the evolutionary history of viruses. In the present study, we estimated the viral diversity in samples of the Baikal endemic gastropod Benedictia baicalensis using metatranscriptomic analysis (total RNA-sequencing); we were able to identify a wide variety of RNA-containing viruses in four samples (pools) of mollusks collected at three stations of Lake Baikal. Most of the identified viral genomes (scaffolds) had only distant similarities to known viruses or (in most cases) to metagenome-assembled viral genomes from various natural samples (mollusks, crustaceans, insects and others) mainly from freshwater ecosystems. We were able to identify viruses similar to those previously identified in mollusks (in particular to the picornaviruses Biomphalaria virus 1 and Biomphalaria virus 3 from the freshwater gastropods); it is possible that picorna-like viruses (as well as a number of other identified viruses) are pathogenic for Baikal gastropods. Our results also suggested that Baikal mollusks, like other species, may bioaccumulate or serve as a reservoir for numerous viruses that infect a variety of organisms (including vertebrates).

Improved high-quality reference genome of red drum facilitates the processes of resistance-related gene exploration

Sciaenops ocellatus is among the most important artificially introduced farmed fish across 11 countries and regions. However, the frequent occurrence of extreme weather events and breeding escapes have placed great pressure on local marine biodiversity and ecosystems. We reported the de novo assembly and annotation with a contig N50 of 28.30 Mb using PacBio HiFi sequencing and Hi-C technologies, which resulted in a 283-fold increase in contig N50 length and improvement in continuity and quality in complex repetitive region for S. ocellatus compared to the previous version. In total, 257.36 Mb of repetitive sequences accounted for 35.48% of the genome, and 22,845 protein-coding genes associated with a BUSCO value of 98.32%, were identified by genome annotation. Moreover, 54 hub genes rapidly responding to hypoosmotic stress were identified by WGCNA. The high-quality chromosome-scale S. ocellatus genome and candidate resistance-related gene sets will not only provide a genomic basis for genetic improvement via molecular breeding, but will also lay an important foundation for investigating the molecular regulation of rapid responses to stress.

A high-quality chromosome-level wild rice genome of Oryza coarctata

Oryza coarctata (2n = 4X = 48, KKLL) is an allotetraploid, undomesticated relative of rice and the only species in the genus Oryza with tolerance to high salinity and submergence. Therefore, it contains important stress and tolerance genes/factors for rice. The initial draft genome published was limited by data and technical restrictions, leading to an incomplete and highly fragmented assembly. This study reports a new, highly contiguous chromosome-level genome assembly and annotation of O. coarctata. PacBio high-quality HiFi reads generated 460 contigs with a total length of 573.4 Mb and an N50 of 23.1 Mb, which were assembled into scaffolds with Hi-C data, anchoring 96.99% of the assembly onto 24 chromosomes. The genome assembly comprises 45,571 genes, and repetitive content contributes 25.5% of the genome. This study provides the novel identification of the KK and LL genome types of the genus Oryza, leading to valuable insights into rice genome evolution. The chromosome-level genome assembly of O. coarctata is a valuable resource for rice research and molecular breeding.

Pangolin Genomes Offer Key Insights and Resources for the World's Most Trafficked Wild Mammals

Pangolins form a group of scaly mammals that are trafficked at record numbers for their meat and purported medicinal properties. Despite their conservation concern, knowledge of their evolution is limited by a paucity of genomic data. We aim to produce exhaustive genomic resources that include 3,238 orthologous genes and whole-genome polymorphisms to assess the evolution of all eight extant pangolin species. Robust orthologous gene-based phylogenies recovered the monophyly of the three genera and highlighted the existence of an undescribed species closely related to Southeast Asian pangolins. Signatures of middle Miocene admixture between an extinct, possibly European, lineage and the ancestor of Southeast Asian pangolins, provide new insights into the early evolutionary history of the group. Demographic trajectories and genome-wide heterozygosity estimates revealed contrasts between continental versus island populations and species lineages, suggesting that conservation planning should consider intraspecific patterns. With the expected loss of genomic diversity from recent, extensive trafficking not yet realized in pangolins, we recommend that populations be genetically surveyed to anticipate any deleterious impact of the illegal trade. Finally, we produce a complete set of genomic resources that will be integral for future conservation management and forensic endeavors for pangolins, including tracing their illegal trade. These comprise the completion of whole-genomes for pangolins through the hybrid assembly of the first reference genome for the giant pangolin (Smutsia gigantea) and new draft genomes (∼43x-77x) for four additional species, as well as a database of orthologous genes with over 3.4 million polymorphic sites.

The pig pangenome provides insights into the roles of coding structural variations in genetic diversity and adaptation

Structural variations have emerged as an important driving force for genome evolution and phenotypic variation in various organisms, yet their contributions to genetic diversity and adaptation in domesticated animals remain largely unknown. Here we constructed a pangenome based on 250 sequenced individuals from 32 pig breeds in Eurasia and systematically characterized coding sequence presence/absence variations (PAVs) within pigs. We identified 308.3-Mb nonreference sequences and 3438 novel genes absent from the current reference genome. Gene PAV analysis showed that 16.8% of the genes in the pangene catalog undergo PAV. A number of newly identified dispensable genes showed close associations with adaptation. For instance, several novel swine leukocyte antigen (SLA) genes discovered in nonreference sequences potentially participate in immune responses to productive and respiratory syndrome virus (PRRSV) infection. We delineated previously unidentified features of the pig mobilome that contained 490,480 transposable element insertion polymorphisms (TIPs) resulting from recent mobilization of 970 TE families, and investigated their population dynamics along with influences on population differentiation and gene expression. In addition, several candidate adaptive TE insertions were detected to be co-opted into genes responsible for responses to hypoxia, skeletal development, regulation of heart contraction, and neuronal cell development, likely contributing to local adaptation of Tibetan wild boars. These findings enhance our understanding on hidden layers of the genetic diversity in pigs and provide novel insights into the role of SVs in the evolutionary adaptation of mammals.

Analyses of 600+ insect genomes reveal repetitive element dynamics and highlight biodiversity-scale repeat annotation challenges

Repetitive elements (REs) are integral to the composition, structure, and function of eukaryotic genomes, yet remain understudied in most taxonomic groups. We investigated REs across 601 insect species and report wide variation in RE dynamics across groups. Analysis of associations between REs and protein-coding genes revealed dynamic evolution at the interface between REs and coding regions across insects, including notably elevated RE-gene associations in lineages with abundant long interspersed nuclear elements (LINEs). We leveraged this large, empirical data set to quantify impacts of long-read technology on RE detection and investigate fundamental challenges to RE annotation in diverse groups. In long-read assemblies, we detected ∼36% more REs than short-read assemblies, with long terminal repeats (LTRs) showing 162% increased detection, whereas DNA transposons and LINEs showed less respective technology-related bias. In most insect lineages, 25%-85% of repetitive sequences were "unclassified" following automated annotation, compared with only ∼13% in Drosophila species. Although the diversity of available insect genomes has rapidly expanded, we show the rate of community contributions to RE databases has not kept pace, preventing efficient annotation and high-resolution study of REs in most groups. We highlight the tremendous opportunity and need for the biodiversity genomics field to embrace REs and suggest collective steps for making progress toward this goal.

Multi-omic insights into the formation and evolution of a novel shell microstructure in oysters

BACKGROUND

Molluscan shell, composed of a diverse range of architectures and microstructures, is a classic model system to study the relationships between molecular evolution and biomineralized structure formation. The shells of oysters differ from those of other molluscs by possessing a novel microstructure, chalky calcite, which facilitates adaptation to the sessile lifestyle. However, the genetic basis and evolutionary origin of this adaptive innovation remain largely unexplored.

RESULTS

We report the first whole-genome assembly and shell proteomes of the Iwagaki oyster Crassostrea nippona. Multi-omic integrative analyses revealed that independently expanded and co-opted tyrosinase, peroxidase, TIMP genes may contribute to the chalky layer formation in oysters. Comparisons with other molluscan shell proteomes imply that von Willebrand factor type A and chitin-binding domains are basic members of molluscan biomineralization toolkit. Genome-wide identification and analyses of these two domains in 19 metazoans enabled us to propose that the well-known Pif may share a common origin in the last common ancestor of Bilateria. Furthermore, Pif and LamG3 genes acquire new genetic function for shell mineralization in bivalves and the chalky calcite formation in oysters likely through a combination of gene duplication and domain reorganization.

CONCLUSIONS

The spatial expression of SMP genes in the mantle and molecular evolution of Pif are potentially involved in regulation of the chalky calcite deposition, thereby shaping the high plasticity of the oyster shell to adapt to a sessile lifestyle. This study further highlights neo-functionalization as a crucial mechanism for the diversification of shell mineralization and microstructures in molluscs, which may be applied more widely for studies on the evolution of metazoan biomineralization.

Chromosome level genome assembly of oriental armyworm Mythimna separata

The oriental armyworm, Mythimna separata, is an extremely destructive polyphagous pest with a broad host range that seriously threatens the safety of agricultural production. Here, a high-quality chromosome-level genome was assembled using Illumina, PacBio HiFi long sequencing, and Hi-C scaffolding technologies. The genome size was 706.30 Mb with a contig N50 of 22.08 Mb, and 99.2% of the assembled sequences were anchored to 31 chromosomes. In addition, 20,375 protein-coding genes and 258.68 Mb transposable elements were identified. The chromosome-level genome assembly of M. separata provides a significant genetic resource for future studies of this insect and contributes to the development of management strategies.

Chromosome-level genome assembly of Chouioia cunea Yang, the parasitic wasp of the fall webworm

Chouioia cunea Yang 1989 is a parasitic wasp of many lepidopteran insects during their pupal stage, and has been successfully used to control pests such as the fall webworm Hyphantria cunea. Here we reported the chromosome-level genome of C. cunea by using short (MGI-SEQ), long (Oxford Nanopore), chromatin-linked (Hi-C) sequencing reads and transcriptomic data, representing the first chromosome-level genome of parasitic wasps of the family Eulophidae. The total assembly length is 171.99 Mb, containing 6 pesudo-chromosomes with a GC content of 36.89% and the scaffold/contig N50 length of 31.70/26.52 Mb. The BUSCO completeness of the assembly was estimated to be 98.7%. A total of 12,258 protein-coding genes (PCGs), 10,547 3'-UTRs, and 10,671 5'-UTRs were annotated. This high-quality genome is an important step toward a better understanding of the genomes of the Eulophidae (Chalcidoidea), and will serve as a valuable resource for analyses of phylogenetic relationships and the evolution of Hymenoptera.

RNA covariation at helix-level resolution for the identification of evolutionarily conserved RNA structure

Many biologically important RNAs fold into specific 3D structures conserved through evolution. Knowing when an RNA sequence includes a conserved RNA structure that could lead to new biology is not trivial and depends on clues left behind by conservation in the form of covariation and variation. For that purpose, the R-scape statistical test was created to identify from alignments of RNA sequences, the base pairs that significantly covary above phylogenetic expectation. R-scape treats base pairs as independent units. However, RNA base pairs do not occur in isolation. The Watson-Crick (WC) base pairs stack together forming helices that constitute the scaffold that facilitates the formation of the non-WC base pairs, and ultimately the complete 3D structure. The helix-forming WC base pairs carry most of the covariation signal in an RNA structure. Here, I introduce a new measure of statistically significant covariation at helix-level by aggregation of the covariation significance and covariation power calculated at base-pair-level resolution. Performance benchmarks show that helix-level aggregated covariation increases sensitivity in the detection of evolutionarily conserved RNA structure without sacrificing specificity. This additional helix-level sensitivity reveals an artifact that results from using covariation to build an alignment for a hypothetical structure and then testing the alignment for whether its covariation significantly supports the structure. Helix-level reanalysis of the evolutionary evidence for a selection of long non-coding RNAs (lncRNAs) reinforces the evidence against these lncRNAs having a conserved secondary structure.

Chromosome genome assembly and annotation of the spiny red gurnard (Chelidonichthys spinosus)

Chelidonichthys spinosus, a secondary economic fish, is increasingly being exploited and valued in China. However, overfishing has led to it being recognized as one of the most depleted marine species in China. In this study, we generated a chromosome-level genome of C. spinosus using PacBio, Illumina, and Hi-C sequencing data. Ultimately, we assembled a 624.7 Mb genome of C. spinosus, with a contig N50 of 13.77 Mb and scaffold N50 of 28.11 Mb. We further anchored and oriented the assembled sequences onto 24 pseudo-chromosomes using Hi-C techniques. In total, 25,358 protein-coding genes were predicted, of which 24,072 (94.93%) genes were functionally annotated. The dot plot reveals a prominent co-linearity between C. spinosus and Cyclopterus lumpus, indicating a remarkably close phylogenetic relationship between these two species. The assembled genome sequences provide valuable information for elucidating the genetic adaptation and potential molecular basis of C. spinosus. They also have the potential to provide insight into the evolutionary investigation of teleost fish and vertebrates.

Towards the Investigation of the Adaptive Divergence in a Species of Exceptional Ecological Plasticity: Chromosome-Scale Genome Assembly of Chouardia litardierei (Hyacinthaceae)

One of the central goals of evolutionary biology is to understand the genomic basis of adaptive divergence. Different aspects of evolutionary processes should be studied through genome-wide approaches, therefore maximizing the investigated genomic space. However, in-depth genome-scale analyses often are restricted to a model or economically important species and their closely related wild congeners with available reference genomes. Here, we present the high-quality chromosome-level genome assembly of Chouardia litardierei, a plant species with exceptional ecological plasticity. By combining PacBio and Hi-C sequencing technologies, we generated a 3.7 Gbp genome with a scaffold N50 size of 210 Mbp. Over 80% of the genome comprised repetitive elements, among which the LTR retrotransposons prevailed. Approximately 86% of the 27,257 predicted genes were functionally annotated using public databases. For the comparative analysis of different ecotypes' genomes, the whole-genome sequencing of two individuals, each from a distinct ecotype, was performed. The detected above-average SNP density within coding regions suggests increased adaptive divergence-related mutation rates, therefore confirming the assumed divergence processes within the group. The constructed genome presents an invaluable resource for future research activities oriented toward the investigation of the genetics underlying the adaptive divergence that is likely unfolding among the studied species' ecotypes.

Multiple and diversified transposon lineages contribute to early and recent bivalve genome evolution

BACKGROUND

Transposable elements (TEs) can represent one of the major sources of genomic variation across eukaryotes, providing novel raw materials for species diversification and innovation. While considerable effort has been made to study their evolutionary dynamics across multiple animal clades, molluscs represent a substantially understudied phylum. Here, we take advantage of the recent increase in mollusc genomic resources and adopt an automated TE annotation pipeline combined with a phylogenetic tree-based classification, as well as extensive manual curation efforts, to characterize TE repertories across 27 bivalve genomes with a particular emphasis on DDE/D class II elements, long interspersed nuclear elements (LINEs), and their evolutionary dynamics.

RESULTS

We found class I elements as highly dominant in bivalve genomes, with LINE elements, despite less represented in terms of copy number per genome, being the most common retroposon group covering up to 10% of their genome. We mined 86,488 reverse transcriptases (RVT) containing LINE coming from 12 clades distributed across all known superfamilies and 14,275 class II DDE/D-containing transposons coming from 16 distinct superfamilies. We uncovered a previously underestimated rich and diverse bivalve ancestral transposon complement that could be traced back to their most recent common ancestor that lived ~ 500 Mya. Moreover, we identified multiple instances of lineage-specific emergence and loss of different LINEs and DDE/D lineages with the interesting cases of CR1- Zenon, Proto2, RTE-X, and Academ elements that underwent a bivalve-specific amplification likely associated with their diversification. Finally, we found that this LINE diversity is maintained in extant species by an equally diverse set of long-living and potentially active elements, as suggested by their evolutionary history and transcription profiles in both male and female gonads.

CONCLUSIONS

We found that bivalves host an exceptional diversity of transposons compared to other molluscs. Their LINE complement could mainly follow a "stealth drivers" model of evolution where multiple and diversified families are able to survive and co-exist for a long period of time in the host genome, potentially shaping both recent and early phases of bivalve genome evolution and diversification. Overall, we provide not only the first comparative study of TE evolutionary dynamics in a large but understudied phylum such as Mollusca, but also a reference library for ORF-containing class II DDE/D and LINE elements, which represents an important genomic resource for their identification and characterization in novel genomes.

Chromosome-level assemblies of cultivated water chestnut Trapa bicornis and its wild relative Trapa incisa

Water chestnut (Trapa L.) is a floating-leaved aquatic plant with high edible and medicinal value. In this study, we presented chromosome-level genome assemblies of cultivated large-seed species Trapa bicornis and its wild small-seed relative Trapa incisa by using PacBio HiFi long reads and Hi-C technology. The T. bicornis and T. incisa assemblies consisted of 479.90 Mb and 463.97 Mb contigs with N50 values of 13.52 Mb and 13.77 Mb, respectively, and repeat contents of 62.88% and 62.49%, respectively. A total of 33,306 and 33,315 protein-coding genes were predicted in T. bicornis and T. incisa assemblies, respectively. There were 159,232 structural variants affecting more than 11 thousand genes detected between the two genomes. The phylogenetic analysis indicated that the lineage leading to Trapa was diverged from the lineage to Sonneratia approximately 23 million years ago. These two assemblies provide valuable resources for future evolutionary and functional genomic research and molecular breeding of water chestnut.

A chromosome-scale genome assembly of Artemisia argyi reveals unbiased subgenome evolution and key contributions of gene duplication to volatile terpenoid diversity

Artemisia argyi Lévl. et Vant., a perennial Artemisia herb with an intense fragrance, is widely used in traditional medicine in China and many other Asian countries. Here, we present a chromosome-scale genome assembly of A. argyi comprising 3.89 Gb assembled into 17 pseudochromosomes. Phylogenetic and comparative genomic analyses revealed that A. argyi underwent a recent lineage-specific whole-genome duplication (WGD) event after divergence from Artemisia annua, resulting in two subgenomes. We deciphered the diploid ancestral genome of A. argyi, and unbiased subgenome evolution was observed. The recent WGD led to a large number of duplicated genes in the A. argyi genome. Expansion of the terpene synthase (TPS) gene family through various types of gene duplication may have greatly contributed to the diversity of volatile terpenoids in A. argyi. In particular, we identified a typical germacrene D synthase gene cluster within the expanded TPS gene family. The entire biosynthetic pathways of germacrenes, (+)-borneol, and (+)-camphor were elucidated in A. argyi. In addition, partial deletion of the amorpha-4,11-diene synthase (ADS) gene and loss of function of ADS homologs may have resulted in the lack of artemisinin production in A. argyi. Our study provides new insights into the genome evolution of Artemisia and lays a foundation for further improvement of the quality of this important medicinal plant.

RNA covariation at helix-level resolution for the identification of evolutionarily conserved RNA structure

Many biologically important RNAs fold into specific 3D structures conserved through evolution. Knowing when an RNA sequence includes a conserved RNA structure that could lead to new biology is not trivial and depends on clues left behind by conservation in the form of covariation and variation. For that purpose, the R-scape statistical test was created to identify from alignments of RNA sequences, the base pairs that significantly covary above phylogenetic expectation. R-scape treats base pairs as independent units. However, RNA base pairs do not occur in isolation. The Watson-Crick (WC) base pairs stack together forming helices that constitute the scaffold that facilitates the formation of the non-WC base pairs, and ultimately the complete 3D structure. The helix-forming WC base pairs carry most of the covariation signal in an RNA structure. Here, I introduce a new measure of statistically significant covariation at helix-level by aggregation of the covariation significance and covariation power calculated at base-pair-level resolution. Performance benchmarks show that helix-level aggregated covariation increases sensitivity in the detection of evolutionarily conserved RNA structure without sacrificing specificity. This additional helix-level sensitivity reveals an artifact that results from using covariation to build an alignment for a hypothetical structure and then testing the alignment for whether its covariation significantly supports the structure. Helix-level reanalysis of the evolutionary evidence for a selection of long non-coding RNAs (lncRNAs) reinforces the evidence against these lncRNAs having a conserved secondary structure.

Availability

Helix aggregated E-values are integrated in the R-scape software package (version 2.0.0.p and higher). The R-scape web server eddylab.org/R-scape includes a link to download the source code.

Contact

elenarivas@fas.harvard.edu.

Supplementary information

Supplementary data and code are provided with this manuscript at rivaslab.org .

Whole-genome identification of transposable elements reveals the equine repetitive element insertion polymorphism in Chinese horses

Transposable elements (TEs) are diverse, abundant, and complicated in genomes. They not only can drive the genome evolution process but can also act as special resources for adaptation. However, little is known about the evolutionary processes that shaped horses. In this work, 126 horse assemblages involved in most horse breeds in China were used to investigate the patterns of TE variation for the first time. By using RepeatMasker and melt software, we found that the horse-specific short interspersed repetitive elements family, equine repetitive elements (ERE1), exhibited polymorphisms in horse genomes. Phylogenetic analysis based on these ERE1 loci (minor allele frequency ≥0.05) revealed three major horse groups, namely, those in northern China, southern China, and Qinghai-Tibetan, which mirrors the result determined by SNPs to some extent. The present ERE1 family emerged ~0.26 to 1.77 Mya ago, with an activity peak at ~0.49 Mya, which matches the early stage of the horse lineage and decreases after the divergence of Equus caballus and Equus ferus przewalskii. To detect the functional ERE1(s) associated with adaptation, locus-specific branch length, genome-wide association study, and absolute allele frequency difference analyses were conducted and resulted in two common protein-coding genes annotated by candidate ERE1s. They were clustered into the vascular smooth muscle contraction (p = 0.01, EDNRA) and apelin signalling pathways (p = 0.02, NRF1). Notably, ERE1 insertion into the EDNRA gene showed a higher association with adaptation among southern China horses and other horses in 15 populations and 451 individuals (p = 4.55 e-8). Our results provide a comprehensive understanding of TE variations to analyse the phylogenetic relationships and traits relevant to adaptive evolution in horses.

An ASO therapy for Angelman syndrome that targets an evolutionarily conserved region at the start of the UBE3A-AS transcript

Angelman syndrome is a devastating neurogenetic disorder for which there is currently no effective treatment. It is caused by mutations or epimutations affecting the expression or function of the maternally inherited allele of the ubiquitin-protein ligase E3A (UBE3A) gene. The paternal UBE3A allele is imprinted in neurons of the central nervous system (CNS) by the UBE3A antisense (UBE3A-AS) transcript, which represents the distal end of the small nucleolar host gene 14 (SNHG14) transcription unit. Reactivating the expression of the paternal UBE3A allele in the CNS has long been pursued as a therapeutic option for Angelman syndrome. Here, we described the development of an antisense oligonucleotide (ASO) therapy for Angelman syndrome that targets an evolutionarily conserved region demarcating the start of the UBE3A-AS transcript. We designed and chemically optimized gapmer ASOs targeting specific sequences at the start of the human UBE3A-AS transcript. We showed that ASOs targeting this region precisely and efficiently repress the transcription of UBE3A-AS, reactivating the expression of the paternal UBE3A allele in neurotypical and Angelman syndrome induced pluripotent stem cell-derived neurons. We further showed that human-targeted ASOs administered to the CNS of cynomolgus macaques by lumbar intrathecal injection repress UBE3A-AS and reactivate the expression of the paternal UBE3A allele throughout the CNS. These findings support the advancement of this investigational molecular therapy for Angelman syndrome into clinical development (ClinicalTrials.gov, NCT04259281).

A reference-grade genome assembly for Astragalus mongholicus and insights into the biosynthesis and high accumulation of triterpenoids and flavonoids in its roots

Astragalus membranaceus var. mongholicus (AMM), a member of the Leguminosae, is one of the most important medicinal plants worldwide. The dried roots of AMM have a wide range of pharmacological effects and are a traditional Chinese medicine. Here, we report the first chromosome-level reference genome of AMM, comprising nine pseudochromosomes with a total size of 1.47 Gb and 27 868 protein-encoding genes. Comparative genomic analysis reveals that AMM has not experienced an independent whole-genome duplication (WGD) event after the WGD event shared by the Papilionoideae species. Analysis of long terminal repeat retrotransposons suggests a recent burst of these elements at approximately 0.13 million years ago, which may explain the large size of the AMM genome. Multiple gene families involved in the biosynthesis of triterpenoids and flavonoids were expanded, and our data indicate that tandem duplication has been the main driver for expansion of these families. Among the expanded families, the phenylalanine ammonia-lyase gene family was primarily expressed in the roots of AMM, suggesting their roles in the biosynthesis of phenylpropanoid compounds. The functional versatility of 2,3-oxidosqualene cyclase genes in cluster III may play a critical role in the diversification of triterpenoids in AMM. Our findings provide novel insights into triterpenoid and flavonoid biosynthesis and can facilitate future research on the genetics and medical applications of AMM.

Chromosome-level genome assembly of the European flat oyster (Ostrea edulis) provides insights into its evolution and adaptation

The European flat oyster (Ostrea edulis) is an endangered and economically important marine bivalve species that plays a critical role in the coastal ecosystem. Here, we report a high-quality chromosome-level genome assembly of O. edulis, generated using PacBio HiFi-CCS long reads and annotated with Nanopore full-length transcriptome. The O. edulis genome covers 946.06 Mb (scaffold N50 94.82 Mb) containing 34,495 protein-coding genes and a high proportion of repeat sequences (58.49 %). The reconstructed demographic histories show that O. edulis population might be shaped by breeding habit (embryo brooding) and historical climatic change. Comparative genomic analysis indicates that transposable elements may drive lineage-specific evolution in oysters. Notably, the O. edulis genome has a Hox gene cluster rearrangement that has never been reported in bivalves, making this species valuable for evolutionary studies of molluscan diversification. Moreover, genome expansion of O. edulis is probably central to its adaptation to filter-feeding and sessile lifestyles, as well as embryo brooding and pathogen resistance, in coastal ecosystems. This chromosome-level genome assembly provides new insights into the genome feature of oysters, and presents an important resource for genetic research, evolutionary studies, and biological conservation of O. edulis.

Genome assembly of wild loquat (Eriobotrya japonica) and resequencing provide new insights into the genomic evolution and fruit domestication in loquat

Wild loquats (Eriobotrya japonica Lindl.) provide remarkable genetic resources for studying domestication and breeding improved varieties. Herein, we generate the first high-quality chromosome-level genome assembly of wild loquat, with 45 791 predicted protein-coding genes. Analysis of comparative genomics indicated that loquat shares a common ancestor with apple and pear, and a recent whole-genome duplication event occurred in loquat prior to its divergence. Genome resequencing showed that the loquat germplasms can be distinctly classified into wild and cultivated groups, and the commercial cultivars have experienced allelic admixture. Compared with cultivated loquats, the wild loquat genome showed very few selected genomic regions and had higher levels of genetic diversity. However, whole-genome scans of selective sweeps were mainly related to fruit quality, size, and flesh color during the domestication process. Large-scale transcriptome and metabolome analyses were further performed to identify differentially expressed genes (DEGs) and differentially accumulated metabolites (DAMs) in wild and cultivated loquats at various fruit development stages. Unlike those in wild loquat, the key DEGs and DAMs involved in carbohydrate metabolism, plant hormone signal transduction, flavonoid biosynthesis, and carotenoid biosynthesis were significantly regulated in cultivated loquats during fruit development. These high-quality reference genome, resequencing, and large-scale transcriptome/metabolome data provide valuable resources for elucidating fruit domestication and molecular breeding in loquat.

TeCD: The eccDNA Collection Database for extrachromosomal circular DNA

BACKGROUND

Extrachromosomal circular DNA (eccDNA) is a kind of DNA that widely exists in eukaryotic cells. Studies in recent years have shown that eccDNA is often enriched during tumors and aging, and participates in the development of cell physiological activities in a special way, so people have paid more and more attention to the eccDNA, and it has also become a critical new topic in modern biological research.

DESCRIPTION

We built a database to collect eccDNA, including animals, plants and fungi, and provide researchers with an eccDNA retrieval platform. The collected eccDNAs were processed in a uniform format and classified according to the species to which it belongs and the chromosome of the source. Each eccDNA record contained sequence length, start and end sites on the corresponding chromosome, order of the bases, genomic elements such as genes and transposons, and other information in the respective sequencing experiment. All the data were stored into the TeCD (The eccDNA Collection Database) and the BLAST (Basic Local Alignment Search Tool) sequence alignment function was also added into the database for analyzing the potential eccDNA sequences.

CONCLUSION

We built TeCD, a platform for users to search and obtain eccDNA data, and analyzed the possible potential functions of eccDNA. These findings may provide a basis and direction for researchers to further explore the biological significance of eccDNA in the future.

Chromosome-scale genomics, metabolomics, and transcriptomics provide insight into the synthesis and regulation of phenols in Vitis adenoclada grapes

Vitis adenoclada is a wild grape unique to China. It exhibits well resistance to heat, humidity, fungal disease, drought, and soil infertility. Here, we report the high-quality, chromosome-level genome assembly of GH6 (V. adenoclada). The 498.27 Mb genome contained 221.78 Mb of transposable elements, 28,660 protein-coding genes, and 481.44 Mb of sequences associated with 19 chromosomes. GH6 shares a common ancestor with PN40024 (Vitis vinifera) from approximately 4.26-9.01 million years ago, whose divergence occurred later than Vitis rotundifolia and Vitis riparia. Widely-targeted metabolome and transcriptome analysis revealed that the profiles and metabolism of phenolic compounds in V. adenoclada varieties significantly were differed from other grape varieties. Specifically, V. adenoclada varieties were rich in phenolic acids and flavonols, whereas the flavan-3-ol and anthocyanin content was lower compared with other varieties that have V. vinifera consanguinity in this study. In addition, ferulic acid and stilbenes content were associated with higher expressions of COMT and STSs in V. adenoclada varieties. Furthermore, MYB2, MYB73-1, and MYB73-2 were presumably responsible for the high expression level of COMT in V. adenoclada berries. MYB12 (MYBF1) was positively correlated with PAL, CHS, FLS and UFGT.Meanwhile, MYB4 and MYBC2-L1 may inhibit the synthesis of flavan-3-ols and anthocyanins in two V. adenoclada varieties (YN2 and GH6). The publication of the V. adenoclada grape genome provides a molecular foundation for further revealing its flavor and quality characteristics, is also important for identifying favorable genes of the East Asian species for future breeding.

Chromosome-level genome assembly of a high-altitude-adapted frog (Rana kukunoris) from the Tibetan plateau provides insight into amphibian genome evolution and adaptation

BACKGROUND

The high-altitude-adapted frog Rana kukunoris, occurring on the Tibetan plateau, is an excellent model to study life history evolution and adaptation to harsh high-altitude environments. However, genomic resources for this species are still underdeveloped constraining attempts to investigate the underpinnings of adaptation.

RESULTS

The R. kukunoris genome was assembled to a size of 4.83 Gb and the contig N50 was 1.80 Mb. The 6555 contigs were clustered and ordered into 12 pseudo-chromosomes covering ~ 93.07% of the assembled genome. In total, 32,304 genes were functionally annotated. Synteny analysis between the genomes of R. kukunoris and a low latitude species Rana temporaria showed a high degree of chromosome level synteny with one fusion event between chr11 and chr13 forming pseudo-chromosome 11 in R. kukunoris. Characterization of features of the R. kukunoris genome identified that 61.5% consisted of transposable elements and expansions of gene families related to cell nucleus structure and taste sense were identified. Ninety-five single-copy orthologous genes were identified as being under positive selection and had functions associated with the positive regulation of proteins in the catabolic process and negative regulation of developmental growth. These gene family expansions and positively selected genes indicate regions for further interrogation to understand adaptation to high altitude.

CONCLUSIONS

Here, we reported a high-quality chromosome-level genome assembly of a high-altitude amphibian species using a combination of Illumina, PacBio and Hi-C sequencing technologies. This genome assembly provides a valuable resource for subsequent research on R. kukunoris genomics and amphibian genome evolution in general.

Endogenous Retrovirus Elements Are Co-Expressed with IFN Stimulation Genes in the JAK-STAT Pathway

Background: Endogenous retrovirus (ERV) elements can act as proximal regulatory elements in promoting interferon (IFN) responses. Previous relevant studies have mainly focused on IFN-stimulated genes (ISGs). However, the role of ERV elements as cis-regulatory motifs in regulating genes of the JAK-STAT pathway remains poorly understood. In our study, we analyzed the changes in ERV elements and genes under both IFN stimulation and blockade of the signaling pathway. Methods: The effects of interferon on cells under normal conditions and knockout of the receptor were compared based on the THP1_IFNAR1_KO and THP1_IFNAR2_mutant cell lines. The correlation between differentially expressed ERVs (DHERVs) and differentially expressed genes (DEGs) as DEHERV-G pairs was explored with construction of gene regulatory networks related to ERV and induced by proinflammatory cytokines. Results: A total of 430 DEHERV loci and 190 DEGs were identified in 842 DEHERV-G pairs that are common to the three groups. More than 87% of DEHERV-G pairs demonstrated a consistent expression pattern. ISGs such as AIM2, IFIT1, IFIT2, IFIT3, STAT1, and IRF were activated via the JAK-STAT pathway in response to interferon stimulation. Thus, STAT1, STAT2, and IRF1 appear to play core roles in regulatory networks and are closely associated with ERVs. Conclusions: The RNA expression of ISGs and ERV elements is correlated, indicating that ERV elements are closely linked to host innate immune responses.

Mobilome of the Rhus Gall Aphid Schlechtendalia chinensis Provides Insight into TE Insertion-Related Inactivation of Functional Genes

Transposable elements (TEs) comprise a considerable proportion of insect genomic DNA; how they contribute to genome structure and organization is still poorly understood. Here, we present an analysis of the TE repertoire in the chromosome-level genome assembly of Rhus gall aphid Schlechtendalia chinensis. The TE fractions are composed of at least 32 different superfamilies and many TEs from different families were transcriptionally active in the S. chinensis genome. Furthermore, different types of transposase-derived proteins were also found in the S. chinensis genome. We also provide insight into the TEs related insertional inactivation, and exogenization of TEs in functional genes. We considered that the presence of TE fragments in the introns of functional genes could impact the activity of functional genes, and a large number of TE fragments in introns could lead to the indirect inactivation of functional genes. The present study will be beneficial in understanding the role and impact of TEs in genomic evolution of their hosts.

Transposable Elements Continuously Remodel the Regulatory Landscape, Transcriptome, and Function of Decidual Stromal Cells

Gene expression evolution underlies the origin, divergence, and conservation of biological characters including cell-types, tissues, and organ systems. Previously we showed that large-scale gene expression changes in decidual stromal cells (DSCs) contributed to the origins of pregnancy in eutherians and the divergence of pregnancy traits in primates and that transposable elements likely contributed to these gene expression changes. Here we show that two large waves of TEs remodeled the transcriptome and regulatory landscape of DSCs, including a major wave in primates. Genes nearby TE-derived regulatory elements are among the most progesterone responsive in the genome and play essential roles in orchestrating progesterone responsiveness and the core function of decidual cells by donating progesterone receptor binding sites to the genome. We tested the regulatory abilities of 89 TE consensus sequences and found that nearly all of them acted as repressors in mammalian cells, but treatment with a histone deacetylase inhibitor unmasked latent enhancer functions. These data indicate that TEs have played an important role in the development, evolution, and function of primate DSCs and suggest a two-step model in which latent enhancer functions of TEs are unmasked after they lose primary repressor functions.

Transposable element-mediated rearrangements are prevalent in human genomes

Transposable elements constitute about half of human genomes, and their role in generating human variation through retrotransposition is broadly studied and appreciated. Structural variants mediated by transposons, which we call transposable element-mediated rearrangements (TEMRs), are less well studied, and the mechanisms leading to their formation as well as their broader impact on human diversity are poorly understood. Here, we identify 493 unique TEMRs across the genomes of three individuals. While homology directed repair is the dominant driver of TEMRs, our sequence-resolved TEMR resource allows us to identify complex inversion breakpoints, triplications or other high copy number polymorphisms, and additional complexities. TEMRs are enriched in genic loci and can create potentially important risk alleles such as a deletion in TRIM65, a known cancer biomarker and therapeutic target. These findings expand our understanding of this important class of structural variation, the mechanisms responsible for their formation, and establish them as an important driver of human diversity.

Chromosome-level assembly and analysis of the Thymus genome provide insights into glandular secretory trichome formation and monoterpenoid biosynthesis in thyme

Thyme has medicinal and aromatic value because of its potent antimicrobial and antioxidant properties. However, the absence of a fully sequenced thyme genome limits functional genomic studies of Chinese native thymes. Thymus quinquecostatus Čelak., which contains large amounts of bioactive monoterpenes such as thymol and carvacrol, is an important wild medicinal and aromatic plant in China. Monoterpenoids are abundant in glandular secretory trichomes. Here, high-fidelity and chromatin conformation capture technologies were used to assemble and annotate the T. quinquecostatus genome at the chromosome level. The 13 chromosomes of T. quinquecostatus had a total length of 528.66 Mb, a contig N50 of 8.06 Mb, and a BUSCO score of 97.34%. We found that T. quinquecostatus had experienced two whole-genome duplications, with the most recent event occurring ∼4.34 million years ago. Deep analyses of the genome, in conjunction with comparative genomic, phylogenetic, transcriptomic, and metabonomic studies, uncovered many regulatory factors and genes related to monoterpenoids and glandular secretory trichome development. Genes encoding terpene synthase (TPS), cytochrome P450 monooxygenases (CYPs), short-chain dehydrogenase/reductase (SDR), R2R3-MYB, and homeodomain-leucine zipper (HD-ZIP) IV were among those present in the T. quinquecostatus genome. Notably, Tq02G002290.1 (TqTPS1) was shown to encode the terpene synthase responsible for catalyzing production of the main monoterpene product γ-terpinene from geranyl diphosphate (GPP). Our study provides significant insight into the mechanisms of glandular secretory trichome formation and monoterpenoid biosynthesis in thyme. This work will facilitate the development of molecular breeding tools to enhance the production of bioactive secondary metabolites in Lamiaceae.