Safety assessment of a novel marine multi-stress-tolerant yeast Meyerozyma guilliermondii GXDK6 according to phenotype and whole genome-sequencing analysis

The application of microorganisms as probiotics is limited due to lack of safety evaluation.Here,a novel multi-stress-tolerant yeast Meyerozyma guilliermondii GXDK6 with aroma-producing properties was identified from marine mangrove microorganisms.Its safety and probiotic properties were assessed in accordance with phenotype and whole-genome sequencing analysis.Results showed that the genes and phenotypic expression of related virulence,antibiotic resistance and retroelement were rarely found.Hyphal morphogenesis genes(SIT4,HOG1,SPA2,ERK1,ICL1,CST20,HSP104,TPS1,and RHO1)and phospholipase secretion gene(VPS4)were annotated.True hyphae and phospholipase were absent.Only one retroelement(Tad1-65_BG)was found.Major biogenic amines(BAs)encoding genes were absent,except for spermidine synthase(JA9_002594),spermine synthase(JA9_004690),and tyrosine decarboxylase(inx).The production of single BAs and total BAs was far below the food-defined thresholds.GXDK6 had no resistance to common antifungal drugs.Virulence enzymes,such as gelatinase,DNase,hemolytic,lecithinase,and thrombin were absent.Acute toxicity test with mice demonstrated that GXDK6 is safe.GXDK6 has a good reproduction ability in the simulation gastrointestinal tract.GXDK6 also has a strong antioxidant ability,β-glucosidase,and inulinase activity.To sum up,GXDK6 is considered as a safe probiotic for human consumption and food fermentation.

To Analyze the Sensitivity of RT-PCR Assays Employing S Gene Target Failure with Whole Genome Sequencing Data during Third Wave by SARS-CoV-2 Omicron Variant

Introduction: Omicron is a highly divergent variant of concern (VOCs) of a severe acute respiratory syndrome SARS-CoV-2. It carries a high number of mutations in its spike protein hence;it is more transmissible in the community by immune evasion mechanisms. Due to mutation within S gene, most Omicron variants have reported S gene target failure (SGTF) with some commercially available PCR kits. Such diagnostic features can be used as markers to screen Omicron. However, Whole Genome Sequencing (WGS) is the only gold standard approach to confirm novel microorganisms at genetically level as similar mutations can also be found in other variants that are circulating at low frequencies worldwide. This Retrospective study is aimed to assess RT-PCR sensitivity in the detection of S gene target failure in comparison with whole genome sequencing to detect variants of Omicron. Methods: We have analysed retrospective data of SARS-CoV-2 positive RT-PCR samples for S gene target failure (SGTF) with TaqPath COVID-19 RT-PCR Combo Kit (ThermoFisher) and combined with sequencing technologies to study the emerged pattern of SARS-CoV-2 variants during third wave at the tertiary care centre, Surat. Results: From the first day of December 2021 till the end of February 2022, a total of 321,803 diagnostic RT-PCR tests for SARS-CoV-2 were performed, of which 20,566 positive cases were reported at our tertiary care centre with an average cumulative positivity of 6.39% over a period of three months. In the month of December 21 samples characterized by the SGTF (70/129) were suggestive of being infected by the Omicron variant and identified as Omicron (B.1.1.529 lineage) when sequence. In the month of January, we analysed a subset of samples (n = 618) with SGTF (24%) and without SGTF (76%) with Ct values Conclusions: During the COVID-19 pandemic, it took almost more than 15 days to diagnose infection and identify pathogen by sequencing technology. In contrast to that molecular assay provided quick identification with the help of SGTF phenomenon within 5 hours of duration. This strategy helps scientists and health policymakers for the quick isolation and identification of clusters. That ultimately results in a decreased transmission of pathogen among the community.

Gene characterization and phylogenetic analysis of four mitochondrial genomes in Caenogastropoda

