Transversions have larger regulatory effects than transitions

dbEBV: A database of Epstein-Barr virus variants and their correlations with human health

Since Epstein-Barr virus (EBV) was discovered in 1964, it has been reported to be associated with various malignancies as well as benign diseases, and the pathogenicity of EBV has been widely studied. Several databases have been established to provide comprehensive information on the virus and its relation to diseases and introduce convenient analysis tools. Although they have greatly facilitated the analysis of EBV at the genome, gene, protein, or epitope level, they did not provide enough insight into the genomic variants of EBV, which have been suggested as relevant to diseases by multiple studies. Here, we introduce dbEBV, a comprehensive database of EBV genomic variation landscape, which contains 942 EBV genomes with 109,893 variants from different tissues or cell lines in 24 countries. The database enables the visualization of information with varying global frequencies and their relationship with the human health of each variant. It also supports phylogenetic analysis at the genome or gene level in subgroups of different characteristics. Information of interest can easily be reached with functions such as searching, browsing, and filtering. In conclusion, dbEBV is a convenient resource for exploring EBV genomic variants, freely available at http://dbebv.omicsbio.info.

RNA-seq transcriptome and pathway analysis of the medicinal mushroom Lignosus tigris (Polyporaceae) offer insights into its bioactive compounds with anticancer and antioxidant potential

ETHNOPHARMACOLOGICAL RELEVANCE

Medicinal mushrooms belonging to the Lignosus spp., colloquially known as Tiger Milk mushrooms (TMMs), are used as traditional medicine by communities across various regions of China and Southeast Asia to enhance immunity and to treat various diseases. At present, three Lignosus species have been identified in Malaysia: L. rhinocerus, L. tigris, and L. cameronensis. Similarities in their macroscopic morphologies and the nearly indistinguishable appearance of their sclerotia often lead to interchangeability between them. Hence, substantiation of their traditional applications via identification of their individual bioactive properties is imperative in ensuring that they are safe for consumption. L. tigris was first identified in 2013. Thus far, studies on L. tigris cultivar sclerotia (Ligno TG-K) have shown that it possesses significant antioxidant activities and has greater antiproliferative action against selected cancer cells in vitro compared to its sister species, L. rhinocerus TM02®. Our previous genomics study also revealed significant genetic dissimilarities between them. Further omics investigations on Ligno TG-K hold immense potential in facilitating the identification of its bioactive compounds and their associated bioactivities.

AIM OF STUDY

The overall aim of this study was to investigate the gene expression profile of Ligno TG-K via de novo RNA-seq and pathway analysis. We also aimed to identify highly expressed genes encoding compounds that contribute to its cytotoxic and antioxidant properties, as well as perform a comparative transcriptomics analysis between Ligno TG-K and its sister species, L. rhinocerus TM02®.

MATERIALS AND METHODS

Total RNA from fresh 3-month-old cultivated L. tigris sclerotia (Ligno TG-K) was extracted and analyzed via de novo RNA sequencing. Expressed genes were analyzed using InterPro and NCBI-Nr databases for domain identification and homology search. Functional categorization based on gene functions and pathways was performed using Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG), and Clusters of Orthologous Genes (COG) databases. Selected genes were subsequently subjected to phylogenetic analysis.

RESULTS

Our transcriptomics analysis of Ligno TG-K revealed that 68.06% of its genes are expressed in the sclerotium; 80.38% of these were coding transcripts. Our analysis identified highly expressed transcripts encoding proteins with prospective medicinal properties. These included serine proteases (FPKM = 7356.68), deoxyribonucleases (FPKM = 3777.98), lectins (FPKM = 3690.87), and fungal immunomodulatory proteins (FPKM = 2337.84), all of which have known associations with anticancer activities. Transcripts linked to proteins with antioxidant activities, such as superoxide dismutase (FPKM = 1161.69) and catalase (FPKM = 1905.83), were also highly expressed. Results of our sequence alignments revealed that these genes and their orthologs can be found in other mushrooms. They exhibit significant sequence similarities, suggesting possible parallels in their anticancer and antioxidant bioactivities.

CONCLUSION

This study is the first to provide a reference transcriptome profile of genes expressed in the sclerotia of L. tigris. The current study also presents distinct COG profiles of highly expressed genes in Ligno TG-K and L. rhinocerus TM02®, highlighting that any distinctions uncovered may be attributed to their interspecies variations and inherent characteristics that are unique to each species. Our findings suggest that Ligno TG-K contains bioactive compounds with prospective medicinal properties that warrant further investigations.

CLASSIFICATION

Systems biology and omics.

Structure and Phylogenetic Relationships of Scolopacidae Mitogenomes (Charadriiformes: Scolopacidae)

The family Scolopacidae presents a valuable subject for evolutionary research; however, molecular studies of Scolopacidae are still relatively understudied, and the phylogenetic relationships of certain species remain unclear. In this study, we sequenced and obtained complete mitochondrial DNA (mtDNA) from Actitis hypoleucos and partial mtDNA from Numenius arquata, Limosa limosa, and Limnodromus semipalmatus. The complete mtDNA contained 13 protein-coding genes (PCGs), two ribosomal RNA genes, 22 tRNA genes, and a control region. Scolopacidae contained three types of start codons and five types of stop codons (including one incomplete stop codon, T--). In 13 protein-coding genes, average uncorrected pairwise distances (Aupd) revealed that ATP8 was the least conserved while COX3 had the lowest evolutionary rate. The ratio of Ka/Ks suggested that all PCGs were under purifying selection. Using two methods (maximum likelihood and Bayesian inference) to analyze the phylogenetic relationships of the family Scolopacidae, it was found that the genera Xenus and Actitis were clustered into another sister group, while the genus Phalaropus is more closely related to the genus Tringa. The genera Limnodromus, Gallinago, and Scolopax form a monophyletic group. This study improves our understanding of the evolutionary patterns and phylogenetic relationships of the family Scolopacidae.

Detection of Transversions and Transitions in HBG2 Cis-Elements Associated with Sickle Cell Allele in Ghanaians

