Domestication has shaped the population structure and agronomic traits of tea plants, yet the complexity of tea population structure and genetic variation that determines these traits remains unclear. We here investigated the resequencing data of 363 diverse tea accessions collected extensively from almost all tea distributions and found that the population structure of tea plants was divided into eight subgroups, which were basically consistent with their geographical distributions. The genetic diversity of tea plants in China decreased from southwest to east as latitude increased. Results also indicated that Camellia sinensis var. assamica (CSA) illustrated divergent selection signatures with Camellia sinensis var. sinensis (CSS). The domesticated genes of CSA were mainly involved in leaf development, flavonoid and alkaloid biosynthesis, while the domesticated genes in CSS mainly participated in amino acid metabolism, aroma compounds biosynthesis, and cold stress. Comparative population genomics further identified ~730 Mb novel sequences, generating 6,058 full-length protein-encoding genes, significantly expanding the gene pool of tea plants. We also discovered 217,376 large-scale structural variations and 56,583 presence and absence variations (PAVs) across diverse tea accessions, some of which were associated with tea quality and stress resistance. Functional experiments demonstrated that two PAV genes (CSS0049975 and CSS0006599) were likely to drive trait diversification in cold tolerance between CSA and CSS tea plants. The overall findings not only revealed the genetic diversity and domestication of tea plants, but also underscored the vital role of structural variations in the diversification of tea plant traits.

BACKGROUND: The genetic diversity of yak, a key domestic animal on the Qinghai-Tibetan Plateau (QTP), is a vital resource for domestication and breeding efforts. This study presents the first yak pangenome obtained through the de novo assembly of 16 yak genomes.
RESULTS: We discovered 290 Mb of nonreference sequences and 504 new genes. Our pangenome-wide presence and absence variation (PAV) analysis revealed 5,120 PAV-related genes, highlighting a wide range of variety-specific genes and genes with varying frequencies across yak populations. Principal component analysis (PCA) based on binary gene PAV data classified yaks into three new groups: wild, domestic, and Jinchuan. Moreover, we proposed a \'two-haplotype genomic hybridization model\' for understanding the hybridization patterns among breeds by integrating gene frequency, heterozygosity, and gene PAV data. A gene PAV-GWAS identified a novel gene (BosGru3G009179) that may be associated with the multirib trait in Jinchuan yaks. Furthermore, an integrated transcriptome and pangenome analysis highlighted the significant differences in the expression of core genes and the mutational burden of differentially expressed genes between yaks from high and low altitudes. Transcriptome analysis across multiple species revealed that yaks have the most unique differentially expressed mRNAs and lncRNAs (between high- and low-altitude regions), especially in the heart and lungs, when comparing high- and low-altitude adaptations.
CONCLUSIONS: The yak pangenome offers a comprehensive resource and new insights for functional genomic studies, supporting future biological research and breeding strategies.

Gummy stem blight (GSB), a widespread disease causing great loss to cucurbit production, has become a major threat to melon cultivation. However, the melon-GSB interaction remains largely unknown. Here, full-length transcriptome and widely targeted metabolome were used to investigate the defence responses of resistant (PI511089) and susceptible (Payzawat) melon accessions to GSB pathogen infection at 24 h. The biosynthesis of secondary metabolites and MAPK signalling pathway were specifically enriched for differentially expressed genes in PI511890, while carbohydrate metabolism and amino acid metabolism were specifically enriched in Payzawat. More than 1000 novel genes were identified and MAPK signalling pathway was specifically enriched for them in PI511890. There were 11 793 alternative splicing events involving in the defence response to GSB. Totally, 910 metabolites were identified in Payzawat and PI511890, and flavonoids were the dominant metabolites. Integrated full-length transcriptome and metabolome analysis showed eriodictyol and oxalic acid were the potential marker metabolites for GSB resistance in melon. Moreover, posttranscription regulation was widely involved in the defence response of melon to GSB pathogen infection. These results not only improve our understanding on the interaction between melon and GSB, but also facilitate the genetic improvement of melon with GSB resistance.