Caenogastropoda is a highly diverse group,containing~60%of all existing gastropods.Species in this subclass predominantly inhabit marine environments and have a high ecological and economic value.Owing to the increase in relevant phylogenetic studies,our understanding of between species relatedness in Caenogastropoda has improved.However,the biodiversity,taxonomic status,and phylogenetic relationships of this group remain unclear.In the present study,we performed next-generation sequencing of four complete mitochondrial genomes from three families(Buccinidae,Columbellidae,and Cypraeidae)and the four mitogenomes were classical circular structures,with a length of 16177 bp in Volutharpa ampullacea,16244 bp in Mitrella albuginosa,16926bp in Mauritia arabica asiatica and 15422 bp in Erronea errones.Base composition analysis indicated that whole sequences were biased toward A and T.Then compared them with 171 complete mitochondrial genomes of Caenogastropoda.The phylogenetic relationship of Caenogastropoda derived from Maximum Likelihood(ML)and Bayesian Inference(BI)trees constructed based on CDS sequences was consistent with the results of traditional morphological analysis,with all three families showing close relationships.This study supported Caenogastropoda at the molecular level as a separate clade of Mollusca.According to our divergence time estimations,Caenogastropoda was formed during the middle Triassic period(~247.2–237 Ma).Our novel mitochondrial genomes provide evidence for the speciation of Caenogastropoda in addition to elucidating the mitochondrial genomic evolution of this subclass.

Action of circulating and infiltrating B cells in the immune microenvironment of colorectal cancer by single-cell sequencing analysis