Short tandem repeats located 5' prime to the β-globin gene, have been observed to be in linkage disequilibrium with the HbS allele, and thought to affect the severity of sickle cell disease. Here, we report on new mutants within the HBG2 region that may impact sickle cell disease. To determine the cis-acting elements microsatellites, indels and single nucleotide polymorphisms (SNPs), within the HBG2 region by sequencing, in subjects with sickle cell disease. The case-control study was located at the Center for Clinical Genetics, Sickle cell unit, Korle-Bu Teaching Hospital. A questionnaire was used for demographic data and clinical information. Hematological profile (red blood cell, white blood cell, platelet, hemoglobin and mean corpuscular volume) were assessed in 83 subjects. A set of 45 samples comprising amplified DNA on the HBG2 gene from HbSS (22), HbSC (17) and 6 controls (HbAA) were sequenced. Differences in the microsatellite region between sickle cell disease (SCD) (HbSS and HbSC) genotypes and control subjects were identified by counting and assessed by Chi-square analysis. Red blood cells, hematocrit, platelets, white blood cells and hemoglobin indices differed in genotypic groups. HbSS subjects were affirmed to have severer hemolytic anemia than HbSC subjects. Two indels (T1824 and C905) were seen in both SS and SC genotypes. Two peculiar SNPs: G:T1860 (transition) and A:G1872 transversions were found within the HBG2 gene that were significantly associated with the HbSS genotype (Fisher's exact test, p = 0.006) and HbS allele respectively (Fisher's exact test, p = 0.006). Cis-acting elements in HbSS and HbSC were different and may contribute to the phenotype seen in the disease state.

Genetic Diversity Analysis and Core Germplasm Construction of Rubus chingii Hu

Rubus chingii Hu is the only species that is used for both edible and medicinal purposes among the 194 species of the genus Rubus in China. It is well known for its sweet and sour fresh fruits that are rich in vitamins and for its dried immature fruits that are used to treat kidney-related ailments. This study aims to evaluate genetic diversity and population structure and build a core germplasm repository of 132 R. chingii accessions from the provinces of Jiangxi and Fujian, using Hyper-seq-derived single-nucleotide polymorphism (SNP) markers. This is the first genetic study of R. chingii based on SNP molecular markers, and a total of 1,303,850 SNPs and 433,159 insertions/deletions (InDels) were identified. Low values for observed heterozygosity, nucleotide diversity (Pi) and fixation indexes (Fis) indicated low genetic diversity within populations, and an analysis of molecular variance (AMOVA) showed that 37.4% and 62.6% of the variations were found between populations and within samples, respectively. Four main clusters were identified by means of neighbor-joining (NJ) trees, the ADMIXTURE program and principal component analysis (PCA). Based on the genetic diversity, we finally constructed 38 representative core collections, representing 50% of the total core germplasm samples and 95.3% of the genotypes. In summary, the results of our study can provide valuable information on the genetic structure of R. chingii germplasm resources, which is helpful for further explorations of potential high-quality genes and for formulating future breeding and conservation strategies.

Effectiveness of DArTseq markers application in genetic diversity and population structure of indigenous chickens in Eastern Province of Rwanda

BACKGROUND

The application of biotechnologies which make use of genetic markers in chicken breeding is developing rapidly. Diversity Array Technology (DArT) is one of the current Genotyping-By-Sequencing techniques allowing the discovery of whole genome sequencing. In livestock, DArT has been applied in cattle, sheep, and horses. Currently, there is no study on the application of DArT markers in chickens. The aim was to study the effectiveness of DArTSeq markers in the genetic diversity and population structure of indigenous chickens (IC) and SASSO in the Eastern Province of Rwanda.

METHODS

In total 87 blood samples were randomly collected from 37 males and 40 females of indigenous chickens and 10 females of SASSO chickens purposively selected from 5 sites located in two districts of the Eastern Province of Rwanda. Genotyping by Sequencing (GBS) using DArTseq technology was employed. This involved the complexity reduction method through digestion of genomic DNA and ligation of barcoded adapters followed by PCR amplification of adapter-ligated fragments.

RESULTS

From 45,677 DArTseq SNPs and 25,444 SilicoDArTs generated, only 8,715 and 6,817 respectively remained for further analysis after quality control. The average call rates observed, 0.99 and 0.98 for DArTseq SNPs and SilicoDArTs respectively were quite similar. The polymorphic information content (PIC) from SilicoDArTs (0.33) was higher than that from DArTseq SNPs (0.22). DArTseq SNPs and SilicoDArTs had 34.4% and 34% of the loci respectively mapped on chromosome 1. DArTseq SNPs revealed distance averages of 0.17 and 0.15 within IC and SASSO chickens respectively while the respective averages observed with SilicoDArTs were 0.42 and 0.36. The average genetic distance between IC and SASSO chickens was moderate for SilicoDArTs (0.120) compared to that of DArTseq SNPs (0.048). The PCoA and population structure clustered the chicken samples into two subpopulations (1 and 2); 1 is composed of IC and 2 by SASSO chickens. An admixture was observed in subpopulation 2 with 12 chickens from subpopulation 1.

CONCLUSIONS

The application of DArTseq markers have been proven to be effective and efficient for genetic relationship between IC and separated IC from exotic breed used which indicate their suitability in genomic studies. However, further studies using all chicken genetic resources available and large big sample sizes are required.

The Construction of lncRNA/circRNA-miRNA-mRNA Networks Reveals Functional Genes Related to Growth Traits in Schima superba

Schima superba is a precious timber and fire-resistant tree species widely distributed in southern China. Currently, there is little knowledge related to its growth traits, especially with respect to molecular breeding. The lack of relevant information has delayed the development of modern breeding. The purpose is to identify probable functional genes involved in S. superba growth through whole transcriptome sequencing. In this study, a total of 32,711 mRNAs, 525 miRNAs, 54,312 lncRNAs, and 1522 circRNAs were identified from 10 S. superba individuals containing different volumes of wood. Four possible regulators, comprising three lncRNAs, one circRNA, and eleven key miRNAs, were identified from the regulatory networks of lncRNA-miRNA-mRNA and circRNA-miRNA-mRNA to supply information on ncRNAs. Several candidate genes involved in phenylpropane and cellulose biosynthesis pathways, including Ss4CL2, SsCSL1, and SsCSL2, and transcription factors, including SsDELLA2 (SsSLR), SsDELLA3 (SsSLN), SsDELLA5 (SsGAI-like2), and SsNAM1, were identified to reveal the molecular regulatory mechanisms regulating the growth traits of S. superba. The results not merely provide candidate functional genes related to S. superba growth trait and will be useful to carry out molecular breeding, but the strategy and method also provide scientists with an effective approach to revealing mechanisms behind important economic traits in other species.

Morphological and molecular phylogenetic characterization of Sarcocystis kani sp. nov. and other novel, closely related Sarcocystis spp. infecting small mammals and colubrid snakes in Asia

We investigated the morphology and phylogenetic relationships of novel and previously recognized Sarcocystis spp. infecting small mammals and colubrid snakes in Asia. The nuclear 18S rRNA and mitochondrial cox1 of Sarcocystis sp.1 from mangrove snakes (Boiga dendrophila) in Thailand and Sarcocystis sp.2 from a ricefield rat (Rattus argentiventer) in Sumatra were partially sequenced. Sporocysts of Sarcocystis sp.1 induced development of sarcocysts in experimentally infected rats, which showed a unique ultrastructure that was observed previously by S.P. Kan in rats from Malaysia; therefore, we describe this species as Sarcocystis kani sp. nov. Its integration into the 18S rRNA phylogeny of Sarcocystis spp. cycling between small mammals and colubrid snakes helped clarify relationships among the so-called S. zuoi-complex of molecularly cryptic species: Sarcocystis kani sp. nov., S. sp.2, S. attenuati, S. scandentiborneensis, and S. zuoi were all included in this clade. Tree topology was resolved into dichotomies congruent with the morphological disparities between the taxa. However, cox1 gene sequencing (including newly sequenced S. singaporensis and S. zamani) revealed that Sarcocystis kani, S. attenuati, and S. scandentiborneensis were identical suggesting a recent, common ancestry. To identify other distinctive features, lineage-specific molecular patterns within both genes were examined revealing that all 18S rRNA sequences of the S. zuoi - complex possess a unique, 7-nt long motif in helix 38 of domain V7 that was different in S. clethrionomyelaphis which branched off basally from the complex. Three-dimensional homology modelling of COX1 protein structure identified amino acid substitutions within the barcode area specific for the S. zuoi-complex and substantial divergence in structurally important amino acids between Sarcocystis species of snakes as definitive hosts and other lineages of the Sarcocystidae. We discuss the utility of selected genes for species delimitation of the Sarcocystis spp. under investigation, which probably evolved during recent radiations of their intermediate and definitive hosts.

Effect of single nucleotide polymorphisms on the structure of long noncoding RNAs and their interaction with RNA binding proteins

Long non-coding RNAs (lncRNA) are emerging as a new class of regulatory RNAs with remarkable potential to be utilized as therapeutic targets against many human diseases. Several genome-wide association studies (GWAS) have catalogued Single Nucleotide Polymorphisms (SNPs) present in the noncoding regions of the genome from where lncRNAs originate. In this study, we have selected 67 lncRNAs with GWAS-tagged SNPs and have also investigated their role in affecting the local secondary structures. Majority of the SNPs lead to changes in the secondary structure of lncRNAs to a different extent by altering the base pairing patterns. These structural changes in lncRNA are also manifested in form of alteration in the binding site for RNA binding proteins (RBPs) along with affecting their binding efficacies. Ultimately, these structural modifications may influence the transcriptional and post-transcriptional pathways of these RNAs, leading to the causation of diseases. Hence, it is important to understand the possible underlying mechanism of RBPs in association with GWAS-tagged SNPs in human diseases.

Association mapping in bambara groundnut [Vigna subterranea (L.) Verdc.] reveals loci associated with agro-morphological traits

BACKGROUND

Genome-wide association studies (GWAS) are important for the acceleration of crop improvement through knowledge of marker-trait association (MTA). This report used DArT SNP markers to successfully perform GWAS on agro-morphological traits using 270 bambara groundnut [Vigna subterranea (L.) Verdc.] landraces sourced from diverse origins. The study aimed to identify marker traits association for nine agronomic traits using GWAS and their candidate genes. The experiment was conducted at two different locations laid out in alpha lattice design. The cowpea [Vigna unguiculata (L.) Walp.] reference genome (i.e. legume genome most closely related to bambara groundnut) assisted in the identification of candidate genes.

RESULTS

The analyses showed that linkage disequilibrium was found to decay rapidly with an average genetic distance of 148 kb. The broadsense heritability was relatively high and ranged from 48.39% (terminal leaf length) to 79.39% (number of pods per plant). The GWAS identified a total of 27 significant marker-trait associations (MTAs) for the nine studied traits explaining 5.27% to 24.86% of phenotypic variations. Among studied traits, the highest number of MTAs was obtained from seed coat colour (6) followed by days to flowering (5), while the least is days to maturity (1), explaining 5.76% to 11.03%, 14.5% to 19.49%, and 11.66% phenotypic variations, respectively. Also, a total of 17 candidate genes were identified, varying in number for different traits; seed coat colour (6), days to flowering (3), terminal leaf length (2), terminal leaf width (2), number of seed per pod (2), pod width (1) and days to maturity (1).

CONCLUSION

These results revealed the prospect of GWAS in identification of SNP variations associated with agronomic traits in bambara groundnut. Also, its present new opportunity to explore GWAS and marker assisted strategies in breeding of bambara groundnut for acceleration of the crop improvement.

Profiling mitochondrial DNA mutations in tumors and circulating extracellular vesicles of triple-negative breast cancer patients for potential biomarker development

Early detection and recurrence prediction are challenging in triple-negative breast cancer (TNBC) patients. We aimed to develop mitochondrial DNA (mtDNA)-based liquid biomarkers to improve TNBC management. Mitochondrial genome (MG) enrichment and next-generation sequencing mapped the entire MG in 73 samples (64 tissues and 9 extracellular vesicles [EV] samples) from 32 metastatic TNBCs. We measured mtDNA and cardiolipin (CL) contents, NDUFB8, and SDHB protein expression in tumors and in corresponding circulating EVs. We identified 168 nonsynonymous mtDNA mutations, with 73% (123/186) coding and 27% (45/168) noncoding in nature. Twenty percent of mutations were nucleotide transversions. Respiratory complex I (RCI) was the key target, which harbored 44% (74/168) of the overall mtDNA mutations. A panel of 11 hotspot mtDNA mutations was identified among 19%-38% TNBCs, which were detectable in the serum-derived EVs with 82% specificity. Overall, 38% of the metastatic tumor-signature mtDNA mutations were traceable in the EVs. An appreciable number of mtDNA mutations were homoplasmic (18%, 31/168), novel (14%, 23/168), and potentially pathogenic (9%, 15/168). The overall and RCI-specific mtDNA mutational load was higher in women with African compared to European ancestry accompanied by an exclusive abundance of respiratory complex (RC) protein NDUFB8 (RCI) and SDHB (RCII) therein. Increased mtDNA (p < 0.0001) content was recorded in both tumors and EVs along with an abundance of CL (p = 0.0001) content in the EVs. Aggressive tumor-signature mtDNA mutation detection and measurement of mtDNA and CL contents in the EVs bear the potential to formulate noninvasive early detection and recurrence prediction strategies.

Loss of a Premature Stop Codon in the Rice Wall-Associated Kinase 91 (WAK91) Gene Is a Candidate for Improving Leaf Sheath Blight Disease Resistance

Leaf sheath blight disease (SB) of rice caused by the soil-borne fungus Rhizoctonia solani results in 10-30% global yield loss annually and can reach 50% under severe outbreaks. Many disease resistance genes and receptor-like kinases (RLKs) are recruited early on by the host plant to respond to pathogens. Wall-associated receptor kinases (WAKs), a subfamily of receptor-like kinases, have been shown to play a role in fungal defense. The rice gene WAK91 (OsWAK91), co-located in the major SB resistance QTL region on chromosome 9, was identified by us as a candidate in defense against rice sheath blight. An SNP mutation T/C in the WAK91 gene was identified in the susceptible rice variety Cocodrie (CCDR) and the resistant line MCR010277 (MCR). The consequence of the resistant allele C is a stop codon loss, resulting in an open reading frame with extra 62 amino acid carrying a longer protein kinase domain and additional phosphorylation sites. Our genotype and phenotype analysis of the parents CCDR and MCR and the top 20 individuals of the double haploid SB population strongly correlate with the SNP. The susceptible allele T is present in the japonica subspecies and most tropical and temperate japonica lines. Multiple US commercial rice varieties with a japonica background carry the susceptible allele and are known for SB susceptibility. This discovery opens the possibility of introducing resistance alleles into high-yielding commercial varieties to reduce yield losses incurred by the sheath blight disease.

In-silico identification and comparison of transcription factor binding sites cluster in anterior-posterior patterning genes in Drosophila melanogaster and Tribolium castaneum

The cis-regulatory data that help in transcriptional regulation is arranged into modular pieces of a few hundred base pairs called CRMs (cis-regulatory modules) and numerous binding sites for multiple transcription factors are prominent characteristics of these cis-regulatory modules. The present study was designed to localize transcription factor binding site (TFBS) clusters on twelve Anterior-posterior (A-P) genes in Tribolium castaneum and compare them to their orthologous gene enhancers in Drosophila melanogaster. Out of the twelve A-P patterning genes, six were gap genes (Kruppel, Knirps, Tailless, Hunchback, Giant, and Caudal) and six were pair rule genes (Hairy, Runt, Even-skipped, Fushi-tarazu, Paired, and Odd-skipped). The genes along with 20 kb upstream and downstream regions were scanned for TFBS clusters using the Motif Cluster Alignment Search Tool (MCAST), a bioinformatics tool that looks for set of nucleotide sequences for statistically significant clusters of non-overlapping occurrence of a given set of motifs. The motifs used in the current study were Hunchback, Caudal, Giant, Kruppel, Knirps, and Even-skipped. The results of the MCAST analysis revealed the maximum number of TFBS for Hunchback, Knirps, Caudal, and Kruppel in both D. melanogaster and T. castaneum, while Bicoid TFBS clusters were found only in D. melanogaster. The size of all the predicted TFBS clusters was less than 1kb in both insect species. These sequences revealed more transversional sites (Tv) than transitional sites (Ti) and the average Ti/Tv ratio was 0.75.

Whole-genome resequencing of common bean elite breeding lines

The expansion of bean genome technologies has prompted new perspectives on generating resources and knowledge essential to research and implementing biotechnological tools for the practical operations of plant breeding programs. This study aimed to resequence the entire genome (whole genome sequencing-WGS) of 40 bean genotypes selected based on their significance in breeding programs worldwide, with the objective of generating an extensive database for the identification of single nucleotide polymorphisms (SNPs). Over 6 million SNPs were identified, distributed across the 11 bean chromosomes. After quality variant filtering, 420,509 high-quality SNPs were established, with an average of 38,228 SNPs per chromosome. These variants were categorized based on their predicted effects, revealing that the majority exerted a modifier impact on non-coding genome regions (94.68%). Notably, a significant proportion of SNPs occurred in intergenic regions (62.89%) and at least one SNP was identified in 58.63% of the genes annotated in the bean genome. Of particular interest, 7841 SNPs were identified in 85% of the putative plant disease defense-related genes, presenting a valuable resource for crop breeding efforts. These findings provide a foundation for the development of innovative and broadly applicable technologies for the routine selection of superior genotypes in global bean improvement and germplasm characterization programs.

Morphological and molecular characteristics of Plasmodium juxtanucleare in layer chicken from three districts of Yogyakarta, Indonesia

Background and Aim

Blood parasite infections in poultry, such as Plasmodium, are a serious threat to the poultry industry due to their potential to cause economic losses. To date, there has been inadequate research on the morphological and molecular detection of the different Plasmodium species that infect poultry in Indonesia. Therefore, this study aimed to analyze the morphological and molecular characteristics of Plasmodium spp. and the several predisposing factors for Plasmodium infection in layer chickens from three districts of Yogyakarta, Indonesia.

Materials and Methods

One hundred and five blood samples from layer chickens were collected from 13 farms located in three districts of Yogyakarta (Sleman, Bantul, and Kulon Progo) between September and November 2022. Blood samples were subjected to microscopic and polymerase chain reaction (PCR) analyses. Sequencing was performed using basic local alignment search tools to identify the nucleotide structure of cytochrome b. Phylogenetic analysis of Plasmodium was performed using the MEGA-X software.

Results

Microscopic examination revealed that 17/105 positives (16.19%) were positive for blood parasite infection. Trophozoites, erythrocytic meronts, and microgametocytes of Plasmodium were found in blood samples. Based on the morphological examination, the species found in the samples was close to Plasmodium juxtanucleare. Polymerase chain reaction examination revealed that 21/60 samples were positive for Plasmodium (35%). The Plasmodium species identified from the sequenced samples were proven to be P. juxtanucleare. The P. juxtanucleare from Thailand was closely related to samples (99.64%-100%) with a genetic distance of 0%-1%. In addition, age, population, and cage type were not significantly associated with Plasmodium infection.

Conclusion

Based on microscopic and PCR examinations, the Plasmodium species found in the three districts of Yogyakarta was P. juxtanucleare. The genetic distance between samples from the three districts of Yogyakarta was closely related (0%-1%) to P. juxtanucleare from Thailand and Japan. There was no correlation between Plasmodium infection and age, cage type, or population.

Genome-Wide SNP and Indel Discovery in Abaca (Musa textilis Née) and among Other Musa spp. for Abaca Genetic Resources Management

Abaca (Musa textilis Née) is an economically important fiber crop in the Philippines. Its economic potential, however, is hampered by biotic and abiotic stresses, which are exacerbated by insufficient genomic resources for varietal identification vital for crop improvement. To address these gaps, this study aimed to discover genome-wide polymorphisms among abaca cultivars and other Musa species and analyze their potential as genetic marker resources. This was achieved through whole-genome Illumina resequencing of abaca cultivars and variant calling using BCFtools, followed by genetic diversity and phylogenetic analyses. A total of 20,590,381 high-quality single-nucleotide polymorphisms (SNP) and DNA insertions/deletions (InDels) were mined across 16 abaca cultivars. Filtering based on linkage disequilibrium (LD) yielded 130,768 SNPs and 13,620 InDels, accounting for 0.396 ± 0.106 and 0.431 ± 0.111 of gene diversity across these cultivars. LD-pruned polymorphisms across abaca, M. troglodytarum, M. acuminata and M. balbisiana enabled genetic differentiation within abaca and across the four Musa spp. Phylogenetic analysis revealed the registered varieties Abuab and Inosa to accumulate a significant number of mutations, eliciting further studies linking mutations to their advantageous phenotypes. Overall, this study pioneered in producing marker resources in abaca based on genome-wide polymorphisms vital for varietal authentication and comparative genotyping with the more studied Musa spp.

Genetic Diversity and Population Structure of Sorghum (Sorghum bicolor (L.) Moench) Landraces Using DArTseq-Derived Single-Nucleotide Polymorphism (SNP) Markers

Genetic integrity of an accession should be preserved in the conservation of germplasm. Characterization of diverse germplasm based on a molecular basis enhances its conservation and use in breeding programs. The aim of this study was to assess the genetic diversity of 169 sorghum accessions using a total of 6977 SNP markers. The polymorphic information content of the markers was 0.31 which is considered to be moderately high. Structure analysis using ADMIXTURE program revealed a total of 10 subpopulations. Neighbor-joining tree revealed the presence of six main clusters among these subpopulations whereas in principal component analysis, seven clusters were identified. Cluster analysis grouped most populations depending on source of collection although other accessions originating from the same source were grouped under different clusters. Analysis of molecular variance (AMOVA) revealed 30% and 70% of the variation occurred within and among accessions, respectively. Gene flow within the populations was, however, limited indicating high differentiation within the subpopulation. Observed heterozygosity among accessions varied from 0.03 to 0.06 with a mean of 0.05 since sorghum is a self-pollinating crop. High genetic diversity among the subpopulations can be further explored for superior genes to develop new sorghum varieties.

Impact of ACE and Endoplasmic Reticulum Aminopeptidases Polymorphisms on COVID-19 Outcome

Background: COVID-19 outcomes display multiple unexpected varieties, ranging from unnoticed symptomless infection to death, without any previous alarm or known aggravating factors. Aim: To appraise the impact of ACErs4291(A/T) and ERAP1rs26618(T/C) human polymorphisms on the outcome of COVID-19. Subjects and methods: In total, 240 individuals were enrolled in the study (80 with severe manifestations, 80 with mild manifestations, and 80 healthy persons). ACErs4291(A/T) and ERAP1rs26618(T/C) genotyping was performed using RT-PCR. Results: The frequency of the ACErs4291AA genotype was higher among the severe COVID-19 group than others (p < 0.001). The ERAP1rs26618TT genotype frequency was higher among the severe COVID-19 group in comparison with the mild group (p < 0.001) and non-infected controls (p = 0.0006). The frequency of the ACErs4291A allele was higher among severe COVID-19 than mild and non-infected groups (64.4% vs. 37.5%, and 34.4%, respectively), and the ERAP1rs26618T allele was also higher in the severe group (67.5% vs. 39.4%, and 49.4%). There was a statistically significant association between severe COVID-19 and ACErs4291A or ERAP1rs26618T alleles. The coexistence of ACErs4291A and ERAP1rs26618T alleles in the same individual increase the severity of the COVID-19 risk by seven times [OR (95%CI) (LL−UL) = 7.058 (3.752−13.277), p < 0.001). A logistic regression analysis revealed that age, male gender, non-vaccination, ACErs4291A, and ERAP1rs26618T alleles are independent risk factors for severe COVID-19. Conclusions: Persons carrying ACErs4291A and/or ERAP1rs26618T alleles are at higher risk of developing severe COVID-19.

Integrative Multi-OMICs Identifies Therapeutic Response Biomarkers and Confirms Fidelity of Clinically Annotated, Serially Passaged Patient-Derived Xenografts Established from Primary and Metastatic Pediatric and AYA Solid Tumors

Establishment of clinically annotated, molecularly characterized, patient-derived xenografts (PDXs) from treatment-naïve and pretreated patients provides a platform to test precision genomics-guided therapies. An integrated multi-OMICS pipeline was developed to identify cancer-associated pathways and evaluate stability of molecular signatures in a panel of pediatric and AYA PDXs following serial passaging in mice. Original solid tumor samples and their corresponding PDXs were evaluated by whole-genome sequencing, RNA-seq, immunoblotting, pathway enrichment analyses, and the drug−gene interaction database to identify as well as cross-validate actionable targets in patients with sarcomas or Wilms tumors. While some divergence between original tumor and the respective PDX was evident, majority of alterations were not functionally impactful, and oncogenic pathway activation was maintained following serial passaging. CDK4/6 and BETs were prioritized as biomarkers of therapeutic response in osteosarcoma PDXs with pertinent molecular signatures. Inhibition of CDK4/6 or BETs decreased osteosarcoma PDX growth (two-way ANOVA, p < 0.05) confirming mechanistic involvement in growth. Linking patient treatment history with molecular and efficacy data in PDX will provide a strong rationale for targeted therapy and improve our understanding of which therapy is most beneficial in patients at diagnosis and in those already exposed to therapy.

Reactive Oxygen Species and Their Consequences on the Structure and Function of Mammalian Spermatozoa

Significance: Among the 200 or so cell types that comprise mammals, spermatozoa have an ambiguous relationship with the reactive oxygen species (ROS) inherent in the consumption of oxygen that supports aerobic metabolism. Recent Advances: In this review, we shall see that spermatozoa need the action of ROS to reach their structural and functional maturity, but that due to intrinsic unique characteristics, they are, perhaps more than any other cell type, susceptible to oxidative damage. Recent studies have improved our knowledge of how oxidative damage affects sperm structures and functions. The focus of this review will be on how genetic and epigenetic oxidative alterations to spermatozoa can have dramatic unintended consequences in terms of both the support and the suppression of sperm function. Critical Issues: Oxidative stress can have dramatic consequences not only for the spermatozoon itself, but also, and above all, on its primary objective, which is to carry out fertilization and to ensure, in part, that the embryonic development program should lead to a healthy progeny. Future Directions: Sperm oxidative DNA damage largely affects the integrity of the paternal genetic material to such an extent that the oocyte may have difficulties in correcting it. Diagnostic and therapeutic actions should be considered more systematically, especially in men with difficulties to conceive. Research is underway to determine whether the epigenetic information carried by spermatozoa is also subject to changes mediated by pro-oxidative situations. Antioxid. Redox Signal. 37, 481-500.

Evolution of a cis-Acting SNP That Controls Type VI Secretion in Vibrio cholerae

Mutations in regulatory mechanisms that control gene expression contribute to phenotypic diversity and thus facilitate the adaptation of microbes and other organisms to new niches. Comparative genomics can be used to infer rewiring of regulatory architecture based on large effect mutations like loss or acquisition of transcription factors but may be insufficient to identify small changes in noncoding, intergenic DNA sequence of regulatory elements that drive phenotypic divergence. In human-derived Vibrio cholerae, the response to distinct chemical cues triggers production of multiple transcription factors that can regulate the type VI secretion system (T6), a broadly distributed weapon for interbacterial competition. However, to date, the signaling network remains poorly understood because no regulatory element has been identified for the major T6 locus. Here we identify a conserved cis-acting single nucleotide polymorphism (SNP) controlling T6 transcription and activity. Sequence alignment of the T6 regulatory region from diverse V. cholerae strains revealed conservation of the SNP that we rewired to interconvert V. cholerae T6 activity between chitin-inducible and constitutive states. This study supports a model of pathogen evolution through a noncoding cis-regulatory mutation and preexisting, active transcription factors that confers a different fitness advantage to tightly regulated strains inside a human host and unfettered strains adapted to environmental niches. IMPORTANCE Organisms sense external cues with regulatory circuits that trigger the production of transcription factors, which bind specific DNA sequences at promoters ("cis" regulatory elements) to activate target genes. Mutations of transcription factors or their regulatory elements create phenotypic diversity, allowing exploitation of new niches. Waterborne pathogen Vibrio cholerae encodes the type VI secretion system "nanoweapon" to kill competitor cells when activated. Despite identification of several transcription factors, no regulatory element has been identified in the promoter of the major type VI locus, to date. Combining phenotypic, genetic, and genomic analysis of diverse V. cholerae strains, we discovered a single nucleotide polymorphism in the type VI promoter that switches its killing activity between a constitutive state beneficial outside hosts and an inducible state for constraint in a host. Our results support a role for noncoding DNA in adaptation of this pathogen.

Genetic Diversity and Population Structure of Schima superba From Southern China

The tree Schima superba is important for afforestation and fire prevention in southern China. The wood of this tree can also be used for furniture and buildings. However, the lack of genetic background and genomic information for this species has lowered wood yield speed and quality improvement. Here, we aimed to discover genome-wide single nucleotide polymorphisms (SNPs) in 302 S. superba germplasms collected from southern China and to use these SNPs to investigate the population structure. Using genotyping by sequencing, a total of 785 high-quality SNP markers (minor allele frequency [MAF] ≥ 0.05) were identified from 302 accessions collected from seven geographical locations. Population structure analyses and principal coordinate analyses (PCoAs) indicated that these germplasm resources can be clearly separated into different populations. The S. superba accessions originating from Yunnan (YN) and Guangxi (GX) fell into the same population, separate from the accessions originating from Guangdong (GD), which indicated that these two regions should be regarded as major provenances of this species. In addition, two independent core germplasm sets with abundant genetic polymorphisms were constructed to support the breeding work. The identification of SNP markers, analyses of population genetics, and construction of core germplasm sets will greatly promote the molecular breeding work of S. superba.

MiRLog and dbmiR: prioritization and functional annotation tools to study human microRNA sequence variants

The recent identification of noncoding variants with pathogenic effects suggests that these variations could underlie a significant number of undiagnosed cases. Several computational methods have been developed to predict the functional impact of noncoding variants, but they exhibit only partial concordance and are not integrated with functional annotation resources, making the interpretation of these variants still challenging. MicroRNAs (miRNAs) are small noncoding RNA molecules that act as fine regulators of gene expression and play crucial functions in several biological processes, such as cell proliferation and differentiation. An increasing number of studies demonstrate a significant impact of miRNA single nucleotide variants (SNVs) both in Mendelian diseases and complex traits. To predict the functional effect of miRNA SNVs, we implemented a new meta‐predictor, MiRLog, and we integrated it into a comprehensive database, dbmiR, which includes a precompiled list of all possible miRNA allelic SNVs, providing their biological annotations at nucleotide and miRNA levels. MiRLog and dbmiR were used to explore the genetic variability of miRNAs in 15,708 human genomes included in the gnomAD project, finding several ultra‐rare SNVs with a potentially deleterious effect on miRNA biogenesis and function representing putative contributors to human phenotypes.

Genetic Diversity and Selection Footprints in the Genome of Brazilian Soybean Cultivars

Although Brazil is currently the largest soybean producer in the world, only a small number of studies have analyzed the genetic diversity of Brazilian soybean. These studies have shown the existence of a narrow genetic base. The objectives of this work were to analyze the population structure and genetic diversity, and to identify selection signatures in the genome of soybean germplasms from different companies in Brazil. A panel consisting of 343 soybean lines from Brazil, North America, and Asia was genotyped using genotyping by sequencing (GBS). Population structure was assessed by Bayesian and multivariate approaches. Genetic diversity was analyzed using metrics such as the fixation index, nucleotide diversity, genetic dissimilarity, and linkage disequilibrium. The software BayeScan was used to detect selection signatures between Brazilian and Asian accessions as well as among Brazilian germplasms. Region of origin, company of origin, and relative maturity group (RMG) all had a significant influence on population structure. Varieties belonging to the same company and especially to the same RMG exhibited a high level of genetic similarity. This result was exacerbated among early maturing accessions. Brazilian soybean showed significantly lower genetic diversity when compared to Asian accessions. This was expected, because the crop's region of origin is its main genetic diversity reserve. We identified 7 genomic regions under selection between the Brazilian and Asian accessions, and 27 among Brazilian varieties developed by different companies. Associated with these genomic regions, we found 96 quantitative trait loci (QTLs) for important soybean breeding traits such as flowering, maturity, plant architecture, productivity components, pathogen resistance, and seed composition. Some of the QTLs associated with the markers under selection have genes of great importance to soybean's regional adaptation. The results reported herein allowed to expand the knowledge about the organization of the genetic variability of the Brazilian soybean germplasm. Furthermore, it was possible to identify genomic regions under selection possibly associated with the adaptation of soybean to Brazilian environments.

Molecular diversity of exotic durian (Durio spp.) germplasm: a case study of Kalimantan, Indonesia

Background

Durian of Indonesia, specifically Durio zibethinus, is a potential agricultural commodity for domestic and international markets. However, its quality is still less competitive or significantly lower to fulfill the export market, compared to a similar one from other countries. This study aimed to determine and analyze the genetic diversity and relationship of the exotic durian (Durio spp.) germplasm originally from Kalimantan, Indonesia, using the rbcL marker.

Results

Based on this marker, the durian germplasm has a low genetic diversity (π%=0.24). It may strongly correspond with the variability sites or mutation present in the region. In this case, the rbcL region of the durian germplasm has generated 23 variable sites with a transition/transversion (Ti/Tv) bias value of 1.00. However, following the phylogenetic and principal component analyses, this germplasm is separated into four main clades and six groups, respectively. In this case, D. zibethinus was very closely related to D. exleyanus. Meanwhile, D. lowianus and D. excelsus were the farthest. In further analysis, 29 durians were very closely related, and the farthest was shown by Durian Burung (D. acutifolius) and Kalih Haliyang (D. kutejensis) as well as Pampaken Burung Kecil (D. kutejensis) and Durian Burung (D. acutifolius) with a divergence coefficient of 0.011. The Pearson correlation analysis confirms that 20 pairs of individual durians have a strong relation, shown by, e.g., Maharawin Hamak and Durian Burung as well as Mantuala Batu Hayam and Durian Burung Besar.

Conclusion

While the durian has a low genetic diversity, the phylogenetic analyses revealed that this germplasm originally from Kalimantan, Indonesia, shows unique relationships. These findings may provide a beneficial task in supporting the durian genetic conservation and breeding practices in the future, locally and globally.

Genome-wide diversity analysis to infer population structure and linkage disequilibrium among Colombian coconut germplasm

Genetic diversity and relatedness of accessions for coconut growing in Colombia was unknown until this study. Here we develop single nucleotide polymorphisms (SNPs) along the coconut genome based on Genotyping by Sequencing (GBS) with the goal of analyze the genetic diversity, population structure, and linkage disequilibrium (LD) of a diverse coconut panel consisting of 112 coconut accessions from the Atlantic and Pacific coasts of Colombia. A comprehensive catalog of approximately 40,000 SNPs with a minor allele frequency (MAF) of > 0.05 is presented. A total of 40,614 SNPs were found but only 19,414 anchored to chromosomes. Of these, 10,338 and 4606 were exclusive to the Atlantic and Pacific gene pools, respectively, and 3432 SNPs could differentiate both gene pools. A filtered subset of unlinked and anchored SNPs (1271) showed a population structure at K = 4, separating accessions from the Pacific and Atlantic coasts that can also be distinguished by palm height, as found in previous studies. The Pacific groups had a slow LD decay, low Fixation Index (Fst) and low nucleotide diversity (π), while the Atlantic group had slightly higher genetic diversity and faster LD decay. Genome-wide diversity analyses are of importance to promote germplasm conservation and breeding programs aimed at developing new cultivars better adapted to the region.

Population structure analysis to explore genetic diversity and geographical distribution characteristics of cultivated-type tea plant in Guizhou Plateau

BACKGROUND

Tea plants originated in southwestern China. Guizhou Plateau is an original center of tea plants, and is rich in germplasm resources. However, the genetic diversity, population structure and distribution characteristics of cultivated-type tea plants in the region are unknown. In this study, we explored the genetic diversity and geographical distribution of cultivated-type tea accessions in Guizhou Plateau.

RESULTS

We used 112,072 high-quality genotyping-by-sequencing to analyze the genetic diversity, principal components, phylogeny, population structure, and linkage disequilibrium, and develop a core collection of 253 cultivated-type tea plant accessions from Guizhou Plateau. The results showed Genetic diversity of the cultivated-type tea accessions of the Pearl River Basin was significantly higher than that of the cultivated-type tea accessions of the Yangtze River Basin. Three inferred pure groups (CG-1, CG-2 and CG-3) and one inferred admixture group (CG-4), were identified by a population structure analysis, and verified by principal component and phylogenetic analyses. The highest genetic distance and differentiation coefficients were determined for CG-2 vs CG-3. The lower genetic distance and differentiation coefficients were determined for CG-4 vs CG-2 and CG-4 vs CG-3, respectively. We developed a core set and a primary set. The primary and core sets contained 77.0 and 33.6% of all individuals in the initial set, respectively. The primary set may serve as the primary population in genome-wide association studies, while the core collection may serve as the core population in multiple treatment setting studies.

CONCLUSIONS

The present study demonstrated the genetic diversity and geographical distribution characteristics of cultivated-type tea plants in Guizhou Plateau. Significant differences in genetic diversity and evolutionary direction were detected between the ancient landraces of the Pearl River Basin and the those of the Yangtze River Basin. Major rivers and ancient hubs were largely responsible for the genetic exchange between the Pearl River Basin and the Yangtze River Basin ancient landraces as well as the formation of the ancient hubs evolutionary group. Genetic diversity, population structure and core collection elucidated by this study will facilitate further genetic studies, germplasm protection, and breeding of tea plants.

Insights on Genetic Diversity, Population Structure, and Linkage Disequilibrium in Globally Diverse Coconut Accessions Using Genotyping-by-Sequencing

Genotyping-by-sequencing (GBS) has emerged as a cost-effective approach for genome-wide discovery of single-nucleotide polymorphism (SNP) markers and high-throughput genotyping. In this study, 96 coconut palms, representing 16 accessions from globally diverse origins, were genotyped using the GBS strategy. A total of 10,835 high-quality SNPs, which were identified after stringent filtering, were utilized to assess genetic diversity, population structure, and linkage disequilibrium (LD) analyses. The polymorphism information content (PIC) values of SNPs ranged from 0.1 to 0.4, with a large proportion of SNPs (8633 nos.; 79.7%) having a higher PIC in the range of 0.3-0.4. The genetic diversity analysis revealed the existence of a high level of variation in coconut accessions, with an average expected heterozygosity (H_e) value of 0.43. Unweighted neighbor-joining phylogenetic tree and Bayesian-based model population structure grouped coconut genotypes into four main clusters. The accessions are generally clustered based on their height (tall or dwarf), with a few accession clusterings based on geographical origins. Investigation of LD pattern in coconut indicated a relatively rapid LD decay with a short range (9 kb). The results obtained in this study will contribute to enhancing the capacity of coconut researchers to utilize genetic diversity for further genetic improvement. In addition, it would open up possibilities for performing genomic studies such as genome-wide association studies and genomic selection to accelerate the efficiency and speed of coconut genetic improvement.

Ciprofloxacin induced antibiotic resistance in Salmonella Typhimurium mutants and genome analysis

Antibiotic resistance of Salmonella species is well reported. Ciprofloxacin is the frontline antibiotic for salmonellosis. The repeated exposure to ciprofloxacin leads to resistant strains. After 20 cycles of antibiotic exposure, resistant bacterial clones were evaluated. The colony size of the mutants was small and had an extended lag phase compared to parent strain. The whole genome sequencing showed 40,513 mutations across the genome. Small percentage (5.2%) of mutations was non-synonymous. Four-fold more transitions were observed than transversions. Ratio of < 1 transition vs transversion showed a positive selection for antibiotic resistant trait. Mutation distribution across the genome was uniform. The native plasmid was an exception and 2 mutations were observed on 90 kb plasmid. The important genes like dnaE, gyrA, iroC, metH and rpoB involved in antibiotic resistance had point mutations. The genome analysis revealed most of the metabolic pathways were affected.

Mutability of mononucleotide repeats, not oxidative stress, explains the discrepancy between laboratory-accumulated mutations and the natural allele-frequency spectrum in C. elegans

Important clues about natural selection can be gleaned from discrepancies between the properties of segregating genetic variants and of mutations accumulated experimentally under minimal selection, provided the mutational process is the same in the laboratory as in nature. The base-substitution spectrum differs between C. elegans laboratory mutation accumulation (MA) experiments and the standing site-frequency spectrum, which has been argued to be in part owing to increased oxidative stress in the laboratory environment. Using genome sequence data from C. elegans MA lines carrying a mutation (mev-1) that increases the cellular titer of reactive oxygen species (ROS), leading to increased oxidative stress, we find the base-substitution spectrum is similar between mev-1, its wild-type progenitor (N2), and another set of MA lines derived from a different wild strain (PB306). Conversely, the rate of short insertions is greater in mev-1, consistent with studies in other organisms in which environmental stress increased the rate of insertion–deletion mutations. Further, the mutational properties of mononucleotide repeats in all strains are different from those of nonmononucleotide sequence, both for indels and base-substitutions, and whereas the nonmononucleotide spectra are fairly similar between MA lines and wild isolates, the mononucleotide spectra are very different, with a greater frequency of A:T → T:A transversions and an increased proportion of ±1-bp indels. The discrepancy in mutational spectra between laboratory MA experiments and natural variation is likely owing to a consistent (but unknown) effect of the laboratory environment that manifests itself via different modes of mutability and/or repair at mononucleotide loci.

DNA barcoding of brackish and marine water fishes and shellfishes of Sundarbans, the world's largest mangrove ecosystem

The present study aims to apply a DNA barcoding tool through amplifying two mitochondrial candidate genes i.e., COI and 16S rRNA for accurate identification of fish, aquatic molluscs and crustaceans of Sundarbans mangrove wetland, to build a reference library of fish and shellfishes of this unique ecosystems. A total of 185 mitochondrial COI barcode sequences and 59 partial sequences of the 16S rRNA gene were obtained from 120 genera, 65 families and 21 orders of fish, crustaceans and molluscs. The collected samples were first identified by examining morphometric characteristics and then assessed by DNA barcoding. The COI and 16S rRNA sequences of fishes and crustaceans were clearly discriminated among genera in their phylogenies. The average Kimura two-parameter (K2P) distances of COI barcode sequences within species, genera, and families of fishes are 1.57±0.06%, 15.16±0.23%, and 17.79±0.02%, respectively, and for 16S rRNA sequences, these values are 1.74±.8%, 0.97±.8%, and 4.29±1.3%, respectively. The minimum and maximum K2P distance based divergences in COI sequences of fishes are 0.19% and 36.27%, respectively. In crustaceans, the K2P distances within genera, families, and orders are 1.4±0.03%, 17.73±0.15%, and 22.81±0.02%, respectively and the minimum and maximum divergences are 0.2% and 33.93%, respectively. Additionally, the present study resolves the misidentification of the mud crab species of the Sundarbans as Scylla olivacea which was previously stated as Scylla serrata. In case of molluscs, values of interspecific divergence ranges from 17.43% to 66.3% in the barcoded species. The present study describes the development of a molecular and morphometric cross-referenced inventory of fish and shellfish of the Sundarbans. This inventory will be useful in future biodiversity studies and in forming future conservation plan.

Placental mitochondrial DNA mutational load and perinatal outcomes: Findings from a multi-ethnic pregnancy cohort

Mitochondria fuel placental activity, with mitochondrial dysfunction implicated in several perinatal complications. We investigated placental mtDNA mutational load using NextGen sequencing in relation to birthweight and gestational length among 358 mother-newborn pairs. We found that higher heteroplasmy, especially in the hypervariable displacement loop region, was associated with shorter gestational length. Results were similar among male and female pregnancies, but stronger in magnitude among females. With regard to growth, we observed that higher mutational load was associated with lower birthweight-for-gestational age (BWGA) among females, but higher BWGA among males. These findings support potential sex-differential fetal biological strategies for coping with increased heteroplasmies.

Development of GBTS and KASP Panels for Genetic Diversity, Population Structure, and Fingerprinting of a Large Collection of Broccoli (Brassica oleracea L. var. italica) in China

Broccoli (Brassica oleracea var. italica) is one of the most important and nutritious vegetables widely cultivated in China. In the recent four decades, several improved varieties were bred and developed by Chinese breeders. However, the efforts for improvement of broccoli are hindered by limited information of genetic diversity and genetic relatedness contained within the available germplasms. This study evaluated the genetic diversity, genetic relationship, population structure, and fingerprinting of 372 accessions of broccoli representing most of the variability of broccoli in China. Millions of SNPs were identified by whole-genome sequencing of 23 representative broccoli genotypes. Through several stringent selection criteria, a total of 1,167 SNPs were selected to characterize genetic diversity and population structure. Of these markers, 1,067 SNPs were genotyped by target sequencing (GBTS), and 100 SNPs were genotyped by kompetitive allele specific PCR (KASP) assay. The average polymorphism information content (PIC) and expected heterozygosity (gene diversity) values were 0.33 and 0.42, respectively. Diversity analysis revealed the prevalence of low to moderate genetic diversity in the broccoli accessions indicating a narrow genetic base. Phylogenetic and principal component analyses revealed that the 372 accessions could be clustered into two main groups but with weak groupings. STRUCTURE analysis also suggested the presence of two subpopulations with weak genetic structure. Analysis of molecular variance (AMOVA) identified 13% variance among populations and 87% within populations revealing very low population differentiation, which could be attributed to massive gene flow and the reproductive biology of the crop. Based on high resolving power, a set of 28 KASP markers was chosen for DNA fingerprinting of the broccoli accessions for seed authentication and varietal identification. To the best of our knowledge, this is the first comprehensive study to measure diversity and population structure of a large collection of broccoli in China and also the first application of GBTS and KASP techniques in genetic characterization of broccoli. This work broadens the understanding of diversity, phylogeny, and population structure of a large collection of broccoli, which may enhance future breeding efforts to achieve higher productivity.

Population dynamics of GC-changing mutations in humans and great apes

The nucleotide composition of the genome is a balance between origin and fixation rates of different mutations. For example, it is well-known that transitions occur more frequently than transversions, particularly at CpG sites. Differences in fixation rates of mutation types are less explored. Specifically, recombination-associated GC-biased gene conversion (gBGC) may differentially impact GC-changing mutations, due to differences in their genomic distributions and efficiency of mismatch repair mechanisms. Given that recombination evolves rapidly across species, we explore gBGC of different mutation types across human populations and great ape species. We report a stronger correlation between segregating GC frequency and recombination for transitions than for transversions. Notably, CpG transitions are most strongly affected by gBGC in humans and chimpanzees. We show that the overall strength of gBGC is generally correlated with effective population sizes in humans, with some notable exceptions, such as a stronger effect of gBGC on non-CpG transitions in populations of European descent. Furthermore, species of the Gorilla and Pongo genus have a greatly reduced gBGC effect on CpG sites. We also study the dependence of gBGC dynamics on flanking nucleotides and show that some mutation types evolve in opposition to the gBGC expectation, likely due to hypermutability of specific nucleotide contexts. Our results highlight the importance of different gBGC dynamics experienced by GC-changing mutations and their impact on nucleotide composition evolution.

Single nucleotide polymorphisms affect miRNA target prediction in bovine

Single nucleotide polymorphisms (SNPs) can have significant effects on phenotypic characteristics in cattle. MicroRNAs (miRNAs) are small, non-coding RNAs that act as post-transcriptional regulators by binding them to target mRNAs. In the present study, we scanned ~56 million SNPs against 1,064 bovine miRNA sequences and analyzed, in silico, their possible effects on target binding prediction, primary miRNA formation, association with QTL regions and the evolutionary conservation for each SNP locus. Following target prediction, we show that 71.6% of miRNA predicted targets were altered as a consequence of SNPs located within the seed region of the mature miRNAs. Next, we identified variations in the Minimum Free Energy (MFE), which represents the capacity to alter molecule stability and, consequently, miRNA maturation. A total of 48.6% of the sequences analyzed showed values within those previously reported as sufficient to alter miRNA maturation. We have also found 131 SNPs in 46 miRNAs, with altered target prediction, occurring in QTL regions. Lastly, analysis of evolutionary conservation scores for each SNP locus suggested that they have a conserved biological function through the evolutionary process. Our results suggest that SNPs in microRNAs have the potential to affect bovine phenotypes and could be of great value for genetic improvement studies, as well as production.

Associations Between Maternal Lifetime Stress and Placental Mitochondrial DNA Mutations in an Urban Multiethnic Cohort

BACKGROUND

Disrupted placental functioning due to stress can have lifelong implications. Cumulative stress and trauma are likely to have lasting impacts on maternal physiological functioning and offspring development, resulting in increased risk for later-life complex disorders for which racial disparities exist.

METHODS

This study examined the association between maternal lifetime stress and placental mitochondrial DNA mutational load in an urban multiethnic cohort. Maternal lifetime exposure to stressful events was assessed using the validated Life Stressor Checklist-Revised. Whole mitochondrial DNA sequencing was performed and mutations were determined for 365 placenta samples with complete exposure and covariate data. Multivariable regression was used to model maternal lifetime stress in relation to placental mitochondrial DNA mutational load. Racial/ethnic differences were examined by cross-product terms and contrast statements. Gene-wise analyses were conducted.

RESULTS

We identified 13,189 heteroplasmies (Phred score > 10,000, minor allele frequency < 0.5, number of mutant reads > 1). Women experiencing increased psychosocial stress over their lifetime exhibited a higher number of total placental mitochondrial mutations (β = .23, 95% confidence interval = .03 to .42) and heteroplasmic mutations (β = .18, 95% confidence interval = .05 to .31) but not homoplasmic mutations (β = -.008, 95% confidence interval = -.03 to .01); the strongest associations were observed among Black women and genes coding for NADH dehydrogenase and cytochrome c oxidase subunits.

CONCLUSIONS

Cumulative maternal lifetime stress is associated with a greater mitochondrial mutational load, particularly among Black women. The impact of racial/ethnic differences in mutational load on placental function directly affecting offspring development and/or leading to chronic disease disparities warrants further investigation.

Virus Prospecting in Crickets-Discovery and Strain Divergence of a Novel Iflavirus in Wild and Cultivated Acheta domesticus

Orthopteran insects have high reproductive rates leading to boom-bust population dynamics with high local densities that are ideal for short, episodic disease epidemics. Viruses are particularly well suited for such host population dynamics, due to their supreme ability to adapt to changing transmission criteria. However, very little is known about the viruses of Orthopteran insects. Since Orthopterans are increasingly reared commercially, for animal feed and human consumption, there is a risk that viruses naturally associated with these insects can adapt to commercial rearing conditions, and cause disease. We therefore explored the virome of the house cricket Acheta domesticus, which is both part of the natural Swedish landscape and reared commercially for the pet feed market. Only 1% of the faecal RNA and DNA from wild-caught A. domesticus consisted of viruses. These included both known and novel viruses associated with crickets/insects, their bacterial-fungal microbiome, or their plant food. Relatively abundant among these viral Operational Taxonomic Units (OTUs) was a novel Iflavirus, tentatively named Acheta domesticus Iflavirus (AdIV). Quantitative analyses showed that AdIV was also abundant in frass and insect samples from commercially reared crickets. Interestingly, the wild and commercial AdIV strains had short, extremely divergent variation hotspots throughout the genome, which may indicate specific adaptation to their hosts’ distinct rearing environments.

Characterizing genomic variants and mutations in SARS-CoV-2 proteins from Indian isolates

SARS-CoV-2 is mutating and creating divergent variants by altering the composition of essential constituent proteins. Pharmacologically, it is crucial to understand the diverse mechanism of mutations for stable vaccine or anti-viral drug design. Our current study concentrates on all the constituent proteins of 469 SARS-CoV-2 genome samples, derived from Indian patients. However, the study may easily be extended to the samples across the globe.

We perform clustering analysis towards identifying unique variants in each of the SARS-CoV-2 proteins. A total of 536 mutated positions within the coding regions of SARS-CoV-2 proteins are detected among the identified variants from Indian isolates. We quantify mutations by focusing on the unique variants of each SARS-CoV-2 protein. We report the average number of mutation per variant, percentage of mutated positions, synonymous and non-synonymous mutations, mutations occurring in three codon positions and so on. Our study reveals the most susceptible six (06) proteins, which are ORF1ab, Spike (S), Nucleocapsid (N), ORF3a, ORF7a, and ORF8. Several non-synonymous substitutions are observed to be unique in different SARS-CoV-2 proteins. A total of 57 possible deleterious amino acid substitutions are predicted, which may impact on the protein functions. Several mutations show a large decrease in protein stability and are observed in putative functional domains of the proteins that might have some role in disease pathogenesis. We observe a good number of physicochemical property change during above deleterious substitutions.

Population structure, genetic diversity and genomic selection signatures among a Brazilian common bean germplasm

Brazil is the world's largest producer of common bean. Knowledge of the genetic diversity and relatedness of accessions adapted to Brazilian conditions is of great importance for the conservation of germplasm and for directing breeding programs aimed at the development of new cultivars. In this context, the objective of this study was to analyze the genetic diversity, population structure, and linkage disequilibrium (LD) of a diversity panel consisting of 219 common bean accessions, most of which belonging to the Mesoamerican gene pool. Genotyping by sequencing (GBS) of these accessions allowed the identification of 49,817 SNPs with minor allele frequency > 0.05. Of these, 17,149 and 12,876 were exclusive to the Mesoamerican and Andean pools, respectively, and 11,805 SNPs could differentiate the two gene pools. Further the separation according to the gene pool, bayesian analysis of the population structure showed a subdivision of the Mesoamerican accessions based on the origin and color of the seed tegument. LD analysis revealed the occurrence of long linkage blocks and low LD decay with physical distance between SNPs (LD half decay in 249 kb, corrected for population structure and relatedness). The GBS technique could effectively characterize the Brazilian common bean germplasms, and the diversity panel used in this study may be of great use in future genome-wide association studies.

Analysis of human mitochondrial genome co-occurrence networks of Asian population at varying altitudes

Networks have been established as an extremely powerful framework to understand and predict the behavior of many large-scale complex systems. We studied network motifs, the basic structural elements of networks, to describe the possible role of co-occurrence of genomic variations behind high altitude adaptation in the Asian human population. Mitochondrial DNA (mtDNA) variations have been acclaimed as one of the key players in understanding the biological mechanisms behind adaptation to extreme conditions. To explore the cumulative effects of variations in the mitochondrial genome with the variation in the altitude, we investigated human mt-DNA sequences from the NCBI database at different altitudes under the co-occurrence motifs framework. Analysis of the co-occurrence motifs using similarity clustering revealed a clear distinction between lower and higher altitude regions. In addition, the previously known high altitude markers 3394 and 7697 (which are definitive sites of haplogroup M9a1a1c1b) were found to co-occur within their own gene complexes indicating the impact of intra-genic constraint on co-evolution of nucleotides. Furthermore, an ancestral 'RSRS50' variant 10,398 was found to co-occur only at higher altitudes supporting the fact that a separate route of colonization at these altitudes might have taken place. Overall, our analysis revealed the presence of co-occurrence interactions specific to high altitude at a whole mitochondrial genome level. This study, combined with the classical haplogroups analysis is useful in understanding the role of co-occurrence of mitochondrial variations in high altitude adaptation.

Identification of nucleotide sequences and some proteins involved in polyadenylation of RNA transcribed by Pol III from SINEs

We have previously reported that not only transcripts of RNA polymerase II (pol II), but also one type of RNA transcribed by RNA polymerase III (pol III), undergo AAUAAA-dependent polyadenylation. Such an unusual feature is inherent in Short Interspersed Elements (SINEs) from genomes of certain mammals. For polyadenylation of its transcript, SINE should contain, besides an AATAAA hexamer and a transcription terminator, two specific regions: β, located downstream of box B of a promoter, and τ, preceding AATAAA. Here, using nucleotide substitutions in SINEs B2 (mouse) and Ves (bat), we identified nucleotides of β regions necessary for polyadenylation of their transcripts. These sequences (β signals) are the following: ACCACATgg in B2 and GGGCATGT in Ves. Using this approach, we identified τ signal of SINE B2 (GCTACagTGTACTTACAT), where TGTA tetramer is most important for polyadenylation. In Ves, τ region is a long polypyrimidine motif which is able to interact with PTB protein in Ves transcripts. We demonstrated by knockdown that B2 and Ves transcript polyadenylation is performed by canonical poly(A) polymerase with the participation of proteins CSPF-160 and Fip1, the known factors of mRNA polyadenylation. We also showed that a factor CFIm partaking in polyadenylation of many mRNAs, is involved only in polyadenylation of B2 transcripts. CFIm seems to interact with τ signal of В2 RNA and thereby facilitates the recruiting of other proteins engaged in polyadenylation. Thus, SINEs utilize at least some proteins involved in polyadenylation of pol II transcripts to polyadenylate their pol III transcripts.

Identification of the immune gene expression signature associated with recurrence of high-grade gliomas

High-grade gliomas (HGGs), the most common and aggressive primary brain tumors in adults, inevitably recur due to incomplete surgery or resistance to therapy. Intratumoral genomic and cellular heterogeneity of HGGs contributes to therapeutic resistance, recurrence, and poor clinical outcomes. Transcriptomic profiles of HGGs at recurrence have not been investigated in detail. Using targeted sequencing of cancer-related genes and transcriptomics, we identified single nucleotide variations, small insertions and deletions, copy number aberrations (CNAs), as well as gene expression changes and pathway deregulation in 16 pairs of primary and recurrent HGGs. Most of the somatic mutations identified in primary HGGs were not detected after relapse, suggesting a subclone substitution during the tumor progression. We found a novel frameshift insertion in the ZNF384 gene which may contribute to extracellular matrix remodeling. An inverse correlation of focal CNAs in EGFR and PTEN genes was detected. Transcriptomic analysis revealed downregulation of genes involved in messenger RNA splicing, cell cycle, and DNA repair, while genes related to interferon signaling and phosphatidylinositol (PI) metabolism are upregulated in secondary HGGs when compared to primary HGGs. In silico analysis of the tumor microenvironment identified M2 macrophages and immature dendritic cells as enriched in recurrent HGGs, suggesting a prominent immunosuppressive signature. Accumulation of those cells in recurrent HGGs was validated by immunostaining. Our findings point to a substantial transcriptomic deregulation and a pronounced infiltration of immature dendritic cells in recurrent HGG, which may impact the effectiveness of frontline immunotherapies in the GBM management. KEY MESSAGES: Most of the somatic mutations identified in primary HGGs were not detected after relapse. Focal CNAs in EGFR and PTEN genes are inversely correlated in primary and recurrent HGGs. Transcriptomic changes and distinct immune-related signatures characterize HGG recurrence. Recurrent HGGs are characterized by a prominent infiltration of immature dendritic and M2 macrophages.

Evolutionary Analysis of Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) Reveals Genomic Divergence with Implications for Universal Vaccine Efficacy