There is an increasing understanding that a reference genome representing an individual cannot capture all the gene repertoire of a species. Here, we conduct a population-scale missing sequences detection of Chinese domestic pigs using whole-genome sequencing data from 534 individuals. We identify 132.41 Mb of sequences absent in the reference assembly, including eight novel genes. In particular, the breeds spread in Chinese high-altitude regions perform significantly different frequencies of new sequences in promoters than other breeds. Furthermore, we dissect the role of non-coding variants and identify a novel sequence inserted in the 3\'UTR of the FMO3 gene, which may be associated with the intramuscular fat phenotype. This novel sequence could be a candidate marker for meat quality. Our study provides a comprehensive overview of the missing sequences in Chinese domestic pigs and indicates that this dataset is a valuable resource for understanding the diversity and biology of pigs.

The ability of the model organism, Caenorhabditis elegans, to distinguish and escape from pathogenic bacteria has been extensively studied; however, studies on the repulsive response of Meloidogyne incognita are still in their infancy. We have recently demonstrated that biocontrol bacteria induce a repulsive response in M. incognita via two classical signaling pathways. The present study aimed to identify the novel genes and signaling molecules of M. incognita that potentially contribute to its defense reaction. Analysis of the transcriptome data of M. incognita with and without a repulsive response against Bacillus nematocida B16 obtained 15 candidate genes, of which the novel genes Minc3s01748g26034 and Minc3s02548g30585 were found to regulate the aversive behavior of M. incognita, and their functions were further validated. To further confirm the neuronal localization of the two novel genes in M. incognita, in situ hybridization was conducted using the digoxin-labeled probes of ten tag genes, and preferentially profiled the localization of amphid sensory neurons of M. incognita. Analysis of the overviewed neuronal map suggested that Minc3s01748g26034 and Minc3s02548g30585 functioned in ASK/ASI and CEPD/V neurons, respectively. During their interactions, the volatile compounds 3-methyl-butyric acid and 2-methyl-butyric acid produced by the biocontrol bacteria were predicted as the primary signaling molecules that promoted the repulsive behavior of M. incognita against biocontrol bacteria. The findings provided novel insights into the mechanisms underlying the repulsive response of M. incognita that are different from the canonical molecular pathways previously found in C. elegans and can aid in developing novel strategies for controlling root-knot nematodes.

BACKGROUND: Oriental river prawn (Macrobrachium nipponense) is one of the most dominant species in shrimp farming in China, which is a rich source of protein and contributes to a significant impact on the quality of human life. Thus, more complete and accurate annotation of gene models are important for the breeding research of oriental river prawn.
RESULTS: A full-length transcriptome of oriental river prawn muscle was obtained using the PacBio Sequel platform. Then, 37.99 Gb of subreads were sequenced, including 584,498 circular consensus sequences, among which 512,216 were full length non-chimeric sequences. After Illumina-based correction of long PacBio reads, 6,599 error-corrected isoforms were identified. Transcriptome structural analysis revealed 2,263 and 2,555 alternative splicing (AS) events and alternative polyadenylation (APA) sites, respectively. In total, 620 novel genes (NGs), 197 putative transcription factors (TFs), and 291 novel long non-coding RNAs (lncRNAs) were identified.
CONCLUSIONS: In summary, this study offers novel insights into the transcriptome complexity and diversity of this prawn species, and provides valuable information for understanding the genomic structure and improving the draft genome annotation of oriental river prawn.

Klebsiella pneumoniae is not only a human and animal opportunistic pathogen, but a food-borne pathogen. Cross-kingdom infection has been focused on since K. pneumoniae was identified as the pathogen of maize, banana, and pomegranate. Although the pathogenicity of K. pneumoniae strains (from ditch water, maize, and human) on plant and mice has been confirmed, there are no reports to explain the molecular mechanisms of the pathogen. This study uncovered the K. pneumoniae KpC4 isolated from maize top rot for the determination of various virulence genes and resistance genes. At least thirteen plant disease-causing genes are found to be involved in the disruption of plant defense. Among them, rcsB is responsible for causing disease in both plants and animals. The novel sequence types provide solid evidence that the pathogen invades plant and has robust ecological adaptability. It is imperative to perform further studies on the verification of these KpC4 genes’ functions to understand the molecular mechanisms involved in plant−pathogen interactions.

