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.

Large-scale across species transcriptomic analysis identifies genetic selection signatures associated with longevity in mammals.

Lifespan varies significantly among mammals, with more than 100-fold difference between the shortest and longest living species. This natural difference may uncover the evolutionary forces and molecular features that define longevity. To understand the relationship between gene expression variation and longevity, we conducted a comparative transcriptomics analysis of liver, kidney, and brain tissues of 103 mammalian species. We found that few genes exhibit common expression patterns with longevity in the three organs analyzed. However, pathways related to translation fidelity, such as nonsense-mediated decay and eukaryotic translation elongation, correlated with longevity across mammals. Analyses of selection pressure found that selection intensity related to the direction of longevity-correlated genes is inconsistent across organs. Furthermore, expression of methionine restriction-related genes correlated with longevity and was under strong selection in long-lived mammals, suggesting that a common strategy is utilized by natural selection and artificial intervention to control lifespan. Our results indicate that lifespan regulation via gene expression is driven through polygenic and indirect natural selection.

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.

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.

Comparative analysis of the superoxide dismutase gene family in Cetartiodactyla.

Cetacea, whales, dolphins and porpoises form an order of mammals adapted to aquatic life. Their transition to an aquatic habitat resulted in exceptional protection against cellular insults, including oxidative and osmotic stress. Here, we considered the structure and molecular evolution of the superoxide dismutase (SOD) gene family, which encodes essential enzymes in the mammalian antioxidant system, in the superorder Cetartiodactyla. To this end, we juxtaposed cetaceans and their closest extant relatives (order Artiodactyla). We identified 94 genes in 23 species, of which 70 are bona fide intact genes. Although the SOD gene family is conserved in Cetartiodactyla, lineage-specific gene duplications and deletions were observed. Phylogenetic analyses show that the SOD2 subfamily diverged from a clade containing SOD1 and SOD3, suggesting that cytoplasmic, extracellular and mitochondrial SODs have started down independent evolutionary paths. Specific-amino acid changes (e.g. K130N in SOD2) that may enhance ROS elimination were identified in cetaceans. In silico analysis suggests that the core transcription factor repertoire of cetartiodactyl SOD genes may include Sp1, NF-κB, Nrf2 and AHR. Putative transcription factors binding sites responding to hypoxia were (e.g. Suppressor of Hairless; Su(H)) found in the cetacean SOD1 gene. We found significant evidence for positive selection in cetaceans using codon models. Cetaceans with different diving abilities also show divergent evolution of SOD1 and SOD2. Our genome-wide analysis of SOD genes helps clarify their relationship and evolutionary trajectory and identify putative functional changes in cetaceans.

Accelerated evolution and diversifying selection drove the adaptation of cetacean bone microstructure.

BACKGROUND: The transition from land to sea by the ancestor of cetaceans approximately 50 million years ago was an incredible evolutionary event that led to a series of morphological, physiological, and behavioral adaptations. During this transition, bone microstructure evolved from the typical terrestrial form to the specialized structure found in modern cetaceans. While the bone microstructure of mammals has been documented before, investigations of its genetic basis lag behind. The increasing number of cetaceans with whole-genome sequences available may shed light on the mechanism underlying bone microstructure evolution as a result of land to water transitions. RESULTS: Cetacean bone microstructure is consistent with their diverse ecological behaviors. Molecular evolution was assessed by correlating bone microstructure and gene substitution rates in terrestrial and aquatic species, and by detecting genes under positive selection along ancestral branches of cetaceans. We found that: 1) Genes involved in osteoclast function are under accelerated evolution in cetaceans, suggestive of important roles in bone remodeling during the adaptation to an aquatic environment; 2) Genes in the Wnt pathway critical for bone development and homeostasis show evidence of divergent evolution in cetaceans; 3) Several genes encoding bone collagens are under selective pressure in cetaceans. CONCLUSIONS: Our results suggest that evolutionary pressures have shaped the bone microstructure of cetaceans, to facilitate life in diverse aquatic environments.

Adaptive aneuploidy protects against thiol peroxidase deficiency by increasing respiration via key mitochondrial proteins.

