Chromosomal-level genome assembly and annotation of the tropical sea cucumber Holothuria scabra.

Holothuria scabra, a commercially valuable yet ecologically vulnerable tropical holothuroid, has experienced a severe decline in its wild populations, especially in China. Genomic resources are crucial for the development of effective genomic breeding projects and stock conservation strategies to restore these natural populations. Until now, a high-quality, chromosome-level reference genome for H. scabra has not been available. Here, we employed Oxford Nanopore and Hi-C sequencing technologies to assemble and annotate a high-quality, chromosome-level reference genome of H. scabra. The final genome comprised 31 scaffolds with a total length of 1.19 Gb and a scaffold N50 length of 53.52 Mb. Remarkably, 1,191.67 Mb (99.95%) of the sequences were anchored to 23 pseudo-chromosomes, with the longest one spanning 79.75 Mb. A total of 34,418 protein-coding genes were annotated in the final genome, with BUSCO analysis revealing 98.01% coverage of metazoa_odb10 genes, marking a significant improvement compared to the previous report. These chromosome-level sequences and annotations will provide an essential genomic basis for further investigation into molecular breeding and conservation management of H. scabra.

A chromosome-level genome assembly of an avivorous bat species (Nyctalus aviator).

Currently, three carnivorous bat species, namely Ia io, Nyctalus lasiopterus, and Nyctalus aviator, are known to actively prey on seasonal migratory birds (hereinafter referred to as "avivorous bats"). However, the absence of reference genomes impedes a thorough comprehension of the molecular adaptations of avivorous bat species. Herein, we present the high-quality chromosome-scale reference genome of N. aviator based on PacBio subreads, DNBSEQ short-reads and Hi-C sequencing data. The genome assembly size of N. aviator is 1.77 Gb, with a scaffold N50 of 102 Mb, of which 99.8% assembly was anchored into 21 pseudo-chromosomes. After masking 635.1 Mb repetitive sequences, a total of 19,412 protein-coding genes were identified, of which 99.3% were functionally annotated. The genome assembly and gene prediction reached 96.1% and 96.1% completeness of Benchmarking Universal Single-Copy Orthologs (BUSCO), respectively. This chromosome-level reference genome of N. aviator fills a gap in the existing information on the genomes of carnivorous bats, especially avivorous ones, and will be valuable for mechanism of adaptations to dietary niche expansion in bat species.

Modeling homologous chromosome recognition via nonspecific interactions.

In many organisms, most notably Drosophila, homologous chromosomes associate in somatic cells, a phenomenon known as somatic pairing, which takes place without double strand breaks or strand invasion, thus requiring some other mechanism for homologs to recognize each other. Several studies have suggested a "specific button" model, in which a series of distinct regions in the genome, known as buttons, can associate with each other, mediated by different proteins that bind to these different regions. Here, we use computational modeling to evaluate an alternative "button barcode" model, in which there is only one type of recognition site or adhesion button, present in many copies in the genome, each of which can associate with any of the others with equal affinity. In this model, buttons are nonuniformly distributed, such that alignment of a chromosome with its correct homolog, compared with a nonhomolog, is energetically favored; since to achieve nonhomologous alignment, chromosomes would be required to mechanically deform in order to bring their buttons into mutual register. By simulating randomly generated nonuniform button distributions, many highly effective button barcodes can be easily found, some of which achieve virtually perfect pairing fidelity. This model is consistent with existing literature on the effect of translocations of different sizes on homolog pairing. We conclude that a button barcode model can attain highly specific homolog recognition, comparable to that seen in actual cells undergoing somatic homolog pairing, without the need for specific interactions. This model may have implications for how meiotic pairing is achieved.

Annotations of four high-quality indigenous chicken genomes identify more than one thousand missing genes in subtelomeric regions and micro-chromosomes with high G/C contents.

