Large-scale snake genome analyses provide insights into vertebrate development.

Snakes are a remarkable squamate lineage with unique morphological adaptations, especially those related to the evolution of vertebrate skeletons, organs, and sensory systems. To clarify the genetic underpinnings of snake phenotypes, we assembled and analyzed 14 de novo genomes from 12 snake families. We also investigated the genetic basis of the morphological characteristics of snakes using functional experiments. We identified genes, regulatory elements, and structural variations that have potentially contributed to the evolution of limb loss, an elongated body plan, asymmetrical lungs, sensory systems, and digestive adaptations in snakes. We identified some of the genes and regulatory elements that might have shaped the evolution of vision, the skeletal system and diet in blind snakes, and thermoreception in infrared-sensitive snakes. Our study provides insights into the evolution and development of snakes and vertebrates.

Genomic adaptations for arboreal locomotion in Asian flying treefrogs.

SignificanceTo adapt to arboreal lifestyles, treefrogs have evolved a suite of complex traits that support vertical movement and gliding, thus presenting a unique case for studying the genetic basis for traits causally linked to vertical niche expansion. Here, based on two de novo-assembled Asian treefrog genomes, we determined that genes involved in limb development and keratin cytoskeleton likely played a role in the evolution of their climbing systems. Behavioral and morphological evaluation and time-ordered gene coexpression network analysis revealed the developmental patterns and regulatory pathways of the webbed feet used for gliding in Rhacophorus kio.

The Genome of Shaw's Sea Snake (Hydrophis curtus) Reveals Secondary Adaptation to Its Marine Environment.

The transition of terrestrial snakes to marine life ∼10 Ma is ideal for exploring adaptive evolution. Sea snakes possess phenotype specializations including laterally compressed bodies, paddle-shaped tails, valvular nostrils, cutaneous respiration, elongated lungs, and salt glands, yet, knowledge on the genetic underpinnings of the transition remains limited. Herein, we report the first genome of Shaw's sea snake (Hydrophis curtus) and use it to investigate sea snake secondary marine adaptation. A hybrid assembly strategy obtains a high-quality genome. Gene family analyses date a pulsed coding-gene expansion to ∼20 Ma, and these genes associate strongly with adaptations to marine environments. Analyses of selection pressure and convergent evolution discover the rapid evolution of protein-coding genes, and some convergent features. Additionally, 108 conserved noncoding elements appear to have evolved quickly, and these may underpin the phenotypic changes. Transposon elements may contribute to adaptive specializations by inserting into genomic regions around functionally related coding genes. The integration of genomic and transcriptomic analyses indicates independent origins and different components in sea snake and terrestrial snake venom; the venom gland of the sea snake harbors the highest PLA2 (17.23%) expression in selected elapids and these genes may organize tandemly in the genome. These analyses provide insights into the genetic mechanisms that underlay the secondary adaptation to marine and venom production of this sea snake.

RNA-binding protein Hfq plays a vital role in cellulose decomposition throughout affecting cellulase gene expression.

OBJECTIVE: To study the function of the RNA-binding protein Hfq in Bacillus subtilis cellulose decomposition. RESULTS: In the medium with sodium carboxymethylcellulose (Na-CMC) as the sole carbon source, the knockout of Hfq resulted in a 38.0% ± 2.1% and 76.6% ± 7.1% decrease in cellulose hydrolysis ability and cellulase activity, respectively. The results of real-time quantitative PCR revealed that several cellulase genes (eglS, bglA, and bglC) were significantly downregulated in the Hfq knockout strain. The isogenic Δhfq complemented strain recovered the cellulose hydrolysis ability, cellulase activity, and expression level of cellulase genes. In addition, the survival of Hfq mutant in stationary phase was significantly affected. CONCLUSION: RNA-binding protein Hfq is involved in the regulation of cellulose hydrolysis ability, cellulase activity, cellulase gene expression, and stationary phase survival.

Comparative genomic investigation of high-elevation adaptation in ectothermic snakes.

