The enormous repetitive Antarctic krill genome reveals environmental adaptations and population insights.

Antarctic krill (Euphausia superba) is Earth's most abundant wild animal, and its enormous biomass is vital to the Southern Ocean ecosystem. Here, we report a 48.01-Gb chromosome-level Antarctic krill genome, whose large genome size appears to have resulted from inter-genic transposable element expansions. Our assembly reveals the molecular architecture of the Antarctic krill circadian clock and uncovers expanded gene families associated with molting and energy metabolism, providing insights into adaptations to the cold and highly seasonal Antarctic environment. Population-level genome re-sequencing from four geographical sites around the Antarctic continent reveals no clear population structure but highlights natural selection associated with environmental variables. An apparent drastic reduction in krill population size 10 mya and a subsequent rebound 100 thousand years ago coincides with climate change events. Our findings uncover the genomic basis of Antarctic krill adaptations to the Southern Ocean and provide valuable resources for future Antarctic research.

African lungfish genome sheds light on the vertebrate water-to-land transition.

Lungfishes are the closest extant relatives of tetrapods and preserve ancestral traits linked with the water-to-land transition. However, their huge genome sizes have hindered understanding of this key transition in evolution. Here, we report a 40-Gb chromosome-level assembly of the African lungfish (Protopterus annectens) genome, which is the largest genome assembly ever reported and has a contig and chromosome N50 of 1.60 Mb and 2.81 Gb, respectively. The large size of the lungfish genome is due mainly to retrotransposons. Genes with ultra-long length show similar expression levels to other genes, indicating that lungfishes have evolved high transcription efficacy to keep gene expression balanced. Together with transcriptome and experimental data, we identified potential genes and regulatory elements related to such terrestrial adaptation traits as pulmonary surfactant, anxiolytic ability, pentadactyl limbs, and pharyngeal remodeling. Our results provide insights and key resources for understanding the evolutionary pathway leading from fishes to humans.

A chickpea genetic variation map based on the sequencing of 3,366 genomes.

Zero hunger and good health could be realized by 2030 through effective conservation, characterization and utilization of germplasm resources1. So far, few chickpea (Cicer arietinum) germplasm accessions have been characterized at the genome sequence level2. Here we present a detailed map of variation in 3,171 cultivated and 195 wild accessions to provide publicly available resources for chickpea genomics research and breeding. We constructed a chickpea pan-genome to describe genomic diversity across cultivated chickpea and its wild progenitor accessions. A divergence tree using genes present in around 80% of individuals in one species allowed us to estimate the divergence of Cicer over the last 21 million years. Our analysis found chromosomal segments and genes that show signatures of selection during domestication, migration and improvement. The chromosomal locations of deleterious mutations responsible for limited genetic diversity and decreased fitness were identified in elite germplasm. We identified superior haplotypes for improvement-related traits in landraces that can be introgressed into elite breeding lines through haplotype-based breeding, and found targets for purging deleterious alleles through genomics-assisted breeding and/or gene editing. Finally, we propose three crop breeding strategies based on genomic prediction to enhance crop productivity for 16 traits while avoiding the erosion of genetic diversity through optimal contribution selection (OCS)-based pre-breeding. The predicted performance for 100-seed weight, an important yield-related trait, increased by up to 23% and 12% with OCS- and haplotype-based genomic approaches, respectively.

A new invertebrate NPY-like polypeptide, ZoaNPY, from the Zoanthus sociatus, as a novel ligand of human NPY Y2 receptor rescues vascular insufficiency via PLC/PKC and Src- FAK-dependent signaling pathways.

Our recent multi-omics studies have revealed rich sources of novel bioactive proteins and polypeptides from marine organisms including cnidarians. In the present study, we initially conducted a transcriptomic analysis to review the composition profile of polypeptides from Zoanthus sociatus. Then, a newly discovered NPY-like polypeptide-ZoaNPY was selected for further in silico structural, binding and virtually pharmacological studies. To evaluate the pro-angiogenic effects of ZoaNPY, we employed an in vitro HUVECs model and an in vivo zebrafish model. Our results indicate that ZoaNPY, at 1-100 pmol, enhances cell survival, migration and tube formation in the endothelial cells. Besides, treatment with ZoaNPY could restore a chemically-induced vascular insufficiency in zebrafish embryos. Western blot results demonstrated the application of ZoaNPY could increase the phosphorylation of proteins related to angiogenesis signaling including PKC, PLC, FAK, Src, Akt, mTOR, MEK, and ERK1/2. Furthermore, through molecular docking and surface plasmon resonance (SPR) verification, ZoaNPY was shown to directly and physically interact with NPY Y2 receptor. In view of this, all evidence showed that the pro-angiogenic effects of ZoaNPY involve the activation of NPY Y2 receptor, thereby activating the Akt/mTOR, PLC/PKC, ERK/MEK and Src- FAK-dependent signaling pathways. Furthermore, in an excision wound model, the treatment with ZoaNPY was shown to accelerate the wound healing process in mice. Our findings provide new insights into the discovery and development of novel pro-angiogenic drugs derived from NPY-like polypeptides in the future.

