Using geneid to Identify Genes

Genomic and transcriptomic analyses of a social hemipteran provide new insights into insect sociality

The origin of sociality represents one of the most important evolutionary transitions. Insect sociality evolved in some hemipteran aphids, which can produce soldiers and normal nymphs with distinct morphology and behaviour through parthenogenesis. The lack of genomic data resources has hindered the investigations into molecular mechanisms underlying their social evolution. Herein, we generated the first chromosomal-level genome of a social hemipteran (Pseudoregma bambucicola) with highly specialized soldiers and performed comparative genomic and transcriptomic analyses to elucidate the molecular signatures and regulatory mechanisms of caste differentiation. P. bambucicola has a larger known aphid genome of 582.2 Mb with an N50 length of 11.24 Mb, and about 99.6% of the assembly was anchored to six chromosomes with a scaffold N50 of 98.27 Mb. A total of 14,027 protein-coding genes were predicted and 37.33% of the assembly were identified as repeat sequences. The social evolution is accompanied by a variety of changes in genome organization, including expansion of gene families related to transcription factors, transposable elements, as well as species-specific expansions of certain sugar transporters and UGPases involved in carbohydrate metabolism. We also characterized large candidate gene sets linked to caste differentiation and found evidence of expression regulation and positive selection acting on energy metabolism and muscle structure, explaining the soldier-specific traits including morphological and behavioural specialization, developmental arrest and infertility. Overall, this study offers new insights into the molecular basis of social aphids and the evolution of insect sociality and also provides valuable data resources for further comparative and functional studies.

Long-term methylome changes after experimental seed demethylation and their interaction with recurrent water stress in Erodium cicutarium (Geraniaceae)

The frequencies and lengths of drought periods are increasing in subtropical and temperate regions worldwide. Epigenetic responses to water stress could be key for plant resilience to these largely unpredictable challenges. Experimental DNA demethylation, together with application of a stress factor is an appropriate strategy to reveal the contribution of epigenetics to plant responses to stress. We analysed leaf cytosine methylation changes in adult plants of the annual Mediterranean herb, Erodium cicutarium, in a greenhouse, after seed demethylation with 5-Azacytidine and/or recurrent water stress. We used bisulfite RADseq (BsRADseq) and a newly reported reference genome for E. cicutarium to characterize methylation changes in a 2 × 2 factorial design, controlling for plant relatedness. In the long term, 5-Azacytidine treatment alone caused both hypo- and hyper-methylation at individual cytosines, with substantial hypomethylation in CG contexts. In control conditions, drought resulted in a decrease in methylation in all but CHH contexts. In contrast, the genome of plants that experienced recurrent water stress and had been treated with 5-Azacytidine increased DNA methylation level by ca. 5%. Seed demethylation and recurrent drought produced a highly significant interaction in terms of global and context-specific cytosine methylation. Most methylation changes occurred around genic regions and within Transposable Elements. The annotation of these Differentially Methylated Regions associated with genes included several with a potential role in stress responses (e.g., PAL, CDKC, and ABCF), confirming an epigenetic contribution in response to stress at the molecular level.

Evolution of the Correlated Genomic Variation Landscape Across a Divergence Continuum in the Genus Castanopsis

The heterogeneous landscape of genomic variation has been well documented in population genomic studies. However, disentangling the intricate interplay of evolutionary forces influencing the genetic variation landscape over time remains challenging. In this study, we assembled a chromosome-level genome for Castanopsis eyrei and sequenced the whole genomes of 276 individuals from 12 Castanopsis species, spanning a broad divergence continuum. We found highly correlated genomic variation landscapes across these species. Furthermore, variations in genetic diversity and differentiation along the genome were strongly associated with recombination rates and gene density. These results suggest that long-term linked selection and conserved genomic features have contributed to the formation of a common genomic variation landscape. By examining how correlations between population summary statistics change throughout the species divergence continuum, we determined that background selection alone does not fully explain the observed patterns of genomic variation; the effects of recurrent selective sweeps must be considered. We further revealed that extensive gene flow has significantly influenced patterns of genomic variation in Castanopsis species. The estimated admixture proportion correlated positively with recombination rate and negatively with gene density, supporting a scenario of selection against gene flow. Additionally, putative introgression regions exhibited strong signals of positive selection, an enrichment of functional genes, and reduced genetic burdens, indicating that adaptive introgression has played a role in shaping the genomes of hybridizing species. This study provides insights into how different evolutionary forces have interacted in driving the evolution of the genomic variation landscape.

Single-cell RNA sequencing facilitates the elucidation of the complete biosynthesis of the antidepressant hyperforin in St. John's wort

Hyperforin is the compound responsible for the effectiveness of St. John's wort (Hypericum perforatum) as an antidepressant, but its complete biosynthetic pathway remains unknown. Gene discovery based on co-expression analysis of bulk RNA-sequencing data or genome mining failed to discover the missing steps in hyperforin biosynthesis. In this study, we sequenced the 1.54-Gb tetraploid H. perforatum genome assembled into 32 chromosomes with the scaffold N50 value of 42.44 Mb. By single-cell RNA sequencing, we identified a type of cell, "Hyper cells", wherein hyperforin biosynthesis de novo takes place in both the leaves and flowers. Through pathway reconstitution in yeast and tobacco, we identified and characterized four transmembrane prenyltransferases (HpPT1-4) that are localized at the plastid envelope and complete the hyperforin biosynthetic pathway. The hyperforin polycyclic scaffold is created by a reaction cascade involving an irregular isoprenoid coupling and a tandem cyclization. Our findings reveal how and where hyperforin is biosynthesized, enabling synthetic-biology reconstitution of the complete pathway. Thus, this study not only deepens our comprehension of specialized metabolism at the cellular level but also provides strategic guidance for elucidation of the biosynthetic pathways of other specializied metabolites in plants.

New insights into decoding the lifestyle of endophytic Fusarium lateritium Fl617 via comparing genomes

Fungal-plant interactions have persisted for 460 million years, and almost all terrestrial plants on Earth have endophytic fungi. However, the mechanism of symbiosis between endophytic fungi and host plants has been inconclusive. In this dissertation, we used a strain of endophytic Fusarium lateritium (Fl617), which was found in the previous stage to promote disease resistance in tomato, and selected the pathogenic Fusarium oxysporum Fo4287 and endophytic Fusarium oxysporum Fo47, which are in the same host and the closest relatives of Fl617, to carry out a comparative genomics analysis of the three systems and to provide a new perspective for the elucidation of the special lifestyle of the fungal endophytes. We found that endophytic F. lateritium has a smaller genome, fewer clusters and genes associated with pathogenicity, and fewer plant cell wall degrading enzymes (PCWDEs). There were also relatively fewer secondary metabolisms and typical Fusarium spp. toxins, and a lack of the key Fusarium spp. pathogenicity factor, secreted in xylem (SIX), but the endophytic fungi may be more sophisticated in their regulation of the colonization process. It is hypothesized that the endophytic fungi may have maintained their symbiosis with plants due to the relatively homogeneous microenvironment in plants for a long period of time, considering only plant interactions and discarding the relevant pathogenicity factors, and that their endophytic evolutionary tendency may tend to be genome streamlining and to enhance the fineness of the regulation of plant interactions, thus maintaining their symbiotic status with plants.

High-quality reference genome of cowpea beetle Callosobruchus maculatus

Callosobruchus maculatus is one of the most competitive stored grain pests, which causes a great loss to agricultural economy. However, due to an inadequacy of high-quality reference genome, the molecular mechanisms for olfactory and hypoxic adaptations to stored environments are unknown and require to be revealed urgently, which will contribute to the detection and prevention of the invasive pests C. maculatus. Here, we presented a high-quality chromosome-level genome of C. maculatus based on Illumina, Nanopore and Hi-C sequencing data. The total size was 1.2 Gb, and 65.17% (797.47 Mb) of it was identified to be repeat sequences. Among assembled chromosomes, chromosome 10 was considered the X chromosome according to the evidence of reads coverage and homologous genes among species. The current version of high-quality genome provides preferable data resources for the adaptive evolution research of C. maculatus.

A high-quality chromosome-level genome assembly of the endangered tree Kmeria septentrionalis

Kmeria septentrionalis is a critically endangered tree endemic to Guangxi, China, and is listed on the International Union for Conservation of Nature's Red List. The lack of genetic information and high-quality genome data has hindered conservation efforts and studies on this species. In this study, we present a chromosome-level genome assembly of K. septentrionalis. The genome was initially assembled to be 2.57 Gb, with a contig N50 of 11.93 Mb. Hi-C guided genome assembly allowed us to anchor 98.83% of the total length of the initial contigs onto 19 pseudochromosomes, resulting in a scaffold N50 of 135.08 Mb. The final chromosome-level genome, spaning 2.54 Gb, achieved a BUSCO completeness of 98.9% and contained 1.67 Gb repetitive elements and 35,927 coding genes. This high-quality genome assembly provides a valuable resource for understanding the genetic basis of conservation-related traits and biological properties of this endangered tree species. Furthermore, it lays a critical foundation for evolutionary studies within the Magnoliaceae family.