Several previous genomic studies have focused on adaptation to high elevations, but these investigations have been largely limited to endotherms. Snakes of the genus Thermophis are endemic to the Tibetan plateau and therefore present an opportunity to study high-elevation adaptations in ectotherms. Here, we report the de novo assembly of the genome of a Tibetan hot-spring snake (Thermophis baileyi) and then compare its genome to the genomes of the other two species of Thermophis, as well as to the genomes of two related species of snakes that occur at lower elevations. We identify 308 putative genes that appear to be under positive selection in Thermophis We also identified genes with shared amino acid replacements in the high-elevation hot-spring snakes compared with snakes and lizards that live at low elevations, including the genes for proteins involved in DNA damage repair (FEN1) and response to hypoxia (EPAS1). Functional assays of the FEN1 alleles reveal that the Thermophis allele is more stable under UV radiation than is the ancestral allele found in low-elevation lizards and snakes. Functional assays of EPAS1 alleles suggest that the Thermophis protein has lower transactivation activity than the low-elevation forms. Our analysis identifies some convergent genetic mechanisms in high-elevation adaptation between endotherms (based on studies of mammals) and ectotherms (based on our studies of Thermophis).

Adaptive evolution to a high purine and fat diet of carnivorans revealed by gut microbiomes and host genomes.

Carnivorous members of the Carnivora reside at the apex of food chains and consume meat-only diets, rich in purine, fats and protein. Here, we aimed to identify potential adaptive evolutionary signatures compatible with high purine and fat metabolism based on analysis of host genomes and symbiotic gut microbial metagenomes. We found that the gut microbiomes of carnivorous Carnivora (e.g., Felidae, Canidae) clustered in the same clade, and other clades comprised omnivorous and herbivorous Carnivora (e.g., badgers, bears and pandas). The relative proportions of genes encoding enzymes involved in uric acid degradation were higher in the gut microbiomes of meat-eating carnivorans than plant-eating species. Adaptive amino acid substitutions in two enzymes, carnitine O-palmitoyltransferase 1 (CPT1A) and lipase F (LIPF), which play a role in fat digestion, were identified in Felidae-Candidae species. Carnivorous carnivorans appear to endure diets high in purines and fats via gut microbiomic and genomic adaptations.

The draft genome of Tibetan hulless barley reveals adaptive patterns to the high stressful Tibetan Plateau.

The Tibetan hulless barley (Hordeum vulgare L. var. nudum), also called "Qingke" in Chinese and "Ne" in Tibetan, is the staple food for Tibetans and an important livestock feed in the Tibetan Plateau. The diploid nature and adaptation to diverse environments of the highland give it unique resources for genetic research and crop improvement. Here we produced a 3.89-Gb draft assembly of Tibetan hulless barley with 36,151 predicted protein-coding genes. Comparative analyses revealed the divergence times and synteny between barley and other representative Poaceae genomes. The expansion of the gene family related to stress responses was found in Tibetan hulless barley. Resequencing of 10 barley accessions uncovered high levels of genetic variation in Tibetan wild barley and genetic divergence between Tibetan and non-Tibetan barley genomes. Selective sweep analyses demonstrate adaptive correlations of genes under selection with extensive environmental variables. Our results not only construct a genomic framework for crop improvement but also provide evolutionary insights of highland adaptation of Tibetan hulless barley.

The genome of the pear (Pyrus bretschneideri Rehd.).

The draft genome of the pear (Pyrus bretschneideri) using a combination of BAC-by-BAC and next-generation sequencing is reported. A 512.0-Mb sequence corresponding to 97.1% of the estimated genome size of this highly heterozygous species is assembled with 194× coverage. High-density genetic maps comprising 2005 SNP markers anchored 75.5% of the sequence to all 17 chromosomes. The pear genome encodes 42,812 protein-coding genes, and of these, ~28.5% encode multiple isoforms. Repetitive sequences of 271.9 Mb in length, accounting for 53.1% of the pear genome, are identified. Simulation of eudicots to the ancestor of Rosaceae has reconstructed nine ancestral chromosomes. Pear and apple diverged from each other ~5.4-21.5 million years ago, and a recent whole-genome duplication (WGD) event must have occurred 30-45 MYA prior to their divergence, but following divergence from strawberry. When compared with the apple genome sequence, size differences between the apple and pear genomes are confirmed mainly due to the presence of repetitive sequences predominantly contributed by transposable elements (TEs), while genic regions are similar in both species. Genes critical for self-incompatibility, lignified stone cells (a unique feature of pear fruit), sorbitol metabolism, and volatile compounds of fruit have also been identified. Multiple candidate SFB genes appear as tandem repeats in the S-locus region of pear; while lignin synthesis-related gene family expansion and highly expressed gene families of HCT, C3'H, and CCOMT contribute to high accumulation of both G-lignin and S-lignin. Moreover, alpha-linolenic acid metabolism is a key pathway for aroma in pear fruit.