Aerobic respiration is a fundamental energy-generating process; however, there is cost associated with living in an oxygen-rich environment, because partially reduced oxygen species can damage cellular components. Organisms evolved enzymes that alleviate this damage and protect the intracellular milieu, most notably thiol peroxidases, which are abundant and conserved enzymes that mediate hydrogen peroxide signaling and act as the first line of defense against oxidants in nearly all living organisms. Deletion of all eight thiol peroxidase genes in yeast (∆8 strain) is not lethal, but results in slow growth and a high mutation rate. Here we characterized mechanisms that allow yeast cells to survive under conditions of thiol peroxidase deficiency. Two independent ∆8 strains increased mitochondrial content, altered mitochondrial distribution, and became dependent on respiration for growth but they were not hypersensitive to H2O2. In addition, both strains independently acquired a second copy of chromosome XI and increased expression of genes encoded by it. Survival of ∆8 cells was dependent on mitochondrial cytochrome-c peroxidase (CCP1) and UTH1, present on chromosome XI. Coexpression of these genes in ∆8 cells led to the elimination of the extra copy of chromosome XI and improved cell growth, whereas deletion of either gene was lethal. Thus, thiol peroxidase deficiency requires dosage compensation of CCP1 and UTH1 via chromosome XI aneuploidy, wherein these proteins support hydroperoxide removal with the reducing equivalents generated by the electron transport chain. To our knowledge, this is the first evidence of adaptive aneuploidy counteracting oxidative stress.

MUC13 overexpression in renal cell carcinoma plays a central role in tumor progression and drug resistance.

Metastatic renal cell carcinoma is a largely incurable disease, and existing treatments targeting angiogenesis and tyrosine kinase receptors are only partially effective. Here we reveal that MUC13, a cell surface mucin glycoprotein, is aberrantly expressed by most renal cell carcinomas, with increasing expression positively correlating with tumor grade. Importantly, we demonstrated that high MUC13 expression was a statistically significant independent predictor of poor survival in two independent cohorts, particularly in stage 1 cancers. In cultured renal cell carcinoma cells MUC13 promoted proliferation and induced the cell cycle regulator, cyclin D1, and inhibited apoptosis by inducing the anti-apoptotic proteins, BCL-xL and survivin. Silencing of MUC13 expression inhibited migration and invasion, and sensitized renal cancer cells to killing by the multi-kinase inhibitors used clinically, sorafenib and sunitinib, and reversed acquired resistance to these drugs. Furthermore, we demonstrated that MUC13 promotion of renal cancer cell growth and survival is mediated by activation of nuclear factor κB, a transcription factor known to regulate the expression of genes that play key roles in the development and progression of cancer. These results show that MUC13 has potential as a prognostic marker for aggressive early stage renal cell cancer and is a plausible target to sensitize these tumors to therapy.

Fusion transcript loci share many genomic features with non-fusion loci.

BACKGROUND: Fusion transcripts are found in many tissues and have the potential to create novel functional products. Here, we investigate the genomic sequences around fusion junctions to better understand the transcriptional mechanisms mediating fusion transcription/splicing. We analyzed data from prostate (cancer) cells as previous studies have shown extensively that these cells readily undergo fusion transcription. RESULTS: We used the FusionMap program to identify high-confidence fusion transcripts from RNAseq data. The RNAseq datasets were from our (N = 8) and other (N = 14) clinical prostate tumors with adjacent non-cancer cells, and from the LNCaP prostate cancer cell line that were mock-, androgen- (DHT), and anti-androgen- (bicalutamide, enzalutamide) treated. In total, 185 fusion transcripts were identified from all RNAseq datasets. The majority (76%) of these fusion transcripts were 'read-through chimeras' derived from adjacent genes in the genome. Characterization of sequences at fusion loci were carried out using a combination of the FusionMap program, custom Perl scripts, and the RNAfold program. Our computational analysis indicated that most fusion junctions (76%) use the consensus GT-AG intron donor-acceptor splice site, and most fusion transcripts (85%) maintained the open reading frame. We assessed whether parental genes of fusion transcripts have the potential to form complementary base pairing between parental genes which might bring them into physical proximity. Our computational analysis of sequences flanking fusion junctions at parental loci indicate that these loci have a similar propensity as non-fusion loci to hybridize. The abundance of repetitive sequences at fusion and non-fusion loci was also investigated given that SINE repeats are involved in aberrant gene transcription. We found few instances of repetitive sequences at both fusion and non-fusion junctions. Finally, RT-qPCR was performed on RNA from both clinical prostate tumors and adjacent non-cancer cells (N = 7), and LNCaP cells treated as above to validate the expression of seven fusion transcripts and their respective parental genes. We reveal that fusion transcript expression is similar to the expression of parental genes. CONCLUSIONS: Fusion transcripts maintain the open reading frame, and likely use the same transcriptional machinery as non-fusion transcripts as they share many genomic features at splice/fusion junctions.