UNASSIGNED: Sugarcane (Saccharum spontaneum L.), the major sugar and biofuel feedstock crop, is cultivated mainly by vegetative propagation worldwide due to the infertility of female reproductive organs resulting in the reduction of quality and output of sugar. Deciphering the gene expression profile during ovule development will improve our understanding of the complications underlying sexual reproduction in sugarcane. Optimal reference genes are essential for elucidating the expression pattern of a given gene by quantitative real-time PCR (qRT-PCR).
UNASSIGNED: In this study, based on transcriptome data obtained from sugarcane ovule, eighteen candidate reference genes were identified, cloned, and their expression levels were evaluated across five developmental stages ovule (AC, MMC, Meiosis, Mitosis, and Mature).
UNASSIGNED: Our results indicated that FAB2 and MOR1 were the most stably expressed genes during sugarcane female gametophyte development. Moreover, two genes, cell cycle-related genes REC8 and CDK, were selected, and their feasibility was validated. This study provides important insights into the female gametophyte development of sugarcane and reports novel reference genes for gene expression research on sugarcane sexual reproduction.

BACKGROUND: Eimeria necatrix is one of the most pathogenic parasites, causing high mortality in chickens. Although its genome sequence has been published, the sequences and complete structures of its mRNA transcripts remain unclear, limiting exploration of novel biomarkers, drug targets and genetic functions in E. necatrix.
METHODS: Second-generation merozoites (MZ-2) of E. necatrix were collected using Percoll density gradients, and high-quality RNA was extracted from them. Single-molecule real-time (SMRT) sequencing and Illumina sequencing were combined to generate the transcripts of MZ-2. Combined with the SMRT sequencing data of sporozoites (SZ) collected in our previous study, the transcriptome and transcript structures of E. necatrix were studied.
RESULTS: SMRT sequencing yielded 21,923 consensus isoforms in MZ-2. A total of 17,151 novel isoforms of known genes and 3918 isoforms of novel genes were successfully identified. We also identified 2752 (SZ) and 3255 (MZ-2) alternative splicing (AS) events, 1705 (SZ) and 1874 (MZ-2) genes with alternative polyadenylation (APA) sites, 4019 (SZ) and 2588 (MZ-2) fusion transcripts, 159 (SZ) and 84 (MZ-2) putative transcription factors (TFs) and 3581 (SZ) and 2039 (MZ-2) long non-coding RNAs (lncRNAs). To validate fusion transcripts, reverse transcription-PCR was performed on 16 candidates, with an accuracy reaching up to 87.5%. Sanger sequencing of the PCR products further confirmed the authenticity of chimeric transcripts. Comparative analysis of transcript structures revealed a total of 3710 consensus isoforms, 815 AS events, 1139 genes with APA sites, 20 putative TFs and 352 lncRNAs in both SZ and MZ-2.
CONCLUSIONS: We obtained many long-read isoforms in E. necatrix SZ and MZ-2, from which a series of lncRNAs, AS events, APA events and fusion transcripts were identified. Information on TFs will improve understanding of transcriptional regulation, and fusion event data will greatly improve draft versions of gene models in E. necatrix. This information offers insights into the mechanisms governing the development of E. necatrix and will aid in the development of novel strategies for coccidiosis control.

The transdifferentiation of adipose-derived mesenchymal stem cells (ADMSCs) into insulin-producing cells (IPCs) is a potential resource for the treatment of diabetes. However, the changes of genes and metabolic pathways on the transdifferentiation of ADMSCs into IPCs are largely unknown. In this study, the transdifferentiation of canine ADMSCs into IPCs was completed using five types of procedures. Absolute Quantitative Transcriptome Sequencing Analysis was performed at different stages of the optimal procedure. A total of 60,151 transcripts were obtained. Differentially expressed genes (DEGs) were divided into five groups: IPC1 vs. ADSC (1169 upregulated genes and 1377 downregulated genes), IPC2 vs. IPC1 (1323 upregulated genes and 803 downregulated genes), IPC3 vs. IPC2 (722 upregulated genes and 680 downregulated genes), IPC4 vs. IPC3 (539 upregulated genes and 1561 downregulated genes), and Beta_cell vs. IPC4 (2816 upregulated genes and 4571 downregulated genes). The gene ontology (GO) and Kyoto Encyclopedia of Genes and Genomes (KEGG) enrichment analysis of DEGs revealed that many genes and signaling pathways that are essential for transdifferentiation. Hnf1B, Dll1, Pbx1, Rfx3, and Foxa1 were screened out, and the functions of five genes were verified further by overexpression and silence. Foxa1, Pbx1, and Rfx3 exhibited significant effects, can be used as specific key regulatory factors in the transdifferentiation of ADMSCs into IPCs. This study provides a foundation for future work to understand the mechanisms of the transdifferentiation of ADMSCs into IPCs and acquire IPCs with high maturity.