Primer premier: program for design of degenerate primers from a protein sequence

Comparative plastomes of eight subgenus Chamaesyce plants and system authentication of Euphorbiae Humifusae Herba

Euphorbiae Humifusae Herba (EHH) was provided with medicinal and edible uses, but frequently was adulterated with its closely related species. Hence, this study sought to identify EHH via an integrated approach comprising data from its morphological evaluation, HPLC analysis, comparative plastomes analysis and allele-specific PCR identification. First, the morphological characteristics of 8 subgenus Chamaesyce plants were summarized. Then, HPLC analysis showed that 18 batches of EHH were adulterated or unqualified. Furthermore, the plastomes of the 8 subg. Chamaesyce species were analyzed. Phylogenetic analysis revealed a sister relationship among the 8 subg. Chamaesyce species. The allele-specific PCR authentication was developed by the nucleotide polymorphisms (SNPs) and insertions or deletions (InDels) analysis. The results of allele-specific PCR showed that 27 batches of EHH were adulterated, indicating that the superior sensitivity of molecular authentication over the other methods used. This study provided a reference for rational use and phylogenetic research of EHH.

A 5.2-kb insertion in the coding sequence of PavSCPL, a serine carboxypeptidase-like enhances fruit firmness in Prunus avium

Fruit firmness is an important trait in sweet cherry breeding because it directly positively influences fruit transportability, storage and shelf life. However, the underlying genes responsible and the molecular mechanisms that control fruit firmness remain unknown. In this study, we identified a candidate gene, PavSCPL, encoding a serine carboxypeptidase-like protein with natural allelic variation, that controls fruit firmness in sweet cherry using map-based cloning and functionally characterized PavSCPL during sweet cherry fruit softening. Genetic analysis revealed that fruit firmness in the 'Rainier' × 'Summit' F1 population was controlled by a single dominant gene. Bulked segregant analysis combined with fine mapping narrowed the candidate gene to a 473-kb region (7418778-7 891 914 bp) on chromosome 6 which included 72 genes. The candidate gene PavSCPL, and a null allele harbouring a 5244-bp insertion in the second exon that completely inactivated PavSCPL expression and resulted in the extra-hard-flesh phenotype, were identified by RNA-sequencing analysis and gene cloning. Quantitative RT-PCR analysis revealed that the PavSCPL expression level was increased with fruit softening. Virus-induced gene silencing of PavSCPL enhanced fruit firmness and suppressed the activities of certain pectin-degrading enzymes in the fruit. In addition, we developed functional molecular markers for PavSCPL and the Pavscpl5.2-k allele that co-segregated with the fruit firmness trait. Overall, this research identified a crucial functional gene for fruit firmness. The results provide insights into the genetic control and molecular mechanism of the fruit firmness trait and present useful molecular markers for molecular-assisted breeding for fruit firmness in sweet cherry.

Whole transcriptome sequencing reveals key genes and ceRNA regulatory networks associated with pimpled eggs in hens

Eggshell is one of the most important indicators of egg quality, and due to low shell strength, pimple eggs (PE) are more susceptible to breakage, thus causing huge economic losses to the egg industry. At the current time, the molecular mechanisms that regulate the formation of pimple eggs are poorly understood. In this study, uterine tissues of PE-laying hens (n = 8) and normal egg (NE) -laying hens (n = 8) were analyzed by whole transcriptome sequencing, and a total of 619 differentially expressed mRNAs (DE mRNAs), 122 differentially expressed lncRNAs (DE lncRNAs) and 21 differentially expressed miRNAs (DE miRNAs) were obtained. Based on the targeting relationship among DE mRNAs, DE lncRNAs and DE miRNAs, we constructed a competitive endogenous RNA (ceRNA) network including 12 DE miRNAs, 19 DE lncRNAs, and 128 DE mRNAs. Considering the large amount of information contained in the network, we constructed a smaller ceRNA network to better understand the complex mechanisms of pimple egg formation. The smaller ceRNA network network contains 7 DE lncRNAs (LOC107056551, LOC121109367, LOC121108909, LOC121108862, LOC112530033, LOC121113165, LOC107054145), 5 DE miRNAs (gga-miR-6568-3p, gga-miR-31-5p, gga-miR-18b-3p, gga-miR-1759-3p, gga-miR-12240-3p) and 7 DE mRNAs (CABP1, DNAJC5, HCN3, HPCA, IBSP, KCNT1, OTOP3), and these differentially expressed genes may play key regulatory roles in the formation of pimpled eggs in hens. This study provides the overall expression profiles of mRNAs, lncRNAs and miRNAs in the uterine tissues of hens, which provides a theoretical basis for further research on the molecular mechanisms of pimpled egg formation, and has potential applications in improving eggshell quality.

The Phylogeography of Deciduous Tree Ulmus macrocarpa (Ulmaceae) in Northern China

Disentangling how climate oscillations and geographical events significantly influence plants' genetic architecture and demographic history is a central topic in phytogeography. The deciduous ancient tree species Ulmus macrocarpa is primarily distributed throughout Northern China and has timber and horticultural value. In the current study, we studied the phylogenic architecture and demographical history of U. macrocarpa using chloroplast DNA with ecological niche modeling. The results indicated that the populations' genetic differentiation coefficient (NST) value was significantly greater than the haplotype frequency (GST) (p < 0.05), suggesting that U. macrocarpa had a clear phylogeographical structure. Phylogenetic inference showed that the putative chloroplast haplotypes could be divided into three groups, in which the group Ⅰ was considered to be ancestral. Despite significant genetic differentiation among these groups, gene flow was detected. The common ancestor of all haplotypes was inferred to originate in the middle-late Miocene, followed by the haplotype overwhelming diversification that occurred in the Quaternary. Combined with demography pattern and ecological niche modeling, we speculated that the surrounding areas of Shanxi and Inner Mongolia were potential refugia for U. macrocarpa during the glacial period in Northern China. Our results illuminated the demography pattern of U. macrocarpa and provided clues and references for further population genetics investigations of precious tree species distributed in Northern China.

Mechanism of Apoptosis in Porcine Ovarian Granulosa Cells Triggered by T-2 Toxin

T-2 toxin (T-2), an A-type mono mycotoxin produced by various Fusarium species, disrupts DNA/RNA and protein synthesis upon entering the body, resulting in pathological conditions in various tissues/organs and posing a significant threat to human and animal health. However, the mechanisms underlying its toxicity remain unclear. With the goal of learning how T-2 affects reproduction in animals, we utilized primary porcine ovarian granulosa cells (pGCs) as a carrier in vitro and constructed concentration models for analyzing cell morphology and RNA-sequencing (RNA-seq). Our findings showed that T-2 could influence pGCs morphology, induce cell cycle arrest, and promote apoptosis in a dose-dependent manner. The results of RNA-seq analyses indicated that a total of 8216 genes exhibited significant differential expression (DEG) following T-2 treatment, of which 4812 were observed to be down-regulated and 3404 were up-regulated. The DEGs following T-2 toxin treatment of pGCs had a notable impact on many metabolic pathways such as PI3K-Akt, Ras, MAPK, and apoptosis, which in turn altered important physiological processes. Gene set enrichment analysis (GSEA) indicated that the differences in the harmful effects of T-2 might be caused by the varying control of cellular processes and the pathway responsible for steroid metabolism. These results present further insights regarding the mechanism of T-2 action on sow reproductive toxicity, enhance our understanding of T-2 reproductive toxicological effects, and lay a theoretical foundation for the judicious prevention of T-2-induced reproductive toxicity.

TMT-Based Quantitative Proteomic Analysis Reveals the Key Role of Cell Proliferation and Apoptosis in Intestine Regeneration of Apostichopus japonicus

Sea cucumbers are widely known for their powerful regenerative abilities, which allow them to regenerate a complete digestive tract within a relatively short time following injury or autotomy. Recently, even though the histological changes and cellular events in the processes of intestinal regeneration have been extensively studied, the molecular machinery behind this faculty remains unclear. In this study, tandem mass tag (TMT)-based quantitation was utilized to investigate protein abundance changes during the process of intestine regeneration. Approximately 538, 445, 397, 1012, and 966 differential proteins (DEPs) were detected (p < 0.05) between the normal and 2, 7, 12, 20, and 28 dpe stages, respectively. These DEPs also mainly focus on pathways of cell proliferation and apoptosis, which were further validated by 5-Ethynyl-2'-deoxyuridine (EdU) or Tunel-based flow cytometry assay. These findings provide a reference for a comprehensive understanding of the regulatory mechanisms of various stages of intestinal regeneration and provide a foundation for subsequent research on changes in cell fate in echinoderms.

Transcriptional dynamics and regulatory function of milRNAs in Ascosphaera apis invading Apis mellifera larvae

In the present study, small RNA (sRNA) data from Ascosphaera apis were filtered from sRNA-seq datasets from the gut tissues of A. apis-infected Apis mellifera ligustica worker larvae, which were combined with the previously gained sRNA-seq data from A. apis spores to screen differentially expressed milRNAs (DEmilRNAs), followed by trend analysis and investigation of the DEmilRNAs in relation to significant trends. Additionally, the interactions between the DEmilRNAs and their target mRNAs were verified using a dual-luciferase reporter assay. In total, 974 A. apis milRNAs were identified. The first base of these milRNAs was biased toward U. The expression of six milRNAs was confirmed by stem-loop RT-PCR, and the sequences of milR-3245-y and milR-10285-y were validated using Sanger sequencing. These miRNAs grouped into four significant trends, with the target mRNAs of DEmilRNAs involving 42 GO terms and 120 KEGG pathways, such as the fungal-type cell wall and biosynthesis of secondary metabolites. Further investigation demonstrated that 299 DEmilRNAs (novel-m0011-3p, milR-10048-y, bantam-y, etc.) potentially targeted nine genes encoding secondary metabolite-associated enzymes, while 258 (milR-25-y, milR-14-y, milR-932-x, etc.) and 419 (milR-4561-y, milR-10125-y, let-7-x, etc.) DEmilRNAs putatively targeted virulence factor-encoded genes and nine genes involved in the MAPK signaling pathway, respectively. Additionally, the interaction between ADM-B and milR-6882-x, as well as between PKIA and milR-7009-x were verified. Together, these results not only offer a basis for clarifying the mechanisms underlying DEmilRNA-regulated pathogenesis of A. apis and a novel insight into the interaction between A. apis and honey bee larvae, but also provide candidate DEmilRNA-gene axis for further investigation.