BACKGROUND The complexity of the immune microenvironment has an impact on the treatment of colorectal cancer(CRC),one of the most prevalent malignancies worldwide.In this study,multi-omics and single-cell sequencing techniques were used to investigate the mechanism of action of circulating and infiltrating B cells in CRC.By revealing the heterogeneity and functional differences of B cells in cancer immunity,we aim to deepen our understanding of immune regulation and provide a scientific basis for the development of more effective cancer treatment strategies.AIM To explore the role of circulating and infiltrating B cell subsets in the immune microenvironment of CRC,explore the potential driving mechanism of B cell development,analyze the interaction between B cells and other immune cells in the immune microenvironment and the functions of communication molecules,and search for possible regulatory pathways to promote the anti-tumor effects of B cells.METHODS A total of 69 paracancer(normal),tumor and peripheral blood samples were collected from 23 patients with CRC from The Cancer Genome Atlas database(https://portal.gdc.cancer.gov/).After the immune cells were sorted by multicolor flow cytometry,the single cell transcriptome and B cell receptor group library were sequenced using the 10X Genomics platform,and the data were analyzed using bioinformatics tools such as Seurat.The differences in the number and function of B cell infiltration between tumor and normal tissue,the interaction between B cell subsets and T cells and myeloid cell subsets,and the transcription factor regulatory network of B cell subsets were explored and analyzed.RESULTS Compared with normal tissue,the infiltrating number of CD20+B cell subsets in tumor tissue increased significantly.Among them,germinal center B cells(GCB)played the most prominent role,with positive clone expansion and heavy chain mutation level increasing,and the trend of differentiation into memory B cells increased.However,the number of plasma cells in the tumor microenvironment decreased significantly,and the plasma cells secreting IgA antibodies decreased most obviously.In addition,compared with the immune microenvironment of normal tissues,GCB cells in tumor tissues became more closely connected with other immune cells such as T cells,and communication molecules that positively regulate immune function were significantly enriched.CONCLUSION The role of GCB in CRC tumor microenvironment is greatly enhanced,and its affinity to tumor antigen is enhanced by its significantly increased heavy chain mutation level.Meanwhile,GCB has enhanced its association with immune cells in the microenvironment,which plays a positive anti-tumor effect.

Safety evaluation and whole genome sequencing for revealing the ability of Penicillium oxalicum WX-209 to safely and effectively degrade citrus segments

The microbial potential of Penicillium has received critical attention.The present research aimed to elucidate the efficacy of crude enzyme secreted from Penicillium oxalicum WX-209 in degrading citrus segments and evaluate the safety of the process.Results showed that citrus segment membranes gradually dissolved after treatment with the crude enzyme solution,indicating good degradation capability.No significant differences in body weight,food ingestion rate,hematology,blood biochemistry,and weight changes of different organs were found between the enzyme intake and control groups.Serial experiments showed that the crude enzyme had high biological safety.Moreover,the whole genome of P.oxalicum WX-209 was sequenced by PacBio and Illumina platforms.Twenty-five scaffolds were assembled to generate 36 Mbp size of genome sequence comprising 11369 predicted genes modeled with a GC content of 48.33%.A total of 592 genes were annotated to encode enzymes related to carbohydrates,and some degradation enzyme genes were identified in strain P.oxalicum WX-209.

Copy number variation sequencing for diagnosis of cytomegalovirus infection based low-depth whole-genome sequencing technology in fetus:Three cases and literature review

Objective:To summarize the application value of copy number variant sequencing(CNV-seq)in the detection of fetal chromosome and cytomegalovirus load.Methods:The study analyzed the clinical basic data,relevant laboratory tests,treatment process,and outcomes of three patients with positive cytomegalovirus load detected by CNV-seq for fetal chromosomes and cytomegalovirus load,and literature review was done simutaneoubly.Results:In all three cases,the amniotic fluid cytomegalovirus load was less than 105 Copies/ml,and there were no significant neurological abnormalities observed during pregnancy or postpartum follow-up.There is no literature review on the application of CNV-seq technology in the detection of cytomegalovirus infection,only literature reports on genome analysis of CMV-DNA in confirmed patients were available.Conclusion:CNV-seq can be used to detect cytomegalovirus load,which may have a certain degree of predictive value for fetal outcome.CNV-seq can simultaneously detect fetal chromosomes and pathogenic microorganisms,which is of great significance for the prevention and control of birth defects.

Complete Genome Sequencing and Genetic Variation Analysis of Two H9N2 Subtype Avian Influenza Virus Strains

[Objective] The study aimed to investigate the genetic variation characters of entire sequences between two H9N2 subtype avian influenza virus strains and other reference strains.[Method] The entire sequences of 8 genes were obtained by using RT-PCR,and these sequences were analyzed with that of six H9N2 subtype avian influenza isolates in homology comparison and genetic evolution relation.[Result] The results showed that the nucleotide sequence of entire gene of the strain shared 91.1%-95.4% homology with other seven reference strains,and PG08 shared the highest homology 91.3% with C/BJ/1/94;ZD06 shared the highest homology 92.3% with D/HK/Y280/97.HA cleavage sites of two H9N2 subtype avian influenza virus isolated strains were PARSSR/GLF,typical of mildly pathogenic avian influenza virus.[Conclusion] Phylogenetic tree for entire gene of eight strains showed that the genetic relationship was the closest between ZD06 and C/Pak/2/99 strains,which belonged to the Eurasian lineage;PG08 shared the highest homology 91.3% with ZD06,it may be the product of gene rearrangements of other sub-lines.

Improving the accuracy of genomic prediction for meat quality traits using whole genome sequence data in pigs

Background Pork quality can directly affect customer purchase tendency and meat quality traits have become valu-able in modern pork production.However,genetic improvement has been slow due to high phenotyping costs.In this study,whole genome sequence(WGS)data was used to evaluate the prediction accuracy of genomic best linear unbiased prediction(GBLUP)for meat quality in large-scale crossbred commercial pigs.Results We produced WGS data(18,695,907 SNPs and 2,106,902 INDELs exceed quality control)from 1,469 sequenced Duroc×(Landrace×Yorkshire)pigs and developed a reference panel for meat quality including meat color score,marbling score,L*(lightness),a*(redness),and b*(yellowness)of genomic prediction.The prediction accuracy was defined as the Pearson correlation coefficient between adjusted phenotypes and genomic estimated breeding values in the validation population.Using different marker density panels derived from WGS data,accuracy differed substantially among meat quality traits,varied from 0.08 to 0.47.Results showed that MultiBLUP outperform GBLUP and yielded accuracy increases ranging from 17.39%to 75%.We optimized the marker density and found medium-and high-density marker panels are beneficial for the estimation of heritability for meat quality.Moreover,we conducted genotype imputation from 50K chip to WGS level in the same population and found average concord-ance rate to exceed 95%and r^(2)=0.81.Conclusions Overall,estimation of heritability for meat quality traits can benefit from the use of WGS data.This study showed the superiority of using WGS data to genetically improve pork quality in genomic prediction.

Breed identification using breed‑informative SNPs and machine learning based on whole genome sequence data and SNP chip data

Background Breed identification is useful in a variety of biological contexts.Breed identification usually involves two stages,i.e.,detection of breed-informative SNPs and breed assignment.For both stages,there are several methods proposed.However,what is the optimal combination of these methods remain unclear.In this study,using the whole genome sequence data available for 13 cattle breeds from Run 8 of the 1,000 Bull Genomes Project,we compared the combinations of three methods(Delta,FST,and In)for breed-informative SNP detection and five machine learning methods(KNN,SVM,RF,NB,and ANN)for breed assignment with respect to different reference population sizes and difference numbers of most breed-informative SNPs.In addition,we evaluated the accuracy of breed identification using SNP chip data of different densities.Results We found that all combinations performed quite well with identification accuracies over 95%in all scenarios.However,there was no combination which performed the best and robust across all scenarios.We proposed to inte-grate the three breed-informative detection methods,named DFI,and integrate the three machine learning methods,KNN,SVM,and RF,named KSR.We found that the combination of these two integrated methods outperformed the other combinations with accuracies over 99%in most cases and was very robust in all scenarios.The accuracies from using SNP chip data were only slightly lower than that from using sequence data in most cases.Conclusions The current study showed that the combination of DFI and KSR was the optimal strategy.Using sequence data resulted in higher accuracies than using chip data in most cases.However,the differences were gener-ally small.In view of the cost of genotyping,using chip data is also a good option for breed identification.

Screening biomarkers for spinal cord injury using weighted gene co-expression network analysis and machine learning

Immune changes and inflammatory responses have been identified as central events in the pathological process of spinal co rd injury.They can greatly affect nerve regeneration and functional recovery.However,there is still limited understanding of the peripheral immune inflammato ry response in spinal cord inju ry.In this study.we obtained microRNA expression profiles from the peripheral blood of patients with spinal co rd injury using high-throughput sequencing.We also obtained the mRNA expression profile of spinal cord injury patients from the Gene Expression Omnibus(GEO)database(GSE151371).We identified 54 differentially expressed microRNAs and 1656 diffe rentially expressed genes using bioinformatics approaches.Functional enrichment analysis revealed that various common immune and inflammation-related signaling pathways,such as neutrophil extracellular trap formation pathway,T cell receptor signaling pathway,and nuclear factor-κB signal pathway,we re abnormally activated or inhibited in spinal cord inju ry patient samples.We applied an integrated strategy that combines weighted gene co-expression network analysis,LASSO logistic regression,and SVM-RFE algorithm and identified three biomarke rs associated with spinal cord injury:ANO10,BST1,and ZFP36L2.We verified the expression levels and diagnostic perfo rmance of these three genes in the original training dataset and clinical samples through the receiver operating characteristic curve.Quantitative polymerase chain reaction results showed that ANO20 and BST1 mRNA levels were increased and ZFP36L2 mRNA was decreased in the peripheral blood of spinal cord injury patients.We also constructed a small RNA-mRNA interaction network using Cytoscape.Additionally,we evaluated the proportion of 22 types of immune cells in the peripheral blood of spinal co rd injury patients using the CIBERSORT tool.The proportions of naive B cells,plasma cells,monocytes,and neutrophils were increased while the proportions of memory B cells,CD8^(+)T cells,resting natural killer cells,resting dendritic cells,and eosinophils were markedly decreased in spinal cord injury patients increased compared with healthy subjects,and ANO10,BST1 and ZFP26L2we re closely related to the proportion of certain immune cell types.The findings from this study provide new directions for the development of treatment strategies related to immune inflammation in spinal co rd inju ry and suggest that ANO10,BST2,and ZFP36L2 are potential biomarkers for spinal cord injury.The study was registe red in the Chinese Clinical Trial Registry(registration No.ChiCTR2200066985,December 12,2022).

Evaluating the potential of(epi)genotype‑by‑low pass nanopore sequencing in dairy cattle:a study on direct genomic value and methylation analysis

Background Genotype-by-sequencing has been proposed as an alternative to SNP genotyping arrays in genomic selection to obtain a high density of markers along the genome.It requires a low sequencing depth to be cost effective,which may increase the error at the genotype assigment.Third generation nanopore sequencing technology offers low cost sequencing and the possibility to detect genome methylation,which provides added value to genotype-by-sequencing.The aim of this study was to evaluate the performance of genotype-by-low pass nanopore sequencing for estimating the direct genomic value in dairy cattle,and the possibility to obtain methylation marks simultaneously.Results Latest nanopore chemistry(LSK14 and Q20)achieved a modal base calling accuracy of 99.55%,whereas previous kit(LSK109)achieved slightly lower accuracy(99.1%).The direct genomic value accuracy from genotype-by-low pass sequencing ranged between 0.79 and 0.99,depending on the trait(milk,fat or protein yield),with a sequencing depth as low as 2×and using the latest chemistry(LSK114).Lower sequencing depth led to biased estimates,yet with high rank correlations.The LSK109 and Q20 achieved lower accuracies(0.57-0.93).More than one million high reliable methylated sites were obtained,even at low sequencing depth,located mainly in distal intergenic(87%)and promoter(5%)regions.Conclusions This study showed that the latest nanopore technology in useful in a LowPass sequencing framework to estimate direct genomic values with high reliability.It may provide advantages in populations with no available SNP chip,or when a large density of markers with a wide range of allele frequencies is needed.In addition,low pass sequencing provided nucleotide methylation status of>1 million nucleotides at≥10×,which is an added value for epigenetic studies.

Surveillance of emerging SARS-CoV-2 variants by nanopore technology-based genome sequencing

Objective:To surveill emerging variants by nanopore technology-based genome sequencing in different COVID-19 waves in Sri Lanka and to examine the association with the sample characteristics,and vaccination status.Methods:The study analyzed 207 RNA positive swab samples received to sequence laboratory during different waves.The N gene cut-off threshold of less than 30 was considered as the major inclusion criteria.Viral RNA was extracted,and elutes were subjected to nanopore sequencing.All the sequencing data were uploaded in the publicly accessible database,GISAID.Results:The Omicron,Delta and Alpha variants accounted for 58%,22%and 4%of the variants throughout the period.Less than 1%were Kappa variant and 16%of the study samples remained unassigned.Omicron variant was circulated among all age groups and in all the provinces.Ct value and variants assigned percentage was 100%in Ct values of 10-15 while only 45%assigned Ct value over 25.Conclusions:The present study examined the emergence,prevalence,and distribution of SARS-CoV-2 variants locally and has shown that nanopore technology-based genome sequencing enables whole genome sequencing in a low resource setting country.

Genome sequencing and phylogenetic analysis of allotetraploid Salix matsudana Koidz

Polyploidy is a common phenomenon among willow species.In this study,genome sequencing was conducted for Salix matsudana Koidz(also named Chinese willow),an important greening and arbor tree species,and the genome of this species was compared with those of four other tree species in Salicaceae.The total genome sequence of S.matsudana was 655.72 Mb in size,with repeated sequences accounting for 45.97%of the total length.In total,531.43 Mb of the genome sequence could be mapped onto 38 chromosomes using the published genetic map as a reference.The genome of S.matsudana could be divided into two groups,the A and B genomes,through homology analysis with the genome of Populus trichocarpa,and the A and B genomes contained 23,985 and 25,107 genes,respectively.4DTv combined transposon analysis predicted that allotetraploidy in S.matsudana appeared~4 million years ago.The results from this study will help reveal the evolutionary history of S.matsudana and lay a genetic basis for its breeding.

Association between childhood obesity and gut microbiota:16S rRNA gene sequencing-based cohort study

BACKGROUND This study aimed to identify characteristic gut genera in obese and normal-weight children(8-12 years old)using 16S rDNA sequencing.The research aimed to provide insights for mechanistic studies and prevention strategies for childhood obesity.Thirty normal-weight and thirty age-and sex-matched obese children were included.Questionnaires and body measurements were collected,and fecal samples underwent 16S rDNA sequencing.Significant differences in body mass index(BMI)and body-fat percentage were observed between the groups.Analysis of gut microbiota diversity revealed lowerα-diversity in obese children.Differences in gut microbiota composition were found between the two groups.Prevotella and Firmicutes were more abundant in the obese group,while Bacteroides and Sanguibacteroides were more prevalent in the control group.AIM To identify the characteristic gut genera in obese and normal-weight children(8-12-year-old)using 16S rDNA sequencing,and provide a basis for subsequent mechanistic studies and prevention strategies for childhood obesity.METHODS Thirty each normal-weight,1:1 matched for age and sex,and obese children,with an obese status from 2020 to 2022,were included in the control and obese groups,respectively.Basic information was collected through questionnaires and body measurements were obtained from both obese and normal-weight children.Fecal samples were collected from both groups and subjected to 16S rDNA sequencing using an Illumina MiSeq sequencing platform for gut microbiota diversity analysis.RESULTS Significant differences in BMI and body-fat percentage were observed between the two groups.The Ace and Chao1 indices were significantly lower in the obese group than those in the control group,whereas differences were not significant in the Shannon and Simpson indices.Kruskal-Wallis tests indicated significant differences in unweighted and weighted UniFrac distances between the gut microbiota of normal-weight and obese children(P<0.01),suggesting substantial disparities in both the species and quantity of gut microbiota between the two groups.Prevotella,Firmicutes,Bacteroides,and Sanguibacteroides were more abundant in the obese and control groups,respectively.Heatmap results demonstrated significant differences in the gut microbiota composition between obese and normal-weight children.CONCLUSION Obese children exhibited lowerα-diversity in their gut microbiota than did the normal-weight children.Significant differences were observed in the composition of gut microbiota between obese and normal-weight children.

Genome sequencing provides potential strategies for drug discovery and synthesis

Medicinal plants are renowned for their abundant production of secondary metabolites,which exhibit notable pharmacological activities and great potential for drug development.The biosynthesis of secondary metabolites is highly intricate and influenced by various intrinsic and extrinsic factors,resulting in substantial species diversity and content variation.Consequently,precise regulation of secondary metabolite synthesis is of utmost importance.In recent years,genome sequencing has emerged as a valuable tool for investigating the synthesis and regulation of secondary metabolites in medicinal plants,facilitated by the widespread use of high-throughput sequencing technologies.This review highlights the latest advancements in genome sequencing within this field and presents several strategies for studying secondary metabolites.Specifically,the article elucidates how genome sequencing can unravel the pathways for secondary metabolite synthesis in medicinal plants,offering insights into the functions and regulatory mechanisms of participating enzymes.Comparative analyses of plant genomes allow identification of shared pathways of metabolite synthesis among species,thereby providing novel avenues for obtaining cost-effective biosynthetic intermediates.By examining individual genomic variations,genes or gene clusters associated with the synthesis of specific compounds can be discovered,indicating potential targets and directions for drug development and the exploration of alternative compound sources.Moreover,the advent of gene-editing technology has enabled the precise modifications of medicinal plant genomes.Optimization of specific secondary metabolite synthesis pathways becomes thus feasible,enabling the precise editing of target genes to regulate secondary metabolite production within cells.These findings serve as valuable references and lessons for future drug development endeavors,conservation of rare resources,and the exploration of new resources.

Characterization of Rheum palmatum mitochondrial genome and comparative analysis among Caryophyllales species

Objective:This work aimed to report the first complete mitochondrial genome(mitogenome)of Rheum palmatum,summarize the features of Caryophyllales mitogenomes,and to reveal the potential of utilizing the mitogenomes of R.palmatum and other Caryophyllales species for inferring phylogenetic relationships and species identification.Methods:Both Illumina short reads and PacBio HiFi reads were utilized to obtain a complete mitogenome of R.palmatum.A variety of bioinformatics tools were employed to characterize the R.palmatum mitogenome,compare the reported mitogenomes in Caryophyllales and conduct phylogenetic analysis.Results:The mitogenome of R.palmatum was assembled into a single master circle of 302,993 bp,encoding 35 known protein-coding genes,18 transfer RNA genes,and three ribosome RNA genes.A total of 249 long repeats and 49 simple sequence repeats were identified in this mitogenome.The sizes of mitogenomes in Caryophyllales varied from 253 kb to 11.3 Mb.Among them,23 mitogenomes were circular molecules,one was linear,and one consisted of relaxed circles,linear molecules,and supercoiled DNA.Out of the total mitogenomes,11 were single-chromosome structure,whereas the remaining 14 were multi-chromosomal organizations.The phylogenetic analysis is consistent with both the Engler system(1964)and the Angiosperm Phylogeny Group III system.Conclusions:We obtained the first mitogenome of R.palmatum,which consists of a master circle.Mitogenomes in Caryophyllales have variable genome sizes and structures even within the same species.Circular molecules are still the dominant pattern in Caryophyllales.Single-chromosome mitogenomes account for nearly a half of all the mitogenomes in Caryophyllales,in contrast to previous studies.It is feasible to utilize mitochondrial genomes for inferring phylogenetic relationships and conducting species identification.

Transcriptomic and bioinformatics analysis of the mechanism by which erythropoietin promotes recovery from traumatic brain injury in mice

Recent studies have found that erythropoietin promotes the recovery of neurological function after traumatic brain injury.However,the precise mechanism of action remains unclea r.In this study,we induced moderate traumatic brain injury in mice by intrape ritoneal injection of erythro poietin for 3 consecutive days.RNA sequencing detected a total of 4065 differentially expressed RNAs,including 1059 mRNAs,92 microRNAs,799 long non-coding RNAs,and 2115circular RNAs.Kyoto Encyclopedia of Genes and Genomes and Gene Ontology analyses revealed that the coding and non-coding RNAs that were differentially expressed after traumatic brain injury and treatment with erythropoietin play roles in the axon guidance pathway,Wnt pathway,and MAPK pathway.Constructing competing endogenous RNA networks showed that regulatory relationship between the differentially expressed non-coding RNAs and mRNAs.Because the axon guidance pathway was repeatedly enriched,the expression of Wnt5a and Ephb6,key factors in the axonal guidance pathway,was assessed.Ephb6 expression decreased and Wnt5a expression increased after traumatic brain injury,and these effects were reversed by treatment with erythro poietin.These findings suggest that erythro poietin can promote recove ry of nerve function after traumatic brain injury through the axon guidance pathway.

Development and parentage analysis of SNP markers for Chestnut-vented Nuthatch(Sitta nagaensis)based on ddRAD-seq data

Extra-pair paternity(EPP)is commonly found in socially monogamous birds,especially in small passerine birds,and there are interspecific and intraspecific variations in the extent of EPP.The Chestnut-vented Nuthatch(Sitta nagaensis)is a socially monogamous passerine bird,and verifying whether this species has EPP relies on parentage testing-S.nagaensis is not known to have EPP.In this study,we developed SNP markers of this species that are informative for parentage analysis from double digest restriction site-associated DNA sequencing(ddRAD-seq)data.A panel consisting of 50 SNP markers,with a mean heterozygosity of 0.343,was used to resolve 95% of nestlings to fathers.The combined exclusion probabilities for the first parent and second parent were 0.991 and 0.9999,respectively.This panel of SNP markers is a powerful tool for parentage assignments in S.nagaensis.In addition,we found that three offspring(7.9%)from three nests(23.1%)were the result of extra-pair fertilization out of 38 offspring in 13 nests.Our study provided information on parentage analysis that has not been reported before in S.nagaensis.It also supplemented the understudied EPP behavior of birds in Asia,contributing to a general understanding of the EPP behaviors of birds.

Fine-mapping and primary analysis of candidate genes associated with seed coat color in mung bean(Vigna radiata L.)

Seed coat color affects the appearance and commodity quality of mung beans(Vigna radiata L.).The substances that affect mung bean seed coat color are mainly flavonoids,which have important medicinal value.Mapping the seed coat color gene in mung beans would facilitate the development of new varieties and improve their value.In this study,an F2 mapping population consisting of 546 plants was constructed using Jilv9(black seed coat)and BIS9805(green seed coat).Using bulk segregated analysis(BSA)sequencing and kompetitive allele-specific PCR(KASP)markers,the candidate region related to seed coat color was finally narrowed to 0.66 Mb on chromosome(Chr.)4 and included eight candidate genes.Combined transcriptome and metabolome analyses showed that three of the eight candidate genes(LOC106758748,LOC106758747,and LOC106759075)were differentially expressed,which may have caused the differences in flavonoid metabolite content between Jilv9 and BIS9805.These findings can provide a research basis for cloning the genes related to seed coat color and accelerate molecular markerassisted selection breeding in mung beans.

Whole Genome Analysis Reveals New Insights into Macrolide Resistance in Mycoplasma pneumoniae

Objective Mutations in 23 S rRNA gene are known to be associated with macrolide resistance in Mycoplasma pneumoniae（M. pneumoniae）. However, these mutations alone do not fully explain the high resistance rates in Asia. The aim of this study was to investigate other possible mutations involved in macrolide resistance in M. pneumoniae. Methods The whole genomes of 10 clinical isolates of M. pneumoniae with macrolide resistance were sequenced by Illumina Hi Seq2000 platform. The role of the macrolide-specific efflux transporter was assessed by efflux-pump inhibition assays with reserpine and carbonyl cyanide m-chlorophenyl-hydrazone（CCCP）. Results A total of 56 single nucleotide polymorphisms（SNPs） were identified in 10 clinical isolates in comparison to the reference strains M129 and FH. Strikingly, 4 of 30 SNPs causing non-synonymous mutations were clustered in macrolide-specific efflux system gene mac B encoding macrolide-specific efflux pump protein of the ATP-binding cassette transporter family. In assays of the minimal inhibitory concentrations（MIC） of macrolide antibiotics in the presence of the efflux pump inhibitors caused a significant decrease of MICs, even under detectable levels in some strains. Conclusion Our study suggests that macrolide efflux pump may contribute to macrolide resistance in M. pneumoniae in addition to the common point mutations in 23 S r RNA gene.