Comparative genomics and genome-wide SNPs of endangered Eld's deer provide breeder selection for inbreeding avoidance.

Eld's deer, a conserved wildlife species of Thailand, is facing inbreeding depression, particularly in the captive Siamese Eld's deer (SED) subspecies. In this study, we constructed genomes of a male SED and a male Burmese Eld's deer (BED), and used genome-wide single nucleotide polymorphisms to evaluate the genetic purity and the inbreeding status of 35 SED and 49 BED with limited pedigree information. The results show that these subspecies diverged approximately 1.26 million years ago. All SED were found to be purebred. A low proportion of admixed SED genetic material was observed in some BED individuals. Six potential breeders from male SED with no genetic relation to any female SED and three purebred male BED with no relation to more than 10 purebred female BED were identified. This study provides valuable insights about Eld's deer populations and appropriate breeder selection in efforts to repopulate this endangered species while avoiding inbreeding.

Resequencing of durian genomes reveals large genetic variations among different cultivars.

Durian (Durio zibethinus), which yields the fruit known as the "King of Fruits," is an important economic crop in Southeast Asia. Several durian cultivars have been developed in this region. In this study, we resequenced the genomes of three popular durian cultivars in Thailand, including Kradumthong (KD), Monthong (MT), and Puangmanee (PM) to investigate genetic diversities of cultivated durians. KD, MT, and PM genome assemblies were 832.7, 762.6, and 821.6 Mb, and their annotations covered 95.7, 92.4, and 92.7% of the embryophyta core proteins, respectively. We constructed the draft durian pangenome and analyzed comparative genomes with related species in Malvales. Long terminal repeat (LTR) sequences and protein families in durian genomes had slower evolution rates than that in cotton genomes. However, protein families with transcriptional regulation function and protein phosphorylation function involved in abiotic and biotic stress responses appeared to evolve faster in durians. The analyses of phylogenetic relationships, copy number variations (CNVs), and presence/absence variations (PAVs) suggested that the genome evolution of Thai durians was different from that of the Malaysian durian, Musang King (MK). Among the three newly sequenced genomes, the PAV and CNV profiles of disease resistance genes and the expressions of methylesterase inhibitor domain containing genes involved in flowering and fruit maturation in MT were different from those in KD and PM. These genome assemblies and their analyses provide valuable resources to gain a better understanding of the genetic diversity of cultivated durians, which may be useful for the future development of new durian cultivars.

The complete chloroplast genome of Sonneratia griffithii Kurz (Lythraceae).

Sonneratia griffithii Kurz is a critically endangered mangrove species that can be found along the western coast of Thailand. In this study, we reported the complete chloroplast genome of S. griffithii. The chloroplast genome is 152,730 bp, consisting of one large single-copy (LSC) region, one small single-copy (SSC) region and a pair of inverted repeats (IRs). The LSC, SSC, and IR lengths are 87,226, 17,764, and 23,870 bp, respectively. The genome contains 113 unique genes, including 79 protein-coding, 30 tRNA, and 4 rRNA genes. The GC content of the chloroplast genome is 37.31%. The phylogenetic analysis based on 76 protein-coding genes showed a monophyletic group of S. griffithii and other Sonneratia species.

Genome assembly of the Pendlebury's roundleaf bat, Hipposideros pendleburyi, revealed the expansion of Tc1/Mariner DNA transposons in Rhinolophoidea.

Bats (Chiroptera) constitute the second largest order of mammals and have several distinctive features, such as true self-powered flight and strong immunity. The Pendlebury's roundleaf bat, Hipposideros pendleburyi, is endemic to Thailand and listed as a vulnerable species. We employed the 10× Genomics linked-read technology to obtain a genome assembly of H. pendleburyi. The assembly size was 2.17 Gb with a scaffold N50 length of 15,398,518 bases. Our phylogenetic analysis placed H. pendleburyi within the rhinolophoid clade of the suborder Yinpterochiroptera. A synteny analysis showed that H. pendleburyi shared conserved chromosome segments (up to 105 Mb) with Rhinolophus ferrumequinum and Phyllostomus discolor albeit having different chromosome numbers and belonging different families. We found positive selection signals in genes involved in inflammation, spermatogenesis and Wnt signalling. The analyses of transposable elements suggested the contraction of short interspersed nuclear elements (SINEs) and the accumulation of young mariner DNA transposons in the analysed hipposiderids. Distinct mariners were likely horizontally transferred to hipposiderid genomes over the evolution of this family. The lineage-specific profiles of SINEs and mariners might involve in the evolution of hipposiderids and be associated with the phylogenetic separations of these bats from other bat families.

Transcriptome sequencing revealed the influence of blue light on the expression levels of light-stress response genes in Centella asiatica.