Transcriptome analysis of Schizothorax oconnori (Cypriniformes: Cyprinidae) oocytes: The role of K+ in promoting yolk globule fusion and regulating oocyte maturation

Schizothorax oconnori (S. oconnori) is an economically important fish in Tibet. Oocyte maturation is a physiological process that is of great significance to reproduction and seed production in S. oconnori, yet little is currently known regarding the molecular mechanisms of oocyte development in this species. To identify candidate genes involved in reproduction of female fish, a combination of PacBio and Illumina HiSeq technologies was employed to provide deep coverage of the oocyte transcriptome. Transcriptome analysis revealed several candidate genes that are potentially involved in the regulation of oocyte maturation in S. oconnori, including GIRK1, CHRM3, NPY2R, GABRA3, GnRH3, mGluR1α, GPER1, GDF9, HSP90, and ESR2. Genes that are significantly expressed during oocyte maturation mainly contribute to the GPCR signaling pathway and the estrogen signaling pathway. Neurotransmitter (Ach, NPY, and GABA) and peptide hormone (GnRH3) binding to G protein-coupled receptors (GPCRs) frees G-protein βγ subunits to interact with the G protein-gated inward rectifier K+ channel 1 (GIRK1). This process helps release K+ from granulosa cells to maturing oocytes, allowing yolk globule fusion. This mechanism may play an important role in oocyte maturation in S. oconnori. In conclusion, this study provides a valuable basis for deciphering the reproductive system in S. oconnori during the oocyte maturation process.

Population genetics reveals new introgression in the nucleus herd of min pigs

OBJECTIVE

Min pigs are a unique genetic resource among local pig breeds in China. They have more excellent characteristics in cold and stress resistance, good meat quality, and a high reproductive rate. However, the genetic structure and driving factors remain unclear in the nucleus herd. In this study, the genetic diversity of Min pigs was studied to reveal the formation mechanism of its unique genetic structure. We hope to protect and develop the genetic resources of Min pigs.

METHODS

We analyzed different types of genes to identify the genetic structure and gene introgression pattern of Min pigs. The nuclear DNA dataset includes information on 21 microsatellite loci and 6 Y-chromosome genes, and the mitochondrial D-loop gene is selected to represent maternal lineages. The above genes are all from the nucleus herd of Min pigs.

RESULTS

The results of genetic structure identification and analysis of potential exogenous gene introgression patterns indicate that the nucleus herd of Min pigs maintains a high level of genetic diversity (polymorphism information content = 0.713, expected heterozygosity = 0.662, observed heterozygosity = 0.612). Compared with other Asian pig breeds, the formation of Min pig breeds is more special. Gene introgression from European pig breeds to Min pigs has occurred, which is characterized by complete introgression of paternal genes and incomplete introgression of maternal genes.

CONCLUSION

Gene introgression caused by cross-breeding is not the main factor leading to the formation of the current genetic structure of Min pigs, but this process has increased the level of genetic diversity in the nucleus herd. Compared with the influence of gene introgression, our research suggest that artificial selection and environmental adaptive evolution make Min pigs form unique genetic characteristics.

Multiple amino acid transporters as carriers load L-valine-phenazine-1-carboxylic acid conjugate into Ricinus sieve tubes for the phloem translocation

Some transporters play important roles in the uptake and acropetal xylem translocation of vectorized agrochemicals. However, it is poorly understood the basipetally phloem-loading functions of transporters toward vectorized agrochemicals. Here, L-Val-PCA (L-valine-phenazine-1-carboxylic acid conjugate) uptake was demonstrated carrier-mediated. RcAAP2, RcANT7, and RcLHT1 showed a similarly up-regulated expression pattern from 62 transporter coding genes in Ricinus at 1 h after L-Val or L-Val-PCA treatment. Subcellular localization revealed that fusion RcAAP2-eGFP, RcANT7-eGFP and RcLHT1-eGFP proteins were expressed in the plasma membrane of mesophyll and phloem cells. Yeast assays found that RcAAP2, RcANT7, and RcLHT1 facilitated L-Val-PCA uptake. To further demonstrate the phloem-loading functions, using vacuum infiltration strategy, an Agrobacterium-mediated RNA interference (RNAi) protocol was constructed in seedlings. HPLC detection indicated that L-Val-PCA phloem sap concentrations were significantly decreased 54.5 %, 27.6 %, and 41.6 % after silencing for 72 h and increased 48.3 %, 52.6 %, and 52.4 % after overexpression, respectively. In conclusion, the plasma membrane-located RcAAP2, RcANT7, and RcLHT1 can loaded L-Val-PCA into Ricinus sieve tubes for the phloem translocation, which may aid in the utilization of transporters and molecular design of phloem-mobile fungicides target root or vascular pathogens.

Genome-Wide Bioinformatics Analysis of SWEET Gene Family and Expression Verification of Candidate PaSWEET Genes in Potentilla anserina

Sugars act as the main energy sources in many fruit and vegetable crops. The biosynthesis and transportation of sugars are crucial and especially contribute to growth and development. SWEET is an important gene family that plays a vital role in plants' growth, development, and adaptation to various types of stresses (biotic and abiotic). Although SWEET genes have been identified in numerous plant species, there is no information on SWEETs in Potentilla anserina. In the present study, we performed a comprehensive genome-wide bioinformatics analysis and identified a total of 23 candidate PaSWEETs genes in the Potentilla anserina genome, which were randomly distributed on ten different chromosomes. The phylogenetic analysis, chromosomal location, gene structure, specific cis-elements, protein interaction network, and physiological characteristics of these genes were systematically examined. The identified results of the phylogenetic relationship with Arabidopsis thaliana revealed that these PaSWEET genes were divided into four clades (I, II, III, and IV). Moreover, tissue-specific gene expression through quantitative real-time polymerase chain reaction (qRT-PCR) validation exposed that the identified PaSWEETs were differentially expressed in various tissues (roots, stems, leaves, and flowers). Mainly, the relative fold gene expression in swollen and unswollen tubers effectively revealed that PaSWEETs (7, 9, and 12) were highly expressed (300-, 120-, and 100-fold) in swollen tubers. To further elucidate the function of PaSWEETs (7, 9, and 12), their subcellular location was confirmed by inserting them into tobacco leaves, and it was noted that these genes were present on the cell membrane. On the basis of the overall results, it is suggested that PaSWEETs (7, 9, and 12) are the candidate genes involved in swollen tuber formation in P. anserina. In crux, we speculated that our study provides a valuable theoretical base for further in-depth function analysis of the PaSWEET gene family and their role in tuber development and further enhancing the molecular breeding of Potentilla anserina.

Development of a SNP-based strain-identified method for Streptococcus thermophilus CICC 6038 and Lactobacillus delbrueckii ssp. bulgaricus CICC 6047 using pan-genomics analysis

The health benefits conferred by probiotics is specific to individual probiotic strains, highlighting the importance of identifying specific strains for research and production purposes. Streptococcus thermophilus CICC 6038 and Lactobacillus delbrueckii ssp. bulgaricus CICC 6047 are exceedingly valuable for commercial use with an excellent mixed-culture fermentation. To differentiate these 2 strains from other S. thermophilus and L. delbrueckii ssp. bulgaricus, a specific, sensitive, accurate, rapid, convenient, and cost-effective method is required. In this study, we conducted a pan-genome analysis of S. thermophilus and L. delbrueckii ssp. bulgaricus to identify species-specific core genes, along with strain-specific single-nucleotide polymorphisms (SNPs). These genes were used to develop suitable PCR primers, and the conformity of sequence length and unique SNPs was confirmed by sequencing for qualitative identification at the strain level. The results demonstrated that SNPs analysis of PCR products derived from these primers could distinguish CICC 6038 and CICC 6047 accurately and reproducibly from the other strains of S. thermophilus and L. delbrueckii ssp. bulgaricus, respectively. The strain-specific PCR method based on SNPs herein is universally applicable for probiotics identification. It offers valuable insights into identifying probiotics at the strain level that is fit-for-purpose in quality control and compliance assessment of commercial dairy products.

A cuticular protein, BgCPLCP1, contributes to insecticide resistance by thickening the cockroach endocuticle

Overuse of insecticides has led to severe environmental problems. Insect cuticle, which consists mainly of chitin, proteins and a thin outer lipid layer, serves multiple functions. Its prominent role is as a physical barrier that impedes the penetration of xenobiotics, including insecticides. Blattella germanica (L.) is a major worldwide indoor pest that causes allergic disease and asthma. Extensive use of pyrethroid insecticides, including β-cypermethrin, has selected for the rapid and independent evolution of resistance in cockroach populations on a global scale. We demonstrated that BgCPLCP1, the first CPLCP (cuticular proteins of low complexity with a highly repetitive proline-rich region) family cuticular protein in order Blattodea, contributes to insecticide penetration resistance. Silencing BgCPLCP1 resulted in 85.0 %-85.7 % and 81.0 %-82.0 % thinner cuticle (and especially thinner endocuticle) in the insecticide-susceptible (S) and β-cypermethrin-resistant (R) strains, respectively. The thinner and more permeable cuticles resulted in 14.4 % and 20.0 % lower survival of β-cypermethrin-treated S- and R-strain cockroaches, respectively. This study advances our understanding of cuticular penetration resistance in insects and opens opportunities for the development of new efficiently and environmentally friendly insecticides targeting the CPLCP family of cuticular proteins.