Comparative genomics analysis reveals sequence characteristics potentially related to host preference in Cryptosporidium xiaoi

Cryptosporidium spp. are important diarrhea-associated pathogens in humans and livestock. Among the known species, Cryptosporidium xiaoi, which causes cryptosporidiosis in sheep and goats, was previously recognized as a genotype of the bovine-specific Cryptosporidium bovis based on their high sequence identity in the ssrRNA gene. However, the lack of genomic data has limited characterization of the genetic differences between the two closely related species. In this study, we sequenced the genomes of two C. xiaoi isolates and performed comparative genomic analysis to identify the sequence uniqueness of this ovine-adapted species compared with other Cryptosporidium spp. Our results showed that C. xiaoi is genetically related to C. bovis as shown by their 95.8% genomic identity and similar gene content. Consistent with this, both C. xiaoi and C. bovis appear to have fewer genes encoding mitochondrial metabolic enzymes and invasion-related protein families. However, they appear to possess several species-specific genes. Further analysis indicates that the sequence differences between these two Cryptosporidium spp. are mainly in 24 highly polymorphic genes, half of which are located in the subtelomeric regions. Some of these subtelomeric genes encode secretory proteins that have undergone positive selection. In addition, the genomes of two C. xiaoi isolates, identified as subtypes XXIIIf and XXIIIh, share 99.9% nucleotide sequence identity, with six highly divergent genes encoding putative secretory proteins. Therefore, these species-specific genes and sequence polymorphism in subtelomeric genes probably contribute to the different host preference of C. xiaoi and C. bovis.

Coping with alpine habitats: genomic insights into the adaptation strategies of Triplostegia glandulifera (Caprifoliaceae)

How plants find a way to thrive in alpine habitats remains largely unknown. Here we present a chromosome-level genome assembly for an alpine medicinal herb, Triplostegia glandulifera (Caprifoliaceae), and 13 transcriptomes from other species of Dipsacales. We detected a whole-genome duplication event in T. glandulifera that occurred prior to the diversification of Dipsacales. Preferential gene retention after whole-genome duplication was found to contribute to increasing cold-related genes in T. glandulifera. A series of genes putatively associated with alpine adaptation (e.g. CBFs, ERF-VIIs, and RAD51C) exhibited higher expression levels in T. glandulifera than in its low-elevation relative, Lonicera japonica. Comparative genomic analysis among five pairs of high- vs low-elevation species, including a comparison of T. glandulifera and L. japonica, indicated that the gene families related to disease resistance experienced a significantly convergent contraction in alpine plants compared with their lowland relatives. The reduction in gene repertory size was largely concentrated in clades of genes for pathogen recognition (e.g. CNLs, prRLPs, and XII RLKs), while the clades for signal transduction and development remained nearly unchanged. This finding reflects an energy-saving strategy for survival in hostile alpine areas, where there is a tradeoff with less challenge from pathogens and limited resources for growth. We also identified candidate genes for alpine adaptation (e.g. RAD1, DMC1, and MSH3) that were under convergent positive selection or that exhibited a convergent acceleration in evolutionary rate in the investigated alpine plants. Overall, our study provides novel insights into the high-elevation adaptation strategies of this and other alpine plants.

New genes helped acorn barnacles adapt to a sessile lifestyle

Barnacles are the only sessile lineages among crustaceans, and their sessile life begins with the settlement of swimming larvae (cyprids) and the formation of protective shells. These processes are crucial for adaptation to a sessile lifestyle, but the underlying molecular mechanisms remain poorly understood. While investigating these mechanisms in the acorn barnacle, Amphibalanus amphitrite, we discovered a new gene, bcs-6, which is involved in the energy metabolism of cyprid settlement and originated from a transposon by acquiring the promoter and cis-regulatory element. Unlike mollusks, the barnacle shell comprises alternate layers of chitin and calcite and requires another new gene, bsf, which generates silk-like fibers that efficiently bind chitin and aggregate calcite in the aquatic environment. Our findings highlight the importance of exploring new genes in unique adaptative scenarios, and the results will provide important insights into gene origin and material development.

Genomic analysis reveals phylogeny of Zygophyllales and mechanism for water retention of a succulent xerophyte

Revealing the genetic basis for stress-resistant traits in extremophile plants will yield important information for crop improvement. Zygophyllum xanthoxylum, an extant species of the ancient Mediterranean, is a succulent xerophyte that can maintain a favorable water status under desert habitats; however, the genetic basis of this adaptive trait is poorly understood. Furthermore, the phylogenetic position of Zygophyllales, to which Z. xanthoxylum belongs, remains controversial. In this study, we sequenced and assembled the chromosome-level genome of Z. xanthoxylum. Phylogenetic analysis showed that Zygophyllales and Myrtales form a separated taxon as a sister to the clade comprising fabids and malvids, clarifying the phylogenetic position of Zygophyllales at whole-genome scale. Analysis of genomic and transcriptomic data revealed multiple critical mechanisms underlying the efficient osmotic adjustment using Na+ and K+ as "cheap" osmolytes that Z. xanthoxylum has evolved through the expansion and synchronized expression of genes encoding key transporters/channels and their regulators involved in Na+/K+ uptake, transport, and compartmentation. It is worth noting that ZxCNGC1;1 (cyclic nucleotide-gated channels) and ZxCNGC1;2 constituted a previously undiscovered energy-saving pathway for Na+ uptake. Meanwhile, the core genes involved in biosynthesis of cuticular wax also featured an expansion and upregulated expression, contributing to the water retention capacity of Z. xanthoxylum under desert environments. Overall, these findings boost the understanding of evolutionary relationships of eudicots, illustrate the unique water retention mechanism in the succulent xerophyte that is distinct from glycophyte, and thus provide valuable genetic resources for the improvement of stress tolerance in crops and insights into the remediation of sodic lands.

Chromosome-scale genome assembly and annotation of Cotoneaster glaucophyllus

Cotoneaster glaucophyllus is a semi-evergreen plant that blossoms in late summer, producing dense, attractive, fragrant white flowers with significant ornamental and ecological value. Here, a chromosome-scale genome assembly was obtained by integrating PacBio and Illumina sequencing data with the aid of Hi-C technology. The genome assembly was 563.3 Mb in length, with contig N50 and scaffold N50 values of ~6 Mb and ~31 Mb, respectively. Most (95.59%) of the sequences were anchored onto 17 pseudochromosomes (538.4 Mb). We predicted 35,856 protein-coding genes, 1,401 miRNAs, 655 tRNAs, 425 rRNAs, and 795 snRNAs. The functions of 34,967 genes (97.52%) were predicted. The availability of this chromosome-level genome will provide valuable resources for molecular studies of this species, facilitating future research on speciation, functional genomics, and comparative genomics within the Rosaceae family.

Genomics of the expanding pine pathogen Lecanosticta acicola reveals patterns of ongoing genetic admixture

Lecanosticta acicola is the causal agent for brown spot needle blight that affects pine trees across the northern hemisphere. Based on marker genes and microsatellite data, two distinct lineages have been identified that were introduced into Europe on two separate occasions. Despite their overall distinct geographic distribution, they have been found to coexist in regions of northern Spain and France. Here, we present the first genome-wide study of Lecanosticta acicola, including assembly of the reference genome and a population genomics analysis of 70 natural isolates from northern Spain. We show that most of the isolates belong to the southern lineage but show signs of introgression with northern lineage isolates, indicating mating between the two lineages. We also identify phenotypic differences between the two lineages based on the activity profiles of 20 enzymes, with introgressed strains being more phenotypically similar to members of the southern lineage. In conclusion, we show undergoing genetic admixture between the two main lineages of L. acicola in a region of recent expansion.

IMPORTANCE

Lecanosticta acicola is a fungal pathogen causing severe defoliation, growth reduction, and even death in more than 70 conifer species. Despite the increasing incidence of this species, little is known about its population dynamics. Two divergent lineages have been described that have now been found together in regions of France and Spain, but it is unknown how these mixed populations evolve. Here we present the first reference genome for this important plant pathogenic fungi and use it to study the population genomics of 70 isolates from an affected forest in the north of Spain. We find signs of introgression between the two main lineages, indicating that active mating is occurring in this region which could propitiate the appearance of novel traits in this species. We also study the phenotypic differences across this population based on enzymatic activities on 20 compounds.

Genomics insights into flowering and floral pattern formation: regional duplication and seasonal pattern of gene expression in Camellia

BACKGROUND

The formation and domestication of ornamental traits are influenced by various aspects, such as the recognition of esthetic values and cultural traditions. Camellia japonica is widely appreciated and domesticated around the world mainly due to its rich variations in ornamental traits. Ornamental camellias have a diverse range of resources, including different bud variations from Camellia spp. as well as inter- and intra- specific hybridization. Despite research on the formation of ornamental traits, a basic understanding of their genetics and genomics is still lacking.

RESULTS

