Search

Antimicrobial resistance in urinary tract infections (UTIs) is a major public health concern. This study aims to characterize the phenotypic and genetic basis of multidrug resistance (MDR) among expanded-spectrum cephalosporin-resistant (ESCR) uropathogenic Escherichia coli (UPEC) causing UTIs in California patient populations. Between February and October 2019, 577 ESCR UPEC isolates were collected from patients at 6 clinical laboratory sites across California. Lineage and antibiotic resistance genes were determined by analysis of whole-genome sequence data. The lineages ST131, ST1193, ST648, and ST69 were predominant, representing 46%, 5.5%, 4.5%, and 4.5% of the collection, respectively. Overall, 527 (91%) isolates had an expanded-spectrum β-lactamase (ESBL) phenotype, with blaCTX-M-15, blaCTX-M-27, blaCTX-M-55, and blaCTX-M-14 being the most prevalent ESBL genes. In the 50 non-ESBL phenotype isolates, 40 (62%) contained blaCMY-2, which was the predominant plasmid-mediated AmpC (pAmpC) gene. Narrow-spectrum β-lactamases, blaTEM-1B and blaOXA-1, were also found in 44.9% and 32.1% of isolates, respectively. Among ESCR UPEC isolates, isolates with an ESBL phenotype had a 1.7-times-greater likelihood of being MDR than non-ESBL phenotype isolates (P < 0.001). The cooccurrence of blaCTX-M-15, blaOXA-1, and aac(6')-Ib-cr within ESCR UPEC isolates was strongly correlated. Cooccurrence of blaCTX-M-15, blaOXA-1, and aac(6')-Ib-cr was associated with an increased risk of nonsusceptibility to piperacillin-tazobactam, cefepime, fluoroquinolones, and amikacin as well as MDR. Multivariate regression revealed the presence of blaCTX-M-55, blaTEM-1B, and the ST131 genotype as predictors of MDR. IMPORTANCE The rising incidence of resistance to expanded-spectrum cephalosporins among Escherichia coli strains, the most common cause of UTIs, is threatening our ability to successfully empirically treat these infections. ESCR E. coli strains are often MDR; therefore, UTI caused by these organisms often leads to treatment failure, increased length of hospital stay, and severe complications (D. G. Mark, Y.-Y. Hung, Z. Salim, N. J. Tarlton, et al., Ann Emerg Med 78:357-369, 2021, https://doi.org/10.1016/j.annemergmed.2021.01.003). Here, we performed an in-depth analysis of genetic factors of ESCR E. coli associated with coresistance and MDR. Such knowledge is critical to advance UTI diagnosis, treatment, and antibiotic stewardship.

HoBi-like pestivirus (HoBiPeV), classified under Pestivirus H species, is an emerging cattle pathogen of high economic impact. However, the origin and evolution of HoBiPeV are not very clear due to a lack of full genomic sequences from diverse clades. This study aimed to determine full-genome sequences of HoBiPeV strains of three novel clades (c, d and e) and perform full-genome-based genetic and evolutionary analyses. Bayesian phylogenetic analyses herein confirmed the existence and independent evolution of four main HoBiPeV clades (a, c, d and e) globally, with genetic divergence ranging from 13.0% to 18.2%. Our Bayesian molecular clock estimates revealed that HoBiPeV most likely originated in India, with a dated tMRCA of 1938 (1762–2000), evidencing a more recent origin of HoBiPeV. The evolution rate of HoBiPeV was estimated to be 2.133 × 10−3 subs/site/year at full-genome level but varied widely among individual genes. Selection pressure analyses identified most of the positively selected sites in E2. Additionally, 21.8% of the ORF codon sites were found under strong episodic diversifying selection, providing first evidence of negative selection in HoBiPeV evolution. No recombination event was evident for HoBiPeV-c, d and e strains. These findings provide new insights into HoBiPeV origin and evolutionary history for better understanding the epidemiology and host–pathogen interactions and stimulate vaccine research.

Salmonella enterica is, globally, an important cause of human illness with beef being a significant attributable source. In the human patient, systemic Salmonella infection requires antibiotic therapy, and when strains are multidrug resistant (MDR), no effective treatment may be available. MDR in bacteria is often associated with the presence of mobile genetic elements (MGE) that mediate horizontal spread of antimicrobial resistance (AMR) genes. In this study, we sought to determine the potential relationship of MDR in bovine Salmonella isolates with MGE. The present study involved 111 bovine Salmonella isolates obtained collectively from specimens derived from healthy cattle or their environments at Midwestern U.S. feedyards (2000–2001, n = 19), or specimens from sick cattle submitted to the Nebraska Veterinary Diagnostic Center (2010–2020, n = 92). Phenotypically, 33/111 isolates (29.7%) were MDR (resistant to ≥3 drug classes). Based on whole-genome sequencing (WGS; n = 41) and PCR (n = 111), a MDR phenotype was strongly associated (OR = 186; p < 0.0001) with carriage of ISVsa3, an IS91-like Family transposase. In all 41 isolates analyzed by WGS ((31 MDR and 10 non-MDR (resistant to 0–2 antibiotic classes)), MDR genes were associated with carriage of ISVsa3, most often on an IncC type plasmid carrying blaCMY-2. The typical arrangement was floR, tet(A), aph(6)-Id, aph(3″)-Ib, and sul2 flanked by ISVsa3. These results suggest that AMR genes in MDR S. enterica isolates of cattle are frequently associated with ISVsa3 and carried on IncC plasmids. Further research is needed to better understand the role of ISVsa3 in dissemination of MDR Salmonella strains.

To prevent invasive pneumococcal disease (IPD), pneumococcal conjugate vaccines (PCVs) have been implemented in many countries; however, many cases of IPD still occur and can be attributable to nonvaccine serotypes of Streptococcus pneumoniae. In Japan, the number of IPD cases attributable to serotype 12F increased from 4.4% in 2015 to 24.6% in 2017 after 13-valent PCV was introduced. To clarify the associated genetic characteristics, we conducted whole-genome sequencing of 75 serotype 12F isolates. We identified 2 sequence types (STs) among the isolates: ST4846, which was the major type, and ST6945. Bayesian analysis suggested that these types diverged in »1942. Among serotype 12F-ST4846, we identified a major cluster, PC-JP12F, whose time of most recent common ancestor was estimated to be »2012. A phylogeographic analysis demonstrated that PC-JP12F isolates spread from the Kanto region, the most populated region in Japan, to other local regions.

Significance Genetic investigations of most structural birth defects, including spina bifida (SB), congenital heart disease, and craniofacial anomalies, have been underpowered for genome-wide association studies because of their rarity, genetic heterogeneity, incomplete penetrance, and environmental influences. Our systems biology strategy to investigate SB predisposition controls for population stratification and avoids much of the bias inherent in candidate gene searches that are pervasive in the field. We examine both protein coding and noncoding regions of whole genomes to analyze sequence variants, collapsed by gene or regulatory region, and apply machine learning, gene enrichment, and pathway analyses to elucidate molecular pathways and genes contributing to human SB., Spina bifida (SB) is a debilitating birth defect caused by multiple gene and environment interactions. Though SB shows non-Mendelian inheritance, genetic factors contribute to an estimated 70% of cases. Nevertheless, identifying human mutations conferring SB risk is challenging due to its relative rarity, genetic heterogeneity, incomplete penetrance, and environmental influences that hamper genome-wide association studies approaches to untargeted discovery. Thus, SB genetic studies may suffer from population substructure and/or selection bias introduced by typical candidate gene searches. We report a population based, ancestry-matched whole-genome sequence analysis of SB genetic predisposition using a systems biology strategy to interrogate 298 case-control subject genomes (149 pairs). Genes that were enriched in likely gene disrupting (LGD), rare protein-coding variants were subjected to machine learning analysis to identify genes in which LGD variants occur with a different frequency in cases versus controls and so discriminate between these groups. Those genes with high discriminatory potential for SB significantly enriched pathways pertaining to carbon metabolism, inflammation, innate immunity, cytoskeletal regulation, and essential transcriptional regulation consistent with their having impact on the pathogenesis of human SB. Additionally, an interrogation of conserved noncoding sequences identified robust variant enrichment in regulatory regions of several transcription factors critical to embryonic development. This genome-wide perspective offers an effective approach to the interrogation of coding and noncoding sequence variant contributions to rare complex genetic disorders.

The true bug, Aspongopus chinensis Dallas, 1851 (Hemiptera: Dinidoridae), is a fascinating insect with prolonged diapause and medicinal properties but also a notorious pest. However, because of the lack of genomic resources, an in-depth understanding of its biological characteristics is lacking. Here, we report the first genome assembly of A. chinensis anchored to 10 pseudochromosomes, which was achieved by combining PacBio long reads and Hi-C sequencing data. This chromosome-level genome assembly was 1.55 Gb in size with a scaffold N50 of 156 Mb. The benchmarking universal single-copy ortholog (BUSCO) analysis of the assembly captured 96.6% of the BUSCO genes. A total of 686,888,052 bp of repeat sequences, 18,511 protein-coding genes, and 1,749 noncoding RNAs were annotated. By comparing the A. chinensis genome with that of 8 homologous insects and 2 model organisms, 213 rapidly evolving gene families were identified, including 83 expanded and 130 contracted gene families. The functional enrichment of Gene Ontology and KEGG pathways showed that the significantly expanded gene families were primarily involved in metabolism, immunity, detoxification, and DNA/RNA replication associated with stress responses. The data reported here shed light on the ecological adaptation of A. chinensis and further expanded our understanding of true bug evolution in general.

Worldwide spread of Enterobacteriaceae resistant to colistin, a polypeptide antibacterial drug for last-resort treatment of carbapenemase-producing Enterobacteriaceae (CPE) infections, is concerning. This study aimed to elucidate colistin MICs and molecular characteristics of mcr-1 to mcr-9 of ESBL-producing Escherichia coli (ESBL-Ec) and CPE in Japan and clarify the genomic structure of strains harboring mcr genes (especially mcr-9). This study included 168 ESBL-Ec and 126 CPE strains isolated at Japanese medical facilities. Colistin susceptibility testing and multiplex PCR targeting mcr-1 to mcr-9 were performed for all strains with S1-nuclease pulsed-field gel electrophoresis, Southern blot hybridization, and whole-genome sequencing (WGS) with hybrid assembly performed for mcr gene-carrying strains. Two CPE strains showed a MIC ≥ 4 μg/ml in colistin susceptibility testing, with no known resistance mechanism detected. However, PCR conducted on all target strains detected three mcr-9-carrying strains showing colistin susceptibility. The blaCTX–M–62-positive E. coli THUN648 strain simultaneously carried blaCTX–M–62 and mcr-9 on a 275-kbp plasmid. Besides, blaIMP–6 + blaCTX–M–2-positive Klebsiella pneumoniae THUN262 and blaGES–24-positive Enterobacter kobei THUN627 had mcr-9 encoded on the chromosome. Only THUN627 encoded qseB/C, which is suggested to be a regulatory gene for mcr-9, downstream of mcr-9. However, this strain showed no increased expression of these genes in mRNA quantitative analysis under colistin exposure. Colistin MICs of ESBL-Ec and CPE in Japan were all below 2 μg/ml, which is below the epidemiological cutoff (ECOFF) value (https://eucast.org/) or clinical breakpoint (CB) (CLSI M100-S30) reported for colistin, indicating neither “microbiological” nor “clinical” resistance. Several colistin-susceptible Enterobacteriaceae carrying silent mcr-9 encoded on plasmids and chromosomes have already spread worldwide along with other antimicrobial resistance genes. However, the mechanism of colistin resistance by mcr-9 remains unclear.

In this study, we investigate the genetic determinants that underlie epilepsy in a captive baboon pedigree and evaluate the potential suitability of this non-human primate model for understanding the genetic etiology of human epilepsy. Archived whole-genome sequence data were analyzed using both a candidate gene approach that targeted variants in baboon homologs of 19 genes (n = 20,881 SNPs) previously implicated in genetic generalized epilepsy (GGE) and a more agnostic approach that examined protein-altering mutations genome-wide as assessed by snpEff (n = 36,169). Measured genotype association tests for baboon cases of epileptic seizure were performed using SOLAR, as well as gene set enrichment analyses (GSEA) and protein–protein interaction (PPI) network construction of top association hits genome-wide (p < 0.01; n = 441 genes). The maximum likelihood estimate of heritability for epileptic seizure in the pedigreed baboon sample is 0.76 (SE = 0.77; p = 0.07). Among candidate genes for GGE, a significant association was detected for an intronic SNP in RBFOX1 (p = 5.92 × 10–6; adjusted p = 0.016). For protein-altering variants, no genome-wide significant results were observed for epilepsy status. However, GSEA revealed significant positive enrichment for genes involved in the extracellular matrix structure (ECM; FDR = 0.0072) and collagen formation (FDR = 0.017), which was reflected in a major PPI network cluster. This preliminary study highlights the potential role of RBFOX1 in the epileptic baboon, a protein involved in transcriptomic regulation of multiple epilepsy candidate genes in humans and itself previously implicated in human epilepsy, both focal and generalized. Moreover, protein-damaging variants from across the genome exhibit a pattern of association that links collagen-containing ECM to epilepsy risk. These findings suggest a shared genetic etiology between baboon and human forms of GGE and lay the foundation for follow-up research.

Coleoptera is the most species-rich insect order, yet is currently underrepresented in genomic databases. An assembly was generated for ca. 1.7-Gb genome of the leaf beetle Gonioctena quinquepunctata by first assembling long-sequence reads (Oxford Nanopore; ± 27-fold coverage) and subsequently polishing the resulting assembly with short sequence reads (Illumina; ± 85-fold coverage). The unusually large size (most Coleoptera species are associated with a reported size below 1 Gb) was at least partially attributed to the presence of a large fraction of repeated elements (73.8%). The final assembly was characterized by an N50 length of 432 kb and a BUSCO score of 95.5%. The heterozygosity rate was ±0.6%. Automated genome annotation informed by RNA-Seq resulted in 40,568 predicted proteins, which is much larger than the typical range 17,000-23,000 predicted for other Coleoptera. However, no evidence of a genome duplication was detected. This new reference genome will contribute to our understanding of genetic variation in the Coleoptera. Among others, it will also allow exploring reproductive barriers between species, investigating introgression in the nuclear genome, and identifying genes involved in resistance to extreme climate conditions., info:eu-repo/semantics/published

Using whole-genome sequence (WGS) data are supposed to be optimal for genome-wide association studies and genomic predictions. However, sequencing thousands of individuals of interest is expensive. Imputation from single nucleotide polymorphisms panels to WGS data is an attractive approach to obtain highly reliable WGS data at low cost. Here, we conducted a genotype imputation study with a combined reference panel in yellow-feather dwarf broiler population. The combined reference panel was assembled by sequencing 24 key individuals of a yellow-feather dwarf broiler population (internal reference panel) and WGS data from 311 chickens in public databases (external reference panel). Three scenarios were investigated to determine how different factors affect the accuracy of imputation from 600 K array data to WGS data, including: genotype imputation with internal, external and combined reference panels; the number of internal reference individuals in the combined reference panel; and different reference sizes and selection strategies of an external reference panel. Results showed that imputation accuracy from 600 K to WGS data were 0.834±0.012, 0.920±0.007 and 0.982±0.003 for the internal, external and combined reference panels, respectively. Increasing the reference size from 50 to 250 improved the accuracy of genotype imputation from 0.848 to 0.974 for the combined reference panel and from 0.647 to 0.917 for the external reference panel. The selection strategies for the external reference panel had no impact on the accuracy of imputation using the combined reference panel. However, if only an external reference panel with reference sizegt;50 was used, the selection strategy of minimizing the average distance to the closest leaf had the greatest imputation accuracy compared with other methods. Generally, using a combined reference panel provided greater imputation accuracy, especially for low-frequency variants. In conclusion, the optimal imputation strategy with a combined reference panel should comprehensively consider genetic diversity of the study population, availability and properties of external reference panels, sequencing and computing costs, and frequency of imputed variants. This work sheds light on how to design and execute genotype imputation with a combined external reference panel in a livestock population.

Here, we report the first imported case of melioidosis from Laos in central China. COMPACT VITEK2 identification system and PCR, as well as sequencing methods confirmed that the patient was infected by Burkholderia pseudomallei, a bacterial species closely related to an isolate detected in Thailand. These findings are highly valuable for an early diagnosis, treatment and to prevent the spread of this emerging infectious disease in central China.

Endophytic fungi are one of the most prolific sources of functional biomolecules with therapeutic potential. Besides playing an important role in serious plant diseases, Fusarium strains possess the powerful capability to produce a diverse array of bioactive secondary metabolites (SMs). In order to in-depth mine gene clusters for SM biosynthesis of the genus Fusarium, an endophytic strain Fusarium sp. R1 isolated from Rumex madaio Makino was extensively investigated by whole-genome sequencing and in-depth bioinformatic analysis, as well as antiSMASH annotation. The results displayed that strain R1 harbors a total of 51.8 Mb genome, which consists of 542 contigs with an N50 scaffold length of 3.21 Mb and 50.4% GC content. Meanwhile, 19,333 functional protein-coding genes, 338 tRNA and 111 rRNA were comprehensively predicted and highly annotated using various BLAST databases including non-redundant (Nr) protein sequence, nucleotide (Nt) sequence, Swiss-Prot, Gene Ontology (GO), Kyoto Encyclopedia of Genes and Genomes (KEGG) and Clusters of Orthologous Groups (COG), as well as Pathogen Host Interactions (PHI) and Carbohydrate-Active enzymes (CAZy) databases. Antibiotics and Secondary Metabolites Analysis Shell (AntiSMASH) results showed that strain R1 has 37 SM biosynthetic gene clusters (BGCs), including 17 nonribosomal peptide synthetases (NRPSs), 13 polyketide synthetases (PKSs), 3 terpene synthases (Ts), 3 hybrid NRPS + PKS and 1 hybrid indole + NRPS. These findings improve our knowledge of the molecular biology of the genus Fusarium and would promote the discovery of new bioactive SMs from strain R1 using gene mining strategies including gene knockout and heteroexpression.

Summary In this article, we analyzed pedigree information on males from 12 bovine breeds born in France between 2015 and 2019. We report an overall small number of paternal lineages with, for example, a minimal number of ancestors accounting for 95% of the Y‐chromosome pool of their breed ranging from only 2 to 15 individuals. Then, we mined whole‐genome sequence data from 811 sires (2 ≤ n ≤ 510 per breed) and built a median‐joining network using 1411 SNPs. Most branches were breed‐specific and in agreement with the geographic and genetic relatedness of these populations. The within‐breed haplotype diversity was lower than expected based on genealogical information, which supports the existence of major male founder effects predating pedigree recording. In addition, we observed de novo mutation events among the descendants of the same ancestors, which are of interest to define paternal sub‐lineages. Our results pave the way to future studies on the estimation of the effects of Y‐chromosome haplotypes on male reproductive performances and on the conservation of Y‐chromosome diversity.

BackgroundThe Chr17q12-21.2 region is the strongest and most consistently associated region with asthma susceptibility. The functional genes or single nucleotide polymorphisms (SNPs) are not obvious due to linkage disequilibrium.ObjectivesWe sought to comprehensively investigate whole-genome sequence and RNA sequence from human bronchial epithelial cells to dissect functional genes/SNPs for asthma severity in the Severe Asthma Research Program.MethodsExpression quantitative trait loci analysis (n= 114), correlation analysis (n= 156) of gene expression and asthma phenotypes, and pathway analysis were performed in bronchialepithelial cells and replicated. Genetic association forasthma severity (426 severe vs 531 nonsevere asthma) andlongitudinal asthma exacerbations (n= 273) was performed.ResultsMultiple SNPs in gasdermin B (GSDMB) associated with asthma severity (odds ratio, >1.25) and longitudinal asthma exacerbations (P< .05). Expression quantitative trait loci analyses identified multiple SNPs associated with expression levels of post-GPI attachment to proteins 3, GSDMB, or gasdermin A(3.1×10-9-4). Higher expression levels of GSDMB correlated with asthma and greater number of exacerbations (P< .05). Expression levels of GSDMB correlated with genes involved in IFN signaling, MHC class I antigen presentation, and immune system pathways (false-discovery rate-adjusted P< .05). rs1031458 and rs3902920 in GSDMB colocalized with IFN regulatory factor binding sites and associated with GSDMB expression, asthma severity, and asthma exacerbations (P< .05).ConclusionsBy using a unique set of gene expression data from lung cells obtained using bronchoscopy from comprehensively characterized subjects with asthma, we show that SNPs in GSDMB associated with asthma severity, exacerbations, and GSDMB expression levels. Furthermore, its expression levels correlated with asthma exacerbations and antiviral pathways. Thus, GSDMB is a functional gene for both asthma susceptibility and severity.

The dark sleeper, Odontobutis potamophila, is a commercially valuable fish that widely distributed in China and Southeast Asia countries. The phenomenon of sexual dimorphism in growth is conspicuous, which the males grow substantially larger and faster than the females. However, the high-quality genome resources for gaining insight into sex-determining mechanisms to develop sex-control breeding are still lacking. Here, a chromosomal-level genome assembly of O. potamophila was generated from a combination of Illumina reads, 10× Genomics sequencing, and Hi-C chromatin interaction sequencing. The assembled genome was 1,134.62 Mb with a contig N50 of 22.25 Mb and a scaffold N50 of 24.85 Mb, representing 94.4% completeness (Benchmarking Universal Single-Copy Orthologs). Using Hi-C data, 96.49% of the total contig bases were anchored to the 22 chromosomes, with a contig N50 of 22.25 Mb and a scaffold N50 of 47.68 Mb. Approximately 54.18% of the genome were identified as repetitive elements, and 23,923 protein-coding genes were annotated in the genome. The assembled genome can be used as a valuable resource for molecular breeding and functional studies of O. potamophila in the future.

Water stress decreases the health and quality of horticulture crops by inhibiting photosynthesis, transpiration, and nutrient uptake. Application of plant growth promoting rhizobacteria (PGPR) can increase the growth, stress tolerance, and overall quality of field and greenhouse grown crops subjected to water stress. Here, we evaluated Serratia plymuthica MBSA-MJ1 for its ability to increase plant growth and quality of Petunia × hybrida (petunia), Impatiens walleriana (impatiens), and Viola × wittrockiana (pansy) plants recovering from severe water stress. Plants were treated weekly with inoculum of MBSA-MJ1, and plant growth and quality were evaluated 2 weeks after recovery from water stress. Application of S. plymuthica MBSA-MJ1 increased the visual quality and shoot biomass of petunia and impatiens and increased the flower number of petunia after recovery from water stress. In addition, in vitro characterizations showed that MBSA-MJ1 is a motile bacterium with moderate levels of antibiotic resistance that can withstand osmotic stress. Further, comprehensive genomic analyses identified genes putatively involved in bacterial osmotic and oxidative stress responses and the synthesis of osmoprotectants and vitamins that could potentially be involved in increasing plant water stress tolerance. This work provides a better understanding of potential mechanisms involved in beneficial plant-microbe interactions under abiotic stress using a novel S. plymuthica strain as a model.

The whole-genome sequence of the type strain Mycoplasma (Mesomycoplasma) hyorhinis DSM 25591 is reported and compared to the available sequences of the corresponding type strains from other strain collections to ascertain conformity. Knowledge of the identity of type strains is of importance for their application in standardized test systems.

The Indian rhinoceros (Rhinoceros unicornis) is a large herbivore found in northern India and southern Nepal. It is a critically endangered species, with an estimated population of approximately 3,600 in the wild. Genetic factors, such as the loss of genetic diversity and the accumulation of deleterious variations, are critical risk factors for the extinction of endangered species, such as the Indian rhinoceros. To support the conservation efforts of the Indian rhinoceros, we assembled its draft genome. The new genomic data will enable the study of functional genes associated with the ecological and physiological characteristics of Indian rhinoceros and help us establish more effective conservation measures. The muscles of an Indian rhinoceros that died from prostration at a zoo were collected, and the samples were stored at the National Institute for Environmental Studies (Tsukuba, Japan). Sequence data were obtained using an Illumina NovaSeq 6000 platform for short reads and an Oxford Nanopore Technologies PromethION for long reads. We generated approximately 235.2 Gbp of data. From these sequences, we assembled a 2,375,051,758 bp genome consisting of 7,615 contigs. The genome data are available from the National Center Biotechnology Information BioProject database under accession number BOSQ00000000.

Devastating fungi are one of the most important biotic factors associated with numerous infectious diseases not only in plants but in animals and humans too. Arthrinium rasikravindrae a devastating fungus is responsible for severe infections in a large number of host plants all over the world. In the present study, we analyzed the whole genome sequence of devastating fungus A. rasikravindrae strain AQZ-20, using Illumina Technology from the Novogene Bio-informatics Co., Ltd. Beijing, China. To identify associated annotation results, various corresponding functional annotations databases were utilized. The genome size was 48.24 MB with an N90 (scaffolds) length of 2,184,859 bp and encoded putative genes were 11,101, respectively. In addition, we evaluated the comparative genomic analyses with 4 fungal strains of Ascomycetes. Two related species showed a strong correlation while others exhibited a weak correlation with the A. rasikravindrae AQZ-20 fungus. This study is a discovery of the genome-scale assembly, as well as annotation for A. rasikravindrae. The results obtained from the whole genome sequencing and genomic resources developed in this study will contribute significantly to genetic improvement applications against diseases caused by A. rasikravindrae. In addition, the phylogenetic tree, followed by genomic RNA, transcriptomic, proteomic, metabolic, as well as pathogenic data reported in current research will provide deep insight for further studies in the future.

Background Methane emission by ruminants has contributed considerably to the global warming and understanding the genomic architecture of methane production may help livestock producers to reduce the methane emission from the livestock production system. The goal of our study was to identify genomic regions affecting the predicted methane emission (PME) from volatile fatty acids (VFAs) indicators and VFA traits using imputed whole-genome sequence data in Iranian Holstein cattle. Results Based on the significant-association threshold (p − 8), 33 single nucleotide polymorphisms (SNPs) were detected for PME per kg milk (n = 2), PME per kg fat (n = 14), and valeric acid (n = 17). Besides, 69 genes were identified for valeric acid (n = 18), PME per kg milk (n = 4) and PME per kg fat (n = 47) that were located within 1 Mb of significant SNPs. Based on the gene ontology (GO) term analysis, six promising candidate genes were significantly clustered in organelle organization (GO:0004984, p = 3.9 × 10− 2) for valeric acid, and 17 candidate genes significantly clustered in olfactory receptors activity (GO:0004984, p = 4 × 10− 10) for PME traits. Annotation results revealed 31 quantitative trait loci (QTLs) for milk yield and its components, body weight, and residual feed intake within 1 Mb of significant SNPs. Conclusions Our results identified 33 SNPs associated with PME and valeric acid traits, as well as 17 olfactory receptors activity genes for PME traits related to feed intake and preference. Identified SNPs were close to 31 QTLs for milk yield and its components, body weight, and residual feed intake traits. In addition, these traits had high correlations with PME trait. Overall, our findings suggest that marker-assisted and genomic selection could be used to improve the difficult and expensive-to-measure phenotypes such as PME. Moreover, prediction of methane emission by VFA indicators could be useful for increasing the size of reference population required in genome-wide association studies and genomic selection.

Streptococcus mutans is regarded as a major pathogen responsible for the onset of dental caries. S. mutans can transmit among people, especially within families. In this study, we established a new epidemiological approach to S. mutans classification. This approach can effectively differentiate among closely related isolates and offers superior reliability relative to that of the traditional MLST molecular typing method. As such, it has the potential to better support effective public health strategies centered around this bacterium that are aimed at preventing and treating dental caries., Streptococcus mutans is one of the primary pathogens responsible for the development of dental caries. Recent whole-genome sequencing (WGS)-based core genome multilocus sequence typing (cgMLST) approaches have been employed in epidemiological studies of specific human pathogens. However, this approach has not been reported in studies of S. mutans. Here, we therefore developed a cgMLST scheme for S. mutans. We surveyed 199 available S. mutans genomes as a means of identifying cgMLST targets, developing a scheme that incorporated 594 targets from the S. mutans UA159 reference genome. Sixty-eight sequence types (STs) were identified in this cgMLST scheme (cgSTs) in 80 S. mutans isolates from 40 children that were sequenced in this study, compared to 35 STs identified by multilocus sequence typing (MLST). Fifty-six cgSTs (82.35%) were associated with a single isolate based on our cgMLST scheme, which is significantly higher than in the MLST scheme (11.43%). In addition, 58.06% of all MLST profiles with ≥2 isolates were further differentiated by our cgMLST scheme. Topological analyses of the maximum likelihood phylogenetic trees revealed that our cgMLST scheme was more reliable than the MLST scheme. A minimum spanning tree of 145 S. mutans isolates from 10 countries developed based upon the cgMLST scheme highlighted the diverse population structure of S. mutans. This cgMLST scheme thus offers a new molecular typing method suitable for evaluating the epidemiological distribution of this pathogen and has the potential to serve as a benchmark for future global studies of the epidemiological nature of dental caries. IMPORTANCE Streptococcus mutans is regarded as a major pathogen responsible for the onset of dental caries. S. mutans can transmit among people, especially within families. In this study, we established a new epidemiological approach to S. mutans classification. This approach can effectively differentiate among closely related isolates and offers superior reliability relative to that of the traditional MLST molecular typing method. As such, it has the potential to better support effective public health strategies centered around this bacterium that are aimed at preventing and treating dental caries.

Staphylococcus aureus is a drug-resistant pathogen, capable of colonizing diverse ecological niches and causing a broad spectrum of infections related to a community and healthcare. In this study, we choose four methicillin-resistant S. aureus (MRSA) clinical isolates from Germany and Hungary based on our previous polyphasic characterization finding. We assumed that the selected strains have a different genetic background in terms of the presence of resistance and virulence genes, prophages, plasmids, and secondary metabolite biosynthesis genes that may play a crucial role in niche adaptation and pathogenesis. To clarify these assumptions, we performed a comparative genome analysis of these strains and observed many differences in their genomic compositions. The Hungarian isolates (SA H27 and SA H32) with ST22-SCCmec type IVa have fewer genes for multiple-drug resistance, virulence, and prophages reported in Germany isolates. Germany isolate, SA G6 acquires aminoglycoside (ant(6)-Ia and aph(3’)-III) and nucleoside (sat-4) resistance genes via phage transduction and may determine its pathogenic potential. The comparative genome study allowed the segregation of isolates of geographical origin and differentiation of the clinical isolates from the commensal isolates. This study suggested that Germany and Hungarian isolates are genetically diverse and showing variation among them due to the gain or loss of mobile genetic elements (MGEs). An interesting finding is the addition of SA G6 genome responsible for the drastic decline of the core/pan-genome ratio curve and causing the pan-genome to open wider. Functional characterizations revealed that S. aureus isolates survival are maintained by the amino acids catabolism and favor adaptation to growing in a protein-rich medium. The dispersible and singleton genes content of S. aureus genomes allows us to understand the genetic variation among the CC5 and CC22 groups. The strains with the same genetic background were clustered together, which suggests that these strains are highly alike; however, comparative genome analysis exposed that the acquisition of phage elements, and plasmids through the events of MGEs transfer contribute to differences in their phenotypic characters. This comparative genome analysis would improve the knowledge about the pathogenic S. aureus strain’s characterization, and responsible for clinically important phenotypic differences among the S. aureus strains.

Background The human plasma proteome is important for many biological processes and targets for diagnostics and therapy. It is therefore of great interest to understand the interplay of genetic and environmental factors to determine the specific protein levels in individuals and to gain a deeper insight of the importance of genetic architecture related to the individual variability of plasma levels of proteins during adult life. Methods We have combined whole-genome sequencing, multiplex plasma protein profiling, and extensive clinical phenotyping in a longitudinal 2-year wellness study of 101 healthy individuals with repeated sampling. Analyses of genetic and non-genetic associations related to the variability of blood levels of proteins in these individuals were performed. Results The analyses showed that each individual has a unique protein profile, and we report on the intra-individual as well as inter-individual variation for 794 plasma proteins. A genome-wide association study (GWAS) using 7.3 million genetic variants identified by whole-genome sequencing revealed 144 independent variants across 107 proteins that showed strong association (P −11) between genetics and the inter-individual variability on protein levels. Many proteins not reported before were identified (67 out of 107) with individual plasma level affected by genetics. Our longitudinal analysis further demonstrates that these levels are stable during the 2-year study period. The variability of protein profiles as a consequence of environmental factors was also analyzed with focus on the effects of weight loss and infections. Conclusions We show that the adult blood levels of many proteins are determined at birth by genetics, which is important for efforts aimed to understand the relationship between plasma proteome profiles and human biology and disease.

Despite the detection of Neorickettsia species DNA sequences in various trematode species and their hosts, only three Neorickettsia species have been cell culture isolated and whole-genome sequenced and are known to infect mammals and/or cause disease. The molecular mechanisms that enable the obligatory intracellular bacterium Neorickettsia to colonize trematodes and to horizontally transmit from trematodes to mammals, as well as the virulence factors associated with specific mammalian hosts, are unknown. Potomac horse fever (PHF) is a severe and acute systemic infectious disease of horses, with clinical signs that include diarrhea. Neorickettsia risticii is the only known bacterial species that causes PHF. Ingestion of insects harboring N. risticii-infected trematodes by horses leads to PHF. Our discovery of a new Neorickettsia species that causes PHF and whole-genome sequence analysis of this bacterium will improve laboratory diagnosis and vaccine development for PHF and will contribute to our understanding of Neorickettsia ecology, pathogenesis, and biology., Potomac horse fever (PHF), a severe and frequently fatal febrile diarrheal disease, has been known to be caused only by Neorickettsia risticii, an endosymbiont of digenean trematodes. Here, we report the cell culture isolation of a new Neorickettsia species found in two locations in eastern Ontario, Canada, in 2016 and 2017 (in addition to 10 variable strains of N. risticii) from N. risticii PCR-negative horses with clinical signs of PHF. Gene sequences of 16S rRNA and the major surface antigen P51 of this new Neorickettsia species were distinct from those of all previously characterized N. risticii strains and Neorickettsia species, except for those from an uncharacterized Neorickettsia species culture isolate from a horse with PHF in northern Ohio in 1991. The new Neorickettsia species nonetheless had the characteristic intramolecular repeats within strain-specific antigen 3 (Ssa3), which were found in all sequenced Ssa3s of N. risticii strains. Experimental inoculation of two naive ponies with the new Neorickettsia species produced severe and subclinical PHF, respectively, and the bacteria were reisolated from both of them, fulfilling Koch’s postulates. Serological assay titers against the new Neorickettsia species were higher than those against N. risticii. Whole-genome sequence analysis of the new Neorickettsia species revealed unique features of this bacterium compared with N. risticii. We propose to classify this new bacterium as Neorickettsia finleia sp. nov. This finding will improve the laboratory diagnosis of and vaccine for PHF, environmental risk assessment of PHF, and understanding of PHF pathogenesis and Neorickettsia biology in general.

The draft genome sequences for eight isolates of Campylobacter helveticus isolated from companion animals are described and compared with that of the type strain. On average, the genomes are 1,825,025 bp long and have a GC content of 34.4% and 1,885 coding DNA sequences (CDSs). CRISPRs were detected in only one isolate and phages in none.

Understanding how genomic variation causes differences in observable phenotypes remains a major challenge in biology. It is difficult to trace the sequence of events originating from genomic variants to changes in transcriptional responses or protein modifications. Ideally, one would conduct experiments with individuals that are at either extreme of the trait of interest, but such resources are often not available. Further, advances in genome editing will enable testing of candidate polymorphisms individually and in combination. Here we have created a resource for the study of sleep with 39 inbred lines of Drosophila—the Sleep Inbred Panel (SIP). SIP lines have stable long- and short-sleeping phenotypes developed from naturally occurring polymorphisms. These lines are fully sequenced, enabling more accurate targeting for genome editing and transgenic constructs. This panel facilitates the study of intermediate transcriptional and proteomic correlates of sleep, and supports genome editing studies to verify polymorphisms associated with sleep duration.

Modern genomic studies are rapidly growing in scale, and the analytical approaches used to analyze genomic data are increasing in complexity. Genomic data management poses logistic and computational challenges, and analyses are increasingly reliant on genomic annotation resources that create their own data management and versioning issues. As a result, genomic datasets are increasingly handled in ways that limit the rigor and reproducibility of 1many analyses. In this work, we examine the use of the Spark infrastructure for the management, access, and analysis of genomic data in comparison to traditional genomic workflows on typical cluster environments. We validate the framework by reproducing previously published results from the Alzheimer’s Disease Sequencing Project. Using the framework and analyses designed using Jupyter notebooks, Spark provides improved workflows, reduces user-driven data partitioning, and enhances the portability and reproducibility of distributed analyses required for large-scale genomic studies.

The present study reports the isolation of antibacterial exhibiting Bacillus pumilus (B. pumilus) SF-4 from soil field. The genome of this strain SF-4 was sequenced and analyzed to acquire in-depth genomic level insight related to functional diversity, evolutionary history, and biosynthetic potential. The genome of the strain SF-4 harbor 12 Biosynthetic Gene Clusters (BGCs) including four Non-ribosomal peptide synthetases (NRPSs), two terpenes, and one each of Type III polyketide synthases (PKSs), hybrid (NRPS/PKS), lipopeptide, β-lactone, and bacteriocin clusters. Plant growth-promoting genes associated with de-nitrification, iron acquisition, phosphate solubilization, and nitrogen metabolism were also observed in the genome. Furthermore, all the available complete genomes of B. pumilus strains were used to highlight species boundaries and diverse niche adaptation strategies. Phylogenetic analyses revealed local diversification and indicate that strain SF-4 is a sister group to SAFR-032 and 150a. Pan-genome analyses of 12 targeted strains showed regions of genome plasticity which regulate function of these strains and proposed direct strain adaptations to specific habitats. The unique genome pool carries genes mostly associated with “biosynthesis of secondary metabolites, transport, and catabolism” (Q), “replication, recombination and repair” (L), and “unknown function” (S) clusters of orthologous groups (COG) categories. Moreover, a total of 952 unique genes and 168 exclusively absent genes were prioritized across the 12 genomes. While newly sequenced B. pumilus SF-4 genome consists of 520 accessory, 59 unique, and seven exclusively absent genes. The current study demonstrates genomic differences among 12 B. pumilus strains and offers comprehensive knowledge of the respective genome architecture which may assist in the agronomic application of this strain in future.

Consistent with a tentative diagnosis of neuronal ceroid lipofuscinosis (NCL), autofluorescent cytoplasmic storage bodies were found in neurons from the brains of 2 related Shiba Inu dogs with a young‐adult onset, progressive neurodegenerative disease. Unexpectedly, no potentially causal NCL‐related variants were identified in a whole‐genome sequence generated with DNA from 1 of the affected dogs. Instead, the whole‐genome sequence contained a homozygous 3 base pair (bp) deletion in a coding region of HEXB. The other affected dog also was homozygous for this 3‐bp deletion. Mutations in the human HEXB ortholog cause Sandhoff disease, a type of GM2 gangliosidosis. Thin‐layer chromatography confirmed that GM2 ganglioside had accumulated in an affected Shiba Inu brain. Enzymatic analysis confirmed that the GM2 gangliosidosis resulted from a deficiency in the HEXB encoded protein and not from a deficiency in products from HEXA or GM2A, which are known alternative causes of GM2 gangliosidosis. We conclude that the homozygous 3‐bp deletion in HEXB is the likely cause of the Shiba Inu neurodegenerative disease and that whole‐genome sequencing can lead to the early identification of potentially disease‐causing DNA variants thereby refocusing subsequent diagnostic analyses toward confirming or refuting candidate variant causality.

International audience; Toxigenic Vibrio cholerae is responsible for 1.4 to 4.3 million cases with about 21,000–143,000 deaths per year. Dominance of O1 and O139 serogroups, classical and El tor biotypes, alterations in CTX phages and the pathogenicity Islands are some of the major features of V. cholerae isolates that are responsible for cholera epidemics. Whole-genome sequencing (WGS) based analyses of single-nucleotide polymorphisms (SNPs) and other infrequent genetic variants provide a robust phylogenetic framework. Recent studies on the global transmission of pandemic V. cholerae O1 strains have shown the existence of eight different phyletic lineages. In these, the classical and El Tor biotype strains were separated as two distinctly evolved lineages. The frequency of SNP accumulation and the temporal and geographical distribution supports the perception that the seventh cholera pandemic (7CP) has spread from the Bay of Bengal region in three independent but overlapping waves. The 2010 Haitian outbreak shared a common ancestor with South-Asian wave-3 strains. In West Africa and East/Southern Africa, cholera epidemics are caused by single expanded lineage, which has been introduced several times since 1970. The Latin American epidemics that occurred in 1991 and 2010 were the result of introductions of two 7CP sublineages. Sublineages representing wave-3 have caused huge outbreaks in Haiti and Yemen. The Ogawa-Inaba serotype switchover in several cholera epidemics are believed to be due to the involvement of certain selection mechanism(s) rather than due to random events. V. cholerae O139 serogroup is phylogenetically related to the 7CP El Tor, and almost all these isolates belonged to the multilocus sequence type-69. Additional phenotypic and genotypic information have been generated to understand the pathogenicity of classical and El Tor vibrios. Presence of integrative conjugative elements (ICE) with antibiotic resistance gene cassettes, clustered regularly interspaced short palindromic repeats-associated protein system and ctxAB promoter based ToxRS expression of cholera toxin (CT) separates classical and El Tor biotypes. With the availability of WGS information, several important applications including, molecular typing, antimicrobial resistance, new diagnostics, and vaccination strategies could be generated.

Pythium irregulare is an oleaginous Oomycete able to accumulate large amounts of lipids, including Eicosapentaenoic acid (EPA). EPA is an important and expensive dietary supplement with a promising and very competitive market, which is dependent on fish-oil extraction. This has prompted several research groups to study biotechnological routes to obtain specific fatty acids rather than a mixture of various lipids. Moreover, microorganisms can use low cost carbon sources for lipid production, thus reducing production costs. Previous studies have highlighted the production of EPA by P. irregulare, exploiting diverse low cost carbon sources that are produced in large amounts, such as vinasse, glycerol, and food wastewater. However, there is still a lack of knowledge about its biosynthetic pathways, because no functional annotation of any Pythium sp. exists yet. The goal of this work was to identify key genes and pathways related to EPA biosynthesis, in P. irregulare CBS 494.86, by sequencing and performing an unprecedented annotation of its genome, considering the possibility of using wastewater as a carbon source., The genome sequencing, strain acquisition and preliminary experiments and data analysis were funded by the São Paulo Research Foundation (FAPESP) and Coordination for the Improvement of Higher Education Personnel (CAPES) (Grante: 2016/10562–4). The experiment performance, data collection, analysis and interpretation of data were supported by the Portuguese Foundation for Science and Technology (FCT) under the scope of the strategic funding of [UID/BIO/04469] unit and COMPETE 2020 [POCI01-0145-FEDER-006684] and BioTecNorte operation [NORTE-01-0145-FEDER000004] funded by the European Regional Development Fund under the scope of Norte2020 - Programa Operacional Regional do Norte. The authors thank the project DD-DeCaF - Bioinformatics Services for Data-Driven Design of Cell Factories and Communities, Ref. H2020-LEIT-BIO-2015-1 686070–1, funded by the European Commission., info:eu-repo/semantics/publishedVersion

Background The complete genome sequences of 44 Bacillus cereus group isolates collected from diverse sources in Japan were analyzed to determine their genetic backgrounds and diversity levels in Japan. Multilocus sequence typing (MLST) and core-genome single-nucleotide polymorphism (SNP) typing data from whole-genome sequences were analyzed to determine genetic diversity levels. Virulence-associated gene profiles were also used to evaluate the genetic backgrounds and relationships among the isolates. Results The 44 B. cereus group isolates, including soil- and animal-derived isolates and isolates recovered from hospitalized patients and food poisoning cases, were genotyped by MLST and core-genome SNP typing. Genetic variation among the isolates was identified by the MLST and core-genome SNP phylogeny comparison against reference strains from countries outside of Japan. Exploratory principal component analysis and nonmetric multidimensional scaling (NMDS) analyses were used to assess the genetic similarities among the isolates using gene presence and absence information and isolate origins as the metadata. A significant correlation was seen between the principal components and the presence of genes encoding hemolysin BL and emetic genetic determinants in B. cereus, and the capsule proteins in B. anthracis. NMDS showed that the cluster of soil isolates overlapped with the cluster comprising animal-derived and clinical isolates. Conclusions Molecular and epidemiological analyses of B. cereus group isolates in Japan suggest that the soil- and animal-derived bacteria from our study are not a significant risk to human health. However, because several of the clinical isolates share close genetic relationships with the environmental isolates, both molecular and epidemiological surveillance studies could be used effectively to estimate virulence in these important pathogens.

Bacteria belonging to the genus Chryseobacterium are ubiquitously distributed in natural environments, plants, and animals. Except C. indologenes and C. gleum, other Chryseobacterium species rarely cause human diseases. This study reported the whole-genome features, comparative genomic analysis, and antimicrobial susceptibility patterns of C. arthrosphaerae ED882-96 isolated in Taiwan. Strain ED882-96 was collected from the blood of a patient who had alcoholic liver cirrhosis and was an intravenous drug abuser. This isolate was initially identified as C. indologenes by using matrix-assisted laser desorption ionization–time of flight mass spectrometry. The analysis of 16S ribosomal RNA gene sequence revealed that ED882-96 shared 100% sequence identity with C. arthrosphaerae type strain CC-VM-7T. The results of whole-genome sequencing of ED882-96 showed two chromosome contigs and one plasmid. The total lengths of the draft genomes of chromosome and plasmid were 4,249,864 bp and 435,667 bp, respectively. The findings of both in silico DNA–DNA hybridization and average nucleotide identity analyses clearly demonstrated that strain ED882-96 was a species of C. arthrosphaerae. A total of 83 potential virulence factor homologs were predicted in the whole-genome sequencing of strain ED882-96. This isolate was resistant to all tested antibiotics, including β-lactams, β-lactam/β-lactamase inhibitor combinations, aminoglycosides, fluoroquinolones, tetracycline, glycylcycline, and trimethoprim-sulfamethoxazole. Only one antibiotic resistance gene was recognized in the plasmid. By contrast, many antibiotic resistance genes were identified in the chromosome. The findings of this study suggest that strain ED882-96 is a highly virulent and multidrug-resistant pathogen. Knowledge regarding genomic characteristics and antimicrobial susceptibility patterns provides valuable insights into this uncommon species.

Copy number variations (CNVs) are gains and losses of large regions of genomic sequence between individuals of a species. Although CNVs have been associated with various phenotypic traits in humans and other species, the extent to which CNVs impact phenotypic variation remains unclear. In swine, as well as many other species, relatively little is understood about the frequency of CNV in the genome, sizes, locations, and other chromosomal properties. In this work, we identified and characterized CNV by utilizing whole-genome sequence from 240 members of an intensely phenotyped experimental swine herd at the U.S. Meat Animal Research Center (USMARC). These animals included all 24 of the purebred founding boars (12 Duroc and 12 Landrace), 48 of the founding Yorkshire-Landrace composite sows, 109 composite animals from generations 4 through 9, 29 composite animals from generation 15, and 30 purebred industry boars (15 Landrace and 15 Yorkshire) used as sires in generations 10 through 15. Using a combination of split reads, paired-end mapping, and read depth approaches, we identified a total of 3,538 copy number variable regions (CNVRs), including 1,820 novel CNVRs not reported in previous studies. The CNVRs covered 0.94% of the porcine genome and overlapped 1,401 genes. Gene ontology analysis identified that CNV-overlapped genes were enriched for functions related to organism development. Additionally, CNVRs overlapped with many known quantitative trait loci (QTL). In particular, analysis of QTL previously identified in the USMARC herd showed that CNVRs were most overlapped with reproductive traits, such as age of puberty and ovulation rate, and CNVRs were significantly enriched for reproductive QTL.

Background Foodborne infections caused by lung flukes of the genus Paragonimus are a significant and widespread public health problem in tropical areas. Approximately 50 Paragonimus species have been reported to infect animals and humans, but Paragonimus westermani is responsible for the bulk of human disease. Despite their medical and economic importance, no genome sequence for any Paragonimus species is available. Results We sequenced and assembled the genome of P. westermani, which is among the largest of the known pathogen genomes with an estimated size of 1.1 Gb. A 922.8 Mb genome assembly was generated from Illumina and Pacific Biosciences (PacBio) sequence data, covering 84% of the estimated genome size. The genome has a high proportion (45%) of repeat-derived DNA, particularly of the long interspersed element and long terminal repeat subtypes, and the expansion of these elements may explain some of the large size. We predicted 12,852 protein coding genes, showing a high level of conservation with related trematode species. The majority of proteins (80%) had homologs in the human liver fluke Opisthorchis viverrini, with an average sequence identity of 64.1%. Assembly of the P. westermani mitochondrial genome from long PacBio reads resulted in a single high-quality circularized 20.6 kb contig. The contig harbored a 6.9 kb region of non-coding repetitive DNA comprised of three distinct repeat units. Our results suggest that the region is highly polymorphic in P. westermani, possibly even within single worm isolates. Conclusions The generated assembly represents the first Paragonimus genome sequence and will facilitate future molecular studies of this important, but neglected, parasite group.

Coronavirus Disease 2019 (COVID-19) warrants comprehensive investigations of publicly available Severe Acute Respiratory Syndrome Coronavirus 2 (SARS-CoV-2) genomes to gain new insight about their epidemiology, mutations, and pathogenesis. Nearly 0.4 million mutations have been identified so far among the ∼60,000 SARS-CoV-2 genomic sequences. In this study, we compared a total of 371 SARS-CoV-2 published whole genomes reported from different parts of Bangladesh with 467 sequences reported globally to understand the origin of viruses, possible patterns of mutations, and availability of unique mutations. Phylogenetic analyses indicated that SARS-CoV-2 viruses might have transmitted through infected travelers from European countries, and the GR clade was found as predominant in Bangladesh. Our analyses revealed 4604 mutations at the RNA level including 2862 missense mutations, 1192 synonymous mutations, 25 insertions and deletions and 525 other types of mutation. In line with the global trend, D614G mutation in spike glycoprotein was predominantly high (98 %) in Bangladeshi isolates. Interestingly, we found the average number of mutations in ORF1ab, S, ORF3a, M, and N were significantly higher (p < 0.001) for sequences containing the G614 variant compared to those having D614. Previously reported frequent mutations, such as R203K, D614G, G204R, P4715L and I300F at protein levels were also prevalent in Bangladeshi isolates. Additionally, 34 unique amino acid changes were revealed and categorized as originating from different cities. These analyses may increase our understanding of variations in SARS-CoV-2 virus genomes, circulating in Bangladesh and elsewhere.

The island inhabitants of Sardinia have long been a focus for studies of complex human traits due to their unique ancestral background and population isolation reflecting geographic and cultural restriction. Population isolates share decreased genomic diversity, increased linkage disequilibrium, and increased inbreeding coefficients. In many regions, dogs and humans have been exposed to the same natural and artificial forces of environment, growth, and migration. Distinct dog breeds have arisen through human-driven selection of characteristics to meet an ideal standard of appearance and function. The Fonni’s Dog, an endemic dog population on Sardinia, has not been subjected to an intensive system of artificial selection, but rather has developed alongside the human population of Sardinia, influenced by geographic isolation and unregulated selection based on its environmental adaptation and aptitude for owner-desired behaviors. Through analysis of 28 dog breeds, represented with whole-genome sequences from 13 dogs and ∼170,000 genome-wide single nucleotide variants from 155 dogs, we have produced a genomic illustration of the Fonni’s Dog. Genomic patterns confirm within-breed similarity, while population and demographic analyses provide spatial identity of Fonni’s Dog to other Mediterranean breeds. Investigation of admixture and fixation indices reveals insights into the involvement of Fonni’s Dogs in breed development throughout the Mediterranean. We describe how characteristics of population isolates are reflected in dog breeds that have undergone artificial selection, and are mirrored in the Fonni’s Dog through traditional isolating factors that affect human populations. Lastly, we show that the genetic history of Fonni’s Dog parallels demographic events in local human populations.

Acetobacter ghanensis LMG 23848(T) and Acetobacter senegalensis 108B are acetic acid bacteria that originate from a spontaneous cocoa bean heap fermentation process and that have been characterised as strains with interesting functionalities through metabolic and kinetic studies. As there is currently little genetic information available for these species, whole-genome sequencing of A. ghanensis LMG 23848(T) and A. senegalensis 108B and subsequent data analysis was performed. This approach not only revealed characteristics such as the metabolic potential and genomic architecture, but also allowed to indicate the genetic adaptations related to the cocoa bean fermentation process. Indeed, evidence was found that both species possessed the genetic ability to be involved in citrate assimilation and displayed adaptations in their respiratory chain that might improve their competitiveness during the cocoa bean fermentation process. In contrast, other properties such as the dependence on glycerol or mannitol and lactate as energy sources or a less efficient acid stress response may explain their low competitiveness. The presence of a gene coding for a proton-translocating transhydrogenase in A. ghanensis LMG 23848(T) and the genes involved in two aromatic compound degradation pathways in A. senegalensis 108B indicate that these strains have an extended functionality compared to Acetobacter species isolated from other ecosystems.

There is an increasing interest in using whole-genome sequence data in genomic selection breeding programmes. Prediction of breeding values is expected to be more accurate when whole-genome sequence is used, because the causal mutations are assumed to be in the data. We performed genomic prediction for the number of eggs in white layers using imputed whole-genome resequence data including ~4.6 million SNPs. The prediction accuracies based on sequence data were compared with the accuracies from the 60 K SNP panel. Predictions were based on genomic best linear unbiased prediction (GBLUP) as well as a Bayesian variable selection model (BayesC). Moreover, the prediction accuracy from using different types of variants (synonymous, non-synonymous and non-coding SNPs) was evaluated. Genomic prediction using the 60 K SNP panel resulted in a prediction accuracy of 0.74 when GBLUP was applied. With sequence data, there was a small increase (~1%) in prediction accuracy over the 60 K genotypes. With both 60 K SNP panel and sequence data, GBLUP slightly outperformed BayesC in predicting the breeding values. Selection of SNPs more likely to affect the phenotype (i.e. non-synonymous SNPs) did not improve the accuracy of genomic prediction. The fact that sequence data were based on imputation from a small number of sequenced animals may have limited the potential to improve the prediction accuracy. A small reference population (n = 1004) and possible exclusion of many causal SNPs during quality control can be other possible reasons for limited benefit of sequence data. We expect, however, that the limited improvement is because the 60 K SNP panel was already sufficiently dense to accurately determine the relationships between animals in our data.

Understanding how genomic variation causes differences in observable phenotypes remains a major challenge in biology. It is difficult to trace the sequence of events originating from genomic variants to changes in transcriptional responses or protein modifications. Ideally, one would conduct experiments with individuals that are at either extreme of the trait of interest, but such resources are often not available. Further, advances in genome editing will enable testing of candidate polymorphisms individually and in combination. Here we have created a resource for the study of sleep with 39 inbred lines of Drosophila—the Sleep Inbred Panel (SIP). SIP lines have stable long- and short-sleeping phenotypes developed from naturally occurring polymorphisms. These lines are fully sequenced, enabling more accurate targeting for genome editing and transgenic constructs. This panel facilitates the study of intermediate transcriptional and proteomic correlates of sleep, and supports genome editing studies to verify polymorphisms associated with sleep duration.

Abyssivirga alkaniphila strain L81T, recently isolated from a black smoker biofilm at the Loki’s Castle hydrothermal vent field, was previously described as a mesophilic, obligately anaerobic heterotroph able to ferment carbohydrates, peptides, and aliphatic hydrocarbons. The strain was classified as a new genus within the family Lachnospiraceae. Herein, its genome is analyzed and A. alkaniphila is reassigned to the genus Vallitalea as a new strain of V. guaymasensis, designated V. guaymasensis strain L81. The 6.4 Mbp genome contained 5651 protein encoding genes, whereof 4043 were given a functional prediction. Pathways for fermentation of mono-saccharides, di-saccharides, peptides, and amino acids were identified whereas a complete pathway for the fermentation of n-alkanes was not found. Growth on carbohydrates and proteinous compounds supported methane production in co-cultures with Methanoplanus limicola. Multiple confurcating hydrogen-producing hydrogenases, a putative bifurcating electron-transferring flavoprotein—butyryl-CoA dehydrogenase complex, and a Rnf-complex form a basis for the observed hydrogen-production and a putative reverse electron-transport in V. guaymasensis strain L81. Combined with the observation that n-alkanes did not support growth in co-cultures with M. limicola, it seemed more plausible that the previously observed degradation patterns of crude-oil in strain L81 are explained by unspecific activation and may represent a detoxification mechanism, representing an interesting ecological function. Genes encoding a capacity for polyketide synthesis, prophages, and resistance to antibiotics shows interactions with the co-occurring microorganisms. This study enlightens the function of the fermentative microorganisms from hydrothermal vents systems and adds valuable information on the bioprospecting potential emerging in deep-sea hydrothermal systems.

The whole-genome sequences of progenies with low-density single-nucleotide polymorphism (SNP) genotypes can be imputed with high accuracy based on the deep-coverage sequences of key ancestors. With this imputation technology, a more powerful genome-wide association study (GWAS) can be carried out using imputed whole-genome variants and the phenotypes of interest to overcome the shortcomings of low-power detection and the large confidence interval derived from low-density SNP markers in classic association studies. In this study, 19 ancestors of a large-scale swine F2 White Duroc × Erhualian population were deeply sequenced for their genome with an average coverage of 25×. Considering 98 pigs from 10 different breeds with high-quality deep sequenced genomes, we imputed the whole genomic variants of 1020 F2 pigs genotyped by the PorcineSNP60 BeadChip with high accuracy and obtained 14,851,440 sequence variants after quality control. Based on this, 87 novel quantitative traits loci (QTLs) for 18 hematological traits at three different physiological stages of the F2 pigs were identified, among which most of the novel QTLs have been repeated in two of the three stages. Literature mining pinpointed that the FGF14 and LCLAT1 genes at SSC11 and SSC3 may affect the MCH at day 240 and MCV at day 18, respectively. The present study shows that combining high-quality imputed genomic variants and correlated phenomic traits into GWAS can improve the capability to detect QTL considerably. The large number of different QTLs for hematological traits identified at multiple growth stages implies the complexity and time specificity of these traits.

Whole-genome sequence (WGS) analysis has revolutionized the food safety industry by enabling high-resolution typing of foodborne bacteria. Higher resolving power allows investigators to identify origins of contamination during illness outbreaks and regulatory activities quickly and accurately. Government agencies and industry stakeholders worldwide are now analyzing WGS data routinely. Although researchers have published many studies that assess the efficacy of WGS data analysis for source attribution, guidance for interpreting WGS analyses is lacking. Here, we provide the framework for interpreting WGS analyses used by the Food and Drug Administration's Center for Food Safety and Applied Nutrition (CFSAN). We based this framework on the experiences of CFSAN investigators, collaborations and interactions with government and industry partners, and evaluation of the published literature. A fundamental question for investigators is whether two or more bacteria arose from the same source of contamination. Analysts often count the numbers of nucleotide differences (single-nucleotide polymorphisms [SNPs]) between two or more genome sequences to measure genetic distances. However, using SNP thresholds in isolation to assess whether bacteria originated from the same source can be misleading. Bacteria that are isolated from food, environmental, or clinical samples are representatives of bacterial populations. These populations are subject to evolutionary forces that can change genome sequences. Therefore, interpreting WGS analyses of foodborne bacteria requires a more sophisticated approach. Here, we present a framework for interpreting WGS analyses that combines SNP counts with phylogenetic tree topologies and bootstrap support. We also clarify the roles of WGS, epidemiological, traceback, and other evidence in forming the conclusions of investigations, making clear that WGS data alone is insufficient for presumptions of relatedness between bacterial isolates. Finally, we present examples that illustrate the application of this framework to real-world situations.

Brucellosis is a worldwide zoonotic disease caused by Brucella spp. In China, brucellosis is recognized as a reemerging disease mainly caused by Brucella melitensis specie. To better understand the currently endemic B. melitensis strains in China, three Brucella genotyping methods were applied to 110 B. melitensis strains obtained in past several years. By MLVA genotyping, five MLVA-8 genotypes were identified, among which genotypes 42 (1-5-3-13-2-2-3-2) was recognized as the predominant genotype, while genotype 63 (1-5-3-13-2-3-3-2) and a novel genotype of 1-5-3-13-2-4-3-2 were second frequently observed. MLVA-16 discerned a total of 57 MLVA-16 genotypes among these Brucella strains, with 41 genotypes being firstly detected and the other 16 genotypes being previously reported. By BruMLSA21 typing, six sequence types (STs) were identified, among them ST8 is the most frequently seen in China while the other five STs were firstly detected and designated as ST137, ST138, ST139, ST140, and ST141 by international multilocus sequence typing database. Whole-genome sequence (WGS)-single-nucleotide polymorphism (SNP)-based typing and phylogenetic analysis resolved Chinese B. melitensis strains into five clusters, reflecting the existence of multiple lineages among these Chinese B. melitensis strains. In phylogeny, Chinese lineages are more closely related to strains collected from East Mediterranean and Middle East countries, such as Turkey, Kuwait, and Iraq. In the next few years, MLVA typing will certainly remain an important epidemiological tool for Brucella infection analysis, as it displays a high discriminatory ability and achieves result largely in agreement with WGS-SNP-based typing. However, WGS-SNP-based typing is found to be the most powerful and reliable method in discerning Brucella strains and will be popular used in the future.

Brucellosis is a worldwide zoonotic disease caused by Brucella spp. In China, brucellosis is recognized as a reemerging disease mainly caused by Brucella melitensis specie. To better understand the currently endemic B. melitensis strains in China, three Brucella genotyping methods were applied to 110 B. melitensis strains obtained in past several years. By MLVA genotyping, five MLVA-8 genotypes were identified, among which genotypes 42 (1-5-3-13-2-2-3-2) was recognized as the predominant genotype, while genotype 63 (1-5-3-13-2-3-3-2) and a novel genotype of 1-5-3-13-2-4-3-2 were second frequently observed. MLVA-16 discerned a total of 57 MLVA-16 genotypes among these Brucella strains, with 41 genotypes being firstly detected and the other 16 genotypes being previously reported. By BruMLSA21 typing, six sequence types (STs) were identified, among them ST8 is the most frequently seen in China while the other five STs were firstly detected and designated as ST137, ST138, ST139, ST140, and ST141 by international multilocus sequence typing database. Whole-genome sequence (WGS)-single-nucleotide polymorphism (SNP)-based typing and phylogenetic analysis resolved Chinese B. melitensis strains into five clusters, reflecting the existence of multiple lineages among these Chinese B. melitensis strains. In phylogeny, Chinese lineages are more closely related to strains collected from East Mediterranean and Middle East countries, such as Turkey, Kuwait, and Iraq. In the next few years, MLVA typing will certainly remain an important epidemiological tool for Brucella infection analysis, as it displays a high discriminatory ability and achieves result largely in agreement with WGS-SNP-based typing. However, WGS-SNP-based typing is found to be the most powerful and reliable method in discerning Brucella strains and will be popular used in the future.

Highly pathogenic avian influenza virus (HPAIV) is a multihost pathogen with lineages that pose health risks for domestic birds, wild birds, and humans. One mechanism of intercontinental HPAIV spread is through wild bird reservoirs, and wild birds were the likely sources of a Eurasian (EA) lineage HPAIV into North America in 2014. The introduction resulted in several reassortment events with North American (NA) lineage low‐pathogenic avian influenza viruses and the reassortant EA/NA H5N2 went on to cause one of the largest HPAIV poultry outbreaks in North America. We evaluated three hypotheses about novel HPAIV introduced into wild and domestic bird hosts: (i) transmission of novel HPAIVs in wild birds was restricted by mechanisms associated with highly pathogenic phenotypes; (ii) the HPAIV poultry outbreak was not self‐sustaining and required viral input from wild birds; and (iii) reassortment of the EA H5N8 generated reassortant EA/NA AIVs with a fitness advantage over fully Eurasian lineages in North American wild birds. We used a time‐rooted phylodynamic model that explicitly incorporated viral population dynamics with evolutionary dynamics to estimate the basic reproductive number (R 0) and viral migration among host types in domestic and wild birds, as well as between the EA H5N8 and EA/NA H5N2 in wild birds. We did not find evidence to support hypothesis (i) or (ii) as our estimates of the transmission parameters suggested that the HPAIV outbreak met or exceeded the threshold for persistence in wild birds (R 0 > 1) and poultry (R 0 ≈ 1) with minimal estimated transmission among host types. There was also no evidence to support hypothesis (iii) because R 0 values were similar among EA H5N8 and EA/NA H5N2 in wild birds. Our results suggest that this novel HPAIV and reassortments did not encounter any transmission barriers sufficient to prevent persistence when introduced to wild or domestic birds.

The objectives of this study are to determine the occurrence of antimicrobial resistance (AMR) genes using whole-genome sequence (WGS) of Streptococcus uberis (S. uberis) and Streptococcus dysgalactiae (S. dysgalactiae) isolates, recovered from dairy cows in the Canadian Maritime Provinces. A secondary objective included the exploration of the association between phenotypic AMR and the genomic characteristics (genome size, guanine–cytosine content, and occurrence of unique gene sequences). Initially, 91 isolates were sequenced, and of these isolates, 89 were assembled. Furthermore, 16 isolates were excluded due to larger than expected genomic sizes (>2.3 bp × 1,000 bp). In the final analysis, 73 were used with complete WGS and minimum inhibitory concentration records, which were part of the previous phenotypic AMR study, representing 18 dairy herds from the Maritime region of Canada (1). A total of 23 unique AMR gene sequences were found in the bacterial genomes, with a mean number of 8.1 (minimum: 5; maximum: 13) per genome. Overall, there were 10 AMR genes [ANT(6), TEM-127, TEM-163, TEM-89, TEM-95, Linb, Lnub, Ermb, Ermc, and TetS] present only in S. uberis genomes and 2 genes unique (EF-TU and TEM-71) to the S. dysgalactiae genomes; 11 AMR genes [APH(3′), TEM-1, TEM-136, TEM-157, TEM-47, TetM, bl2b, gyrA, parE, phoP, and rpoB] were found in both bacterial species. Two-way tabulations showed association between the phenotypic susceptibility to lincosamides and the presence of linB (P = 0.002) and lnuB (P 11 AMR genes present in the genome, compared with 250,000 cells/mL, a trend toward higher odds of resistance compared with the baseline category of

Zika virus (ZIKV) (genus Flavivirus, family Flaviviridae) is an emerging pathogen associated with microcephaly and Guillain-Barre syndrome. The rapid spread of ZIKV disease in over 60 countries and the large numbers of travel-associated cases have caused worldwide concern. Thus, intensified surveillance of cases among immigrants and tourists from ZIKV-endemic areas is important for disease control and prevention. In this study, using Next Generation Sequencing, we reported the first whole-genome sequence of ZIKV strain AFMC-U, amplified from the urine of a traveler returning to Korea from the Philippines. Phylogenetic analysis showed geographic-specific clustering. Our results underscore the importance of examining urine in the diagnosis of ZIKV infection.

Shrimp is extensively cultured worldwide. Shrimp farming is suffering from a variety of diseases. Probiotics are considered to be one of the effective methods to prevent and cure shrimp diseases. Exiguobacterium arabatum W‐01, a gram‐positive and orange‐pigmented bacterium, was isolated from the intestine of a healthy Penaeus vannamei specimen. Whole‐genome sequencing revealed a genome of 2,914,854 bp, with 48.02% GC content. In total, 3,083 open reading frames (ORFs) were identified, with an average length of 843.98 bp and a mean GC content of 48.11%, accounting for 89.27% of the genome. Among these ORFs, 2,884 (93.5%) genes were classified into Clusters of Orthologous Groups (COG) families comprising 21 functional categories, and 1,650 ORFs were classified into 83 functional Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways. A total of 27 rRNA operons and 68 tRNAs were identified, with all 20 amino acids represented. In addition, 91 genomic islands, 68 potential prophages, and 33 tandem repeats, but no clustered regularly interspaced short palindromic repeats (CRISPRs), were found. No resistance genes and only one virulence gene were identified. Among the 150 secreted proteins of E. arabatum W‐01, a variety of transport system substrate‐binding proteins, enzymes, and biosynthetic proteins, which play important roles in the uptake and metabolism of nutrients, were found. Two adherence‐related protein genes and 31 flagellum‐related protein genes were also identified. Taken together, these results indicate potential probiotic functions for E. arabatum W‐01.