Genome- and transcriptome-wide identification of trehalose-6-phosphate phosphatases (TPP) gene family and their expression patterns under abiotic stress and exogenous trehalose in soybean

BACKGROUND

Trehalose-6-phosphate phosphatase (TPP) is an essential enzyme catalyzing trehalose synthesis, an important regulatory factor for plant development and stress response in higher plants. However, the TPP gene family in soybean has not been reported.

RESULTS

A comprehensive analysis of the TPP gene family identified 18 GmTPPs classified into eight groups based on the phylogenetic relationships and the conservation of protein in six monocot and eudicot plants. The closely linked subfamilies had similar motifs and intron/exon numbers. Segmental duplication was the main driving force of soybean GmTPPs expansion. In addition, analysis of the cis-regulatory elements and promoter regions of GmTPPs revealed that GmTPPs regulated the response to several abiotic stresses. Moreover, RNA-seq and qRT-PCR analysis of the tissue-specific GmTPPs under different abiotic stresses revealed that most GmTPPs were associated with response to different stresses, including cold, drought, saline-alkali, and exogenous trehalose. Notably, exogenous trehalose treatment up-regulated the expression of most TPP genes under saline-alkali conditions while increasing the carbohydrate and trehalose levels and reducing reactive oxygen species (ROS) accumulation in soybean sprouts, especially in the saline-alkali tolerant genotype. Furthermore, the interaction network and miRNA target prediction revealed that GmTPPs interacted with abiotic stress response-related transcription factors.

CONCLUSIONS

The findings in this study lay a foundation for further functional studies on TPP-based breeding to improve soybean development and stress tolerance.

Stripe rust effector Pst21674 compromises wheat resistance by targeting transcription factor TaASR3

Pathogens compromise host defense responses by strategically secreting effector proteins. However, the molecular mechanisms by which effectors manipulate disease-resistance factors to evade host surveillance remain poorly understood. In this study, we characterized a Puccinia striiformis f. sp. tritici (Pst) effector Pst21674 with a signal peptide. Pst21674 was significantly upregulated during Pst infections in wheat (Triticum aestivum L.) and knocking down Pst21674 by host-induced gene silencing led to reduced Pst pathogenicity and restricted hyphal spread in wheat. Pst21674 interaction with the abscisic acid-, stress-, and ripening-induced protein TaASR3 was validated mainly in the nucleus. Size exclusion chromatography, bimolecular fluorescence complementation, and luciferase complementation imaging assays confirmed that TaASR3 could form a functional tetramer. Virus-induced gene silencing and overexpression demonstrated that TaASR3 contributes to wheat resistance to stripe rust by promoting accumulation of reactive oxygen species and cell death. Additionally, transcriptome analysis revealed that the expression of defense-related genes was regulated in transgenic wheat plants overexpressing TaASR3. Interaction between Pst21674 and TaASR3 interfered with the polymerization of TaASR3 and suppressed TaASR3-mediated transcriptional activation of defense-related genes. These results indicate that Pst21674 serves as an important virulence factor secreted into the host nucleus to impede wheat resistance to Pst, possibly by targeting and preventing polymerization of TaASR3.

Comparative transcriptome analysis reveals candidate genes related to the sex differentiation of Schisandra chinensis

Schisandra chinensis is a monoecious plant with unisex flowers. The fruit of S. chinensis is of high medical with economic value. The yield of S. chinensis fruit is related to the ratio of its female and male flowers. However, there is little research on its floral development and sex differentiation. To elucidate the possible mechanism for the sex differentiation of S. chinensis, we collected 18 samples of female and male flowers from three developmental stages and performed a comparative RNA-seq analysis aimed at identifying differentially expressed genes (DEGs) that may be related to sex differentiation. The results showed 936, 7179, and 6890 differentially expressed genes between female and male flowers at three developmental stages, respectively, and 466 candidate genes may play roles in sex differentiation. KEGG analysis showed genes involved in the flavonoid biosynthesis pathway and DNA replication pathway were essential for the development of female flowers. 51 MADS-box genes and 10 YABBY genes were identified in S. chinensis. The DEGs analysis indicated that MADS-box and YABBY genes were strongly related to the sex determination of S. chinensis. RT-qPCR confirmed the RNA-seq results of 20 differentially expressed genes, including three male-biased genes and 17 female-biased genes. A possible regulatory model of sex differentiation in S. chinensis was proposed according to our results. This study helps reveal the sex-differentiation mechanism of S. chinensis and lays the foundation for regulating the male-female ratio of S. chinensis in the future.

The novel developed microsatellite markers revealed potential hybridization among Cymbidium species and the interspecies sub-division of C. goeringii and C. ensifolium

BACKGROUND

Orchids (Cymbidium spp.) exhibit significant variations in floral morphology, pollinator relations, and ecological habitats. Due to their exceptional economic and ornamental value, Cymbidium spp. have been commercially cultivated for centuries. SSR markers are extensively used genetic tools for biology identification and population genetics analysis.

RESULT

In this study, nine polymorphic EST-SSR loci were isolated from Cymbidium goeringii using RNA-Seq technology. All nine SSR loci showed transferability in seven other congeneric species, including 51 cultivars. The novel SSR markers detected inter-species gene flow among the Cymbidium species and intra-species sub-division of C. goeringii and C. ensifolium, as revealed by neighborhood-joining and Structure clustering analyses.

CONCLUSION

In this study, we developed nine microsatellites using RNA-Seq technology. These SSR markers aided in detecting potential gene flow among Cymbidium species and identified the intra-species sub-division of C. goeringii and C. ensifolium.

Developmental validation of a high-resolution panel genotyping 639 Y-chromosome SNP and InDel markers and its evolutionary features in Chinese populations

Uniparental-inherited haploid genetic marker of Y-chromosome single nucleotide polymorphisms (Y-SNP) have the power to provide a deep understanding of the human evolutionary past, forensic pedigree, and bio-geographical ancestry information. Several international cross-continental or regional Y-panels instead of Y-whole sequencing have recently been developed to promote Y-tools in forensic practice. However, panels based on next-generation sequencing (NGS) explicitly developed for Chinese populations are insufficient to represent the Chinese Y-chromosome genetic diversity and complex population structures, especially for Chinese-predominant haplogroup O. We developed and validated a 639-plex panel including 633 Y-SNPs and 6 Y-Insertion/deletions, which covered 573 Y haplogroups on the Y-DNA haplogroup tree. In this panel, subgroups from haplogroup O accounted for 64.4% of total inferable haplogroups. We reported the sequencing metrics of 354 libraries sequenced with this panel, with the average sequencing depth among 226 individuals being 3,741×. We illuminated the high level of concordance, accuracy, reproducibility, and specificity of the 639-plex panel and found that 610 loci were genotyped with as little as 0.03 ng of genomic DNA in the sensitivity test. 94.05% of the 639 loci were detectable in male-female mixed DNA samples with a mix ratio of 1:500. Nearly all of the loci were genotyped correctly when no more than 25 ng/μL tannic acid, 20 ng/μL humic acid, or 37.5 μM hematin was added to the amplification mixture. More than 80% of genotypes were obtained from degraded DNA samples with a degradation index of 11.76. Individuals from the same pedigree shared identical genotypes in 11 male pedigrees. Finally, we presented the complex evolutionary history of 183 northern Chinese Hans and six other Chinese populations, and found multiple founding lineages that contributed to the northern Han Chinese gene pool. The 639-plex panel proved an efficient tool for Chinese paternal studies and forensic applications.

Mutation to a cytochrome P450 -like gene alters the leaf color by affecting the heme and chlorophyll biosynthesis pathways in Brassica napus

The regulated biosynthesis of chlorophyll is important because of its effects on plant photosynthesis and dry biomass production. In this study, a map-based cloning approach was used to isolate the cytochrome P450 -like gene BnaC08g34840D (BnCDE1) from a chlorophyll-deficient mutant (cde1) of Brassica napus obtained by ethyl methanesulfonate (EMS) mutagenization. Sequence analyses revealed that BnaC08g34840D in the cde1 mutant (BnCDE1I320T ) encodes a substitution at amino acid 320 (Ile320Thr) in the conserved region. The over-expression of BnCDE1I320T in ZS11 (i.e., gene-mapping parent with green leaves) recapitulated a yellow-green leaf phenotype. The CRISPR/Cas9 genome-editing system was used to design two single-guide RNAs (sgRNAs) targeting BnCDE1I320T in the cde1 mutant. The knockout of BnCDE1I320T in the cde1 mutant via a gene-editing method restored normal leaf coloration (i.e., green leaves). These results indicate that the substitution in BnaC08g34840D alters the leaf color. Physiological analyses showed that the over-expression of BnCDE1I320T leads to decreases in the number of chloroplasts per mesophyll cell and in the contents of the intermediates of the chlorophyll biosynthesis pathway in leaves, while it increases heme biosynthesis, thereby lowering the photosynthetic efficiency of the cde1 mutant. The Ile320Thr mutation in the highly conserved region of BnaC08g34840D inhibited chlorophyll biosynthesis and disrupted the balance between heme and chlorophyll biosynthesis. Our findings may further reveal how the proper balance between the chlorophyll and heme biosynthesis pathways is maintained.