Coronavirus disease (COVID-19), caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2), is one of the pressing contemporary public health challenges. Investigations into the genomic structure of SARS-CoV-2 may inform ongoing vaccine development efforts and/or provide insights into vaccine efficacy to fight against COVID-19. Evolutionary analysis of 540 genomes spanning 20 different countries/territories was conducted and revealed an increase in the genomic divergence across successive generations. The ancestor of the phylogeny was found to be the isolate from the 2019/2020 Wuhan outbreak. Its transmission was outlined across 20 countries/territories as per genomic similarity. Our results demonstrate faster evolving variations in the genomic structure of SARS-CoV-2 when compared to the isolates from early stages of the pandemic. Genomic alterations were predominantly located and mapped onto the reported vaccine candidates of structural genes, which are the main targets for vaccine candidates. S protein showed 34, N protein 25, E protein 2, and M protein 3 amino acid variations in 246 genomes among 540. Among identified mutations, 23 in S protein, 1 in E, 2 from M, and 7 from N protein were mapped with the reported vaccine candidates explaining the possible implications on universal vaccines. Hence, potential target regions for vaccines would be ideally chosen from the structural regions of the genome that lack high variation. The increasing variations in the genome of SARS-CoV-2 together with our observations in structural genes have important implications for the efficacy of a successful universal vaccine against SARS-CoV-2.