Comparative analysis reveals loss of the appetite-regulating peptide hormone ghrelin in falcons.

Ghrelin and leptin are key peripherally secreted appetite-regulating hormones in vertebrates. Here we consider the ghrelin gene (GHRL) of birds (class Aves), where it has been reported that ghrelin inhibits rather than augments feeding. Thirty-one bird species were compared, revealing that most species harbour a functional copy of GHRL and the coding region for its derived peptides ghrelin and obestatin. We provide evidence for loss of GHRL in saker and peregrine falcons, and this is likely to result from the insertion of an ERVK retrotransposon in intron 0. We hypothesise that the loss of anorexigenic ghrelin is a predatory adaptation that results in increased food-seeking behaviour and feeding in falcons.

Associations between ghrelin and ghrelin receptor polymorphisms and cancer in Caucasian populations: a meta-analysis.

BACKGROUND: There is growing evidence that the ghrelin axis, including ghrelin (GHRL) and its receptor, the growth hormone secretagogue receptor (GHSR), play a role in cancer progression. Ghrelin gene and ghrelin receptor gene polymorphisms have been reported to have a range of effects in cancer, from increased risk, to protection from cancer, or having no association. In this study we aimed to clarify the role of ghrelin and ghrelin receptor polymorphisms in cancer by performing a meta-analysis of published case-control studies. RESULTS: In the overall analysis, homozygous and recessive associations indicated that the minor alleles of rs696217 and rs2075356 GHRL polymorphisms conferred reduced cancer risk (odds ratio [OR] 0.61-0.78). The risk was unchanged for breast cancer patients when analysed separately (OR 0.73-0.83). In contrast, the rs4684677 GHRL and the rs572169 GHSR polymorphisms conferred increased breast cancer risk (OR 1.97-1.98, p = 0.08 and OR 1.42-1.43, p = 0.08, respectively). All dominant and co-dominant effects showed null effects (OR 0.96-1.05), except for the rs572169 co-dominant effect, with borderline increased risk (OR 1.08, p = 0.05). CONCLUSIONS: This study suggests that the rs696217 and rs2075356 ghrelin gene (GHRL) polymorphisms may protect carriers against breast cancer, and the rs4684677 GHRL and rs572169 GHSR polymorphisms may increase the risk among carriers. In addition, larger studies are required to confirm these findings.

Cloning and tissue distribution of novel splice variants of the ovine ghrelin gene.

BACKGROUND: The ghrelin axis is involved in the regulation of metabolism, energy balance, and the immune, cardiovascular and reproductive systems. The manipulation of this axis has potential for improving economically valuable traits in production animals, and polymorphisms in the ghrelin (GHRL) and ghrelin receptor (GHSR) genes have been associated with growth and carcass traits. Here we investigate the structure and expression of the ghrelin gene (GHRL) in sheep, Ovis aries. RESULTS: We identify two ghrelin mRNA isoforms, which we have designated Δex2 preproghrelin and Δex2,3 preproghrelin. Expression of Δex2,3 preproghrelin is likely to be restricted to ruminants, and would encode truncated ghrelin and a novel C-terminal peptide. Both Δex2 preproghrelin and canonical preproghrelin mRNA isoforms were expressed in a range of tissues. Expression of the Δex2,3 preproghrelin isoform, however, was restricted to white blood cells (WBC; where the wild-type preproghrelin isoform is not co-expressed), and gastrointestinal tissues. Expression of Δex2 preproghrelin and Δex2,3 preproghrelin mRNA was elevated in white blood cells in response to parasitic worm (helminth) infection in genetically susceptible sheep, but not in resistant sheep. CONCLUSIONS: The restricted expression of the novel preproghrelin variants and their distinct WBC expression pattern during parasite infection may indicate a novel link between the ghrelin axis and metabolic and immune function in ruminants.