Identification of a new allele of BraA09g066480.3C controlling the wax-less phenotype of Chinese cabbage

BACKGROUND

Epidermal wax covers the surfaces of terrestrial plants to resist biotic and abiotic stresses. Wax-less flowering Chinese cabbage (Brassica campestris L. ssp. chinesis var. utilis tsen et lee) has the charateristics of lustrous green leaves and flower stalks, which are of high commercial value.

RESULTS

To clarify the mechanism of the wax deficiency, the wax-less flowering Chinese cabbage doubled-haploid (DH) line 'CX001' and Chinese cabbage DH line 'FT', obtained from isolated microspore culture, were used in the experiments. Genetic analysis showed that the wax-less phenotype of 'CX001' was controlled by a recessive nuclear gene, named wlm1 (wax-less mutation 1), which was fine-mapped on chromosome A09 by bulked segregant analysis sequencing (BSA-seq) of B.rapa genome V3.0. There was only one gene (BraA09g066480.3C) present in the mapping region. The homologous gene in Arabidopsis thaliana is AT1G02205 (CER1) that encodes an aldehyde decarboxylase in the epidermal wax metabolism pathway. Semi-quantitative reverse transcription PCR and transcriptome analysis indicated that BraA09g066480.3C was expressed in 'FT' but not in 'CX001'. BraA09g066480.3C was lost in the CXA genome to which 'CX001' belonged.

CONCLUSION

The work presented herein demonstrated that BraA09g066480.3C was the causal gene for wax-less flowering Chinese cabbage 'CX001'. This study will lay a foundation for further research on the molecular mechanism of epidermal wax synthesis in flowering Chinese cabbage.

Analysis of the RNA Editing Sites and Orthologous Gene Function of Transcriptome and Chloroplast Genomes in the Evolution of Five Deutzia Species

In this study, the chloroplast genomes and transcriptomes of five Deutzia genus species were sequenced, characterized, combined, and analyzed. A phylogenetic tree was constructed, including 32 other chloroplast genome sequences of Hydrangeoideae species. The results showed that the five Deutzia chloroplast genomes were typical circular genomes 156,860-157,025 bp in length, with 37.58-37.6% GC content. Repeat analysis showed that the Deutzia species had 41-45 scattered repeats and 199-201 simple sequence repeats. Comparative genomic and pi analyses indicated that the genomes are conservative and that the gene structures are stable. According to the phylogenetic tree, Deutzia species appear to be closely related to Kirengeshoma palmata and Philadelphus. By combining chloroplast genomic and transcriptomic analyses, 29-31 RNA editing events and 163-194 orthologous genes were identified. The ndh, rpo, rps, and atp genes had the most editing sites, and all RNA editing events were of the C-to-U type. Most of the orthologous genes were annotated to the chloroplast, mitochondria, and nucleus, with functions including energy production and conversion, translation, and protein transport. Genes related to the biosynthesis of monoterpenoids and flavonoids were also identified from the transcriptome of Deutzia spp. Our results will contribute to further studies of the genomic information and potential uses of the Deutzia spp.

Multiple-Genome-Based Simple Sequence Repeat Is an Efficient and Successful Method in Genotyping and Classifying Different Jujube Germplasm Resources

Jujube (Ziziphus jujuba Mill.) is a commercially important tree native to China, known for its high nutritional value and widespread distribution, as well as its diverse germplasm resources. Being resilient to harsh climatic conditions, the cultivation of jujube could provide a solution to food insecurity and income for people of arid and semi-arid regions in and outside of China. The evaluation of germplasm resources and genetic diversity in jujube necessitates the use of Simple Sequence Repeat (SSR) markers. SSR markers are highly polymorphic and can be used to evaluate the genetic diversity within and between cultivars of Chinese jujube, and are important for conservation biology, breeding programs, and the discovery of important traits for Chinese jujube improvement in China and abroad. However, traditional methods of SSR development are time-consuming and inadequate to meet the growing research demands. To address this issue, we developed a novel approach called Multiple-Genome-Based SSR identification (MGB-SSR), which utilizes the genomes of three jujube cultivars to rapidly screen for polymorphic SSRs in the jujube genome. Through the screening process, we identified 12 pairs of SSR primers, which were then used to successfully classify 249 jujube genotypes. Based on the genotyping results, a digital ID card was established, enabling the complete identification of all 249 jujube plants. The MGB-SSR approach proved efficient in rapidly detecting polymorphic SSRs within the jujube genome. Notably, this study represents the first successful differentiation of jujube germplasm resources using 12 SSR markers, with 4 markers successfully identifying triploid jujube genotypes. These findings offer valuable information for the classification of Chinese jujube germplasm, thereby providing significant assistance to jujube researchers and breeders in identifying unknown jujube germplasm.

DNA-Binding One Finger Transcription Factor PhDof28 Regulates Petal Size in Petunia

Petal size is a key indicator of the ornamental value of plants, such as Petunia hybrida L., which is a popular ornamental species worldwide. Our previous study identified a flower-specific expression pattern of a DNA-binding one finger (Dof)-type transcription factor (TF) PhDof28, in the semi-flowering and full-flowering stages of petunia. In this study, subcellular localization and activation assays showed that PhDof28 was localized in the cell nucleus and could undergo in vitro self-activation. The expression levels of PhDof28 tended to be significantly up-regulated at the top parts of petals during petunia flower opening. Transgenic petunia 'W115' and tobacco plants overexpressing PhDof28 showed similar larger petal phenotypes. The cell sizes at the middle and top parts of transgenic petunia petals were significantly increased, along with higher levels of endogenous indole-3-acetic acid (IAA) hormone. Interestingly, the expression levels of two TFs, PhNAC100 and PhBPEp, which were reported as negative regulators for flower development, were dramatically increased, while the accumulation of jasmonic acid (JA), which induces PhBPEp expression, was also significantly enhanced in the transgenic petals. These results indicated that PhDof28 overexpression could increase petal size by enhancing the synthesis of endogenous IAA in petunias. Moreover, a JA-related feedback regulation mechanism was potentially activated to prevent overgrowth of petals in transgenic plants. This study will not only enhance our knowledge of the Dof TF family, but also provide crucial genetic resources for future improvements of plant ornamental traits.

bHLH transcription factor family identification, phylogeny, and its response to abiotic stress in Chenopodium quinoa

The second-largest transcription factor superfamily in plants is that of the basic helix-loop-helix (bHLH) family, which plays an important complex physiological role in plant growth, tissue development, and environmental adaptation. Systematic research on the Chenopodium quinoa bHLH family will enable a better understanding of this species. Herein, authors used a variety of bioinformatics methods and quantitative Real-Time Polymerase Chain Reaction (qRT-PCR) to explore the evolution and function of the 218 CqbHLH genes identified. A total of 218 CqbHLH transcription factor genes were identified in the whole genome, located on 18 chromosomes. A phylogenetic tree was constructed using the CqbHLH and AtbHLH proteins to determine their homology, and the members were divided into 20 subgroups and one unclustered gene. Authors also analyzed 218 CqbHLH genes, conservative motifs, chromosome diffusion, and gene replication. The author constructed one Neighbor-Joining (NJ) tree and a collinearity analysis map of the bHLH family in C. quinoa and six other plant species to study the evolutionary relationship and homology among multiple species. In addition, the expression levels of 20 CqbHLH members from different subgroups in various tissues, different fruit developmental stages, and six abiotic stresses were analyzed. Authors identified 218 CqbHLH genes and studied their biological functions, providing a basis for better understanding and further studying the bHLH family in quinoa.

Genome-Wide Identification and Characterization of the TLP Gene Family in Phyllostachys edulis and Association with Witches' Broom Disease Resistance in Bamboo

Thaumatin-like proteins (TLPs) are pathogenesis-related proteins with pivotal roles in plant defense mechanisms. In this study, various bioinformatics and RNA-seq methods were used to analyze the biotic and abiotic stress responses of the TLP family in Phyllostachys edulis. Overall, 81 TLP genes were identified in P. edulis; 166 TLPs from four plant species were divided into three groups and ten subclasses, with genetic covariance observed between these species. Subcellular localization in silico studies indicated that TLPs were primarily distributed in the extracellular. Analysis of the upstream sequences of TLPs demonstrated the presence of cis-acting elements related to disease defense, environmental stress, and hormonal responses. Multiple sequence alignment demonstrated that most TLPs possessed five conserved REDDD amino acid sequences with only a few amino acid residue differences. RNA-seq analysis of P. edulis responses to Aciculosporium take, the pathogenic fungus that causes witches' broom disease, showed that P. edulis TLPs (PeTLPs) were expressed in different organs, with the highest expression in buds. PeTLPs responded to both abscisic acid and salicylic acid stress. These PeTLP expression patterns were consistent with their gene and protein structures. Collectively, our findings provide a basis for further comprehensive analyses of the genes related to witches' broom in P. edulis.

Amino acid permease RcAAP1 increases the uptake and phloem translocation of an L-valine-phenazine-1-carboxylic acid conjugate