Mutation bias interacts with composition bias to influence adaptive evolution

Mutation is a biased stochastic process, with some types of mutations occurring more frequently than others. Previous work has used synthetic genotype-phenotype landscapes to study how such mutation bias affects adaptive evolution. Here, we consider 746 empirical genotype-phenotype landscapes, each of which describes the binding affinity of target DNA sequences to a transcription factor, to study the influence of mutation bias on adaptive evolution of increased binding affinity. By using empirical genotype-phenotype landscapes, we need to make only few assumptions about landscape topography and about the DNA sequences that each landscape contains. The latter is particularly important because the set of sequences that a landscape contains determines the types of mutations that can occur along a mutational path to an adaptive peak. That is, landscapes can exhibit a composition bias—a statistical enrichment of a particular type of mutation relative to a null expectation, throughout an entire landscape or along particular mutational paths—that is independent of any bias in the mutation process. Our results reveal the way in which composition bias interacts with biases in the mutation process under different population genetic conditions, and how such interaction impacts fundamental properties of adaptive evolution, such as its predictability, as well as the evolution of genetic diversity and mutational robustness.

Mutation is often depicted as a random process due its unpredictable nature. However, such randomness does not imply uniformly distributed outcomes, because some DNA sequence changes happen more frequently than others. Mutation bias can be an orienting factor in adaptive evolution, influencing the mutational trajectories populations follow toward higher-fitness genotypes. Because these trajectories are typically just a small subset of all possible mutational trajectories, they can exhibit composition bias—an enrichment of a particular kind of DNA sequence change, such as transition or transversion mutations. Here, we use empirical data from eukaryotic transcriptional regulation to study how mutation bias and composition bias interact to influence adaptive evolution.