Here, we report the chromosomal-level reference genome of C. japonica through combining multiple DNA-sequencing technologies and obtain a high-density genetic linkage map of 4255 markers by sequencing 98 interspecific F1 hybrids between C. japonica and C. chekiangoleosa. We identify two whole-genome duplication events in C. japonica: one is a shared ancient γ event, and the other is revealed to be specific to genus Camellia. Based on the micro-collinearity analysis, we find large-scale segmental duplication of chromosome 8, resulting to two copies of the AGAMOUS loci, which may play a key role in the domestication of floral shapes. To explore the regulatory mechanisms of seasonal flowering, we have analyzed year-round gene expression patterns of C. japonica and C. azalea-a sister plant of continuous flowering that has been widely used for cross breeding. Through comparative analyses of gene co-expression networks and annual gene expression patterns, we show that annual expression rhythms of some important regulators of seasonal growth and development, including GIGANTEA and CONSTANS of the photoperiod pathway, have been disrupted in C. azalea. Furthermore, we reveal that the distinctive expression patterns of FLOWERING LOCUS T can be correlated with the seasonal activities of flowering and flushing. We demonstrate that the regulatory module involved in GIGANTEA, CONSTANS, and FLOWERING LOCUS T is central to achieve seasonality.

CONCLUSIONS

Through the genomic and comparative genomics characterizations of ornamental Camellia spp., we propose that duplication of chromosomal segments as well as the establishment of gene expression patterns has played a key role in the formation of ornamental traits (e.g., flower shape, flowering time). This work provides a valuable genomic platform for understanding the molecular basis of ornamental traits.

A chromosome-level genome assembly of Agave hybrid NO.11648 provides insights into the CAM photosynthesis

The subfamily Agavoideae comprises crassulacean acid metabolism (CAM), C3, and C4 plants with a young age of speciation and slower mutation accumulation, making it a model crop for studying CAM evolution. However, the genetic mechanism underlying CAM evolution remains unclear because of lacking genomic information. This study assembled the genome of Agave hybrid NO.11648, a constitutive CAM plant belonging to subfamily Agavoideae, at the chromosome level using data generated from high-throughput chromosome conformation capture, Nanopore, and Illumina techniques, resulting in 30 pseudo-chromosomes with a size of 4.87 Gb and scaffold N50 of 186.42 Mb. The genome annotation revealed 58 841 protein-coding genes and 76.91% repetitive sequences, with the dominant repetitive sequences being the I-type repeats (Copia and Gypsy accounting for 18.34% and 13.5% of the genome, respectively). Our findings also provide support for a whole genome duplication event in the lineage leading to A. hybrid, which occurred after its divergence from subfamily Asparagoideae. Moreover, we identified a gene duplication event in the phosphoenolpyruvate carboxylase kinase (PEPCK) gene family and revealed that three PEPCK genes (PEPCK3, PEPCK5, and PEPCK12) were involved in the CAM pathway. More importantly, we identified transcription factors enriched in the circadian rhythm, MAPK signaling, and plant hormone signal pathway that regulate the PEPCK3 expression by analysing the transcriptome and using yeast one-hybrid assays. Our results shed light on CAM evolution and offer an essential resource for the molecular breeding program of Agave spp.

A high-quality chromosomal genome assembly of the sea cucumber Chiridota heheva and its hydrothermal adaptation

BACKGROUND

Chiridota heheva is a cosmopolitan holothurian well adapted to diverse deep-sea ecosystems, especially chemosynthetic environments. Besides high hydrostatic pressure and limited light, high concentrations of metal ions also represent harsh conditions in hydrothermal environments. Few holothurian species can live in such extreme conditions. Therefore, it is valuable to elucidate the adaptive genetic mechanisms of C. heheva in hydrothermal environments.

FINDINGS

Herein, we report a high-quality reference genome assembly of C. heheva from the Kairei vent, which is the first chromosome-level genome of Apodida. The chromosome-level genome size was 1.43 Gb, with a scaffold N50 of 53.24 Mb and BUSCO completeness score of 94.5%. Contig sequences were clustered, ordered, and assembled into 19 natural chromosomes. Comparative genome analysis found that the expanded gene families and positively selected genes of C. heheva were involved in the DNA damage repair process. The expanded gene families and the unique genes contributed to maintaining iron homeostasis in an iron-enriched environment. The positively selected gene RFC2 with 10 positively selected sites played an essential role in DNA repair under extreme environments.

CONCLUSIONS

This first chromosome-level genome assembly of C. heheva reveals the hydrothermal adaptation of holothurians. As the first chromosome-level genome of order Apodida, this genome will provide the resource for investigating the evolution of class Holothuroidea.

A chromosome-level genome assembly of the Asian giant softshell turtle Pelochelys cantorii

The Asian giant softshell turtle Pelochelys cantorii is one of the largest aquatic turtles in China and has been designated a First Grade Protected Animal in China. To advance conservation research, a combination of Illumina short-read, PacBio long-read, and Hi-C scaffolding technologies was used to develop a high-quality chromosome-level genome assembly for P. cantorii. A total of 262.77 Gb of clean data were produced (121.6 × depth) and then the genome was assembled into 2.16 Gb with a contig N50 of 41.44 Mb and scaffold N50 length of 120.17 Mb, respectively. Moreover, about 99.98% assembly genome sequences were clustered and ordered onto 33 pseudochromosomes. Genome annotation revealed that 21,833 protein-coding genes were predicted, and 96.40% of them were annotated. This new chromosome-level assembly will be an enabling resource for genetic and genomic studies to support fundamental insight into P. cantorii biology.

The hologenome of Daphnia magna reveals possible DNA methylation and microbiome-mediated evolution of the host genome

Properties that make organisms ideal laboratory models in developmental and medical research are often the ones that also make them less representative of wild relatives. The waterflea Daphnia magna is an exception, by both sharing many properties with established laboratory models and being a keystone species, a sentinel species for assessing water quality, an indicator of environmental change and an established ecotoxicology model. Yet, Daphnia's full potential has not been fully exploited because of the challenges associated with assembling and annotating its gene-rich genome. Here, we present the first hologenome of Daphnia magna, consisting of a chromosomal-level assembly of the D. magna genome and the draft assembly of its metagenome. By sequencing and mapping transcriptomes from exposures to environmental conditions and from developmental morphological landmarks, we expand the previously annotates gene set for this species. We also provide evidence for the potential role of gene-body DNA-methylation as a mutagen mediating genome evolution. For the first time, our study shows that the gut microbes provide resistance to commonly used antibiotics and virulence factors, potentially mediating Daphnia's environmental-driven rapid evolution. Key findings in this study improve our understanding of the contribution of DNA methylation and gut microbiota to genome evolution in response to rapidly changing environments.

Chromosome-level assembly and annotation of the Xyrichtys novacula (Linnaeus, 1758) genome

The pearly razorfish (Xyrichtys novacula), commonly known as raor in the Balearic Islands, is a wrasse within the family Labridae. This fish species has particular biological and socio-cultural characteristics making it an ideal model organism in the fields of behavioural ecology, molecular ecology and conservation biology. In this study, we present the first annotated chromosome-level assembly for this species. Sequencing involved a combination of long reads with Oxford Nanopore Technologies, Illumina paired-end short reads (2 × 151 bp), Hi-C and RNA-seq from different tissues. The nuclear genome assembly has a scaffold N50 of 34.33 Mb, a total assembly span of 775.53 Mb and 99.63% of the sequence assembled into 24 superscaffolds, consistent with its known karyotype. Quality metrics revealed a consensus accuracy (QV) of 42.92 and gene completeness > 98%. The genome annotation resulted in 26,690 protein-coding genes and 12,737 non-coding transcripts. The coding regions encoded 39,613 unique protein products, 93% of them with assigned function. Overall, the publication of the X. novacula's reference genome will broaden the scope and impact of genomic research conducted on this iconic and colourful species.

Chromosome-level assembly of Dictyophora rubrovolvata genome using third-generation DNA sequencing and Hi-C analysis

Dictyophora rubrovolvata, a rare edible mushroom with both nutritional and medicinal values, was regarded as the "queen of the mushroom" for its attractive appearance. Dictyophora rubrovolvata has been widely cultivated in China in recent years, and many researchers were focusing on its nutrition, culture condition, and artificial cultivation. Due to a lack of genomic information, research on bioactive substances, cross breeding, lignocellulose degradation, and molecular biology is limited. In this study, we report a chromosome-level reference genome of D. rubrovolvata using the PacBio single-molecule real-time-sequencing technique and high-throughput chromosome conformation capture (Hi-C) technologies. A total of 1.83 Gb circular consensus sequencing reads representing ∼983.34 coverage of the D. rubrovolvata genome were generated. The final genome was assembled into 136 contigs with a total length of 32.89 Mb. The scaffold and contig N50 length were 2.71 and 2.48 Mb, respectively. After chromosome-level scaffolding, 11 chromosomes with a total length of 28.24 Mb were constructed. Genome annotation further revealed that 9.86% of the genome was composed of repetitive sequences, and a total of 508 noncoding RNA (rRNA: 329, tRNA: 150, ncRNA: 29) were annotated. In addition, 9,725 protein-coding genes were predicted, among which 8,830 (90.79%) genes were predicted using homology or RNA-seq. Benchmarking Universal Single-Copy Orthologs results further revealed that there were 80.34% complete single-copy fungal orthologs. In this study, a total of 360 genes were annotated as belonging to the carbohydrate-active enzymes family. Further analysis also predicted 425 cytochromes P450 genes, which can be classified into 41 families. This highly accurate, chromosome-level reference genome of D. rubrovolvata will provide essential genomic information for understanding the molecular mechanism in its fruiting body formation during morphological development and facilitate the exploitation of medicinal compounds produced by this mushroom.