Amino acid conjugates of pesticides can promote the phloem translocation of parent ingredients, allowing for the reduction of usage, and decreased environmental pollution. Plant transporters play important roles in the uptake and phloem translocation of such amino acid-pesticide conjugates such as L-Val-PCA (L-valine-phenazine-1-carboxylic acid conjugate). However, the effects of an amino acid permease, RcAAP1, on the uptake and phloem mobility of L-Val-PCA are still unclear. Here, the relative expression levels of RcAAP1 were found to be up-regulated 2.7-fold and 2.2-fold by the qRT-PCR after L-Val-PCA treatments of Ricinus cotyledons for 1 h and 3 h, respectively. Subsequently, expression of RcAAP1 in yeast cells increased the L-Val-PCA uptake (0.36 μmol/10⁷ cells), which was 2.1-fold higher than the control (0.17 μmol/10⁷ cells). Pfam analysis suggested RcAAP1 with its 11 transmembrane domains belongs to the amino acid transporter family. Phylogenetic analysis found RcAAP1 to be strongly similar to AAP3 in nine other species. Subcellular localization showed that fusion RcAAP1-eGFP proteins were observed in the plasma membrane of mesophyll cells and phloem cells. Furthermore, overexpression of RcAAP1 for 72 h significantly increased the phloem mobility of L-Val-PCA in Ricinus seedlings, and phloem sap concentration of the conjugate was 1.8-fold higher than the control. Our study suggested that RcAAP1 as carrier was involved in the uptake and phloem translocation of L-Val-PCA, which could lay foundation for the utilization of amino acids and further development of vectorized agrochemicals.

Insight into the metabolic potential and ecological function of a novel Magnetotactic Nitrospirota in coral reef habitat

Magnetotactic bacteria (MTB) within the Nitrospirota phylum play important roles in biogeochemical cycles due to their outstanding ability to biomineralize large amounts of magnetite magnetosomes and intracellular sulfur globules. For several decades, Nitrospirota MTB were believed to only live in freshwater or low-salinity environments. While this group have recently been found in marine sediments, their physiological features and ecological roles have remained unclear. In this study, we combine electron microscopy with genomics to characterize a novel population of Nitrospirota MTB in a coral reef area of the South China Sea. Both phylogenetic and genomic analyses revealed it as representative of a novel genus, named as Candidatus Magnetocorallium paracelense XS-1. The cells of XS-1 are small and vibrioid-shaped, and have bundled chains of bullet-shaped magnetite magnetosomes, sulfur globules, and cytoplasmic vacuole-like structures. Genomic analysis revealed that XS-1 has the potential to respire sulfate and nitrate, and utilize the Wood-Ljungdahl pathway for carbon fixation. XS-1 has versatile metabolic traits that make it different from freshwater Nitrospirota MTB, including Pta-ackA pathway, anaerobic sulfite reduction, and thiosulfate disproportionation. XS-1 also encodes both the cbb₃-type and the aa₃-type cytochrome c oxidases, which may function as respiratory energy-transducing enzymes under high oxygen conditions and anaerobic or microaerophilic conditions, respectively. XS-1 has multiple copies of circadian related genes in response to variability in coral reef habitat. Our results implied that XS-1 has a remarkable plasticity to adapt the environment and can play a beneficial role in coral reef ecosystems.

Expression Profile, Regulatory Network, and Putative Role of microRNAs in the Developmental Process of Asian Honey Bee Larval Guts

MiRNAs, as a kind of key regulators in gene expression, play vital roles in numerous life activities from cellular proliferation and differentiation to development and immunity. However, little is known about the regulatory manner of miRNAs in the development of Asian honey bee (Apis cerana) guts. Here, on basis of our previously gained high-quality transcriptome data, transcriptome-wide identification of miRNAs in the larval guts of Apis cerana cerana was conducted, followed by investigation of the miRNAs' differential expression profile during the gut development. In addition to the regulatory network, the potential function of differentially expressed miRNAs (DEmiRNAs) was further analyzed. In total, 330, 351, and 321 miRNAs were identified in the 4-, 5-, and 6-day-old larval guts, respectively; among these, 257 miRNAs were shared, while 38, 51, and 36 ones were specifically expressed. Sequences of six miRNAs were confirmed by stem-loop RT-PCR and Sanger sequencing. Additionally, in the "Ac4 vs. Ac5" comparison group, there were seven up-regulated and eight down-regulated miRNAs; these DEmiRNAs could target 5041 mRNAs, involving a series of GO terms and KEGG pathways associated with growth and development, such as cellular process, cell part, Wnt, and Hippo. Comparatively, four up-regulated and six down-regulated miRNAs detected in the "Ac5 vs. Ac6" comparison group and the targets were associated with diverse development-related terms and pathways, including cell, organelle, Notch and Wnt. Intriguingly, it was noticed that miR-6001-y presented a continuous up-regulation trend across the developmental process of larval guts, implying that miR-6001-y may be a potential essential modulator in the development process of larval guts. Further investigation indicated that 43 targets in the "Ac4 vs. Ac5" comparison group and 31 targets in the "Ac5 vs. Ac6" comparison group were engaged in several crucial development-associated signaling pathways such as Wnt, Hippo, and Notch. Ultimately, the expression trends of five randomly selected DEmiRNAs were verified using RT-qPCR. These results demonstrated that dynamic expression and structural alteration of miRNAs were accompanied by the development of A. c. cerana larval guts, and DEmiRNAs were likely to participate in the modulation of growth as well as development of larval guts by affecting several critical pathways via regulation of the expression of target genes. Our data offer a basis for elucidating the developmental mechanism underlying Asian honey bee larval guts.

Induction of Hibernation and Changes in Physiological and Metabolic Indices in Pelodiscus sinensis

Pelodiscus sinensis (P. sinensis) is a commonly cultivated turtle species with a habit of hibernation. To study the changes in histone expression and methylation of P. sinensis during hibernation induction, a model was established by artificial induction. Physiological and metabolic indices were measured, and the expression and localization of histone (H1, H2A, H2B, H3, and H4) and methylation-related genes (ASH2L, KMT2A, KMT2E, KDM1A, KDM1B, and KDM5A) were measured by quantitative PCR, immunohistochemistry, and Western blot analysis. The results indicated that the metabolism, antioxidation index, and relative expression of histone methyltransferase were significantly decreased (p < 0.05), whereas the activity and expression of histone demethyltransferase were significantly increased (p < 0.05). Although our results showed significant changes in physiological and gene expression after hibernation induction, we could not confirm that P. sinensis entered deep hibernation. Therefore, for the state after cooling-induced hibernation, cold torpor might be a more accurate description. The results indicate that the P. sinensis can enter cold torpor through artificial induction, and the expression of histones may promote gene transcription. Unlike histones expressed under normal conditions, histone methylation may activate gene transcription during hibernation induction. Western blot analysis revealed that the ASH2L and KDM5A proteins were differentially expressed in the testis at different months (p < 0.05), which may perform a role in regulating gene transcription. The immunohistochemical localization of ASH2L and KDM5A in spermatogonia and spermatozoa suggests that ASH2L and KDM5A may perform a role in mitosis and meiosis. In conclusion, this study is the first to report changes in histone-related genes in reptiles, which provides insight for further studies on the physiological metabolism and histone methylation regulation of P. sinensis during the hibernation induction and hibernation period.

Protein disulfide-isomerase A4 confers glioblastoma angiogenesis promotion capacity and resistance to anti-angiogenic therapy

Introduction

Increasing evidence has revealed the key activity of protein disulfide isomerase A4 (PDIA4) in the endoplasmic reticulum stress (ERS) response. However, the role of PDIA4 in regulating glioblastoma (GBM)-specific pro-angiogenesis is still unknown.

Methods

The expression and prognostic role of PDIA4 were analyzed using a bioinformatics approach and were validated in 32 clinical samples and follow-up data. RNA-sequencing was used to search for PDIA4-associated biological processes in GBM cells, and proteomic mass spectrum (MS) analysis was used to screen for potential PDIA4 substrates. Western blotting, real-time quantitative polymerase chain reaction (RT-qPCR), and enzyme-linked immunosorbent assays (ELISA) were used to measure the levels of the involved factors. Cell migration and tube formation assays determined the pro-angiogenesis activity of PDIA4 in vitro. An intracranial U87 xenograft GBM animal model was constructed to evaluate the pro-angiogenesis role of PDIA4 in vivo.

Results

Aberrant overexpression of PDIA4 was associated with a poor prognosis in patients with GBM, although PDIA4 could also functionally regulate intrinsic GBM secretion of vascular endothelial growth factor-A (VEGF-A) through its active domains of Cys-X-X-Cys (CXXC) oxidoreductase. Functionally, PDIA4 exhibits pro-angiogenesis activity both in vitro and in vivo, and can be upregulated by ERS through transcriptional regulation of X-box binding protein 1 (XBP1). The XBP1/PDIA4/VEGFA axis partially supports the mechanism underlying GBM cell survival under ER stress. Further, GBM cells with higher expression of PDIA4 showed resistance to antiangiogenic therapy in vivo.

Conclusions

Our findings revealed the pro-angiogenesis role of PDIA4 in GBM progression and its potential impact on GBM survival under a harsh microenvironment. Targeting PDIA4 might help to improve the efficacy of antiangiogenic therapy in patients with GBM.

A Putative TRAPα Protein of Microsporidia Nosema bombycis Exhibits Non-Canonical Alternative Polyadenylation in Transcripts

Microsporidia are obligate intracellular eukaryotic parasites that have significantly reduced genomes and that have lost most of their introns. In the current study, we characterized a gene in microsporidia Nosema bombycis, annotated as TRAPα (HNbTRAPα). The homologous of TRAPα are a functional component of ER translocon and facilitates the initiation of protein translocation in a substrate-specific manner, which is conserved in animals but absent from most fungi. The coding sequence of HNbTRAPα consists of 2226 nucleotides, longer than the majority of homologs in microsporidia. A 3′ RACE analysis indicated that there were two mRNA isoforms resulting from non-canonical alternative polyadenylation (APA), and the polyadenylate tail was synthesized after the C951 or C1167 nucleotide, respectively. Indirect immunofluorescence analysis showed two different localization characteristics of HNbTRAPα, which are mainly located around the nuclear throughout the proliferation stage and co-localized with the nuclear in mature spores. This study demonstrated that the post-transcriptional regulation mechanism exists in Microsporidia and expands the mRNA isoform repertoire.