Centella asiatica is rich in medical and cosmetic properties. While physiological responses of C. asiatica to light have been widely reported, the knowledge of the effects of light on its gene expression is sparse. In this study, we used RNA sequencing (RNA-seq) to investigate the expression of the C. asiatica genes in response to monochromatic red and blue light. Most of the differentially expressed genes (DEGs) under blue light were up-regulated but those under red light were down-regulated. The DEGs encoded for CRY-DASH and UVR3 were among up-regulated genes that play significant roles in responses under blue light. The DEGs involved in the response to photosystem II photodamages and in the biosynthesis of photoprotective xanthophylls were also up-regulated. The expression of flavonoid biosynthetic DEGs under blue light was up-regulated but that under red light was down-regulated. Correspondingly, total flavonoid content under blue light was higher than that under red light. The ABI5, MYB4, and HYH transcription factors appeared as hub nodes in the protein-protein interaction network of the DEGs under blue light while ERF38 was a hub node among the DEGs under red light. In summary, stress-responsive genes were predominantly up-regulated under blue light to respond to stresses that could be induced under high energy light. The information obtained from this study can be useful to better understand the responses of C. asiatica to different light qualities.

Establishment of Human-Induced Pluripotent Stem Cell-Derived Neurons-A Promising In Vitro Model for a Molecular Study of Rabies Virus and Host Interaction.

Rabies is a deadly viral disease caused by the rabies virus (RABV), transmitted through a bite of an infected host, resulting in irreversible neurological symptoms and a 100% fatality rate in humans. Despite many aspects describing rabies neuropathogenesis, numerous hypotheses remain unanswered and concealed. Observations obtained from infected primary neurons or mouse brain samples are more relevant to human clinical rabies than permissive cell lines; however, limitations regarding the ethical issue and sample accessibility become a hurdle for discovering new insights into virus-host interplays. To better understand RABV pathogenesis in humans, we generated human-induced pluripotent stem cell (hiPSC)-derived neurons to offer the opportunity for an inimitable study of RABV infection at a molecular level in a pathologically relevant cell type. This study describes the characteristics and detailed proteomic changes of hiPSC-derived neurons in response to RABV infection using LC-MS/MS quantitative analysis. Gene ontology (GO) enrichment of differentially expressed proteins (DEPs) reveals temporal changes of proteins related to metabolic process, immune response, neurotransmitter transport/synaptic vesicle cycle, cytoskeleton organization, and cell stress response, demonstrating fundamental underlying mechanisms of neuropathogenesis in a time-course dependence. Lastly, we highlighted plausible functions of heat shock cognate protein 70 (HSC70 or HSPA8) that might play a pivotal role in regulating RABV replication and pathogenesis. Our findings acquired from this hiPSC-derived neuron platform help to define novel cellular mechanisms during RABV infection, which could be applicable to further studies to widen views of RABV-host interaction.

A chromosome-scale assembly of the black gram (Vigna mungo) genome.

Black gram (Vigna mungo) is an important short duration grain legume crop. Black gram seeds provide an inexpensive source of dietary protein. Here, we applied the 10X Genomics linked-read technology to obtain a de novo whole genome assembly of V. mungo cultivated variety Chai Nat 80 (CN80). The preliminary assembly contained 12,228 contigs and had an N50 length of 5.2 Mb. Subsequent scaffolding using the long-range Chicago and HiC techniques yielded the first high-quality, chromosome-level assembly of 499 Mb comprising 11 pseudomolecules. Comparative genomics analyses based on sequence information from single-copy orthologous genes revealed that black gram and mungbean (Vigna radiata) diverged about 2.7 million years ago . The transversion rate (4DTv) analysis in V. mungo revealed no evidence supporting a recent genome-wide duplication event observed in the tetraploid créole bean (Vigna reflexo-pilosa). The proportion of repetitive elements in the black gram genome is slightly lower than the numbers reported for related Vigna species. The majority of long terminal repeat retrotransposons appeared to integrate into the genome within the last five million years. We also examined alternative splicing events in V. mungo using full-length transcript sequences. While intron retention was the most prevalent mode of alternative splicing in several plant species, alternative 3' acceptor site selection represented the majority of events in black gram. Our high-quality genome assembly along with the genomic variation information from the germplasm provides valuable resources for accelerating the development of elite varieties through marker-assisted breeding and for future comparative genomics and phylogenetic studies in legume species.

De novo assemblies of Luffa acutangula and Luffa cylindrica genomes reveal an expansion associated with substantial accumulation of transposable elements.