The Complete Chloroplast Genome Sequences of the Medicinal Plant Pogostemon cablin.

Pogostemon cablin, the natural source of patchouli alcohol, is an important herb in the Lamiaceae family. Here, we present the entire chloroplast genome of P. cablin. This genome, with 38.24% GC content, is 152,460 bp in length. The genome presents a typical quadripartite structure with two inverted repeats (each 25,417 bp in length), separated by one small and one large single-copy region (17,652 and 83,974 bp in length, respectively). The chloroplast genome encodes 127 genes, of which 107 genes are single-copy, including 79 protein-coding genes, four rRNA genes, and 24 tRNA genes. The genome structure, GC content, and codon usage of this chloroplast genome are similar to those of other species in the family, except that it encodes less protein-coding genes and tRNA genes. Phylogenetic analysis reveals that P. cablin diverged from the Scutellarioideae clade about 29.45 million years ago (Mya). Furthermore, most of the simple sequence repeats (SSRs) are short polyadenine or polythymine repeats that contribute to high AT content in the chloroplast genome. Complete sequences and annotation of P. cablin chloroplast genome will facilitate phylogenic, population and genetic engineering research investigations involving this particular species.

The WRKY transcription factor genes in eggplant (Solanum melongena L.) and Turkey Berry (Solanum torvum Sw.).

WRKY transcription factors, which play critical roles in stress responses, have not been characterized in eggplant or its wild relative, turkey berry. The recent availability of RNA-sequencing data provides the opportunity to examine WRKY genes from a global perspective. We identified 50 and 62 WRKY genes in eggplant (SmelWRKYs) and turkey berry (StorWRKYs), respectively, all of which could be classified into three groups (I-III) based on the WRKY protein structure. The SmelWRKYs and StorWRKYs contain ~76% and ~95% of the number of WRKYs found in other sequenced asterid species, respectively. Positive selection analysis revealed that different selection constraints could have affected the evolution of these groups. Positively-selected sites were found in Groups IIc and III. Branch-specific selection pressure analysis indicated that most WRKY domains from SmelWRKYs and StorWRKYs are conserved and have evolved at low rates since their divergence. Comparison to homologous WRKY genes in Arabidopsis revealed several potential pathogen resistance-related SmelWRKYs and StorWRKYs, providing possible candidate genetic resources for improving stress tolerance in eggplant and probably other Solanaceae plants. To our knowledge, this is the first report of a genome-wide analyses of the SmelWRKYs and StorWRKYs.

Comparative transcriptome analysis of eggplant (Solanum melongena L.) and turkey berry (Solanum torvum Sw.): phylogenomics and disease resistance analysis.

BACKGROUND: Eggplant (Solanum melongena L.) and turkey berry (S. torvum Sw.), a wild ally of eggplant with promising multi-disease resistance traits, are of great economic, medicinal and genetic importance, but genomic resources for these species are lacking. In the present study, we sequenced the transcriptomes of eggplant and turkey berry to accelerate research on these two non-model species. RESULTS: We built comprehensive, high-quality de novo transcriptome assemblies of the two Leptostemonum clade Solanum species from short-read RNA-Sequencing data. We obtained 34,174 unigenes for eggplant and 38,185 unigenes for turkey berry. Functional annotations based on sequence similarity to known plant datasets revealed a distribution of functional categories for both species very similar to that of tomato. Comparison of eggplant, turkey berry and another 11 plant proteomes resulted in 276 high-confidence single-copy orthologous groups, reasonable phylogenetic tree inferences and reliable divergence time estimations. From these data, it appears that eggplant and its wild Leptostemonum clade relative turkey berry split from each other in the late Miocene, ~6.66 million years ago, and that Leptostemonum split from the Potatoe clade in the middle Miocene, ~15.75 million years ago. Furthermore, 621 and 815 plant resistance genes were identified in eggplant and turkey berry respectively, indicating the variation of disease resistance genes between them. CONCLUSIONS: This study provides a comprehensive transcriptome resource for two Leptostemonum clade Solanum species and insight into their evolutionary history and biological characteristics. These resources establish a foundation for further investigations of eggplant biology and for agricultural improvement of this important vegetable. More generally, we show that RNA-Seq is a fast, reliable and cost-effective method for assessing genome evolution in non-model species.