Assembly of Genome and Resequencing Provide Insights into Genetic Differentiation between Parents of Hulong Hybrid Grouper (Epinephelus fuscoguttatus ♀ × E. lanceolatus ♂)

The Hulong hybrid grouper was bred from the brown-marbled grouper (Epinephelus fuscoguttatus) ♀ and the giant grouper (E. lanceolatus) ♂, combining the advantageous traits of both parents. Possessing an excellent performance, this hybrid's cultivation promotes the development of the grouper industry. Its male parent, the giant grouper, possesses the fastest growth and the largest body size among all coral-reef-dwelling fish. This species is not only an economically important species in marine aquaculture, but it is also an ideal male parent in the interspecific crossing of grouper species. In the present study, a high-quality chromosome-level genome of the giant grouper was constructed with a total length of 1.06 Gb, consisting of 24 chromosomes and 69 scaffolds. To analyze the genetic differences between the parents of the Hulong hybrid grouper, the structural variations (SVs) between both parental genomes were detected, and a total of 46,643 SVs were obtained. High-quality SNPs were identified from resequencing data. There were significant differences between the two genomes, and the average FST reached 0.685. A total of 234 highly differentiated regions were detected with an FST > 0.9. The protein-coding genes involved in SVs and highly differentiated regions were significantly enriched in metabolic pathways, including fatty metabolism, carbohydrate metabolism, amino acid metabolism and the TCA cycle. These genes may be related to the differences in feeding preferences and the ability to digest carbohydrates between the two grouper species under natural conditions. In addition, protein-coding genes related to the cell cycle and p53-signaling pathway were also detected. These genes may play important roles in the regulation of body size and growth performance. This research provides genomic resources for further breeding works and evolutionary analyses.

High subtelomeric GC content in the genome of a zoonotic Cryptosporidium species

Cryptosporidium canis is a zoonotic species causing cryptosporidiosis in humans in addition to its natural hosts dogs and other fur animals. To understand the genetic basis for host adaptation, we sequenced the genomes of C. canis from dogs, minks, and foxes and conducted a comparative genomics analysis. While the genomes of C. canis have similar gene contents and organisations, they (~41.0 %) and C. felis (39.6 %) have GC content much higher than other Cryptosporidium spp. (24.3-32.9 %) sequenced to date. The high GC content is mostly restricted to subtelomeric regions of the eight chromosomes. Most of these GC-balanced genes encode Cryptosporidium-specific proteins that have intrinsically disordered regions and are involved in host-parasite interactions. Natural selection appears to play a more important role in the evolution of codon usage in GC-balanced C. canis, and most of the GC-balanced genes have undergone positive selection. While the identity in whole genome sequences between the mink- and dog-derived isolates is 99.9 % (9365 SNVs), it is only 96.0 % (362 894 SNVs) between them and the fox-derived isolate. In agreement with this, the fox-derived isolate possesses more subtelomeric genes encoding invasion-related protein families. Therefore, the change in subtelomeric GC content appears to be responsible for the more GC-balanced C. canis genomes, and the fox-derived isolate could represent a new Cryptosporidium species.

Genome Sequence Resource of a Colletotrichum graminicola Field Strain from China

The maize anthracnose stalk rot and leaf blight diseases caused by the fungal pathogen Colletotrichum graminicola is emerging as an important threat to corn production worldwide. In this work, we provide an improved genome assembly of a C. graminicola strain (TZ-3) by using the PacBio Sequel II and Illumina high-throughput sequencing technologies. The genome of TZ-3 consists of 36 contigs with a length of 59.3 Mb. After correction and evaluation with the Illumina sequencing data and BUSCO, this genome showed a high assembly quality and integrity. Gene annotation of this genome predicted 11,911 protein-coding genes, among which 983 secreted protein-coding genes and 332 effector genes were predicted. Compared with previous genomes of C. graminicola strains, TZ-3 genome is superior in nearly all parameters. The genome assembly and annotation will enhance our knowledge of the genetic makeup of the pathogen and molecular mechanisms underlying its pathogenicity and will provide valuable insights into genome variation across different regions. [Formula: see text] Copyright © 2023 The Author(s). This is an open access article distributed under the CC BY-NC-ND 4.0 International license.

Chromosome-Level Genome Assembly Provides Insights into the Evolution of the Special Morphology and Behaviour of Lepturacanthus savala

Savalani hairtail Lepturacanthus savala is a widely distributed fish along the Indo-Western Pacific coast, and contributes substantially to trichiurid fishery resources worldwide. In this study, the first chromosome-level genome assembly of L. savala was obtained by PacBio SMRT-Seq, Illumina HiSeq, and Hi-C technologies. The final assembled L. savala genome was 790.02 Mb with contig N50 and scaffold N50 values of 19.01 Mb and 32.77 Mb, respectively. The assembled sequences were anchored to 24 chromosomes by using Hi-C data. Combined with RNA sequencing data, 23,625 protein-coding genes were predicted, of which 96.0% were successfully annotated. In total, 67 gene family expansions and 93 gene family contractions were detected in the L. savala genome. Additionally, 1825 positively selected genes were identified. Based on a comparative genomic analysis, we screened a number of candidate genes associated with the specific morphology, behaviour-related immune system, and DNA repair mechanisms in L. savala. Our results preliminarily revealed mechanisms underlying the special morphological and behavioural characteristics of L. savala from a genomic perspective. Furthermore, this study provides valuable reference data for subsequent molecular ecology studies of L. savala and whole-genome analyses of other trichiurid fishes.

The Sapria himalayana genome provides new insights into the lifestyle of endoparasitic plants

BACKGROUND

Sapria himalayana (Rafflesiaceae) is an endoparasitic plant characterized by a greatly reduced vegetative body and giant flowers; however, the mechanisms underlying its special lifestyle and greatly altered plant form remain unknown. To illustrate the evolution and adaptation of S. himalayasna, we report its de novo assembled genome and key insights into the molecular basis of its floral development, flowering time, fatty acid biosynthesis, and defense responses.

RESULTS

The genome of S. himalayana is ~ 1.92 Gb with 13,670 protein-coding genes, indicating remarkable gene loss (~ 54%), especially genes involved in photosynthesis, plant body, nutrients, and defense response. Genes specifying floral organ identity and controlling organ size were identified in S. himalayana and Rafflesia cantleyi, and showed analogous spatiotemporal expression patterns in both plant species. Although the plastid genome had been lost, plastids likely biosynthesize essential fatty acids and amino acids (aromatic amino acids and lysine). A set of credible and functional horizontal gene transfer (HGT) events (involving genes and mRNAs) were identified in the nuclear and mitochondrial genomes of S. himalayana, most of which were under purifying selection. Convergent HGTs in Cuscuta, Orobanchaceae, and S. himalayana were mainly expressed at the parasite-host interface. Together, these results suggest that HGTs act as a bridge between the parasite and host, assisting the parasite in acquiring nutrients from the host.

CONCLUSIONS

Our results provide new insights into the flower development process and endoparasitic lifestyle of Rafflesiaceae plants. The amount of gene loss in S. himalayana is consistent with the degree of reduction in its body plan. HGT events are common among endoparasites and play an important role in their lifestyle adaptation.

Chromosome-level genome assembly of Lilford's wall lizard, Podarcis lilfordi (Günther, 1874) from the Balearic Islands (Spain)

The Mediterranean lizard Podarcis lilfordi is an emblematic species of the Balearic Islands. The extensive phenotypic diversity among extant isolated populations makes the species a great insular model system for eco-evolutionary studies, as well as a challenging target for conservation management plans. Here we report the first high-quality chromosome-level assembly and annotation of the P. lilfordi genome, along with its mitogenome, based on a mixed sequencing strategy (10X Genomics linked reads, Oxford Nanopore Technologies long reads and Hi-C scaffolding) coupled with extensive transcriptomic data (Illumina and PacBio). The genome assembly (1.5 Gb) is highly contiguous (N50 = 90 Mb) and complete, with 99% of the sequence assigned to candidate chromosomal sequences and >97% gene completeness. We annotated a total of 25,663 protein-coding genes translating into 38,615 proteins. Comparison to the genome of the related species Podarcis muralis revealed substantial similarity in genome size, annotation metrics, repeat content, and a strong collinearity, despite their evolutionary distance (~18-20 MYA). This genome expands the repertoire of available reptilian genomes and will facilitate the exploration of the molecular and evolutionary processes underlying the extraordinary phenotypic diversity of this insular species, while providing a critical resource for conservation genomics.