BACKGROUND: Although multiple chicken genomes have been assembled and annotated, the numbers of protein-coding genes in chicken genomes and their variation among breeds are still uncertain due to the low quality of these genome assemblies and limited resources used in their gene annotations. To fill these gaps, we recently assembled genomes of four indigenous chicken breeds with distinct traits at chromosome-level. In this study, we annotated genes in each of these assembled genomes using a combination of RNA-seq- and homology-based approaches. RESULTS: We identified varying numbers (17,497-17,718) of protein-coding genes in the four indigenous chicken genomes, while recovering 51 of the 274 "missing" genes in birds in general, and 36 of the 174 "missing" genes in chickens in particular. Intriguingly, based on deeply sequenced RNA-seq data collected in multiple tissues in the four breeds, we found 571 ~ 627 protein-coding genes in each genome, which were missing in the annotations of the reference chicken genomes (GRCg6a and GRCg7b/w). After removing redundancy, we ended up with a total of 1,420 newly annotated genes (NAGs). The NAGs tend to be found in subtelomeric regions of macro-chromosomes (chr1 to chr5, plus chrZ) and middle chromosomes (chr6 to chr13, plus chrW), as well as in micro-chromosomes (chr14 to chr39) and unplaced contigs, where G/C contents are high. Moreover, the NAGs have elevated quadruplexes G frequencies, while both G/C contents and quadruplexes G frequencies in their surrounding regions are also high. The NAGs showed tissue-specific expression, and we were able to verify 39 (92.9%) of 42 randomly selected ones in various tissues of the four chicken breeds using RT-qPCR experiments. Most of the NAGs were also encoded in the reference chicken genomes, thus, these genomes might harbor more genes than previously thought. CONCLUSION: The NAGs are widely distributed in wild, indigenous and commercial chickens, and they might play critical roles in chicken physiology. Counting these new genes, chicken genomes harbor more genes than originally thought.

An improved chromosome-level genome assembly and annotation of Echeneis naucrates.

Echeneis naucrates, as known as live sharksucker, is famous for the behavior of attaching to hosts using a highly modified dorsal fin with oval-shaped sucking disc. Here, we generated an improved high-quality chromosome-level genome assembly of E. naucrates using Illumina short reads, PacBio long reads and Hi-C data. Our assembled genome spans 572.85 Mb with a contig N50 of 23.19 Mb and is positioned to 24 pseudo-chromosomes. Additionally, at least one telomere was identified for 23 out of 24 chromosomes. Furthermore, we identified a total of 22,161 protein-coding genes, of which 21,402 genes (96.9%) were annotated successfully with functions. The combination of ab initio predictions and Repbase-based searches revealed that 15.57% of the assembled E. naucrates genome was identified as repetitive sequences. The completeness of the genome assembly and the gene annotation were estimated to be 97.5% and 95.4% with BUSCO analyses. This work enhances the utility of the live sharksucker genome and provides a valuable groundwork for the future study of genomics, biology and adaptive evolution in this species.

A chromosome-level genome assembly of the spider mite Tetranychus piercei McGregor.

Despite the rapid advances in sequencing technology, limited genomic resources are currently available for phytophagous spider mites, which include many important agricultural pests. One of these pests is Tetranychus piercei (McGregor), a serious banana pest in East Asia exhibiting remarkable tolerance to high temperature. In this study, we assembled a high-quality genome of T. piercei using a combination of PacBio long reads and Illumina short reads sequencing. With the assistance of chromatin conformation capture technology, 99.9% of the contigs were anchored into three pseudochromosomes with a total size of 86.02 Mb. Repetitive elements, accounting for 14.16% of this genome (12.20 Mb), are predominantly composed of long-terminal repeats (30.7%). By combining evidence of ab initio prediction, transcripts, and homologous proteins, we annotated 11,881 protein-coding genes. Both the genome and proteins have high BUSCO completeness scores (>94%). This high-quality genome, along with reliable annotation, provides a valuable resource for investigating the high-temperature tolerance of this species and exploring the genomic basis that underlies the host range evolution of spider mites.

Elucidation of Spatial Positioning of Ribosomes around Chromosome in Escherichia coli Cytoplasm via a Data-Informed Polymer-Based Model.