2,5-Hexanedione Affects Ovarian Granulosa Cells in Swine by Regulating the CDKN1A Gene: A Transcriptome Analysis

N-hexane, a common industrial organic solvent, causes multiple organ damage owing to its metabolite, 2,5-hexanedione (2,5-HD). To identify and evaluate the effects of 2,5-HD on sows’ reproductive performance, we used porcine ovarian granulosa cells (pGCs) as a vehicle and carried out cell morphology and transcriptome analyses. 2,5-HD has the potential to inhibit the proliferation of pGCs and induce morphological changes and apoptosis depending on the dose. RNA-seq analyses identified 4817 differentially expressed genes (DEGs), with 2394 down-regulated and 2423 up-regulated following 2,5-HD exposure treatment. The DEG, cyclin-dependent kinase inhibitor 1A (CDKN1A), according to the Kyoto Encyclopedia of Genes and Genomes enrichment analysis, was significantly enriched in the p53 signaling pathway. Thus, we evaluated its function in pGC apoptosis in vitro. Then, we knocked down the CDKN1A gene in the pGCs to identify its effects on pGCs. Its knockdown decreased pGC apoptosis, with significantly fewer cells in the G1 phase (p < 0.05) and very significantly more cells in the S phase (p < 0.01). Herein, we revealed novel candidate genes that influence pGCs apoptosis and cell cycle and provided new insights into the role of CDKN1A in pGCs during apoptosis and cell cycle arrest.

Identifying Genes Associated with Female Flower Development of Phellodendron amurense Rupr. Using a Transcriptomics Approach

Phellodendron amurense Rupr., a species of Rutaceae, is a nationally protected and valuable medicinal plant. It is generally considered to be dioecious. With the discovery of monoecious P. amurense, the phenomenon that its sex development is regulated by epigenetics has been revealed, but the way epigenetics affects the sex differentiation of P. amurense is still unclear. In this study, we investigated the effect of DNA methylation on the sexual development of P. amurense. The young inflorescences of male plants were treated with the demethylation agent 5-azaC, and the induced female flowers were obtained. The induced female flowers’ morphological functions and transcriptome levels were close to those of normally developed plants. Genes associated with the development of female flowers were studied by comparing the differences in transcriptome levels between the male and female flowers. Referring to sex-related genes reported in other plants, 188 candidate genes related to the development of female flowers were obtained, including sex-regulating genes, genes related to the formation and development of sexual organs, genes related to biochemical pathways, and hormone-related genes. RPP0W, PAL3, MCM2, MCM6, SUP, PIN1, AINTEGUMENTA, AINTEGUMENTA-LIKE6, AGL11, SEUSS, SHI-RELATED SEQUENCE 5, and ESR2 were preliminarily considered the key genes for female flower development. This study has demonstrated that epigenetics was involved in the sex regulation of P. amurense, with DNA methylation as one of its regulatory modes. Moreover, some candidate genes related to the sexual differentiation of P. amurense were obtained with analysis. These results are of great significance for further exploring the mechanism of sex differentiation of P. amurense and studying of sex differentiation of plants.

Complete Mitochondrial Genomes of Four Pelodiscus sinensis Strains and Comparison with Other Trionychidae Species

The Chinese soft-shelled turtle (Pelodiscus sinensis) is an important aquaculture reptile with rich nutritional and medicinal values. In recent decades, the wild resources of P. sinensis have been depleting due to natural and artificial factors. Herein, we report the complete mitochondrial genome of four P. sinensis strains, including the Japanese (RB) strain, Qingxi Huabie (HB) strain, Jiangxi (JB) strain, and Qingxi Wubie (WB) strain. The nucleotide composition within the complete mitogenomes was biased towards A + T with a variable frequency ranging from 59.28% (cox3) to 70.31% (atp8). The mitogenomes of all four strains contained 13 protein-coding genes (PCGs), 22 tRNAs, 2 rRNAs, 1 control region, and a replication origin region of the L-strand replication (OL), which was consistent with most vertebrates. Additionally, the atp8, nad4l, nad6, and nad3 genes possessed high genetic variation and can be used as potential markers for the identification of these P. sinensis strains. Additionally, all PCGs genes were evolving primarily under purifying selection. Through comparative analysis, it was revealed that most of the tRNAs were structurally different in the TψC stem, DHU stem, and acceptor stem. The length of the tandem repeats in the control region was variable in the four P. sinensis strains, ranging from 2 bp to 50 bp. Phylogenetic analysis indicated that all P. sinensis strains clustered into one branch and were closely related to other Trionychinae species. Overall, this study provides mitochondrial genome information for different P. sinensis strains to support further species identification and germplasm resource conservation.

Identification of Trehalose-6-Phosphate Synthase (TPS) Genes Associated with Both Source-/Sink-Related Yield Traits and Drought Response in Rapeseed (Brassica napus L.)

Trehalose-6-phosphate synthase (TPS) is an important enzyme for the synthesis of Trehalose-6-phosphate (T6P). In addition to being a signaling regulator of carbon allocation that improves crop yields, T6P also plays essential roles in desiccation tolerance. However, comprehensive studies, such as evolutionary analysis, expression analysis, and functional classification of the TPS family in rapeseed (Brassica napus L.) are lacking. Here, we identified 35 BnTPSs, 14 BoTPSs, and 17 BrTPSs in cruciferous plants, which were classified into three subfamilies. Phylogenetic and syntenic analysis of TPS genes in four cruciferous species indicated that only gene elimination occurred during their evolution. Combined phylogenetic, protein property, and expression analysis of the 35 BnTPSs suggested that changes in gene structures might have led to changes in their expression profiles and further functional differentiation during their evolution. In addition, we analyzed one set of transcriptome data from Zhongshuang11 (ZS11) and two sets of data from extreme materials associated with source-/sink-related yield traits and the drought response. The expression levels of four BnTPSs (BnTPS6, BnTPS8, BnTPS9, and BnTPS11) increased sharply after drought stress, and three differentially expressed genes (BnTPS1, BnTPS5, and BnTPS9) exhibited variable expression patterns among source and sink tissues between yield-related materials. Our findings provide a reference for fundamental studies of TPSs in rapeseed and a framework for future functional research of the roles of BnTPSs in both yield and drought resistance.

Fine Mapping of a Pleiotropic Locus (BnUD1) Responsible for the Up-Curling Leaves and Downward-Pointing Siliques in Brassica napus

Leaves and siliques are important organs associated with dry matter biosynthesis and vegetable oil accumulation in plants. We identified and characterized a novel locus controlling leaf and silique development using the Brassica napus mutant Bnud1, which has downward-pointing siliques and up-curling leaves. The inheritance analysis showed that the up-curling leaf and downward-pointing silique traits are controlled by one dominant locus (BnUD1) in populations derived from NJAU5773 and Zhongshuang 11. The BnUD1 locus was initially mapped to a 3.99 Mb interval on the A05 chromosome with a BC₆F₂ population by a bulked segregant analysis-sequencing approach. To more precisely map BnUD1, 103 InDel primer pairs uniformly covering the mapping interval and the BC₅F₃ and BC₆F₂ populations consisting of 1042 individuals were used to narrow the mapping interval to a 54.84 kb region. The mapping interval included 11 annotated genes. The bioinformatic analysis and gene sequencing data suggested that BnaA05G0157900ZS and BnaA05G0158100ZS may be responsible for the mutant traits. Protein sequence analyses showed that the mutations in the candidate gene BnaA05G0157900ZS altered the encoded PME in the trans-membrane region (G45A), the PMEI domain (G122S), and the pectinesterase domain (G394D). In addition, a 573 bp insertion was detected in the pectinesterase domain of the BnaA05G0157900ZS gene in the Bnud1 mutant. Other primary experiments indicated that the locus responsible for the downward-pointing siliques and up-curling leaves negatively affected the plant height and 1000-seed weight, but it significantly increased the seeds per silique and positively affected photosynthetic efficiency to some extent. Furthermore, plants carrying the BnUD1 locus were compact, implying they may be useful for increasing B. napus planting density. The findings of this study provide an important foundation for future research on the genetic mechanism regulating the dicotyledonous plant growth status, and the Bnud1 plants can be used directly in breeding.

Genome-Wide Identification and Posttranscriptional Regulation Analyses Elucidate Roles of Key Argonautes and Their miRNA Triggers in Regulating Complex Yield Traits in Rapeseed

Argonautes (AGOs) interact with microRNAs (miRNAs) to form the RNA-induced silencing complex (RISC), which can posttranscriptionally regulate the expression of targeted genes. To date, however, the AGOs and their miRNA triggers remain elusive in rapeseed (Brassica napus). Here, we systematically performed a phylogenetic analysis and examined the collinear relationships of the AGOs among four Brassicaceae species. Their physicochemical properties, gene structures, and expression patterns among 81 tissues from multiple materials and developmental stages were further analyzed. Additionally, their posttranscriptional regulation was analyzed using psRNATarget prediction, miRNA-/mRNA-Seq analyses, and a qRT-PCR verification. We finally identified 10 AtAGOs, 13 BolAGOs, 11 BraAGOs, and 24 BnaAGOs. An expression analysis of the BnaAGOs in the B. napus cultivar ZS11, as well as genotypes with extreme phenotypes in various yield-related traits, revealed the conservation and diversity of these genes. Furthermore, we speculated the posttranscriptional regulation of the B. napus miR168a-AGO1s and miR403-AGO2s modules. Combining miRNA-Seq and mRNA-Seq analyses, we found that the B. napus miR168a-AGO1s module may play an essential role in negatively regulating yield traits, whereas the miR403-AGO2s module positively impacts yield. This is the first attempt to comprehensively analyze the AGOs and their miRNA triggers in B. napus and provides a theoretical basis for breeding high-yielding varieties through the manipulation of the miRNA-AGOs modules.