Draft Genome of White-blotched River Stingray Provides Novel Clues for Niche Adaptation and Skeleton Formation

The white-blotched river stingray (Potamotrygon leopoldi) is a cartilaginous fish native to the Xingu River, a tributary of the Amazon River system. As a rare freshwater-dwelling cartilaginous fish in the Potamotrygonidae family in which no member has the genome sequencing information available, P. leopoldi provides the evolutionary details in fish phylogeny, niche adaptation, and skeleton formation. In this study, we present its draft genome of 4.11 Gb comprising 16,227 contigs and 13,238 scaffolds, with contig N50 of 3937 kb and scaffold N50 of 5675 kb in size. Our analysis shows that P. leopoldi is a slow-evolving fish that diverged from elephant sharks about 96 million years ago. Moreover, two gene families related to the immune system (immunoglobulin heavy constant delta genes and T-cell receptor alpha/delta variable genes) exhibit expansion in P. leopoldi only. We also identified the Hox gene clusters in P. leopoldi and discovered that seven Hox genes shared by five representative fish species are missing in P. leopoldi. The RNA sequencing data from P. leopoldi and other three fish species demonstrate that fishes have a more diversified tissue expression spectrum when compared to mammals. Our functional studies suggest that lack of the gc gene encoding vitamin D-binding protein in cartilaginous fishes (both P. leopoldi and Callorhinchus milii) could partly explain the absence of hard bone in their endoskeleton. Overall, this genome resource provides new insights into the niche adaptation, body plan, and skeleton formation of P. leopoldi, as well as the genome evolution in cartilaginous fishes.

Spatio-temporal transcriptome dynamics coordinate rapid transition of core crop functions in 'lactating' pigeon

Pigeons (Columba livia) are among a select few avian species that have developed a specialized reproductive mode wherein the parents produce a 'milk' in their crop to feed newborn squabs. Nonetheless, the transcriptomic dynamics and role in the rapid transition of core crop functions during 'lactation' remain largely unexplored. Here, we generated a de novo pigeon genome assembly to construct a high resolution spatio-temporal transcriptomic landscape of the crop epithelium across the entire breeding stage. This multi-omics analysis identified a set of 'lactation'-related genes involved in lipid and protein metabolism, which contribute to the rapid functional transitions in the crop. Analysis of in situ high-throughput chromatin conformation capture (Hi-C) sequencing revealed extensive reorganization of promoter-enhancer interactions linked to the dynamic expression of these 'lactation'-related genes between stages. Moreover, their expression is spatially localized in specific epithelial layers, and can be correlated with phenotypic changes in the crop. These results illustrate the preferential de novo synthesis of 'milk' lipids and proteins in the crop, and provides candidate enhancer loci for further investigation of the regulatory elements controlling pigeon 'lactation'.

A chromosome-level genome assembly of radish (Raphanus sativus L.) reveals insights into genome adaptation and differential bolting regulation

High-quality radish (Raphanus sativus) genome represents a valuable resource for agronomical trait improvements and understanding genome evolution among Brassicaceae species. However, existing radish genome assembly remains fragmentary, which greatly hampered functional genomics research and genome-assisted breeding. Here, using a NAU-LB radish inbred line, we generated a reference genome of 476.32 Mb with a scaffold N50 of 56.88 Mb by incorporating Illumina, PacBio and BioNano optical mapping techniques. Utilizing Hi-C data, 448.12 Mb (94.08%) of the assembled sequences were anchored to nine radish chromosomes with 40 306 protein-coding genes annotated. In total, 249.14 Mb (52.31%) comprised the repetitive sequences, among which long terminal repeats (LTRs, 30.31%) were the most abundant class. Beyond confirming the whole-genome triplication (WGT) event in R. sativus lineage, we found several tandem arrayed genes were involved in stress response process, which may account for the distinctive phenotype of high disease resistance in R. sativus. By comparing against the existing Xin-li-mei radish genome, a total of 2 108 573 SNPs, 7740 large insertions, 7757 deletions and 84 inversions were identified. Interestingly, a 647-bp insertion in the promoter of RsVRN1 gene can be directly bound by the DOF transcription repressor RsCDF3, resulting into its low promoter activity and late-bolting phenotype of NAU-LB cultivar. Importantly, introgression of this 647-bp insertion allele, RsVRN1^In-536 , into early-bolting genotype could contribute to delayed bolting time, indicating that it is a potential genetic resource for radish late-bolting breeding. Together, this genome resource provides valuable information to facilitate comparative genomic analysis and accelerate genome-guided breeding and improvement in radish.

A desert lncRNA HIDEN regulates human endoderm differentiation via interacting with IMP1 and stabilizing FZD5 mRNA

BACKGROUND

Extensive studies have revealed the function and mechanism of lncRNAs in development and differentiation, but the majority have focused on those lncRNAs adjacent to protein-coding genes. In contrast, lncRNAs located in gene deserts are rarely explored. Here, we utilize multiple differentiation systems to dissect the role of a desert lncRNA, HIDEN (human IMP1-associated "desert" definitive endoderm lncRNA), in definitive endoderm differentiation from human pluripotent stem cells.

RESULTS

We show that desert lncRNAs are highly expressed with cell-stage-specific patterns and conserved subcellular localization during stem cell differentiation. We then focus on the desert lncRNA HIDEN which is upregulated and plays a vital role during human endoderm differentiation. We find depletion of HIDEN by either shRNA or promoter deletion significantly impairs human endoderm differentiation. HIDEN functionally interacts with RNA-binding protein IMP1 (IGF2BP1), which is also required for endoderm differentiation. Loss of HIDEN or IMP1 results in reduced WNT activity, and WNT agonist rescues endoderm differentiation deficiency caused by the depletion of HIDEN or IMP1. Moreover, HIDEN depletion reduces the interaction between IMP1 protein and FZD5 mRNA and causes the destabilization of FZD5 mRNA, which is a WNT receptor and necessary for definitive endoderm differentiation.

CONCLUSIONS

These data suggest that desert lncRNA HIDEN facilitates the interaction between IMP1 and FZD5 mRNA, stabilizing FZD5 mRNA which activates WNT signaling and promotes human definitive endoderm differentiation.

Chromosome-level genome assembly of Dastarcus helophoroides provides insights into CYP450 genes expression upon insecticide exposure

BACKGROUND

Dastarcus helophoroides is an important natural enemy of cerambycids, and is wildly used in biological control of pests. Nevertheless, the absence of complete genomic information limits the investigation of the underlying molecular mechanisms. Here, a chromosome-level of Dastarcus helophoroides genome is assembled using a combination strategy of Illumina, PacBio, 10x™ Genomics, and Hi-C.

RESULTS

The final assembly is 609.09 Mb with contig N50, scaffold N50 and GC content of 5.46 Mb, 42.56 Mb and 31.50%, respectively, and 95.25% of the contigs anchor into 13 chromosomes. In total 14 890 protein-coding genes and 65.37% repeat sequences are predicted in the assembly genome. The phylogenetic analysis of single-copy gene families shared among 20 insect species indicates that Dastarcus helophoroides is placed as the sister species to clade (Nitidulidae+Curculionoidea+Chrysomeloidea) + Tenebrionoidea, and diverges from the related species ~242.9 Mya. In total 36 expanded gene families are identified in Dastarcus helophoroides genome, and are functionally related to drug metabolism and metabolism of xenobiotics by cytochrome P450. Some members of CYP4 Clade and CYP6 Clade are up-regulated in Dastarcus helophoroides adults upon insecticide exposure, of which expressions of DhCYP4Q, DhCYP6A14X1 and DhCYP4C1 are significantly up-regulated. The silencing of the three genes leads to adults more sensitive to insecticide and increased knocked-down rate, which may indicate their critical roles in stress resistance and detoxication.

CONCLUSION

Our study systematically integrated the chromosome-level genome, transcriptome and gene expression of Dastarcus helophoroides, which will provide valuable resources for understanding mechanisms of pesticide metabolism, growth and development, and utilization of the natural enemy in integrated control. © 2022 Society of Chemical Industry.

Comparative genomics reveals the diversification of triterpenoid biosynthesis and origin of ocotillol-type triterpenes in Panax

Gene duplication is assumed to be the major force driving the evolution of metabolite biosynthesis in plants. Freed from functional burdens, duplicated genes can mutate toward novelties until fixed due to selective fitness. However, the extent to which this mechanism has driven the diversification of metabolite biosynthesis remains to be tested. Here we performed comparative genomics analysis and functional characterization to evaluate the impact of gene duplication on the evolution of triterpenoid biosynthesis using Panax species as models. We found that whole-genome duplications (WGDs) occurred independently in Araliaceae and Apiaceae lineages. Comparative genomics revealed the evolutionary trajectories of triterpenoid biosynthesis in plants, which was mainly promoted by WGDs and tandem duplication. Lanosterol synthase (LAS) was likely derived from a tandem duplicate of cycloartenol synthase that predated the emergence of Nymphaeales. Under episodic diversifying selection, the LAS gene duplicates produced by γ whole-genome triplication have given rise to triterpene biosynthesis in core eudicots through neofunctionalization. Moreover, functional characterization revealed that oxidosqualene cyclases (OSCs) responsible for synthesizing dammarane-type triterpenes in Panax species were also capable of producing ocotillol-type triterpenes. Genomic and biochemical evidence suggested that Panax genes encoding the above OSCs originated from the specialization of one OSC gene duplicate produced from a recent WGD shared by Araliaceae (Pg-β). Our results reveal the crucial role of gene duplication in diversification of triterpenoid biosynthesis in plants and provide insight into the origin of ocotillol-type triterpenes in Panax species.