The spatial arrangement of ribosomes and chromosome in Escherichia coli's cytoplasm challenges conventional wisdom. Contrary to the notion of ribosomes acting as inert crowders to the chromosome in the cytoplasm, here we propose a nuanced view by integrating a wide array of experimental data sets into a polymer-based computer model. A set of data-informed computer simulations determines that a delicate balance of attractive and repulsive interactions between ribosomes and the chromosome is required in order to reproduce experimentally obtained linear densities and brings forth the view that ribosomes are not mere inert crowders in the cytoplasm. The model finds that the ribosomes represent themselves as a poor solvent for the chromosome with a 50 nm mesh size, consistent with previous experimental analysis. Our multidimensional analysis of ribosome distribution, both free (30S and 50S) and bound (70S polysome), uncovers a relatively less pronounced segregation pattern than previously thought. Notably, we identify a ribosome-rich central region within the innermost core of the nucleoid. Moreover, our exploration of the chromosome mesh size and the conformation of bound ribosomes suggests that these ribosomes maintain elongated shapes, enabling them to navigate through the chromosome mesh and access the central core. This dynamic localization challenges the static segregation model and underscores the pivotal role of ribosome-chromosome interactions in cellular media.

The chromosome-level genome assembly of the giant dobsonfly Acanthacorydalis orientalis (McLachlan, 1899).

Acanthacorydalis orientalis (McLachlan, 1899) (Megaloptera: Corydalidae) is an important freshwater-benthic invertebrate species that serves as an indicator for water-quality biomonitoring and is valuable for conservation from East Asia. Here, a high-quality reference genome for A. orientalis was constructed using Oxford Nanopore sequencing and High throughput Chromosome Conformation Capture (Hi-C) technology. The final genome size is 547.98 Mb, with the N50 values of contig and scaffold being 7.77 Mb and 50.53 Mb, respectively. The longest contig and scaffold are 20.57 Mb and 62.26 Mb in length, respectively. There are 99.75% contigs anchored onto 13 pseudo-chromosomes. Benchmarking Universal Single-Copy Orthologs (BUSCO) analysis showed that the completeness of the genome assembly is 99.01%. There are 10,977 protein-coding genes identified, of which 84.00% are functionally annotated. The genome contains 44.86% repeat sequences. This high-quality genome provides substantial data for future studies on population genetics, aquatic adaptation, and evolution of Megaloptera and other related insect groups.

Identification and functional analysis of a chromosome 2D fragment harboring TaFPF1 gene with the potential for yield improvement using a late heading wheat mutant.

KEY MESSAGE: A chromosome fragment influencing wheat heading and grain size was identified using mapping of m406 mutant. The study of TaFPF1 in this fragment provides more insights into wheat yield improvement. In recent years, wheat production has faced formidable challenges driven by rapid population growth and climate change, emphasizing the importance of improving specific agronomic traits such as heading date, spike length, and grain size. To identify potential genes for improving these traits, we screened a wheat EMS mutant library and identified a mutant, designated m406, which exhibited a significantly delayed heading date compared to the wild-type. Intriguingly, the mutant also displayed significantly longer spike and larger grain size. Genetic analysis revealed that a single recessive gene was responsible for the delayed heading. Surprisingly, a large 46.58 Mb deletion at the terminal region of chromosome arm 2DS in the mutant was identified through fine mapping and fluorescence in situ hybridization. Thus, the phenotypes of the mutant m406 are controlled by a group of linked genes. This deletion encompassed 917 annotated high-confidence genes, including the previously studied wheat genes Ppd1 and TaDA1, which could affect heading date and grain size. Multiple genes in this region probably contribute to the phenotypes of m406. We further investigated the function of TaFPF1 using gene editing. TaFPF1 knockout mutants showed delayed heading and increased grain size. Moreover, we identified the direct upstream gene of TaFPF1 and investigated its relationship with other important flowering genes. Our study not only identified more genes affecting heading and grain development within this deleted region but also highlighted the potential of combining these genes for improvement of wheat traits.

Chromatin looping links gene expression to the assembly of transcription factories (Review).