Unveiling the circRNA-Mediated Immune Responses of Western Honey Bee Larvae to Ascosphaera apis Invasion

Western honey bee (Apis mellifera), a eusocial insect with a superior economic and ecological value, is widely used in the beekeeping industry throughout the world. As a new class of non-coding RNAs (ncRNAs), circular RNAs (circRNAs) participate in the modulation of considerable biological processes, such as the immune response via diverse manners. Here, the identification, characteristic investigation, and molecular verification of circRNAs in the Apis mellifera ligustica larval guts were conducted, and the expression pattern of larval circRNAs during the Ascosphaera apis infection was analyzed, followed by the exploration of the potential regulatory part of differentially expressed circRNAs (DEcircRNAs) in host immune responses. A total of 2083 circRNAs in the larval guts of A. m. ligustcia were identified, with a length distribution ranging from 106 nt to 92,798 nt. Among these, exonic circRNAs were the most abundant type and LG1 was the most distributed chromosome. Additionally, 25, 14, and 30 up-regulated circRNAs as well as 26, 25, and 62 down-regulated ones were identified in the A. apis-inoculated 4-, 5-, and 6-day-old larval guts in comparison with the corresponding un-inoculated larval guts. These DEcircRNAs were predicted to target 35, 70, and 129 source genes, which were relative to 12, 23, and 20 GO terms as well as 11, 10, and 27 KEGG pathways, including 5 cellular and humoral immune pathways containing apoptosis, autophagy, endocytosis, MAPK, Toll, and Imd signaling pathways. Furthermore, complex competing endogenous RNA (ceRNA) regulatory networks were detected to be formed among DEcircRNAs, DEmiRNAs, and DEmRNAs. The Target DEmRNAs were engaged in 24, 20, and 25 functional terms as well as 62, 80, and 159 pathways, including several vital immune defense-associated pathways, namely the lysosome, endocytosis, phagosome, autophagy, apoptosis, MAPK, Jak-STAT, Toll, and Imd signaling pathways. Finally, back-splicing sites within 15 circRNAs and the difference in the 9 DEcircRNAs' expression between un-inoculated and A. apis-inoculated larval guts were confirmed utilizing molecular methods. These findings not only enrich our understanding of bee host-fungal pathogen interactions, but also lay a foundation for illuminating the mechanism underlying the DEcircRNA-mediated immune defense of A. m. ligustica larvae against A. apis invasion.

Identification of two tandem genes associated with primary rosette branching in flowering Chinese cabbage

Branching is an important agronomic trait determining plant architecture and yield; however, the molecular mechanisms underlying branching in the stalk vegetable, flowering Chinese cabbage, remain unclear. The present study identified two tandem genes responsible for primary rosette branching in flowering Chinese cabbage by GradedPool-Seq (GPS) combined with Kompetitive Allele Specific PCR (KASP) genotyping. A 900 kb candidate region was mapped in the 28.0-28.9 Mb interval of chromosome A07 through whole-genome sequencing of three graded-pool samples from the F₂ population derived by crossing the branching and non-branching lines. KASP genotyping narrowed the candidate region to 24.6 kb. Two tandem genes, BraA07g041560.3C and BraA07g041570.3C, homologous to AT1G78440 encoding GA2ox1 oxidase, were identified as the candidate genes. The BraA07g041560.3C sequence was identical between the branching and non-branching lines, but BraA07g041570.3C had a synonymous single nucleotide polymorphic (SNP) mutation in the first exon (290^th bp, A to G). In addition, an ERE cis-regulatory element was absent in the promoter of BraA07g041560.3C, and an MYB cis-regulatory element in the promoter of BraA07g041570.3C in the branching line. Gibberellic acid (GA₃) treatment decreased the primary rosette branch number in the branching line, indicating the significant role of GA in regulating branching in flowering Chinese cabbage. These results provide valuable information for revealing the regulatory mechanisms of branching and contributing to the breeding programs of developing high-yielding species in flowering Chinese cabbage.

Comprehensive investigation and regulatory function of lncRNAs engaged in western honey bee larval immune response to Ascosphaera apis invasion

Ascosphaera apis is a fungal pathogen that exclusively infects bee larvae, causing chalkbrood disease, which results in severe damage for beekeeping industry. Long non-coding RNAs (lncRNAs) are versatile regulators in various biological processes such as immune defense and host-pathogen interaction. However, expression pattern and regulatory role of lncRNAs involved in immune response of bee host to A. apis invasion is still very limited. Here, the gut tissues of Apis mellifera ligustica 4-, 5-, and 6-day-old larvae inoculated by A. apis spores (AmT1, AmT2, and AmT3 groups) and corresponding un-inoculated larval guts (AmCK1, AmCK2, and AmCK3 groups) were prepared and subjected to deep sequencing, followed by identification of lncRNAs, analysis of differentially expressed lncRNAs (DElncRNAs), and investigation of competing endogenous RNA (ceRNA) network. In total, 3,746 A. m. ligustica lncRNAs were identified, including 78 sense lncRNAs, 891 antisense lncRNAs, 1,893 intergenic lncRNAs, 346 bidirectional lncRNAs, and 210 intronic lncRNAs. In the 4-, 5-, and 6- comparison groups, 357, 236, and 505 DElncRNAs were discovered. Additionally, 217, 129, and 272 DElncRNAs were respectively predicted to regulate neighboring genes via cis-acting manner, and these targets were associated with a series of GO terms and KEGG pathways of great importance, such as response to stimulus and Jak-STAT signaling pathway. Moreover, 197, 95, and 356 DElncRNAs were observed to target 10, eight, and 21 DEmiRNAs and further target 147, 79, and 315 DEmRNAs, forming complex regulatory networks. Further investigation suggested that these targets were engaged in several key cellular and humoral immune pathways, such as phagosome and MAPK signaling pathway. Ultimately, the expression trends of nine randomly selected DElncRNAs were verified by RT-qPCR, confirming the authenticity and reliability of our transcriptome data. Findings in this current work not only provide candidate DElncRNAs for functional study, but also lay a foundation for unclosing the mechanism underlying DElncRNA-regulated larval immune responses to A. apis invasion.

Comparative analysis of mitochondrial genomes of Broussonetia spp. (Moraceae) reveals heterogeneity in structure, synteny, intercellular gene transfer, and RNA editing

The genus Broussonetia (Moraceae) is comprised of three non-hybrid recognized species that all produce high quality fiber essential in the development of papermaking and barkcloth-making technology. In addition, these species also have medicinal value in several countries. Despite their important economical, medicinal, and ecological values, the complete mitogenome of Broussonetia has not been reported and investigated, which would greatly facilitate molecular phylogenetics, species identification and understanding evolutionary processes. Here, we assembled the first-reported three complete Broussonetia (B. papyrifera, B. kaempferi, and B. monoica) mitochondrial genomes (mitogenome) based on a hybrid strategy using Illumina and Oxford Nanopore Technology sequencing data, and performed comprehensive comparisons in terms of their structure, gene content, synteny, intercellular gene transfer, phylogeny, and RNA editing. Our results showed their huge heterogeneities among the three species. Interestingly, the mitogenomes of B. monoica and B. papyrifera consisted of a single circular structure, whereas the B. kaempferi mitogenome was unique and consisted of a double circular structure. Gene content was consistent except for a few transfer RNA (tRNA) genes. The Broussonetia spp. mitogenomes had high sequence conservation but B. monoica with B. kaempferi contained more synteny blocks and were more related, a finding that was well-supported in organellar phylogeny. Fragments that had been transferred between mitogenomes were detected at plastome hotspots that had integrated under potential mediation of tRNA genes. In addition, RNA editing sites showed great differences in abundance, type, location and efficiency among species and tissues. The availability of these complete gap-free mitogenomes of Broussonetia spp. will provide a valuable genetic resource for evolutionary research and understanding the communications between the two organelle genomes.

Comprehensive Molecular Characterization of the Mitochondrial Genome of the Takin Lungworm Varestrongylus eleguneniensis (Strongylida: Protostrongylidae)

The takin lungworm Varestrongylus eleguneniensis (Strongylida: Protostrongylidae) causes lethal bronchopneumonia and represents severe threats to captive and wild populations. However, until now there has been very limited information available concerning the molecular epidemiology and evolutionary biology of V. eleguneniensis. Mitochondrial genomes (mtDNAs) can provide resources for investigations in these areas and, therefore, can assist with the surveillance and control of this lungworm. Herein, the complete mtDNA of V. eleguneniensis was sequenced and characterized with Illumina pipeline analyses. This circular genome (13,625 bp) encoded twelve protein-coding genes (PCGs), two rRNAs, and twenty-two tRNAs, with notable levels of AT and GC skews. Comparative genomics revealed a purifying selection among PCGs, with cox1 and nad6 having the lowest and the highest evolutionary rate, respectively. Genome-wide phylogenies showed a close relationship between V. eleguneniensis and Protostrongylus rufescens in Strongylida. Single gene (PCGs or rRNAs)-based phylogenies indicated that cox1 and nad5 genes shared the same family-level topology with that inferred from genomic datasets, suggesting that both genes could be suitable genetic markers for evolutionary and phylogenetic studies of Strongylida species. This was the first mtDNA of any member of the genus Varestrongylus, and its comprehensive molecular characterization represents a new resource for systematic, population genetic and evolutionary biological studies of Varestrongylus lungworms in wildlife.

A small chromosomal inversion mediated by MITE transposons confers cleistogamy in Brassica napus