Polyploidization of Indotyphlops braminus: evidence from isoform-sequencing.

BACKGROUND: Indotyphlops braminus, the only known triploid parthenogenetic snake, is a compelling species for revealing the mechanism of polyploid emergence in vertebrates. METHODS: In this study, we applied PacBio isoform sequencing technology to generate the first full-length transcriptome of I. braminus, aiming to improve the understanding of the molecular characteristics of this species. RESULTS: A total of 51,849 nonredundant full-length transcript assemblies (with an N50 length of 2980 bp) from I. braminus were generated and fully annotated using various gene function databases. Our analysis provides preliminary evidence supporting a recent genome duplication event in I. braminus. Phylogenetic analysis indicated that the divergence of I. braminus subgenomes occurred approximately 11.5 ~ 15 million years ago (Mya). The full-length transcript resource generated as part of this research will facilitate transcriptome analysis and genomic evolution studies in the future.

Chromosome-level assembly of the synthetic hexaploid wheat-derived cultivar Chuanmai 104.

Synthetic hexaploid wheats (SHWs) are effective genetic resources for transferring agronomically important genes from wild relatives to common wheat (Triticum aestivum L.). Dozens of reference-quality pseudomolecule assemblies of hexaploid wheat have been generated, but none is reported for SHW-derived cultivars. Here, we generated a chromosome-scale assembly for the SHW-derived cultivar 'Chuanmai 104' based on PacBio HiFi reads and chromosome conformation capture sequencing. The total assembly size was 14.81 Gb with a contig N50 length of 58.25 Mb. A BUSCO analysis yielded a completeness score of 99.30%. In total, repetitive elements comprised 81.36% of the genome and 122,554 high-confidence protein-coding gene models were predicted. In summary, the first chromosome-level assembly for a SHW-derived cultivar presents a promising outlook for the study and utilization of SHWs in wheat improvement, which is essential to meet the global food demand.

The Spatio-Temporal Expression Profiles of Silkworm Pseudogenes Provide Valuable Insights into Their Biological Roles.

Background: Pseudogenes are sequences that have lost the ability to transcribe RNA molecules or encode truncated but possibly functional proteins. While they were once considered to be meaningless remnants of evolution, recent researches have shown that pseudogenes play important roles in various biological processes. However, the studies of pseudogenes in the silkworm, an important model organism, are limited and have focused on single or only a few specific genes. Objective: To fill these gaps, we present a systematic genome-wide studies of pseudogenes in the silkworm. Methods: We identified the pseudogenes in the silkworm using the silkworm genome assemblies, transcriptome, protein sequences from silkworm and its related species. Then we used transcriptome datasets from 832 RNA-seq analyses to construct spatio-temporal expression profiles for these pseudogenes. Additionally, we identified tissue-specifically expressed and differentially expressed pseudogenes to further understand their characteristics. Finally, the functional roles of pseudogenes as lncRNAs were systematically analyzed. Results: We identified a total of 4410 pseudogenes, which were grouped into 4 groups, including duplications (DUPs), unitary pseudogenes (Unitary), processed pseudogenes (retropseudogenes, RETs), and fragments (FRAGs). The most of pseudogenes in the domestic silkworm were generated before the divergence of wild and domestic silkworm, however, the domestication may also involve in the accumulation of pseudogenes. These pseudogenes were clearly divided into 2 cluster, a highly expressed and a lowly expressed, and the posterior silk gland was the tissue with the most tissue-specific pseudogenes (199), implying these pseudogenes may be involved in the development and function of silkgland. We identified 3299 lncRNAs in these pseudogenes, and the target genes of these lncRNAs in silkworm pseudogenes were enriched in the egg formation and olfactory function. Conclusions: This study replenishes the genome annotations for silkworm, provide valuable insights into the biological roles of pseudogenes. It will also contribute to our understanding of the complex gene regulatory networks in the silkworm and will potentially have implications for other organisms as well.

Systematic identification of smORFs in domestic silkworm (Bombyx mori).

The silkworm (Bombyx mori) is not only an excellent model species, but also an important agricultural economic insect. Taking it as the research object, its advantages of low maintenance cost and no biohazard risks are considered. Small open reading frames (smORFs) are an important class of genomic elements that can produce bioactive peptides. However, the smORFs in silkworm had been poorly identified and studied. To further study the smORFs in silkworm, systematic genome-wide identification is essential. Here, we identified and analyzed smORFs in the silkworm using comprehensive methods. Our results showed that at least 738 highly reliable smORFs were found in B. mori and that 34,401 possible smORFs were partially supported. We also identified some differentially expressed and tissue-specific-expressed smORFs, which may be closely related to the characteristics and functions of the tissues. This article provides a basis for subsequent research on smORFs in silkworm, and also hopes to provide a reference point for future research methods for smORFs in other species.