Single-cell transcriptome landscape of zebrafish liver reveals hepatocytes and immune cell interactions in understanding nonalcoholic fatty liver disease.

Nonalcoholic fatty liver disease (NAFLD) is becoming the most common chronic liver disease in the world. Immunity is the major contributing factor in NAFLD; however, the interaction of immune cells and hepatocytes in disease progression has not been fully elucidated. As a popular species for studying NAFLD, zebrafish, whose liver is a complex immune system mediated by immune cells and non-immune cells in maintaining immune tolerance and homeostasis. Understanding the cellular composition and immune environment of zebrafish liver is of great significance for its application in NAFLD. Here, we established a liver atlas that consists of 10 cell types using single-cell RNA sequencing (scRNA-seq). By examining the heterogeneity of hepatocytes and analyzing the expression of NAFLD-associated genes in the specific cluster, we provide a potential target cell model to study NAFLD. Additionally, our analysis identified two subtypes of distinct resident macrophages with inflammatory and non-inflammatory functions and characterized the successive stepwise development of T cell subclusters in the liver. Importantly, we uncovered the possible regulation of macrophages and T cells on target cells of fatty liver by analyzing the cellular interaction between hepatocytes and immune cells. Our data provide valuable information for an in-depth study of immune cells targeting hepatocytes to regulate the immune balance in NAFLD.

Comparative genomics provides insights into the aquatic adaptations of mammals.

The ancestors of marine mammals once roamed the land and independently committed to an aquatic lifestyle. These macroevolutionary transitions have intrigued scientists for centuries. Here, we generated high-quality genome assemblies of 17 marine mammals (11 cetaceans and six pinnipeds), including eight assemblies at the chromosome level. Incorporating previously published data, we reconstructed the marine mammal phylogeny and population histories and identified numerous idiosyncratic and convergent genomic variations that possibly contributed to the transition from land to water in marine mammal lineages. Genes associated with the formation of blubber (NFIA), vascular development (SEMA3E), and heat production by brown adipose tissue (UCP1) had unique changes that may contribute to marine mammal thermoregulation. We also observed many lineage-specific changes in the marine mammals, including genes associated with deep diving and navigation. Our study advances understanding of the timing, pattern, and molecular changes associated with the evolution of mammalian lineages adapting to aquatic life.

Hologenome analysis reveals independent evolution to chemosymbiosis by deep-sea bivalves.

BACKGROUND: Bivalves have independently evolved a variety of symbiotic relationships with chemosynthetic bacteria. These relationships range from endo- to extracellular interactions, making them ideal for studies on symbiosis-related evolution. It is still unclear whether there are universal patterns to symbiosis across bivalves. Here, we investigate the hologenome of an extracellular symbiotic thyasirid clam that represents the early stages of symbiosis evolution. RESULTS: We present a hologenome of Conchocele bisecta (Bivalvia: Thyasiridae) collected from deep-sea hydrothermal vents with extracellular symbionts, along with related ultrastructural evidence and expression data. Based on ultrastructural and sequencing evidence, only one dominant Thioglobaceae bacteria was densely aggregated in the large bacterial chambers of C. bisecta, and the bacterial genome shows nutritional complementarity and immune interactions with the host. Overall, gene family expansions may contribute to the symbiosis-related phenotypic variations in different bivalves. For instance, convergent expansions of gaseous substrate transport families in the endosymbiotic bivalves are absent in C. bisecta. Compared to endosymbiotic relatives, the thyasirid genome exhibits large-scale expansion in phagocytosis, which may facilitate symbiont digestion and account for extracellular symbiotic phenotypes. We also reveal that distinct immune system evolution, including expansion in lipopolysaccharide scavenging and contraction of IAP (inhibitor of apoptosis protein), may contribute to the different manners of bacterial virulence resistance in C. bisecta. CONCLUSIONS: Thus, bivalves employ different pathways to adapt to the long-term co-existence with their bacterial symbionts, further highlighting the contribution of stochastic evolution to the independent gain of a symbiotic lifestyle in the lineage.

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.