The Torreya grandis genome illuminates the origin and evolution of gymnosperm-specific sciadonic acid biosynthesis

Torreya plants produce dry fruits with assorted functions. Here, we report the 19-Gb chromosome-level genome assembly of T. grandis. The genome is shaped by ancient whole-genome duplications and recurrent LTR retrotransposon bursts. Comparative genomic analyses reveal key genes involved in reproductive organ development, cell wall biosynthesis and seed storage. Two genes encoding a C₁₈ Δ⁹-elongase and a C₂₀ Δ⁵-desaturase are identified to be responsible for sciadonic acid biosynthesis and both are present in diverse plant lineages except angiosperms. We demonstrate that the histidine-rich boxes of the Δ⁵-desaturase are crucial for its catalytic activity. Methylome analysis reveals that methylation valleys of the T. grandis seed genome harbor genes associated with important seed activities, including cell wall and lipid biosynthesis. Moreover, seed development is accompanied by DNA methylation changes that possibly fuel energy production. This study provides important genomic resources and elucidates the evolutionary mechanism of sciadonic acid biosynthesis in land plants.

A chromosome-level genome assembly of an early matured aromatic Japonica rice variety Qigeng10 to accelerate rice breeding for high grain quality in Northeast China

Early-matured aromatic japonica rice from the Northeast is the most popular rice commodity in the Chinese market. The Qigeng10 (QG10) was one of the varieties with the largest planting area in this region in recent years. It was an early-matured japonica rice variety with a lot of superior traits such as semi-dwarf, lodging resistance, long grain, aromatic and good quality. Therefore, a high-quality assembly of Qigeng10 genome is critical and useful for japonica research and breeding. In this study, we produced a high-precision QG10 chromosome-level genome by using a combination of Nanopore and Hi-C platforms. Finally, we assembled the QG10 genome into 77 contigs with an N50 length of 11.80 Mb in 27 scaffolds with an N50 length of 30.55 Mb. The assembled genome size was 378.31Mb with 65 contigs and constituted approximately 99.59% of the 12 chromosomes. We identified a total of 1,080,819 SNPs and 682,392 InDels between QG10 and Nipponbare. We also annotated 57,599 genes by the Ab initio method, homology-based technique, and RNA-seq. Based on the assembled genome sequence, we detected the sequence variation in a total of 63 cloned genes involved in grain yield, grain size, disease tolerance, lodging resistance, fragrance, and many other important traits. Finally, we identified five elite alleles (qTGW2^Nipponbare , qTGW3^Nanyangzhan , GW5^IR24 , GW6^Suyunuo , and qGW8^Basmati385 ) controlling long grain size, four elite alleles (COLD1^Nipponbare , bZIP73^Nipponbare , CTB4a^{Kunmingxiaobaigu} , and CTB2^{Kunmingxiaobaigu} ) controlling cold tolerance, three non-functional alleles (DTH7^Kitaake , Ghd7^Hejiang19 , and Hd1^Longgeng31 ) for early heading, two resistant alleles (Pia^Akihikari and Pid4^Digu ) for rice blast, a resistant allele STV11^Kasalath for rice stripe virus, an NRT1.1B^IR24 allele for higher nitrate absorption activity, an elite allele SCM3^Chugoku117 for stronger culms, and the typical aromatic gene badh2-E2 for fragrance in QG10. These results not only help us to better elucidate the genetic mechanisms underlying excellent agronomic traits in QG10 but also have wide-ranging implications for genomics-assisted breeding in early-matured fragrant japonica rice.

Genomic, transcriptomic, and metabolomic analysis of Oldenlandia corymbosa reveals the biosynthesis and mode of action of anti-cancer metabolites

Plants accumulate a vast array of secondary metabolites, which constitute a natural resource for pharmaceuticals. Oldenlandia corymbosa belongs to the Rubiaceae family, and has been used in traditional medicine to treat different diseases, including cancer. However, the active metabolites of the plant, their biosynthetic pathway and mode of action in cancer are unknown. To fill these gaps, we exposed this plant to eight different stress conditions and combined different omics data capturing gene expression, metabolic profiles, and anti-cancer activity. Our results show that O. corymbosa extracts are active against breast cancer cell lines and that ursolic acid is responsible for this activity. Moreover, we assembled a high-quality genome and uncovered two genes involved in the biosynthesis of ursolic acid. Finally, we also revealed that ursolic acid causes mitotic catastrophe in cancer cells and identified three high-confidence protein binding targets by Cellular Thermal Shift Assay (CETSA) and reverse docking. Altogether, these results constitute a valuable resource to further characterize the biosynthesis of active metabolites in the Oldenlandia group, while the mode of action of ursolic acid will allow us to further develop this valuable compound.

Citrus genomic resources unravel putative genetic determinants of Huanglongbing pathogenicity

Citrus HLB caused by Candidatus Liberibacter asiaticus is a pathogen-triggered immune disease. Here, we identified putative genetic determinants of HLB pathogenicity by integrating citrus genomic resources to characterize the pan-genome of accessions that differ in their response to HLB. Genome-wide association mapping and analysis of allele-specific expression between susceptible, tolerant, and resistant accessions further refined candidates underlying the response to HLB. We first developed a phased diploid assembly of Citrus sinensis 'Newhall' genome and produced resequencing data for 91 citrus accessions that differ in their response to HLB. These data were combined with previous resequencing data from 356 accessions for genome-wide association mapping of the HLB response. Genes determinants for HLB pathogenicity were associated with host immune response, ROS production, and antioxidants. Overall, this study has provided a significant resource of citrus genomic data and identified candidate genes to be further explored to understand the genetic determinants of HLB pathogenicity.

Large-Scale Chromosomal Changes Lead to Genome-Level Expression Alterations, Environmental Adaptation, and Speciation in the Gayal (Bos frontalis)

Determining the functional consequences of karyotypic changes is invariably challenging because evolution tends to obscure many of its own footprints, such as accumulated mutations, recombination events, and demographic perturbations. Here, we describe the assembly of a chromosome-level reference genome of the gayal (Bos frontalis) thereby revealing the structure, at base-pair-level resolution, of a telo/acrocentric-to-telo/acrocentric Robertsonian translocation (2;28) (T/A-to-T/A rob[2;28]). The absence of any reduction in the recombination rate or genetic introgression within the fusion region of gayal served to challenge the long-standing view of a role for fusion-induced meiotic dysfunction in speciation. The disproportionate increase noted in the distant interactions across pro-chr2 and pro-chr28, and the change in open-chromatin accessibility following rob(2;28), may, however, have led to the various gene expression irregularities observed in the gayal. Indeed, we found that many muscle-related genes, located synthetically on pro-chr2 and pro-chr28, exhibited significant changes in expression. This, combined with genome-scale structural variants and expression alterations in genes involved in myofibril composition, may have driven the rapid sarcomere adaptation of gayal to its rugged mountain habitat. Our findings not only suggest that large-scale chromosomal changes can lead to alterations in genome-level expression, thereby promoting both adaptation and speciation, but also illuminate novel avenues for studying the relationship between karyotype evolution and speciation.

The gold-ringed octopus (Amphioctopus fangsiao) genome and cerebral single-nucleus transcriptomes provide insights into the evolution of karyotype and neural novelties

BACKGROUND

Coleoid cephalopods have distinctive neural and morphological characteristics compared to other invertebrates. Early studies reported massive genomic rearrangements occurred before the split of octopus and squid lineages (Proc Natl Acad Sci U S A 116:3030-5, 2019), which might be related to the neural innovations of their brain, yet the details remain elusive. Here we combine genomic and single-nucleus transcriptome analyses to investigate the octopod chromosome evolution and cerebral characteristics.

RESULTS

We present a chromosome-level genome assembly of a gold-ringed octopus, Amphioctopus fangsiao, and a single-nucleus transcriptome of its supra-esophageal brain. Chromosome-level synteny analyses estimate that the chromosomes of the ancestral octopods experienced multiple chromosome fission/fusion and loss/gain events by comparing with the nautilus genome as outgroup, and that a conserved genome organization was detected during the evolutionary process from the last common octopod ancestor to their descendants. Besides, protocadherin, GPCR, and C2H2 ZNF genes are thought to be highly related to the neural innovations in cephalopods (Nature 524:220-4, 2015), and the chromosome analyses pinpointed several collinear modes of these genes on the octopod chromosomes, such as the collinearity between PCDH and C2H2 ZNF, as well as between GPCR and C2H2 ZNF. Phylogenetic analyses show that the expansion of the octopod protocadherin genes is driven by a tandem-duplication mechanism on one single chromosome, including two separate expansions at 65 million years ago (Ma) and 8-14 Ma, respectively. Furthermore, we identify eight cell types (i.e., cholinergic and glutamatergic neurons) in the supra-esophageal brain of A. fangsiao, and the single-cell expression analyses reveal the co-expression of protocadherin and GPCR in specific neural cells, which may contribute to the neural development and signal transductions in the octopod brain.