Luffa spp. (sponge gourd or ridge gourd) is an economically important vegetable crop widely cultivated in China, India and Southeast Asia. Here, we employed PacBio long-read single-molecule real-time (SMRT) sequencing to perform de novo genome assemblies of two commonly cultivated Luffa species, L. acutangula and L. cylindrica. We obtained preliminary draft genomes of 734.6 Mb and 689.8 Mb with scaffold N50 of 786,130 and 578,616 bases for L. acutangula and L. cylindrica, respectively. We also applied long-range Chicago and HiC techniques to obtain the first chromosome-scale whole-genome assembly of L. acutangula. The final assembly contained 13 pseudomolecules, corresponding to the haploid chromosome number in Luffa spp. (1n = 13, 2n = 26). The sizes of the assembled Luffa genomes are approximately twice as large as the genome assemblies of related Cucurbitaceae. A large proportion of L. acutangula (62.17%; 456.69 Mb) and L. cylindrica (56.78%; 391.65 Mb) genome assemblies contained repetitive elements. Phylogenetic analyses revealed that the substantial accumulation of transposable elements likely contributed to the expansion of the Luffa genomes. We also investigated alternative splicing events in Luffa using full-length transcript sequences obtained from PacBio Isoform Sequencing (Iso-seq). While the predominant form of alternative splicing in most plant species examined was intron retention, alternative 3' acceptor site selection appeared to be a major event observed in Luffa. High-quality genome assemblies for L. acutangula and L. cylindrica reported here provide valuable resources for Luffa breeding and future genetics and comparative genomics studies in Cucurbitaceae.

Differential expression between drought-tolerant and drought-sensitive sugarcane under mild and moderate water stress as revealed by a comparative analysis of leaf transcriptome.

Sugarcane contributes 80% of global sugar production and to bioethanol generation for the bioenergy industry. Its productivity is threatened by drought that can cause up to 60% yield loss. This study used RNA-Seq to gain a better understanding of the underlying mechanism by which drought-tolerant sugarcane copes with water stress. We compared gene expression in KPS01-12 (drought-tolerant genotype) and UT12 (drought-sensitive genotype) that have significantly different yield loss rates under drought conditions. We treated KPS01-12 and UT12 with mild and moderate water stress and found differentially expressed genes in various biological processes. KPS01-12 had higher expression of genes that were involved in water retention, antioxidant secondary metabolite biosynthesis, and oxidative and osmotic stress response than UT12. In contrast, the sensitive genotype had more down-regulated genes that were involved in photosynthesis, carbon fixation and Calvin cycle than the tolerant genotype. Our obtained expression profiles suggest that the tolerant sugarcane has a more effective genetic response than the sensitive genotype at the initiation of drought stress. The knowledge gained from this study may be applied in breeding programs to improve sugarcane production in drought conditions.

Chloroplast genome data of Luffa acutangula and Luffa aegyptiaca and their phylogenetic relationships.

Luffa acutangula and Luffa aegyptiaca are domesticated plants in the family Cucurbitaceae. They are mainly cultivated in the tropical and subtropical regions of Asia. The chloroplast genomes of many Cucurbitaceae species were sequenced to examine gene content and evolution. However, the chloroplast genome sequences of L. acutangula and L. aegyptiaca have not been reported. We report the first complete sequences of L. acutangula and L. aegyptiaca chloroplast genomes obtained from Pacific Biosciences sequencing and use them to infer evolutionary relationships. The chloroplast genomes of L. acutangula and L. aegyptiaca are 157,202 and 157,275 bp, respectively. Both genomes possessed the typical quadripartite structure and contained 131 genes, including 87 coding genes, 36 tRNA genes and 8 rRNA genes. We identified simple sequence repeats (SSR) and single nucleotide polymorphisms (SNP) from both chloroplast genomes. Polycistronic mRNA was examined in L. acutangula and L. aegyptiaca using RNA sequences from Isoform sequencing to identify co-transcribed genes. IR size and locations were compared to other species and found to be relatively unchanged. Phylogenetic analysis confirmed the close relationship between L. acutangula and L. aegyptiaca in the Cucurbitaceae lineage and showed separation of the Luffa monophyletic clade from other species in the subtribe Sicyocae. The results obtained from this study can be useful for studying the evolution of Cucurbitaceae plants.

The Genome and Transcriptome Analysis of the Vigna mungo Chloroplast.