Placental mitochondrial DNA mutations and copy numbers in intrauterine growth restricted (IUGR) pregnancy

Intrauterine Growth Restriction (IUGR) is a common and significant complication that arises during pregnancy wherein the fetus fails to attain its full growth potential. Mitochondria being one of the primary sources of energy, plays an important role in placentation and fetal development. In IUGR pregnancy, increased oxidative stress due to inadequate oxygen and nutrient supply could possibly alter mitochondrial functions and homeostasis. In this study, we evaluated the biochemical and molecular changes in mitochondria as biosignature for early and better characterization of IUGR pregnancies. We identified significant increase in mtDNA copy number in both IUGR (p = 0.0001) and Small for Gestational Age (SGA) but healthy (p = 0.0005) placental samples when compared to control. Whole mitochondrial genome sequencing identified novel mutations in both coding and non-coding regions of mtDNA in multiple IUGR placental samples. Sirtuin-3 (Sirt3) protein expression was significantly downregulated (p = 0.027) in IUGR placenta but there was no significant difference in Nrf1 expression in IUGR when compared to control group. Our study provides an evidence for altered mitochondrial homeostasis and paves a way towards interrogating mitochondrial abnormalities in IUGR pregnancies.

Genetic diversity analysis of a flax (Linum usitatissimum L.) global collection