Genes are not randomly dispersed within the nuclear space, instead they occupy precise sites either with respect to the nuclear lamina as well as to each other. This observation stands at the basis of the today well accepted concept of nuclear territories where any chromosome shows reproducible spatial connections with a selection of others in a general picture that meets a functional criterion where genes that answer the same stimuli are grouped in the same sites. In fact, transcription is not visible widely dispersed throughout the nucleus but is gathered in several 'granules', called transcription factories that accommodates ~10 genes concurrently transcribed. This dynamic behavior of chromosomes is allowed by changes in chromatin plasticity that are governed by several classes of proteins that either modify its building or induce post­translational modifications in the protein component of nucleosomes, triggering formation of chromosome loops that modify the location of specific sites along the DNA strand. For example, transcription associated to nuclear receptors benefits of the generation of nuclear ROS that induce nicks following activation of the DNA repair apparatus that enhance helix unfolding and chromosome bridging. In the present review, the role that protocols facing elucidation of chromosome architecture are playing and will play in the near future were highlighted in order to investigate composition of the transcription factories assembled in response of a specific trigger: The estrogen­sensitive transcription was cited but the authors are convinced that the same portrait will be observed with a multitude of (if not all) other stimuli.

Chromosome-level genome assembly of Hippophae tibetana provides insights into high-altitude adaptation and flavonoid biosynthesis.

BACKGROUND: As an endemic shrub of the Qinghai-Tibetan Plateau (QTP), the distribution of Hippophae tibetana Schlecht. ranges between 2800 and 5200 m above sea level. As the most basal branch of the Hippophae genus, H. tibetana has an extensive evolutionary history. The H. tibetana is a valuable tree for studying the ecological evolution of species under extreme conditions. RESULTS: Here, we generated a high-quality chromosome-level genome of H. tibetana. The total size of the assembly genome is 917 Mb. The phylogenomic analysis of 1064 single-copy genes showed a divergence between 3.4 and 12.8 Mya for H. tibetana. Multiple gene families associated with DNA repair and disease resistance were significantly expanded in H. tibetana. We also identified many genes related to DNA repair with signs of positive selection. These results showed expansion and positive selection likely play important roles in H. tibetana's adaptation to comprehensive extreme environments in the QTP. A comprehensive genomic and transcriptomic analysis identified 49 genes involved in the flavonoid biosynthesis pathway in H. tibetana. We generated transgenic sea buckthorn hairy root producing high levels of flavonoid. CONCLUSIONS: Taken together, this H. tibetana high-quality genome provides insights into the plant adaptation mechanisms of plant under extreme environments and lay foundation for the functional genomic research and molecular breeding of H. tibetana.

Machine-learning analysis reveals an important role for negative selection in shaping cancer aneuploidy landscapes.

BACKGROUND: Aneuploidy, an abnormal number of chromosomes within a cell, is a hallmark of cancer. Patterns of aneuploidy differ across cancers, yet are similar in cancers affecting closely related tissues. The selection pressures underlying aneuploidy patterns are not fully understood, hindering our understanding of cancer development and progression. RESULTS: Here, we apply interpretable machine learning methods to study tissue-selective aneuploidy patterns. We define 20 types of features corresponding to genomic attributes of chromosome-arms, normal tissues, primary tumors, and cancer cell lines (CCLs), and use them to model gains and losses of chromosome arms in 24 cancer types. To reveal the factors that shape the tissue-specific cancer aneuploidy landscapes, we interpret the machine learning models by estimating the relative contribution of each feature to the models. While confirming known drivers of positive selection, our quantitative analysis highlights the importance of negative selection for shaping aneuploidy landscapes. This is exemplified by tumor suppressor gene density being a better predictor of gain patterns than oncogene density, and vice versa for loss patterns. We also identify the importance of tissue-selective features and demonstrate them experimentally, revealing KLF5 as an important driver for chr13q gain in colon cancer. Further supporting an important role for negative selection in shaping the aneuploidy landscapes, we find compensation by paralogs to be among the top predictors of chromosome arm loss prevalence and demonstrate this relationship for one paralog interaction. Similar factors shape aneuploidy patterns in human CCLs, demonstrating their relevance for aneuploidy research. CONCLUSIONS: Our quantitative, interpretable machine learning models improve the understanding of the genomic properties that shape cancer aneuploidy landscapes.

Entangling and disentangling mitotic chromosomes.

The DNA topological challenges generated by cellular manipulation of extremely long DNA fibers remain poorly understood. In this issue of Molecular Cell, Hildebrand et al.1 describe how mitotic chromosomes are self entangled and that disentanglement requires TOP2 activity in late mitosis.