Vigna mungo is cultivated in approximately 5 million hectares worldwide. The chloroplast genome of this species has not been previously reported. In this study, we sequenced the genome and transcriptome of the V. mungo chloroplast. We identified many positively selected genes in the photosynthetic pathway (e.g., rbcL, ndhF, and atpF) and RNA polymerase genes (e.g., rpoC2) from the comparison of the chloroplast genome of V. mungo, temperate legume species, and tropical legume species. Our transcriptome data from PacBio isoform sequencing showed that the 51-kb DNA inversion could affect the transcriptional regulation of accD polycistronic. Using Illumina deep RNA sequencing, we found RNA editing of clpP in the leaf, shoot, flower, fruit, and root tissues of V. mungo. We also found three G-to-A RNA editing events that change guanine to adenine in the transcripts transcribed from the adenine-rich regions of the ycf4 gene. The edited guanine bases were found particularly in the chloroplast genome of the Vigna species. These G-to-A RNA editing events were likely to provide a mechanism for correcting DNA base mutations. The V. mungo chloroplast genome sequence and the analysis results obtained in this study can apply to phylogenetic studies and chloroplast genome engineering.

Molecular dynamics exploration of poration and leaking caused by Kalata B1 in HIV-infected cell membrane compared to host and HIV membranes.

The membrane disruption activities of kalata B1 (kB1) were investigated using molecular dynamics simulations with membrane models. The models were constructed to mimic the lipid microdomain formation in membranes of HIV particle, HIV-infected cell, and host cell. The differences in the lipid ratios of these membranes caused the formation of liquid ordered (lo) domains of different sizes, which affected the binding and activity of kB1. Stronger kB1 disruptive activity was observed for the membrane with small sized lo domain. Our results show that kB1 causes membrane leaking without bilayer penetration. The membrane poration mechanism involved in the disorganization of the lo domain and in cholesterol inter-leaflet translocation is described. This study enhances our understanding of the membrane activity of kB1, which may be useful for designing novel and potentially therapeutic peptides based on the kB1 framework.

Computational Design of Hypothetical New Peptides Based on a Cyclotide Scaffold as HIV gp120 Inhibitor.

Cyclotides are a family of triple disulfide cyclic peptides with exceptional resistance to thermal/chemical denaturation and enzymatic degradation. Several cyclotides have been shown to possess anti-HIV activity, including kalata B1 (KB1). However, the use of cyclotides as anti-HIV therapies remains limited due to the high toxicity in normal cells. Therefore, grafting anti-HIV epitopes onto a cyclotide might be a promising approach for reducing toxicity and simultaneously improving anti-HIV activity. Viral envelope glycoprotein gp120 is required for entry of HIV into CD4+ T cells. However, due to a high degree of variability and physical shielding, the design of drugs targeting gp120 remains challenging. We created a computational protocol in which molecular modeling techniques were combined with a genetic algorithm (GA) to automate the design of new cyclotides with improved binding to HIV gp120. We found that the group of modified cyclotides has better binding scores (23.1%) compared to the KB1. By using molecular dynamic (MD) simulation as a post filter for the final candidates, we identified two novel cyclotides, GA763 and GA190, which exhibited better interaction energies (36.6% and 22.8%, respectively) when binding to gp120 compared to KB1. This computational design represents an alternative tool for modifying peptides, including cyclotides and other stable peptides, as therapeutic agents before the synthesis process.

Dynamic scenario of membrane binding process of kalata b1.

Kalata B1 (kB1), a cyclotide that has been used in medical applications, displays cytotoxicity related to membrane binding and oligomerization. Our molecular dynamics simulation results demonstrate that Trp19 in loop 5 of both monomeric and tetrameric kB1 is a key residue for initial anchoring in the membrane binding process. This residue also facilitates the formation of kB1 tetramers. Additionally, we elucidate that kB1 preferentially binds to the membrane interfacial zone and is unable to penetrate into the membrane. In particular, significant roles of amino acid residues in loop 5 and loop 6 on the localization of kB1 to this membrane-water interface zone are found. This study reveals the roles of amino acid residues in the bioactivity of kB1, which is information that can be useful for designing new therapeutic cyclotides with less toxicity.

Defining the membrane disruption mechanism of kalata B1 via coarse-grained molecular dynamics simulations.

Kalata B1 has been demonstrated to have bioactivity relating to membrane disruption. In this study, we conducted coarse-grained molecular dynamics simulations to gain further insight into kB1 bioactivity. The simulations were performed at various concentrations of kB1 to capture the overall progression of its activity. Two configurations of kB1 oligomers, termed tower-like and wall-like clusters, were detected. The conjugation between the wall-like oligomers resulted in the formation of a ring-like hollow in the kB1 cluster on the membrane surface. Our results indicated that the molecules of kB1 were trapped at the membrane-water interface. The interfacial membrane binding of kB1 induced a positive membrane curvature, and the lipids were eventually extracted from the membrane through the kB1 ring-like hollow into the space inside the kB1 cluster. These findings provide an alternative view of the mechanism of kB1 bioactivity that corresponds with the concept of an interfacial bioactivity model.