Background

A sustainable breeding program requires a minimum level of germplasm diversity to provide varied options for the selection of new breeding lines. To maximize genetic gain of the North Dakota State University (NDSU) flax breeding program, we aimed to increase the genetic diversity of its parental stocks by incorporating diverse genotypes. For this purpose, we analyzed the genetic diversity, linkage disequilibrium, and population sub-structure of 350 globally-distributed flax genotypes with 6200 SNP markers.

Results

All the genotypes tested clustered into seven sub-populations (P1 to P7) based on the admixture model and the output of neighbor-joining (NJ) tree analysis and principal coordinate analysis were in line with that of structure analysis. The largest sub-population separation arose from a cluster of NDSU/American genotypes with Turkish and Asian genotypes. All sub-populations showed moderate genetic diversity (average H = 0.22 and I = 0.34). The pairwise F_st comparison revealed a great degree of divergence (F_st > 0.25) between most of the combinations. A whole collection mantel test showed significant positive correlation (r = 0.30 and p < 0.01) between genetic and geographic distances, whereas it was non-significant for all sub-populations except P4 and P5 (r = 0.251, 0.349 respectively and p < 0.05). In the entire collection, the mean linkage disequilibrium was 0.03 and it decayed to its half maximum within < 21 kb distance.

Conclusions