Reconstruction of the Origin of a Neo-Y Sex Chromosome and Its Evolution in the Spotted Knifejaw, Oplegnathus punctatus.

Sex chromosomes are a peculiar constituent of the genome because the evolutionary forces that fix the primary sex-determining gene cause genic degeneration and accumulation of junk DNA in the heterogametic partner. One of the most spectacular phenomena in sex chromosome evolution is the occurrence of neo-Y chromosomes, which lead to X1X2Y sex-determining systems. Such neo-sex chromosomes are critical for understanding the processes of sex chromosome evolution because they rejuvenate their total gene content. We assembled the male and female genomes at the chromosome level of the spotted knifejaw (Oplegnathus punctatus), which has a cytogenetically recognized neo-Y chromosome. The full assembly and annotation of all three sex chromosomes allowed us to reconstruct their evolutionary history. Contrary to other neo-Y chromosomes, the fusion to X2 is quite ancient, estimated at 48 Ma. Despite its old age and being even older in the X1 homologous region which carries a huge inversion that occurred as early as 55-48 Ma, genetic degeneration of the neo-Y appears to be only moderate. Transcriptomic analysis showed that sex chromosomes harbor 87 genes, which may serve important functions in the testis. The accumulation of such male-beneficial genes, a large inversion on the X1 homologous region and fusion to X2 appear to be the main drivers of neo-Y evolution in the spotted knifejaw. The availability of high-quality assemblies of the neo-Y and both X chromosomes make this fish an ideal model for a better understanding of the variability of sex determination mechanisms and of sex chromosome evolution.

A chromosome-level assembly of the Atlantic herring genome-detection of a supergene and other signals of selection.

The Atlantic herring is a model species for exploring the genetic basis for ecological adaptation, due to its huge population size and extremely low genetic differentiation at selectively neutral loci. However, such studies have so far been hampered because of a highly fragmented genome assembly. Here, we deliver a chromosome-level genome assembly based on a hybrid approach combining a de novo Pacific Biosciences (PacBio) assembly with Hi-C-supported scaffolding. The assembly comprises 26 autosomes with sizes ranging from 12.4 to 33.1 Mb and a total size, in chromosomes, of 726 Mb, which has been corroborated by a high-resolution linkage map. A comparison between the herring genome assembly with other high-quality assemblies from bony fishes revealed few inter-chromosomal but frequent intra-chromosomal rearrangements. The improved assembly facilitates analysis of previously intractable large-scale structural variation, allowing, for example, the detection of a 7.8-Mb inversion on Chromosome 12 underlying ecological adaptation. This supergene shows strong genetic differentiation between populations. The chromosome-based assembly also markedly improves the interpretation of previously detected signals of selection, allowing us to reveal hundreds of independent loci associated with ecological adaptation.

Accurate haplotype-resolved assembly reveals the origin of structural variants for human trios.

MOTIVATION: Achieving a near complete understanding of how the genome of an individual affects the phenotypes of that individual requires deciphering the order of variations along homologous chromosomes in species with diploid genomes. However, true diploid assembly of long-range haplotypes remains challenging. RESULTS: To address this, we have developed Haplotype-resolved Assembly for Synthetic long reads using a Trio-binning strategy, or HAST, which uses parental information to classify reads into maternal or paternal. Once sorted, these reads are used to independently de novo assemble the parent-specific haplotypes. We applied HAST to cobarcoded second-generation sequencing data from an Asian individual, resulting in a haplotype assembly covering 94.7% of the reference genome with a scaffold N50 longer than 11 Mb. The high haplotyping precision (∼99.7%) and recall (∼95.9%) represents a substantial improvement over the commonly used tool for assembling cobarcoded reads (Supernova), and is comparable to a trio-binning-based third generation long-read-based assembly method (TrioCanu) but with a significantly higher single-base accuracy [up to 99.99997% (Q65)]. This makes HAST a superior tool for accurate haplotyping and future haplotype-based studies. AVAILABILITY AND IMPLEMENTATION: The code of the analysis is available at https://github.com/BGI-Qingdao/HAST. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

SLR-superscaffolder: a de novo scaffolding tool for synthetic long reads using a top-to-bottom scheme.