A chromosome-level genome of electric catfish (Malapterurus electricus) provided new insights into order Siluriformes evolution.

The electric catfish (Malapterurus electricus), belonging to the family Malapteruridae, order Siluriformes (Actinopterygii: Ostariophysi), is one of the six branches that has independently evolved electrical organs. We assembled a 796.75 Mb M. electricus genome and anchored 88.72% sequences into 28 chromosomes. Gene family analysis revealed 295 expanded gene families that were enriched on functions related to glutamate receptors. Convergent evolutionary analyses of electric organs among different lineage of electric fishes further revealed that the coding gene of rho guanine nucleotide exchange factor 4-like (arhgef4), which is associated with G-protein coupled receptor (GPCR) signaling pathway, underwent adaptive parallel evolution. Gene identification suggests visual degradation in catfishes, and an important role for taste in environmental adaptation. Our findings fill in the genomic data for a branch of electric fish and provide a relevant genetic basis for the adaptive evolution of Siluriformes. Supplementary Information: The online version contains supplementary material available at 10.1007/s42995-023-00197-8.

Genomics of post-bottleneck recovery in the northern elephant seal.

Populations and species are threatened by human pressure, but their fate is variable. Some depleted populations, such as that of the northern elephant seal (Mirounga angustirostris), recover rapidly even when the surviving population was small. The northern elephant seal was hunted extensively and taken by collectors between the early 1800s and 1892, suffering an extreme population bottleneck as a consequence. Recovery was rapid and now there are over 200,000 individuals. We sequenced 260 modern and 8 historical northern elephant seal nuclear genomes to assess the impact of the population bottleneck on individual northern elephant seals and to better understand their recovery. Here we show that inbreeding, an increase in the frequency of alleles compromised by lost function, and allele frequency distortion, reduced the fitness of breeding males and females, as well as the performance of adult females on foraging migrations. We provide a detailed investigation of the impact of a severe bottleneck on fitness at the genomic level and report on the role of specific gene systems.

Silencing of ghrelin receptor expression inhibits endometrial cancer cell growth in vitro and in vivo.

Ghrelin is a 28-amino acid peptide hormone produced predominantly in the stomach but also in a range of normal cell types and tumors, where it has endocrine, paracrine, and autocrine roles. Previously, we have demonstrated that ghrelin has proliferative and antiapoptotic effects in endometrial cancer cell lines, suggesting a potential role in promoting tumor growth. In the present study, we investigated the effect of ghrelin receptor, GHSR, and gene silencing in vitro and in vivo and characterized ghrelin and GHSR1a protein expression in human endometrial tumors. GHSR gene silencing was achieved in the Ishikawa and KLE endometrial cancer cell lines, using a lentiviral short-hairpin RNA targeting GHSR. The effects of GHSR1a knockdown were further analyzed in vivo using the Ishikawa cell line in a NOD/SCID xenograft model. Cell proliferation was reduced in cultured GHSR1a knockdown Ishikawa and KLE cells compared with scrambled controls in the absence of exogenously applied ghrelin and in response to exogenous ghrelin (1,000 nM). The tumor volumes were reduced significantly in GHSR1a knockdown Ishikawa mouse xenograft tumors compared with scrambled control tumours. Using immunohistochemistry, we demonstrated that ghrelin and GHSR1a are expressed in benign and cancerous glands in human endometrial tissue specimens, although there was no correlation between the intensity of staining and cancer grade. These data indicate that downregulation of GHSR expression significantly inhibits endometrial cancer cell line and mouse xenograft tumour growth. This is the first preclinical evidence that downregulation of GHSR may be therapeutic in endometrial cancer.