To maximize genetic gain, hybridization between NDSU stock (P5) and Asian individuals (P6) are potentially the best option as genetic differentiation between them is highest (F_st > 0.50). In contrast, low genetic differentiation between P5 and P2 may enhance the accumulation of favorable alleles for oil and fiber upon crossing to develop dual purpose varieties. As each sub-population consists of many genotypes, a Neighbor-Joining tree and kinship matrix assist to identify distantly related genotypes. These results also inform genotyping decisions for future association mapping studies to ensure the identification of a sufficient number of molecular markers to tag all linkage blocks.

ycf1-ndhF genes, the most promising plastid genomic barcode, sheds light on phylogeny at low taxonomic levels in Prunus persica

Background

Chloroplast genome sequencing is becoming a valuable process for developing several DNA barcodes. At present, plastid DNA barcode for systematics and evolution in flowering plant rely heavily on the use of non-coding genes. The present study was performed to verify the novelty and suitability of the two hotspot barcode plastid coding gene ycf1 and ndhF, to estimate the rate of molecular evolution in the Prunus genus at low taxonomic levels.

Results

Here, 25 chloroplast genomes of Prunus genus were selected for sequences annotation to search for the highly variable coding DNA barcode regions. Among them, 5 genera were of our own data, including the ornamental, cultivated, and wild haplotype, while 20 genera have been downloaded from the GenBank database. The results indicated that the two hotspot plastid gene ycf1 and ndhF were the most variable regions within the coding genes in Prunus with an average of 3268 to 3416 bp in length, which have been predicted to have the highest nucleotide diversity, with the overall transition/transversion bias (R = 1.06). The ycf1-ndhF structural domains showed a positive trend evident in structure variation among the 25 specimens tested, due to the variant overlap^’s gene annotation and insertion or deletion with a broad trend of the full form of IGS sequence. As a result, the principal component analysis (PCA) and the ML tree data drew an accurate monophyletic annotations cluster in Prunus species, offering unambiguous identification without overlapping groups between peach, almond, and cherry.