CONCLUSIONS

The octopod genome analyses reveal the dynamic evolutionary history of octopod chromosomes and neural-related gene families. The single-nucleus transcriptomes of the supra-esophageal brain indicate their cellular heterogeneities and functional interactions with other tissues (i.e., gill), which provides a foundation for further octopod cerebral studies.

Chromosome-scale assemblies of the male and female Populus euphratica genomes reveal the molecular basis of sex determination and sexual dimorphism

Reference-quality genomes of both sexes are essential for studying sex determination and sex-chromosome evolution, as their gene contents and expression profiles differ. Here, we present independent chromosome-level genome assemblies for the female (XX) and male (XY) genomes of desert poplar (Populus euphratica), resolving a 22.7-Mb X and 24.8-Mb Y chromosome. We also identified a relatively complete 761-kb sex-linked region (SLR) in the peritelomeric region on chromosome 14 (Y). Within the SLR, recombination around the partial repeats for the feminizing factor ARR17 (ARABIDOPSIS RESPONSE REGULATOR 17) was potentially suppressed by flanking palindromic arms and the dense accumulation of retrotransposons. The inverted small segments S1 and S2 of ARR17 exhibited relaxed selective pressure and triggered sex determination by generating 24-nt small interfering RNAs that induce male-specific hyper-methylation at the promoter of the autosomal targeted ARR17. We also detected two male-specific fusion genes encoding proteins with NB-ARC domains at the breakpoint region of an inversion in the SLR that may be responsible for the observed sexual dimorphism in immune responses. Our results show that the SLR appears to follow proposed evolutionary dynamics for sex chromosomes and advance our understanding of sex determination and the evolution of sex chromosomes in Populus.

Reference-quality genomes of both sexes of the dioecious tree species, Populus euphratica, provide further insight into the evolution of Populus sex chromosomes and highlight male-specific fusion genes that may contribute to sexual dimorphism.

An improved, chromosome-level genome of the giant panda (Ailuropoda melanoleuca)

BACKGROUND

The iconic giant panda (Ailuropoda melanoleuca), as both a flagship and umbrella species endemic to China, is a world famous symbol for wildlife conservation. The giant panda has several specific biological traits and holds a relatively small place in evolution. A high-quality genome of the giant panda is key to understanding the biology of this vulnerable species.

FINDINGS

We generated a 2.48-Gb chromosome-level genome (GPv1) of the giant panda named "Jing Jing" with a contig N50 of 28.56 Mb and scaffold N50 of 134.17 Mb, respectively. The total length of chromosomes (n = 21) was 2.39-Gb, accounting for 96.4% of the whole genome. Compared with the previously published four genomes of the giant panda, our genome is characterized by the highest completeness and the correct sequence orientation. A gap-free and 850 kb length of immunoglobulin heavy-chain gene cluster was manually annotated in close proximity to the telomere of chromosome 14. Additionally, we developed an algorithm to predict the centromere position of each chromosome. We also constructed a complete chromatin structure for "Jing Jing", which includes inter-chromosome interaction pattern, A/B compartment, topologically associated domain (TAD), TAD-clique and promoter-enhancer interaction (PEI).

CONCLUSIONS

We presented an improved chromosome-level genome and complete chromatin structure for the giant panda. This is a valuable resource for the future genetic and genomic studies on giant panda.

Identification of Novel Genes and Proteoforms in Angiostrongylus costaricensis through a Proteogenomic Approach

RNA sequencing (RNA-Seq) and mass-spectrometry-based proteomics data are often integrated in proteogenomic studies to assist in the prediction of eukaryote genome features, such as genes, splicing, single-nucleotide (SNVs), and single-amino-acid variants (SAAVs). Most genomes of parasite nematodes are draft versions that lack transcript- and protein-level information and whose gene annotations rely only on computational predictions. Angiostrongylus costaricensis is a roundworm species that causes an intestinal inflammatory disease, known as abdominal angiostrongyliasis (AA). Currently, there is no drug available that acts directly on this parasite, mostly due to the sparse understanding of its molecular characteristics. The available genome of A. costaricensis, specific to the Costa Rica strain, is a draft version that is not supported by transcript- or protein-level evidence. This study used RNA-Seq and MS/MS data to perform an in-depth annotation of the A. costaricensis genome. Our prediction improved the reference annotation with (a) novel coding and non-coding genes; (b) pieces of evidence of alternative splicing generating new proteoforms; and (c) a list of SNVs between the Brazilian (Crissiumal) and the Costa Rica strain. To the best of our knowledge, this is the first time that a multi-omics approach has been used to improve the genome annotation of A. costaricensis. We hope this improved genome annotation can assist in the future development of drugs, kits, and vaccines to treat, diagnose, and prevent AA caused by either the Brazil strain (Crissiumal) or the Costa Rica strain.

Genomes of Two Flying Squid Species Provide Novel Sights into Adaptations of Cephalopods to Pelagic Life

Pelagic cephalopods have evolved a series of fascinating traits, such as excellent visual acuity, high-speed agility, and photophores for adaptation to open pelagic oceans. However, the genetic mechanisms underpinning these traits are not well understood. Thus, in this study, we obtained high-quality genomes of two purpleback flying squid species (Sthenoteuthis oualaniensis and Sthenoteuthis sp.), with sizes of 5450 Mb and 5651 Mb, respectively. Comparative genomic analyses revealed that the S-crystallin subfamily SL20-1 associated with visual acuity in the purpleback flying squid lineage was significantly expanded, and the evolution of high-speed agility for the species was accompanied by significant positive selection pressure on genes related to energy metabolism. These molecular signals might have contributed to the evolution of their adaptative predatory and anti-predatory traits. In addition, the transcriptomic analysis provided clear indications of the evolution of the photophores of purpleback flying squids, especially the recruitment of new genes and energy metabolism-related genes which may have played key functional roles in the process.

Long read genome assemblies complemented by single cell RNA-sequencing reveal genetic and cellular mechanisms underlying the adaptive evolution of yak

Wild yak (Bos mutus) and domestic yak (Bos grunniens) are adapted to high altitude environment and have ecological, economic, and cultural significances on the Qinghai-Tibetan Plateau (QTP). Currently, the genetic and cellular bases underlying adaptations of yak to extreme conditions remains elusive. In the present study, we assembled two chromosome-level genomes, one each for wild yak and domestic yak, and screened structural variants (SVs) through the long-read data of yak and taurine cattle. The results revealed that 6733 genes contained high-FST SVs. 127 genes carrying special type of SVs were differentially expressed in lungs of the taurine cattle and yak. We then constructed the first single-cell gene expression atlas of yak and taurine cattle lung tissues and identified a yak-specific endothelial cell subtype. By integrating SVs and single-cell transcriptome data, we revealed that the endothelial cells expressed the highest proportion of marker genes carrying high-FST SVs in taurine cattle lungs. Furthermore, we identified pathways which were related to the medial thickness and formation of elastic fibers in yak lungs. These findings provide new insights into the high-altitude adaptation of yak and have important implications for understanding the physiological and pathological responses of large mammals and humans to hypoxia.

The genetic bases of yak adaptations to extreme conditions remains elusive. This study compares yak and cattle at a genomic and transcriptomic level, revealing a new type of endothelial cell and candidate genes related with elastic fiber formation in yak lungs that might contribute to high altitude adaptation.

De novo assembly of two chromosome-level rice genomes and bin-based QTL mapping reveal genetic diversity of grain weight trait in rice

Following the "green revolution," indica and japonica hybrid breeding has been recognized as a new breakthrough in further improving rice yields. However, heterosis-related grain weight QTLs and the basis of yield advantage among subspecies has not been well elucidated. We herein de novo assembled the chromosome level genomes of an indica/xian rice (Luohui 9) and a japonica/geng rice (RPY geng) and found that gene number differences and structural variations between these two genomes contribute to the differences in agronomic traits and also provide two different favorable allele pools to produce better derived recombinant inbred lines (RILs). In addition, we generated a high-generation (> F₁₅) population of 272 RILs from the cross between Luohui 9 and RPY geng and two testcross hybrid populations derived from the crosses of RILs and two cytoplasmic male sterile lines (YTA, indica and Z7A, japonica). Based on three derived populations, we totally identified eight 1,000-grain weight (KGW) QTLs and eight KGW heterosis loci. Of QTLs, qKGW-6.1 and qKGW-8.1 were accepted as novel KGW QTLs that have not been reported previously. Interestingly, allele genotyping results revealed that heading date related gene (Ghd8) in qKGW-8.1 and qLH-KGW-8.1, can affect grain weight in RILs and rice core accessions and may also play an important role in grain weight heterosis. Our results provided two high-quality genomes and novel gene editing targets for grain weight for future rice yield improvement project.