BACKGROUND: Synthetic long reads (SLR) with long-range co-barcoding information are now widely applied in genomics research. Although several tools have been developed for each specific SLR technique, a robust standalone scaffolder with high efficiency is warranted for hybrid genome assembly. RESULTS: In this work, we developed a standalone scaffolding tool, SLR-superscaffolder, to link together contigs in draft assemblies using co-barcoding and paired-end read information. Our top-to-bottom scheme first builds a global scaffold graph based on Jaccard Similarity to determine the order and orientation of contigs, and then locally improves the scaffolds with the aid of paired-end information. We also exploited a screening algorithm to reduce the negative effect of misassembled contigs in the input assembly. We applied SLR-superscaffolder to a human single tube long fragment read sequencing dataset and increased the scaffold NG50 of its corresponding draft assembly 1349 fold. Moreover, benchmarking on different input contigs showed that this approach overall outperformed existing SLR scaffolders, providing longer contiguity and fewer misassemblies, especially for short contigs assembled by next-generation sequencing data. The open-source code of SLR-superscaffolder is available at https://github.com/BGI-Qingdao/SLR-superscaffolder . CONCLUSIONS: SLR-superscaffolder can dramatically improve the contiguity of a draft assembly by integrating a hybrid assembly strategy.

Insights into amphicarpy from the compact genome of the legume Amphicarpaea edgeworthii.

Amphicarpy (seed heteromorphy) is a unique and fascinating reproductive strategy wherein a single plant produces both aerial and subterranean fruits. This strategy is believed to be an adaptation to life under stressful or uncertain environments. Here, we sequenced and de novo assembled a chromosome-level genome assembly of the legume Amphicarpaea edgeworthii Benth. The 299-Mb A. edgeworthii genome encodes 27 899 protein-coding genes and is the most compact sequenced legume genome reported until date. Its reduced genome size may be attributed to the reduced long-terminal repeat retrotransposon content, which stems from the unequal homologous recombination. Gene families related to immunity and stress resistance have been contracted in A. edgeworthii, which is consistent with the notion that the amphicarpic reproductive strategy may be a complementary mechanism for its weak environmental-adaptation ability. We demonstrated the 'ABCE' model for the differentiation of chasmogamous and cleistogamous flowers. In addition, the characteristics of aerial and subterranean seeds in hard-seededness were explored. Thus, we suggest that the A. edgeworthii genome, which is the first of an amphicarpic plant, offers significant insights into its unusual reproductive strategy that is a key resource towards comprehending the evolution of angiosperms.

Runs of homozygosity in killer whale genomes provide a global record of demographic histories.

Runs of homozygosity (ROH) occur when offspring inherit haplotypes that are identical by descent from each parent. Length distributions of ROH are informative about population history; specifically, the probability of inbreeding mediated by mating system and/or population demography. Here, we investigated whether variation in killer whale (Orcinus orca) demographic history is reflected in genome-wide heterozygosity and ROH length distributions, using a global data set of 26 genomes representative of geographic and ecotypic variation in this species, and two F1 admixed individuals with Pacific-Atlantic parentage. We first reconstructed demographic history for each population as changes in effective population size through time using the pairwise sequential Markovian coalescent (PSMC) method. We found a subset of populations declined in effective population size during the Late Pleistocene, while others had more stable demography. Genomes inferred to have undergone ancestral declines in effective population size, were autozygous at hundreds of short ROH (<1 Mb), reflecting high background relatedness due to coalescence of haplotypes deep within the pedigree. In contrast, longer and therefore younger ROH (>1.5 Mb) were found in low latitude populations, and populations of known conservation concern. These include a Scottish killer whale, for which 37.8% of the autosomes were comprised of ROH >1.5 Mb in length. The fate of this population, in which only two adult males have been sighted in the past five years, and zero fecundity over the last two decades, may be inextricably linked to its demographic history and consequential inbreeding depression.

Performance of copy number variants detection based on whole-genome sequencing by DNBSEQ platforms.

BACKGROUND: DNBSEQ™ platforms are new massively parallel sequencing (MPS) platforms that use DNA nanoball technology. Use of data generated from DNBSEQ™ platforms to detect single nucleotide variants (SNVs) and small insertions and deletions (indels) has proven to be quite effective, while the feasibility of copy number variants (CNVs) detection is unclear. RESULTS: Here, we first benchmarked different CNV detection tools based on Illumina whole-genome sequencing (WGS) data of NA12878 and then assessed these tools in CNV detection based on DNBSEQ™ sequencing data from the same sample. When the same tool was used, the CNVs detected based on DNBSEQ™ and Illumina data were similar in quantity, length and distribution, while great differences existed within results from different tools and even based on data from a single platform. We further estimated the CNV detection power based on available CNV benchmarks of NA12878 and found similar precision and sensitivity between the DNBSEQ™ and Illumina platforms. We also found higher precision of CNVs shorter than 1 kbp based on DNBSEQ™ platforms than those based on Illumina platforms by using Pindel, DELLY and LUMPY. We carefully compared these two available benchmarks and found a large proportion of specific CNVs between them. Thus, we constructed a more complete CNV benchmark of NA12878 containing 3512 CNV regions. CONCLUSIONS: We assessed and benchmarked CNV detections based on WGS with DNBSEQ™ platforms and provide guidelines for future studies.