Conclusion

To this end, we put forward the domain of the two-locus ycf1-ndhF genes as the most promising coding plastid DNA barcode in P. persica at low taxonomic levels. We believe that the discovering of further variable loci with high evolutionary rates is extremely useful and potential uses as a DNA barcode in P. persica for further phylogeny study and species identification.

Genetic diversity and population structure of Quercus fabri Hance in China revealed by genotyping-by-sequencing

Analysis of genetic diversity and population structure among Quercus fabri populations is essential for the conservation and utilization of Q. fabri resources. Here, the genetic diversity and structure of 158 individuals from 13 natural populations of Quercus fabri in China were analyzed using genotyping-by-sequencing (GBS). A total of 459,564 high-quality single nucleotide polymorphisms (SNPs) were obtained after filtration for subsequent analysis. Genetic structure analysis revealed that these individuals can be clustered into two groups and the structure can be explained mainly by the geographic barrier, showed gene introgression from coastal to inland areas and high mountains could significantly hinder the mutual introgression of genes. Genetic diversity analysis indicated that the individual differences within groups are greater than the differences between the two groups. These results will help us better understand the genetic backgrounds of Q. fabri.

Description of Sarcocystis scandentiborneensis sp. nov. from treeshrews (Tupaia minor, T. tana) in northern Borneo with annotations on the utility of COI and 18S rDNA sequences for species delineation

Sarcocystis scandentiborneensis sp. nov. was discovered in histological sections of striated musculature of treeshrews (Tupaia minor, T. tana) from Northern Borneo. Sarcocysts were cigar-shaped, 102 μm–545 μm long, and on average 53 μm in diameter. The striated cyst wall varied in thickness (2–10 μm), depending on whether the finger-like, villous protrusions (VP) were bent. Ultrastructurally, sarcocysts were similar to wall type 12 but basal microtubules extended into VPs that tapered off with a unique U-shaped, electron-dense apical structure. In phylogenetic trees of the nuclear 18S rRNA gene, S. scandentiborneensis formed a distinct branch within a monophyletic subclade of Sarcocystis spp. with (colubrid) snake-rodent life cycle. We mapped all intraspecific (two haplotypes) and interspecific nucleotide substitutions to the secondary structure of the 18S rRNA gene: in both cases, the highest variability occurred within helices V2 and V4 but intraspecific variability mostly related to transitions, while transition/transversion ratios between S. scandentiborneensis, S. zuoi, and S. clethrionomyelaphis were skewed towards transversions. Lack of relevant sequences restricted phylogenetic analysis of the mitochondrial Cytochrome C oxidase subunit I (COI) gene to include only one species of Sarcocystis recovered from a snake host (S. pantherophisi) with which the new species formed a sister relationship. We confirm the presence of the functionally important elements of the COI barcode amino acid sequence of S. scandentiborneensis, whereby the frequency of functionally important amino acids (Alanine, Serine) was markedly different to other taxa of the Sarcocystidae. We regard S. scandentiborneensis a new species, highlighting that structurally or functionally important aspects of the 18S rRNA and COI could expand their utility for delineation of species. We also address the question why treeshrews, believed to be close to primates, carry a parasite that is genetically close to a Sarcocystis lineage preferably developing in the Rodentia as intermediate hosts.

•

Sarcocystis scandentiborneensis sp. nov. identified from Tupaia tana and T. minor.

•

Putative snake-treeshrew life cycle.

•

Unique brush-like cyst wall with an electron-dense apical structure.

•

Transition/transversion ratios of 18S rDNA improve species discrimination.

•

COI: Interspecific differences of functionally important amino acids of barcode area.

CRISPR C-to-G base editors for inducing targeted DNA transversions in human cells

CRISPR-guided DNA cytosine and adenine base editors (CBEs and ABEs) are widely used for many applications^1–4 but primarily create DNA base transitions (i.e., pyrimidine-to-pyrimidine, or purine-to-purine). Here we describe the engineering of two base editor architectures that can efficiently induce targeted C-to-G base transversions, with reduced levels of unwanted C-to-W (W = A or T) and indel mutations. One of these C-to-G base editors (CGBE1), consists of an RNA-guided Cas9 nickase, an E. coli-derived uracil DNA N-glycosylase (eUNG), and a rat APOBEC1 cytidine deaminase variant (R33A) previously shown to have reduced off-target RNA and DNA editing activities^{5, 6}. We show that CGBE1 can efficiently induce C-to-G edits, particularly in AT-rich sequence contexts in human cells. We also removed the eUNG domain to yield miniCGBE1, which reduced indel frequencies but only modestly decreased editing efficiency. CGBE1 and miniCGBE1 enable C-to-G edits and will serve as a basis for optimizing C-to-G base editors for research and therapeutic applications.

A new base editor enables the creation of C-to-G base changes in human cells.