A chromosome-level genome reveals genome evolution and molecular basis of anthraquinone biosynthesis in Rheum palmatum.

BACKGROUND: Rhubarb is one of common traditional Chinese medicine with a diverse array of therapeutic efficacies. Despite its widespread use, molecular research into rhubarb remains limited, constraining our comprehension of the geoherbalism. RESULTS: We assembled the genome of Rheum palmatum L., one of the source plants of rhubarb, to elucidate its genome evolution and unpack the biosynthetic pathways of its bioactive compounds using a combination of PacBio HiFi, Oxford Nanopore, Illumina, and Hi-C scaffolding approaches. Around 2.8 Gb genome was obtained after assembly with more than 99.9% sequences anchored to 11 pseudochromosomes (scaffold N50 = 259.19 Mb). Transposable elements (TE) with a continuous expansion of long terminal repeat retrotransposons (LTRs) is predominant in genome size, contributing to the genome expansion of R. palmatum. Totally 30,480 genes were predicted to be protein-coding genes with 473 significantly expanded gene families enriched in diverse pathways associated with high-altitude adaptation for this species. Two successive rounds of whole genome duplication event (WGD) shared by Fagopyrum tataricum and R. palmatum were confirmed. We also identified 54 genes involved in anthraquinone biosynthesis and other 97 genes entangled in flavonoid biosynthesis. Notably, RpALS emerged as a compelling candidate gene for the octaketide biosynthesis after the key residual screening. CONCLUSION: Overall, our findings offer not only an enhanced understanding of this remarkable medicinal plant but also pave the way for future innovations in its genetic breeding, molecular design, and functional genomic studies.

A chromosome-level genome assembly of the soybean pod borer: insights into larval transcriptional response to transgenic soybean expressing the pesticidal Cry1Ac protein.

BACKGROUND: Genetically modified (GM) crop plants with transgenic expression of Bacillus thuringiensis (Bt) pesticidal proteins are used to manage feeding damage by pest insects. The durability of this technology is threatened by the selection for resistance in pest populations. The molecular mechanism(s) involved in insect physiological response or evolution of resistance to Bt is not fully understood. RESULTS: To investigate the response of a susceptible target insect to Bt, the soybean pod borer, Leguminivora glycinivorella (Lepidoptera: Tortricidae), was exposed to soybean, Glycine max, expressing Cry1Ac pesticidal protein or the non-transgenic parental cultivar. Assessment of larval changes in gene expression was facilitated by a third-generation sequenced and scaffolded chromosome-level assembly of the L. glycinivorella genome (657.4 Mb; 27 autosomes + Z chromosome), and subsequent structural annotation of 18,197 RefSeq gene models encoding 23,735 putative mRNA transcripts. Exposure of L. glycinivorella larvae to transgenic Cry1Ac G. max resulted in prediction of significant differential gene expression for 204 gene models (64 up- and 140 down-regulated) and differential splicing among isoforms for 10 genes compared to unexposed cohorts. Differentially expressed genes (DEGs) included putative peritrophic membrane constituents, orthologs of Bt receptor-encoding genes previously linked or associated with Bt resistance, and those involved in stress responses. Putative functional Gene Ontology (GO) annotations assigned to DEGs were significantly enriched for 36 categories at GO level 2, respectively. Most significantly enriched cellular component (CC), biological process (BP), and molecular function (MF) categories corresponded to vacuolar and microbody, transport and metabolic processes, and binding and reductase activities. The DEGs in enriched GO categories were biased for those that were down-regulated (≥ 0.783), with only MF categories GTPase and iron binding activities were bias for up-regulation genes. CONCLUSIONS: This study provides insights into pathways and processes involved larval response to Bt intoxication, which may inform future unbiased investigations into mechanisms of resistance that show no evidence of alteration in midgut receptors.

The scales, mechanisms, and dynamics of the genome architecture.

Even when split into several chromosomes, DNA molecules that make up our genome are too long to fit into the cell nuclei unless massively folded. Such folding must accommodate the need for timely access to selected parts of the genome by transcription factors, RNA polymerases, and DNA replication machinery. Here, we review our current understanding of the genome folding inside the interphase nuclei. We consider the resulting genome architecture at three scales with a particular focus on the intermediate (meso) scale and summarize the insights gained from recent experimental observations and diverse computational models.