The chromosome-level holly (Ilex latifolia) genome reveals key enzymes in triterpenoid saponin biosynthesis and fruit color change

The Ilex L. (hollies) genus of Aquifoliaceae shows high species diversity in tropical and subtropical regions of Asia and South America. Throughout the range of the genus, Ilex species have been widely used in beverage and medicine production and as ornamentals. Here, we assembled a high-quality, chromosome-level genome of Ilex latifolia, which has extremely high economic value because of its useful secondary metabolite production and the high ornamental value of its decorative red berries. The 99.8% genome sequence was anchored to 20 pseudochromosomes, with a total length of 766.02 Mb and a scaffold N50 of 33.45 Mb. Based on the comparative genomic analysis of 14 angiosperm species, we recovered I. latifolia as the sister group to all other campanulids. Two whole-genome duplication (WGD) events were identified in hollies: one shared ancient WGD in the ancestor of all eudicots and a recent and independent WGD in hollies. We performed a genome-wide search to screen candidate genes involved in the biosynthesis of pentacyclic triterpenoid saponins in I. latifolia. Three subfamilies of CYP450 (CYP71A, CYP72A, and CYP716A) appear to have expanded. The transcriptomic analysis of I. latifolia leaves at five developmental stages revealed that two CYP716A genes and one CYP72A gene probably play important roles in this biosynthetic pathway. In addition, we totally identified 12 genes in the biosynthesis pathways of pelargonidin and cyanidin and observed their differential expression in green and red fruit pericarps, suggesting an association between pelargonidin and cyanidin biosynthesis and fruit pericarp color change. The accumulation of pelargonidin and cyanidin is expected to play an important role in the ornamental value of I. latifolia. Altogether, this study elucidated the molecular basis of the medicinal and ornamental value of I. latifolia, providing a data basis and promising clues for further applications.

Assembly of high-quality genomes of the locoweed Oxytropis ochrocephala and its endophyte Alternaria oxytropis provides new evidence for their symbiotic relationship and swainsonine biosynthesis

Locoweeds are perennial forbs poisonous to livestock and cause extreme losses to animal husbandry. Locoweed toxicity is attributed to the symbiotic endophytes in Alternaria sect. Undifilum, which produce a mycotoxin swainsonine (SW). We performed a de novo whole genome sequencing of the most common locoweed in China, Oxytropis ochrocephala (2n = 16), and assembled a high-quality, chromosome-level reference genome. Its genome size is 958.83 Mb with 930.94 Mb (97.09 %) anchored and oriented onto 8 chromosomes, and 31,700 protein-coding genes were annotated. Phylogenetic and collinearity analysis showed it is closely related to Medicago truncatula with a pair of large interchromosomal rearrangements, and both species underwent a whole-genome duplication event. We also derived the genome of A. oxytropis at 74.48 Mb with a contig N50 of 8.87 Mb and 10,657 protein-coding genes, and refined the genes of SW biosynthesis. Multiple Alternaria species containing the swnK gene were grouped into a single clade, but in other genera, swnK's homologues are diverse. Resequencing of 41 A. oxytropis strains revealed one SNP in the SWN cluster causing changes in SW concentration. Comparing the transcriptomes of symbiotic and non-symbiotic interactions identified differentially expressed genes (DEG) linked to defense and secondary metabolism in the host. Within the endophyte DEGs were linked to cell wall degradation, fatty acids and nitrogen metabolism. Symbiosis induced the up-regulation of most of the SW biosynthetic genes. These two genomes and relevant sequencing data should provide valuable genetic resources for the study of the evolution, interaction, and SW biosynthesis in the symbiont.

Chromosome-level genome assembly of largemouth bass (Micropterus salmoides) using PacBio and Hi-C technologies

The largemouth bass (Micropterus salmoides) has become a cosmopolitan species due to its widespread introduction as game or domesticated fish. Here a high-quality chromosome-level reference genome of M. salmoides was produced by combining Illumina paired-end sequencing, PacBio single molecule sequencing technique (SMRT) and High-through chromosome conformation capture (Hi-C) technologies. Ultimately, the genome was assembled into 844.88 Mb with a contig N50 of 15.68 Mb and scaffold N50 length of 35.77 Mb. About 99.9% assembly genome sequences (844.00 Mb) could be anchored to 23 chromosomes, and 98.03% assembly genome sequences could be ordered and directed. The genome contained 38.19% repeat sequences and 2693 noncoding RNAs. A total of 26,370 protein-coding genes from 3415 gene families were predicted, of which 97.69% were functionally annotated. The high-quality genome assembly will be a fundamental resource to study and understand how M. salmoides adapt to novel and changing environments around the world, and also be expected to contribute to the genetic breeding and other research.

Chromosome-level genome assembly defines female-biased genes associated with sex determination and differentiation in the human blood fluke Schistosoma japonicum

Schistosomiasis is a neglected tropical disease of humans caused by blood flukes of the genus Schistosoma, the only dioecious parasitic flatworm. Although aspects of sex determination, differentiation and reproduction have been studied in some Schistosoma species, almost nothing is known for Schistosoma japonicum, the causative agent of schistosomiasis japonica. This mainly reflects the lack of high-quality genomic and transcriptomic resources for this species. As current genomes for S. japonicum are highly fragmented, we assembled and report a chromosome-level reference genome (seven autosomes, the Z-chromosome and partial W-chromosome), achieving a substantially enhanced gene annotation. Utilizing this genome, we discovered that the sex chromosomes of S. japonicum and its congener S. mansoni independently suppressed recombination during evolution, forming five and two evolutionary strata, respectively. By exploring the W-chromosome and sex-specific transcriptomes, we identified 35 W-linked genes and 257 female-preferentially transcribed genes (FTGs) from our chromosomal assembly and uncovered a signature for sex determination and differentiation in S. japonicum. These FTGs clustering within autosomes or the Z-chromosome exhibit a highly dynamic transcription profile during the pairing of female and male schistosomula, thereby representing a critical phase for the maturation of the female worms and suggesting distinct layers of regulatory control of gene transcription at this development stage. Collectively, these data provide a valuable resource for further functional genomic characterization of S. japonicum, shed light on the evolution of sex chromosomes in this highly virulent human blood fluke, and provide a pathway to identify novel targets for development of intervention tools against schistosomiasis.

Chromosome‐level genome assembly for the horned‐gall aphid provides insights into interactions between gall‐making insect and its host plant

The aphid Schlechtendalia chinensis is an economically important insect that can induce horned galls, which are valuable for the medicinal and chemical industries. Up to now, more than twenty aphid genomes have been reported. Most of the sequenced genomes are derived from free‐living aphids. Here, we generated a high‐quality genome assembly from a galling aphid. The final genome assembly is 271.52 Mb, representing one of the smallest sequenced genomes of aphids. The genome assembly is based on contig and scaffold N50 values of the genome sequence are 3.77 Mb and 20.41 Mb, respectively. Nine‐seven percent of the assembled sequences was anchored onto 13 chromosomes. Based on BUSCO analysis, the assembly involved 96.9% of conserved arthropod and 98.5% of the conserved Hemiptera single‐copy orthologous genes. A total of 14,089 protein‐coding genes were predicted. Phylogenetic analysis revealed that S. chinensis diverged from the common ancestor of Eriosoma lanigerum approximately 57 million years ago (MYA). In addition, 35 genes encoding salivary gland proteins showed differentially when S. chinensis forms a gall, suggesting they have potential roles in gall formation and plant defense suppression. Taken together, this high‐quality S. chinensis genome assembly and annotation provide a solid genetic foundation for future research to reveal the mechanism of gall formation and to explore the interaction between aphids and their host plants.

The Gastrodia menghaiensis (Orchidaceae) genome provides new insights of orchid mycorrhizal interactions

BACKGROUND

To illustrate the molecular mechanism of mycoheterotrophic interactions between orchids and fungi, we assembled chromosome-level reference genome of Gastrodia menghaiensis (Orchidaceae) and analyzed the genomes of two species of Gastrodia.

RESULTS

Our analyses indicated that the genomes of Gastrodia are globally diminished in comparison to autotrophic orchids, even compared to Cuscuta (a plant parasite). Genes involved in arbuscular mycorrhizae colonization were found in genomes of Gastrodia, and many of the genes involved biological interaction between Gatrodia and symbiotic microbionts are more numerous than in photosynthetic orchids. The highly expressed genes for fatty acid and ammonium root transporters suggest that fungi receive material from orchids, although most raw materials flow from the fungi. Many nuclear genes (e.g. biosynthesis of aromatic amino acid L-tryptophan) supporting plastid functions are expanded compared to photosynthetic orchids, an indication of the importance of plastids even in totally mycoheterotrophic species.

CONCLUSION

Gastrodia menghaiensis has the smallest proteome thus far among angiosperms. Many of the genes involved biological interaction between Gatrodia and symbiotic microbionts are more numerous than in photosynthetic orchids.