Reference genome assemblies reveal the origin and evolution of allohexaploid oat.

Common oat (Avena sativa) is an important cereal crop serving as a valuable source of forage and human food. Although reference genomes of many important crops have been generated, such work in oat has lagged behind, primarily owing to its large, repeat-rich polyploid genome. Here, using Oxford Nanopore ultralong sequencing and Hi-C technologies, we have generated a reference-quality genome assembly of hulless common oat, comprising 21 pseudomolecules with a total length of 10.76 Gb and contig N50 of 75.27 Mb. We also produced genome assemblies for diploid and tetraploid Avena ancestors, which enabled the identification of oat subgenomes and provided insights into oat chromosomal evolution. The origin of hexaploid oat is inferred from whole-genome sequencing, chloroplast genomes and transcriptome assemblies of different Avena species. These findings and the high-quality reference genomes presented here will facilitate the full use of crop genetic resources to accelerate oat improvement.

Proteomics revisits the cancer metabolome.

Many cancer cells utilize aerobic glycolysis (also known as the 'Warburg effect'), instead of mitochondrial oxidative phosphorylation, to generate the energy necessary for diverse cellular processes. In tumor cells, mitochondria play more important roles in anabolism, for instance, de novo lipid biosynthesis and glutamine-dependent anaplerosis to fuel robust cell growth and proliferation. Proteomic analysis of tumor-related alterations of metabolism-associated proteins clearly indicates that such metabolic reprogramming contributes to cancer cell survival and cancer progression. Moreover, proteomics-based systems biology provides a powerful tool to re-evaluate the metabolic phenotype and regulatory mechanism associated with malignant cancer cells, and underscores their implications for cancer diagnosis and therapy. This article will address recent exciting advances in the understanding of cancer cell metabolism using proteomics-based systems biology approaches.

Cytochrome P450 monooxygenase genes in the wild silkworm, Bombyx mandarina.

Wild (Bombyx mandarina) and domestic silkworms (B. mori) are good models for investigating insect domestication, as 5000 years of artificial breeding and selection have resulted in significant differences between B. mandarina and B. mori. In this study, we improved the genome assemblies to the chromosome level and updated the protein-coding gene annotations for B. mandarina. Based on this updated genome, we identified 68 cytochrome P450 genes in B. mandarina. The cytochrome P450 repository in B. mandarina is smaller than in B. mori. Certain currently unknown key genes, rather than gene number, are critical for insecticide resistance in B. mandarina, which shows greater resistance to insecticides than B. mori. Based on the physical maps of B. mandarina, we located 66 cytochrome P450s on 18 different chromosomes, and 27 of the cytochrome P450 genes were concentrated into seven clusters. KEGG enrichment analysis of the P450 genes revealed the involvement of cytochrome P450 genes in hormone biosynthesis. Analyses of the silk gland transcriptome identified candidate cytochrome P450 genes (CYP306A) involved in ecdysteroidogenesis and insecticide metabolism in B. mandarina.

Genome-wide signatures of mammalian skin covering evolution.

Animal body coverings provide protection and allow for adaptation to environmental pressures such as heat, ultraviolet radiation, water loss, and mechanical forces. Here, using a comparative genomics analysis of 39 mammal species spanning three skin covering types (hairless, scaly and spiny), we found some genes (e.g., UVRAG, POLH, and XPC) involved in skin inflammation, skin innate immunity, and ultraviolet radiation damage repair were under selection in hairless ocean mammals (e.g., whales and manatees). These signatures might be associated with a high risk of skin diseases from pathogens and ultraviolet radiation. Moreover, the genomes from three spiny mammal species shared convergent genomic regions (EPHB2, EPHA4, and NIN) and unique positively selected genes (FZD6, INVS, and CDC42) involved in skin cell polarity, which might be related to the development of spines. In scaly mammals, the shared convergent genomic regions (e.g., FREM2) were associated with the integrity of the skin epithelium and epidermal adhesion. This study identifies potential convergent genomic features among distantly related mammals with the same skin covering type.