Chromosome Transplantation: Opportunities and Limitations.

There are thousands of rare genetic diseases that could be treated with classical gene therapy strategies such as the addition of the defective gene via viral or non-viral delivery or by direct gene editing. However, several genetic defects are too complex for these approaches. These "genomic mutations" include aneuploidies, intra and inter chromosomal rearrangements, large deletions, or inversion and copy number variations. Chromosome transplantation (CT) refers to the precise substitution of an endogenous chromosome with an exogenous one. By the addition of an exogenous chromosome and the concomitant elimination of the endogenous one, every genetic defect, irrespective of its nature, could be resolved. In the current review, we analyze the state of the art of this technique and discuss its possible application to human pathology. CT might not be limited to the treatment of human diseases. By working on sex chromosomes, we showed that female cells can be obtained from male cells, since chromosome-transplanted cells can lose either sex chromosome, giving rise to 46,XY or 46,XX diploid cells, a modification that could be exploited to obtain female gametes from male cells. Moreover, CT could be used in veterinary biology, since entire chromosomes containing an advantageous locus could be transferred to animals of zootechnical interest without altering their specific genetic background and the need for long and complex interbreeding. CT could also be useful to rescue extinct species if only male cells were available. Finally, the generation of "synthetic" cells could be achieved by repeated CT into a recipient cell. CT is an additional tool for genetic modification of mammalian cells.

Chromosome-scale assembly of the streamlined picoeukaryote Picochlorum sp. SENEW3 genome reveals Rabl-like chromatin structure and potential for C4 photosynthesis.

Genome sequencing and assembly of the photosynthetic picoeukaryotic Picochlorum sp. SENEW3 revealed a compact genome with a reduced gene set, few repetitive sequences, and an organized Rabl-like chromatin structure. Hi-C chromosome conformation capture revealed evidence of possible chromosomal translocations, as well as putative centromere locations. Maintenance of a relatively few selenoproteins, as compared to similarly sized marine picoprasinophytes Mamiellales, and broad halotolerance compared to others in Trebouxiophyceae, suggests evolutionary adaptation to variable salinity environments. Such adaptation may have driven size and genome minimization and have been enabled by the retention of a high number of membrane transporters. Identification of required pathway genes for both CAM and C4 photosynthetic carbon fixation, known to exist in the marine mamiellale pico-prasinophytes and seaweed Ulva, but few other chlorophyte species, further highlights the unique adaptations of this robust alga. This high-quality assembly provides a significant advance in the resources available for genomic investigations of this and other photosynthetic picoeukaryotes.

Insight into the adaptive role of arachnid genome-wide duplication through chromosome-level genome assembly of the Western black widow spider.

Although spiders are one of the most diverse groups of arthropods, the genetic architecture of their evolutionary adaptations is largely unknown. Specifically, ancient genome-wide duplication occurring during arachnid evolution ~450 mya resulted in a vast assembly of gene families, yet the extent to which selection has shaped this variation is understudied. To aid in comparative genome sequence analyses, we provide a chromosome-level genome of the Western black widow spider (Latrodectus hesperus)-a focus due to its silk properties, venom applications, and as a model for urban adaptation. We used long-read and Hi-C sequencing data, combined with transcriptomes, to assemble 14 chromosomes in a 1.46 Gb genome, with 38,393 genes annotated, and a BUSCO score of 95.3%. Our analyses identified high repetitive gene content and heterozygosity, consistent with other spider genomes, which has led to challenges in genome characterization. Our comparative evolutionary analyses of eight genomes available for species within the Araneoidea group (orb weavers and their descendants) identified 1,827 single-copy orthologs. Of these, 155 exhibit significant positive selection primarily associated with developmental genes, and with traits linked to sensory perception. These results support the hypothesis that several traits unique to spiders emerged from the adaptive evolution of ohnologs-or retained ancestrally duplicated genes-from ancient genome-wide duplication. These comparative spider genome analyses can serve as a model to understand how positive selection continually shapes ancestral duplications in generating novel traits today within and between diverse taxonomic groups.