Cleistogamy, self-pollination within closed flowers, can help maintain seed purity, accelerate breeding speed, and aid in the development of ornamental flowers. However, the mechanism underlying petal closing/opening behavior remains elusive. Here, we found that a Brassica napus petal closing/opening behavior was inherited in a Mendelian manner. Fine mapping and positional cloning experiments revealed that the Mendelian factor originated from a short (29.8 kb) inversion mediated by BnDTH9 miniature inverted-repeat transposable elements (MITEs) on chromosome C03. This inversion led to tissue-specific gene promoter exchange between BnaC03.FBA (BnaC03G0156800ZS encoding an F-Box-associated domain-containing protein) and BnaC03.EFO1 (BnaC03G0157400ZS encoding an EARLY FLOWERING BY OVEREXPRESSION 1 protein) positioned near the respective inversion breakpoints. Our genetic transformation work demonstrated that the cleistogamy originated from high tissue-specific expression of the BnaC03.FBA gene caused by promoter changes due to the MITE-mediated inversion. BnaC03.FBA is involved in the formation of an SCF (Skp1-Cullin-F-box) complex, which participates in ubiquitin-mediated protein targeting for degradation through the ubiquitin 26S-proteasome system. Our results shed light on a molecular model of petal-closing behavior.

Development of InDel markers for Oryza sativa ssp. javanica based on whole-genome resequencing

Oryza sativa ssp. javanica rice varieties exhibit a wide variation in the phenotypes of several important agronomic traits, including grain quality, grain shape, plant architecture, disease resistance, and high adaption to an unfavorable environment, indicating a great potential for rice improvement. DNA molecular markers are basic and critical tools in genetic analysis and gene mining. However, only a few whole-genome variation analyses have been performed in Oryza sativa ssp. Javanica (tropical japonica rice), and this has hampered the utilization of such an important resource. In this study, the length of insertions/deletions variation greater larger than 10 bp from 10 Oryza sativa ssp. indica rice and 10 Oryza sativa ssp. tropical japonica rice were extracted by using the Nipponbare genome as a reference. A total of 118 primer pairs which were almost evenly distributed on each chromosome corresponding to the loci of InDels were designed by the Primer 5 program. We confirmed 85 InDel markers from 60 rice varieties, including indica and tropical japonica, by running polyacrylamide gels. The InDel markers function like SSRs in identifying hybrids, calculating genetic distance, constructing the genetic linkage map, and gene mining. The InDel markers developed in this study might help in genetic studies and to investigate the tropical japonica rice varieties.

Functional characterization of the Serine acetyltransferase family genes uncovers the diversification and conservation of cysteine biosynthesis in tomato

Sulfur-containing compounds are essential for plant development and environmental adaptation, and closely related to the flavor and nutrition of the agricultural products. Cysteine, the first organic sulfur-containing molecule generated in plants, is the precursor for most of these active substances. Serine acetyltransferase (SERAT) catalyzes the rate-limiting step of its formation. However, despite their importance, systematic analyses of these enzymes in individual species, especially in economically important crops, are still limited. Here, The SERAT members (SlSERATs, four in total) were identified and characterized in tomato. Phylogenetically, the four SlSERAT proteins were classified into three subgroups with distinct genomic structures and subcellular localizations. On the function, it was interesting to find that SlSERAT3;1, possessed a high ability to catalyze the formation of OAS, even though it contained a long C-terminus. However, it retained the essential C-terminal Ile, which seems to be a characteristic feature of SERAT3 subfamily members in Solanaceae. Besides, SlSERAT1;1 and SlSERAT2;2 also had high activity levels and their catalyzing abilities were significantly improved by the addition of an OAS-(thiol)-lyase protein. At the transcriptional level, the four SlSERAT genes had distinct expression patterns during tomato plant development. Under abiotic stress conditions, the chloroplast-localized SlSERATs were the main responders, and the SlSERATs adopted different strategies to cope with osmotic, ion toxicity and other stresses. Finally, analyses in the loss-of-function and overexpression lines of SlSERAT1;1 suggested that function redundancy existed in the tomato SERAT members, and the tomato SERAT member was ideal target for S-assimilation manipulating in molecular breeding.

Molecular beacon based real-time PCR p1 gene genotyping, macrolide resistance mutation detection and clinical characteristics analysis of Mycoplasma pneumoniae infections in children

Background

Mycoplasma pneumoniae can be divided into different subtypes on the basis of the sequence differences of adhesive protein P1, but the relationship between different subtypes, macrolide resistance and clinical manifestations are still unclear. In the present study, we established a molecular beacon based real-time polymerase chain reaction (real-time PCR) p1 gene genotyping method, analyzed the macrolide resistance gene mutations and the relationship of clinical characteristics with the genotypes.

Methods

A molecular beacon based real-time PCR p1 gene genotyping method was established, the mutation sites of macrolide resistance genes were analyzed by PCR and sequenced, and the relationship of clinical characteristics with the genotypes was analyzed.

Results

The detection limit was 1–100 copies/reaction. No cross-reactivity was observed in the two subtypes. In total, samples from 100 patients with positive M. pneumoniae detection results in 2019 and 2021 were genotyped using the beacon based real-time PCR method and P1-1 M. pneumoniae accounted for 69.0%. All the patients had the A2063G mutation in the macrolide resistance related 23S rRNA gene. Novel mutations were also found, which were C2622T, C2150A, C2202G and C2443A mutations. The relationship between p1 gene genotyping and the clinical characteristics were not statistically related.

Conclusion

A rapid and easy clinical application molecular beacon based real-time PCR genotyping method targeting the p1 gene was established. A shift from type 1 to type 2 was found and 100.0% macrolide resistance was detected. Our study provided an efficient method for genotyping M. pneumoniae, valuable epidemiological monitoring information and clinical treatment guidance to control high macrolide resistance.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12879-022-07715-6.

Development of a Real-Time Quantitative PCR Assay for the Specific Detection of Bacillus velezensis and Its Application in the Study of Colonization Ability

Bacillus velezensis is a widely used biocontrol agent closely related to B. amyloliquefaciens, and the two species cannot be distinguished by universal primers that are currently available. The study aimed to establish a rapid, specific detection approach for B. velezensis. Many unique gene sequences of B. velezensis were selected through whole genome sequence alignment of B. velezensis strains and were used to design a series of forward and reverse primers, which were then screened by PCR and qPCR using different Bacillus samples as templates. The colonization ability of B. velezensis ZF2 in different soils and different soil environmental conditions was measured by qPCR and a 10-fold dilution plating assay. A specific primer pair targeting the sequence of the D3N19_RS13500 gene of B. velezensis ZF2 was screened and could successfully distinguish B. velezensis from B. amyloliquefaciens. A rapid specific real-time qPCR detection system for B. velezensis was established. B. velezensis ZF2 had a very strong colonization ability in desert soil, and the optimal soil pH was 7-8. Moreover, the colonization ability of strain ZF2 was significantly enhanced when organic matter from different nitrogen sources was added to the substrate. This study will provide assistance for rapid specificity detection and biocontrol application of B. velezensis strains.

Accelerated evolution of Vkorc1 in desert rodent species reveals genetic preadaptation to anticoagulant rodenticides

BACKGROUND

Some rodent species living in arid areas show elevated physiological tolerance to anti-vitamin K rodenticides (AVKs), which seems to be due to some unknown selective pressures that rodents may experience in desert habitats. Genes involved in the ϒ-carboxylation of blood coagulation, including vitamin K epoxide reductase complex, subunit 1 (Vkorc1), ϒ-glutamyl-carboxylase (Ggcx) and NAD(P)H quinone one dehydrogenase (Nqo1) are associated with anticoagulant resistance, or some levels of elevated tolerance, in rodents. To detect whether the DNA sequences of the three genes are also under natural selection in the desert rodent species, we analyzed the Vkorc1, Ggcx and Nqo1 genes of the desert rodents and compared them with other rodent species.

RESULTS

We found an accelerated evolutionary rate in Vkorc1 of desert rodents, especially in Mus spretus, Nannospalax galili and Psammomys obesus. By contrast, signals of positive selection were absent for Ggcx and Nqo1 in all species. Mapping the amino acid variations on the VKORC1 protein three-dimensional model suggested most interspecific amino acid variations occur on the outer surface of the VKORC1 pocket, whereas most intraspecific amino acid changes and known AVK resistance mutations occurred on the inner surface and endoplasmic reticulum luminal loop regions. Some desert-species-specific amino acid variations were found on the positions where known resistance mutations occurred, indicating these variations might be related to the elevated physical tolerance to AVKs in desert rodents.

CONCLUSION

The evolution of Vkorc1 has been accelerated in some desert rodent species, indicating genetic preadaptation to anticoagulant rodenticides. Positive selection and relaxed selection have been detected in Psammomys obesus and Nannospalax galili, indicating the two rodent species might also show tolerance to AVKs, which needs further verification. © 2022 Society of Chemical Industry.

A Method for Degenerate Primer Design Based on Artificial Bee Colony Algorithm

Aiming to address the complex degenerate primer design problem in the biological field, in this paper, we design a degenerate primer optimization model considering primer coverage and degeneracy that allows a small number of base mismatches, and propose a global optimization method based on the artificial bee colony algorithm. The proposed algorithm combines the idea of the ant colony algorithm with the optimization process of the artificial bee colony algorithm, overcomes the disadvantage of the uncertain candidate solution length of the artificial bee colony algorithm in solving discrete optimization problems, designs the search space model according to the construction process of candidate solution in ant colony optimization algorithm, and redesigns various bee foraging strategies according to the optimization process information. In the comparative experiments on DNA template sequences of different scales, the degenerate primer designed by the proposed algorithm is superior to the existing methods in terms of stability, specificity, coverage and degeneracy.