New insights from Opisthorchis felineus genome: update on genomics of the epidemiologically important liver flukes.

BACKGROUND: The three epidemiologically important Opisthorchiidae liver flukes Opisthorchis felineus, O. viverrini, and Clonorchis sinensis, are believed to harbour similar potencies to provoke hepatobiliary diseases in their definitive hosts, although their populations have substantially different ecogeographical aspects including habitat, preferred hosts, population structure. Lack of O. felineus genomic data is an obstacle to the development of comparative molecular biological approaches necessary to obtain new knowledge about the biology of Opisthorchiidae trematodes, to identify essential pathways linked to parasite-host interaction, to predict genes that contribute to liver fluke pathogenesis and for the effective prevention and control of the disease. RESULTS: Here we present the first draft genome assembly of O. felineus and its gene repertoire accompanied by a comparative analysis with that of O. viverrini and Clonorchis sinensis. We observed both noticeably high heterozygosity of the sequenced individual and substantial genetic diversity in a pooled sample. This indicates that potency of O. felineus population for rapid adaptive response to control and preventive measures of opisthorchiasis is higher than in O. viverrini and C. sinensis. We also have found that all three species are characterized by more intensive involvement of trans-splicing in RNA processing compared to other trematodes. CONCLUSION: All revealed peculiarities of structural organization of genomes are of extreme importance for a proper description of genes and their products in these parasitic species. This should be taken into account both in academic and applied research of epidemiologically important liver flukes. Further comparative genomics studies of liver flukes and non-carcinogenic flatworms allow for generation of well-grounded hypotheses on the mechanisms underlying development of cholangiocarcinoma associated with opisthorchiasis and clonorchiasis as well as species-specific mechanisms of these diseases.

Chromosome level comparative analysis of Brassica genomes.

KEY MESSAGE: We provided a chromosome-length assembly of B. nigra and show the comprehensive chromosome-scale variations among Brassica genomes. Chromosome-level assembly of the Brassica species, which include many important crops, is essential for the agricultural and evolutionary studies. While the present B. nigra chromosomes was connected with genetic map of B. juncea, hindering the comparative analysis of the B chromosomes. Here we present a chromosome-length B. nigra assembly constructed with Hi-C connections and its variations on chromosome level compared with other Brassica species. We produced an assembly of 484 Mb annotated with 51,829 genes, of which 393 Mb were anchored onto 8 chromosomes, taking 81.26% of the assembly. Comparison of the B chromosomes shows high concordance of the two B. nigra assemblies and reveals comprehensive variations of the B chromosomes after polyploidization and gene loss in syntenic regions. Chromosome blocks with variations have lower gene density and higher TE content. Furthermore, we compared the chromosomes of the three major Brassica diploids, which showed that most of the variations between B and A/C had completed before A/C divergence and there are more variations on C chromosomes after their divergence. In summary, our work presents a chromosome-length assembly of B. nigra and comprehensive comparative analysis of the Brassica chromosomes, which provides a useful reference for other studies and comprehensive information of Brassica chromosome evolution.

Proteogenomic analysis reveals alternative splicing and translation as part of the abscisic acid response in Arabidopsis seedlings.

In eukaryotes, mechanisms such as alternative splicing (AS) and alternative translation initiation (ATI) contribute to organismal protein diversity. Specifically, splicing factors play crucial roles in responses to environment and development cues; however, the underlying mechanisms are not well investigated in plants. Here, we report the parallel employment of short-read RNA sequencing, single molecule long-read sequencing and proteomic identification to unravel AS isoforms and previously unannotated proteins in response to abscisic acid (ABA) treatment. Combining the data from the two sequencing methods, approximately 83.4% of intron-containing genes were alternatively spliced. Two AS types, which are referred to as alternative first exon (AFE) and alternative last exon (ALE), were more abundant than intron retention (IR); however, by contrast to AS events detected under normal conditions, differentially expressed AS isoforms were more likely to be translated. ABA extensively affects the AS pattern, indicated by the increasing number of non-conventional splicing sites. This work also identified thousands of unannotated peptides and proteins by ATI based on mass spectrometry and a virtual peptide library deduced from both strands of coding regions within the Arabidopsis genome. The results enhance our understanding of AS and alternative translation mechanisms under normal conditions, and in response to ABA treatment.