Unique genetic architecture of prolificacy in 'Sikkim Primitive' maize unraveled through whole-genome resequencing-based DNA polymorphism

KEY MESSAGE

'Sikkim Primitive' maize landrace, unique for prolificacy (7-9 ears per plant) possesses unique genomic architecture in branching and inflorescence-related gene(s), and locus Zm00001eb365210 encoding glycosyltransferases was identified as the putative candidate gene underlying QTL (qProl-SP-8.05) for prolificacy. The genotype possesses immense usage in breeding high-yielding baby-corn genotypes. 'Sikkim Primitive' is a native landrace of North Eastern Himalayas, and is characterized by having 7-9 ears per plant compared to 1-2 ears in normal maize. Though 'Sikkim Primitive' was identified in the 1960s, it has not been characterized at a whole-genome scale. Here, we sequenced the entire genome of an inbred (MGUSP101) derived from 'Sikkim Primitive' along with three non-prolific (HKI1128, UMI1200, and HKI1105) and three prolific (CM150Q, CM151Q and HKI323) inbreds. A total of 942,417 SNPs, 24,160 insertions, and 27,600 deletions were identified in 'Sikkim Primitive'. The gene-specific functional mutations in 'Sikkim Primitive' were classified as 10,847 missense (54.36%), 402 non-sense (2.015%), and 8,705 silent (43.625%) mutations. The number of transitions and transversions specific to 'Sikkim Primitive' were 666,021 and 279,950, respectively. Among all base changes, (G to A) was the most frequent (215,772), while (C to G) was the rarest (22,520). Polygalacturonate 4-α-galacturonosyltransferase enzyme involved in pectin biosynthesis, cell-wall organization, nucleotide sugar, and amino-sugar metabolism was found to have unique alleles in 'Sikkim Primitive'. The analysis further revealed the Zm00001eb365210 gene encoding glycosyltransferases as the putative candidate underlying QTL (qProl-SP-8.05) for prolificacy in 'Sikkim Primitive'. High-impact nucleotide variations were found in ramosa3 (Zm00001eb327910) and zeaxanthin epoxidase1 (Zm00001eb081460) genes having a role in branching and inflorescence development in 'Sikkim Primitive'. The information generated unraveled the genetic architecture and identified key genes/alleles unique to the 'Sikkim Primitive' genome. This is the first report of whole-genome characterization of the 'Sikkim Primitive' landrace unique for its high prolificacy.

Genomic characterization of Rocahepevirus ratti hepatitis E virus genotype C1 in Yunnan province of China

The Rocahepevirus ratti hepatitis E virus genotype C1 (HEV-C1) has been documented to infect humans. However, the understanding of HEV-C1 remains constrained. This study aims to determine the prevalence and genomic characteristics of HEV-C1 in small animals in Yunnan province of southwestern China. A total of 444 liver tissues were collected from animals covering the orders Rodentia, Soricomorpha, Scandentia and Erinaceomorpha in three regions in Yunnan. Then Paslahepevirus balayani and Rocahepevirus were examined using RT-qPCR. The detection rate of Rocahepevirus was 12.95 % (36/278) in animals of order Rodentia, with 14.77 % (35/237) in Rattus tanezumi and 33.33 % (1/3) in Niviventer fulvescens. No Paslahepevirus balayani was detected. Additionally, two full-length Rocahepevirus sequences (MSE-17 and LHK-54) and thirty-three partial ORF1 sequences were amplified and determined to be HEV-C1. MSE-17 and LHK-54 shared moderate nucleotide identity (78.9 %-80.3 %) with HEV-C1 isolated in rats and humans. The HEV-C1 isolated from Niviventer fulvescens demonstrated a 100 % nucleotide identity with that from Rattus tanezumi. The rat HEV-C1 sequences isolated in our study and other Asian HEV-C1 sequences were phylogenetically distant from those isolated in North America and Europe. Furthermore, the two full-length sequences isolated in our study had less amino acid substitutions in the motifs of RNA-dependent RNA polymerase domain (F204L and L238F), compared with other Asian sequences. In summary, HEV-C1 commonly spreads in rats in Yunnan province of China. Our findings suggest a spatially associated phylogeny, and potential cross-species transmission of HEV-C1.

A person-to-person transmission cluster of severe fever with thrombocytopenia syndrome characterized by mixed viral infections with familial and nosocomial clustering

Severe fever with thrombocytopenia syndrome (SFTS) is an emerging tick-borne infectious disease with sporadic occurrence and high mortality. Herein, we report an example of the in-hospital transmission of SFTS virus (SFTSV) infections with familial and nosocomial clustering in Zhejiang Province, eastern China, from March to April 2023. The epidemiological investigation and genomic analysis revealed that at least eight suspected cases of SFTS occurred in this cluster, including one death and one asymptomatic case. Our report reemphasizes the risk of familial and nosocomial SFTSV infections in healthcare settings and the urgent need for the long-term systematic surveillance of SFTSV evolution in humans and animals in the eastern coastal regions of China.

Isolation and genomic characterization of a cypovirus from Clanis bilineata

Clanis bilineata Walker, soybean hawkmoth, belongs to the subfamily Ambulicinae (Sphingidae, Lepidoptera) and is an edible insect that usually grows on soybean leaves. In this study, we isolated a new cypovirus from naturally diseased Clanis bilineata larvae (named CbCPV), scanned its structure, sequenced its genome, and studied its phylogenetic relationship to other cypoviruses. Microscopy showed that CbCPV polyhedral occlusion bodies were about 1.878 μm on average and contained many virions in the ultrathin sections. The complete genome sequence of CbCPV is 22,812 bp comprising 10 segmented double-stranded RNAs. Apart from segment 1 containing one open reading frame (ORF) and one sub-ORF, the other nine segments all contain one open reading frame and encoded one putative protein. The non-coding regions contained conserved sequences at 5' termini (AGUCAAA) and 3' termini (AGC), except segment 4 containing a different 5' termini (AUGUUUA). The whole sequence of the polyhedrin gene in CbCPV contained 892 nucleotides, encoding a protein of 246 amino acids. Based on amino acid sequences of polyhedrin or RNA dependent RNA polymerase (RdRp), the phylogenetic analysis indicated that CbCPV was closely related to DnCPV-23. The putative function of all segments differed from each other, but the most closely related species of segments were DnCPV-23 with 98.2-99.8% nucleotide identity. Overall, the evidence of morphology, protein analysis and nucleic acids (genomic pattern) showed that CbCPV is a new isolate in the cypovirus-23 type and can be termed Clanis bilineata cypovirus type 23 (CbCPV-23).

Characterization and genomic analysis of a novel bacteriophage BUCT_49532 lysing Klebsiella pneumoniae

Bacteriophages are a type of virus widely distributed in nature that demonstrates a remarkable aptitude for selectively recognizing and infecting bacteria. In particular, Klebsiella pneumoniae is acknowledged as a clinical pathogen responsible for nosocomial infections and frequently develops multidrug resistance. Considering the increasing prevalence of antibiotic-resistant bacteria, bacteriophages have emerged as a compelling alternative therapeutic approach. In this study, a novel phage named BUCT_49532 was isolated from sewage using K. pneumoniae K1119 as the host. Electron microscopy revealed that BUCT_49532 belongs to the Caudoviricetes class. Further analysis through whole genome sequencing demonstrated that BUCT_49532 is a Jedunavirus comprised of linear double-stranded DNA with a length of 49,532 bp. Comparative genomics analysis based on average nucleotide identity (ANI) values revealed that BUCT_49532 should be identified as a novel species. Characterized by a good safety profile, high environmental stability, and strong lytic performance, phage BUCT_49532 presents an interesting case for consideration. Although its host range is relatively narrow, its application potential can be expanded by utilizing phage cocktails, making it a promising candidate for biocontrol approaches.

Epidemiology and Genomic characteristics of arenavirus in rodents from the southeast coast of P.R. China

BACKGROUND

Wenzhou virus (WENV), a member of the Mammarenavirus genus in the Arenaviridae family, has been detected in wild rodents from eight provinces in China, including Zhejiang, Shandong, Hainan, Xinjiang, Hunan, Guangdong, Yunnan, and Jiangxi provinces, and some countries from Southeast Asia. The IgG-antibodies of WENV have been detected in both healthy populations and patients with unknown fever and respiratory symptoms. However, the potential harmfulness of WENV to humans has been underestimated due to mild symptoms after infection, similar to respiratory diseases. Thus, it is imperative to enhance the surveillance of WENV in wild rodents, particularly Rattus norvegicus, and continuously monitor its prevalence.

RESULTS

From 2017 to 2021, a total of 390 wild rodents were collected from six provinces in the eastern and southern coastal areas, containing nine species of rats. Samples of each tissue were collected, and PCR amplified for identification. Four R. norvegicus samples were detected to be WENV-positive. No genomic sequence of WENV was detected in Rattus flavipectus, Rattus losea, Suncus murinus, Apodemus agrarius, Mus musculus, Microtus fortis, Micromys minutus, and Niviventer niviventer from Jiangsu, Zhejiang, Fujian, Hainan, Guangdong and Guangxi provinces. Three genomic sequences were identified to be WENV by phylogenetic analysis. The full-length sequences of HAIKOU-40 were amplified in R. norvegicus from Hainan, which showed a close relationship to Wufeng/ WFS, sharing 84.5-89.4% homology at the nucleotide level and 91.6-98.9% homology at the amino acid level. Phylogenetic analysis revealed that HAIKOU-40 formed an Asia-specific cluster with all WENVs and Loie River mammarenavirus (LORV), provisionally named Asian ancestry. This cluster has diverged earlier from the remaining mammarenavirus. The sequences obtained in Xiamen, Fujian province showed more than 90% nucleotide identities with WENV, which may be a strain of WENV. Additionally, the sequence of Wuxi-87 which was a positive sequence detected in Wuxi, Jiangsu province exhibited 83% nucleotide identity with Lassa virus (LASV). Further efforts will be made to isolate and identify this virus strain, verify the relationship between Wuxi-87 and LASV, and confirm whether R. norvegicus is a new host of LASV.

CONCLUSIONS

In this study, we conducted a systematic examination of the prevalence of WENV among rodents on the southeast coast of China. Additionally, we characterized the genome of a newly discovered WENV strain, that confirmed the role of R. norvegicus in the transmission of WENV. This highlights the importance of investigating the prevalence of WENV in both wild rodents and humans.

Genomic characterization and phylogenetic analysis of Aleutian mink disease virus identified in a sudden death mink case

Aleutian mink disease (AMD) is one of the most serious diseases in minks worldwide, it brings tremendous financial losses in mink farming. AMD virus (AMDV) has unusually high genetic diversity, its genomic structure remains unclear. In 2014, sudden death of breeding minks was occurred in northeast China. After clinical signs evaluation and virus isolation, AMDV was identified in all sudden death minks, we investigated the complete genomic sequence of AMDV-LM isolated from the sudden death case. The full-genome sequence of AMDV-LM was 7 nucleotides (nts) or 8 nts longer than isolates AMDV-BJ and AMDV-G. AMDV-LM contained two unique nucleotide changes in VP2 (G79T, T710C), which led to two amino acid changes G27W and L237S. For NS1, some unique point mutations, such as A374C, A428C, A463C, and T476A were found and resulted in four unique amino acid mutations at N24V, H125P, V143P, K155Q, and V159N, respectively. The predicted secondary structure of the 5' terminal of AMDV-LM formed a large bubble formation near the 5' end, which affected the stability of the U-shaped hairpin. Phylogenetic analysis demonstrated that AMDV-LM was closely related to Chinese isolates and confirmed that AMDV strains circulating in China had different origins of ancestors. This study was first to investigate the association of sudden death of adult breeding minks with AMDV infection. Our findings provide useful suggestions for evaluation of the pathogenic potential of AMDV, additional details on AMDV genome characterization were also presented. Future work should focus on the importance of AMDV-LM strain in mink infection.

Genetic characterization of Sus scrofa papillomavirus type 1 from domestic pigs in Guangxi Province, China

Sus scrofa papillomatosis (SsP) is a tumour caused by Sus scrofa papillomaviruses (SsPVs). To investigate the presence of SsPVs in China, 354 domestic pig skin samples collected from Guangxi Province were examined for SsPV DNA by PCR. Three SsPV1s (GX12, GX14, and GX18) were identified with a prevalence of 0.847% (3/354). Sequence analysis showed that L1 of SsPV1/GX12 and SsPV1/GX14 had 99.7% and 99.6% nucleotide identify with the reference SsPV1a, respectively. Phylogenetic and evolutionary analyses showed that SsPV1/GX12 and SsPV1/14 clustered into SsPV1a and that SsPV1/GX18 clustered into SsPV1b. Compared with other SsPV L1 and L2 proteins, we found that the SsPV1/GX18 and SsPV1b strains shared the same unique substitutions, and SsPV1/GX12, SsPV1/GX14, and SsPV1a shared almost identical amino acid sequences. This study reports the first detection of SsPV DNA in China based on whole genome information and provides a scientific basis for the development of SsPV pathogenic biology, epidemiology, and prevention, as well as control technology research.

Diverse single-stranded DNA viruses from viral metagenomics on a cynopterus bat in China

Bats serve as reservoirs for many emerging viruses. Cressdnaviruses can infect a wide range of animals, including agricultural species, such as pigs, in which porcine circoviruses cause severe gastroenteritis. New cressdnaviruses have also attracted considerable attention recently, due to their involvement with infectious diseases. However, little is known about their host range and many cressdnaviruses remain poorly characterized. We identified and characterized 11 contigs consisting of previously unknown cressdnaviruses from a rectal swab sample of a Cynopterus bat collected in Yunnan Province, China, in 2011. Full genomes of two cressdnaviruses (OQ267680, 2069 nt; OQ351951, 2382 nt), and a nearly complete genome for a third (OQ267683, 2361 nt) were obtained. Phylogenetic analyses and the characteristics of these viral genomes suggest a high degree of ssDNA virus diversity. These results shed light on cressdnavirus diversity and the probable role of Cynopterus bats as their hosts.

SARS-CoV-2 genomic characterization and evolution in China

The pandemic caused by severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) affected global health worldwide due to its high contagiousness. During the viral spread, many mutations occurred within the virus genome. China has adopted nonpharmaceutical intervention (NPI) to contain COVID-19 outbreaks. In order to understand the evolution and genomic variation of SARS-CoV-2 in China under this policy, a total of 524 sequences downloaded from Global Initiative on Sharing All Influenza Data (GISAID) between 2019 and 2022 were included in this study. The time-scaled evolutionary analysis showed that these sequences clustered in three groups (Group A-C). Group B and C accounted for the majority of the sequences whose divergence times were around 2020 and distributed in multiple regions. Group A was mainly composed of G variants, which were mainly isolated from several regions. Moreover, we found that 191 sites had mutations with no less than 3 times, including 30 amino acids that were deleted. Finally, we found that spike and nucleocapsid genes underwent positive selection evolution, indicating that the mutations within spike and nucleocapsid genes increased the SARS-CoV-2 contagiousness. Hence, this study preliminarily elucidates the evolutionary characteristics and genomic mutations of SARS-CoV-2 under the implementation of the NPI policy in China, providing scientific basis for further understanding the control effect of the NPI policy on the epidemic.

Evaluating the mouse neural precursor line, SN4741, as a suitable proxy for midbrain dopaminergic neurons

To overcome the ethical and technical limitations of in vivo human disease models, the broader scientific community frequently employs model organism-derived cell lines to investigate disease mechanisms, pathways, and therapeutic strategies. Despite the widespread use of certain in vitro models, many still lack contemporary genomic analysis supporting their use as a proxy for the affected human cells and tissues. Consequently, it is imperative to determine how accurately and effectively any proposed biological surrogate may reflect the biological processes it is assumed to model. One such cellular surrogate of human disease is the established mouse neural precursor cell line, SN4741, which has been used to elucidate mechanisms of neurotoxicity in Parkinson disease for over 25 years. Here, we are using a combination of classic and contemporary genomic techniques - karyotyping, RT-qPCR, single cell RNA-seq, bulk RNA-seq, and ATAC-seq - to characterize the transcriptional landscape, chromatin landscape, and genomic architecture of this cell line, and evaluate its suitability as a proxy for midbrain dopaminergic neurons in the study of Parkinson disease. We find that SN4741 cells possess an unstable triploidy and consistently exhibits low expression of dopaminergic neuron markers across assays, even when the cell line is shifted to the non-permissive temperature that drives differentiation. The transcriptional signatures of SN4741 cells suggest that they are maintained in an undifferentiated state at the permissive temperature and differentiate into immature neurons at the non-permissive temperature; however, they may not be dopaminergic neuron precursors, as previously suggested. Additionally, the chromatin landscapes of SN4741 cells, in both the differentiated and undifferentiated states, are not concordant with the open chromatin profiles of ex vivo, mouse E15.5 forebrain- or midbrain-derived dopaminergic neurons. Overall, our data suggest that SN4741 cells may reflect early aspects of neuronal differentiation but are likely not a suitable proxy for dopaminergic neurons as previously thought. The implications of this study extend broadly, illuminating the need for robust biological and genomic rationale underpinning the use of in vitro models of molecular processes.

Functional and comparative analysis of THI1 gene in grasses with a focus on sugarcane

De novo synthesis of thiamine (vitamin B1) in plants depends on the action of thiamine thiazole synthase, which synthesizes the thiazole ring, and is encoded by the THI1 gene. Here, we investigated the evolution and diversity of THI1 in Poaceae, where C4 and C3 photosynthetic plants co-evolved. An ancestral duplication of THI1 is observed in Panicoideae that remains in many modern monocots, including sugarcane. In addition to the two sugarcane copies (ScTHI1-1 and ScTHI1-2), we identified ScTHI1-2 alleles showing differences in their sequence, indicating divergence between ScTHI1-2a and ScTHI1-2b. Such variations are observed only in the Saccharum complex, corroborating the phylogeny. At least five THI1 genomic environments were found in Poaceae, two in sugarcane, M. sinensis, and S. bicolor. The THI1 promoter in Poaceae is highly conserved at 300 bp upstream of the start codon ATG and has cis-regulatory elements that putatively bind to transcription factors associated with development, growth, development and biological rhythms. An experiment set to compare gene expression levels in different tissues across the sugarcane R570 life cycle showed that ScTHI1-1 was expressed mainly in leaves regardless of age. Furthermore, ScTHI1 displayed relatively high expression levels in meristem and culm, which varied with the plant age. Finally, yeast complementation studies with THI4-defective strain demonstrate that only ScTHI1-1 and ScTHI1-2b isoforms can partially restore thiamine auxotrophy, albeit at a low frequency. Taken together, the present work supports the existence of multiple origins of THI1 harboring genomic regions in Poaceae with predicted functional redundancy. In addition, it questions the contribution of the levels of the thiazole ring in C4 photosynthetic plant tissues or potentially the relevance of the THI1 protein activity.

Genomic characterization of ESBL-producing Salmonella Thompson isolates harboring mcr-9 from dead chick embryos in China

The emergence and dissemination of the extended spectrum β-lactamase (ESBL)-producing Enterobacteriaceae harbouring antimicrobial resistance (AMR) genes has diminished the potential options for treating multidrug-resistant (MDR) bacterial infections. Until now, numerous studies reported the spreading of critical plasmid-borne AMR genes from different sources worldwide. While the knowledge on the occurrence of the plasmid-borne AMR genes, especially mcr genes in the dead chick embryos, remains obscure. A retrospective study was conducted to detect the presence of the mcr genes in forty-five Salmonella enterica isolates recovered from 2139 dead chick embryo samples, from breeding chicken hatcheries in Henan, China. Using multiplex PCR, we found only four isolates out of the forty-five were mcr-9-positive. These four isolates were found to be MDR, ESBL- producing and showed resistance to 10 antimicrobial drugs. Additionally, mcr-9 harbouring plasmids were successfully transferred into Escherichia coli (E. coli) J53 by conjugation and the mcr-9 gene was confirmed by PCR. We also found that the transconjugants exhibited higher MICs for ampicillin, gentamycin and colistin than the recipient. Whole-genome sequence analysis showed that the four isolates belonged to Salmonella Thompson ST26 and harboured IncHI2 plasmid replicon. Furthermore, the mcr-9 harbouring plasmids were reconstructed using in silico tools and found to be carried other AMR genes (bla_DHA-1 and qnrB4). The studied isolates carried the typical virulence factors from Salmonella pathogenicity islands 1 and 2 (SPI-1 and SPI-2), in addition to pef and csg operons which are important in host adhesion and biofilm formation. The mgtC gene, which is involved in phagocytosis, has also been identified. Together, the increase in the phenotypic resistance of the transconjugants and the plasmid in silico reconstruction analysis confirmed that the corresponding resistance genes might be located together on the same plasmid. To track the potential phylogenomic relations of our detected ESBL S. Thompson isolates, we constructed a phylogenomic tree with available ESBL S. Thompson genomes (n = 26) that were reported worldwide. The studied isolates were independently clustered together with four other Chinese isolates of food origin in one clade, providing strong evidence of a potential recent and wide dissemination of ESBL S. Thompson across the food chain in China. In conclusion, we report the detection of four highly virulent ESBL-producing S. Thompson ST26 isolates harbouring mcr-9 gene obtained from dead chick embryos in Henan, China.

Colonization by ceftazidime/avibactam-resistant KPC-producing Klebsiella pneumoniae following therapy in critically ill patients

OBJECTIVES

Ceftazidime-avibactam (CAZ-AVI)-based treatments have been associated with the emergence of resistance in KPC-producing Klebsiella pneumoniae (KPC-Kp) isolates after antimicrobial exposure. Here, we evaluated the CAZ-AVI resistance development in KPC-Kp isolated from patients treated with CAZ-AVI-based therapy.

METHODS

We enrolled adult patients treated with CAZ-AVI-based regimens between January 2020 and January 2021. Carbapenemase-producing isolates collected from clinical samples and rectal swabs were evaluated for CAZ-AVI resistance development after antimicrobial exposure. KPC-Kp developing CAZ-AVI resistance and parental susceptible strains were genomically characterized. Whole genome sequencing was performed by using the Illumina iSeq100 platform and genomes were analyzed for antimicrobial-resistance genes, plasmid and porins sequences.

RESULTS

We enrolled 90 patients treated with CAZ-AVI-based therapy and 62.2% (56/90) of them were colonized by KPC-producers before CAZ-AVI-based treatment and 6.6% acquired colonization during therapy. Six (6.6%) patients developed infections because of resistant KPC-Kp after CAZ-AVI exposure and 3 (3.3%) of them developed CAZ-AVI resistance in the rectum. Development of resistance among KPC in the rectum occurred after 32 (IQR, 9-35) days of therapy and after 30 (IQR, 22-40) days in clinical specimens. Genetic analysis demonstrated that the development of CAZ-AVI resistance was associated with mutated bla_KPC-3 (bla_KPC-31, bla_KPC-53, bla_KPC-89, and bla_KPC-130) and phylogenetic analysis demonstrated a close genomic relationship between KCP-Kp collected from rectum and clinical samples of the same patient.

DISCUSSION

Antimicrobial exposure induce a higher incidence of CAZ-AVI resistance development in the blood and respiratory tract than in the rectum (6.7% vs. 3.3%) of CAZ-AVI-treated patients and genome analysis showed that resistance was associated with mutated bla_KPC-3 variants.

An ultrastructural and genomic study on the SARS-CoV-2 variant B.1.210 circulating during the first wave of COVID-19 pandemic in India

PURPOSE

The study aims to isolate and understand cytopathogenesis, ultrastructure, genomic characteristics and phylogenetic analysis of SARS-CoV-2 virus of B.1.210 lineage, that circulated in India during first wave of the pandemic.

METHODS

Clinical specimen from an interstate traveller from Maharashtra to Karnataka, in May 2020, who was positive by RT PCR for SARS-CoV-2 infection was subjected to virus isolation and Whole Genome Sequencing. Vero cells were used to study cytopathogenesis and ultrastructural features by Transmission Electron Microscopy (TEM). Phylogenetic analysis of the whole genome sequences of several SARS-CoV-2 variants downloaded from GISAID was performed in comparison with the B.1.210 variant identified in this study.

RESULTS

The virus was isolated in Vero cells and identified by immunofluorescence assay and RT PCR. The growth kinetics in infected Vero cells revealed a peak viral titre at 24 ​h post-infection. Ultrastructural studies revealed distinct morphological changes with accumulation of membrane-bound vesicles containing pleomorphic virions in the cytoplasm, with single or multiple intranuclear filamentous inclusions and dilated rough endoplasmic reticulum with viral particles. Whole genome sequence of the clinical specimen as well as the isolated virus revealed the virus to be of lineage B.1.210 with the D614G mutation in the spike protein. Phylogenetic analysis of the whole genome sequence in comparison with other variants reported globally revealed that the isolated SARS-CoV-2 virus of lineage B.1.210 is closely related to the original Wuhan virus reference sequence.

CONCLUSIONS

The SARS-CoV-2 variant B.1.210 virus isolated here showed ultrastructural features and cytopathogenesis similar to that of the virus reported during early phase of pandemic. Phylogenetic analysis showed that the isolated virus is closely related to the original Wuhan virus, thereby suggesting that the SARS-CoV-2 lineage B.1.210 that was circulating in India during the early phase of pandemic is likely to have evolved from the original Wuhan strain.

Genomic characterization reveals distinct mutation landscapes and therapeutic implications in neuroendocrine carcinomas of the gastrointestinal tract

BACKGROUND

Neuroendocrine carcinomas of the gastrointestinal tract (GI-NECs) remain a disease of grim prognosis with limited therapeutic options. Their molecular characteristics are still undefined. This study aimed to explore the underlying genetic basis and heterogeneity of GI-NECs.

METHODS

Comprehensive genomic analysis using whole-exome sequencing was performed on 143 formalin-fixed, paraffin-embedded samples of surgically resected GI-NEC with a thorough histological evaluation. Mutational signatures, somatic mutations, and copy number aberrations were analyzed and compared across anatomic locations and histological subtypes. Survival analysis was conducted to identify the independent factors.

RESULTS

In total, 143 GI-NECs were examined: the stomach, 87 cases (60.8%); the esophagus, 29 cases (20.3%); the colorectum, 20 cases (14.0%); and the small intestine, 7 cases (4.9%). Eighty-three (58.0%) and 60 (42.0%) cases were subclassified into small cell and large cell subtypes, respectively. GI-NECs showed distinct genetic alterations from their lung counterparts and non-neuroendocrine carcinomas in the same locations. Obvious heterogeneity of mutational signatures, somatic mutations, and copy number variations was revealed across anatomic locations rather than histological subtypes. Except for tumor protein p53 (TP53) and retinoblastoma 1 (RB1), the most frequently mutated genes in the stomach, esophagus, colorectum, and small intestine were low-density lipoprotein receptor-related protein 1B (LRP1B), notch receptor 1 (NOTCH1), adenomatosis polyposis coli (APC), catenin beta 1 (CTNNB1), respectively. Mutations in the WNT-β-catenin, NOTCH and erythroblastic leukemia viral oncogene B (ERBB) pathways were prevalently identified in gastric, esophageal, and colorectal NECs, respectively. Importantly, 104 (72.7%) GI-NECs harbored putative clinically relevant alterations, and non-gastric location and RB1 bi-allelic inactivation with copy number alterations were identified as two independent poor prognostic factors. Furthermore, we found that tumor cells in GI-NECs first gain clonal mutations in TP53, RB1, NOTCH1 and APC, followed by subsequent whole-genome doubling (WGD) and post-WGD clonal mutations in LRP1B, CUB and Sushi multiple domains 3 (CSMD3), FAT tumor suppressor homolog 4 (FAT4) and erb-b2 receptor tyrosine kinase 4 (ERBB4), and finally develop subclonal mutations.

CONCLUSIONS

GI-NECs harbor distinct genomic landscapes and demonstrate significant genetic heterogeneity across different anatomic locations. Moreover, potentially actionable alterations and prognostic factors were revealed for GI-NECs.

Comparative genomics analysis of two Helcococcus kunzii strains co-isolated with Staphylococcus aureus from diabetic foot ulcers

Helcococcus kunzii is a commensal Gram-positive bacterial species recovered from the human skin microbiota and considered as an opportunistic pathogen. Although little is known about its clinical significance, its increased abundance has been reported in infected wounds, particularly in foot ulcers in persons with diabetes. This species is usually detected in mixed cultures from human specimens and frequently isolated with Staphylococcus aureus. Modulation of staphylococci virulence by H. kunzii has been shown in an infection model of Caenorhabditis elegans. The aim of this study was to compare the genomes of two H. kunzii strains isolated from foot ulcers -isolate H13 and H10 showing high or low impact on S. aureus virulence, respectively- and the H. kunzii ATCC51366 strain. Whole genome analyses revealed some differences between the two strains: length (2.06 Mb (H13) and 2.05 Mb (H10) bp), GC content (29.3% (H13) and 29.5% (H10)) and gene content (1,884 (H13) and 1,786 (H10) predicted genes). The core-proteome phylogenies within the genus characterised H. kunzii H13 and H10 as genetically similar to their ancestor. The main differences between the strains were mainly in sugar-associated transporters and various hypothetical proteins. Five targets were identified as potentially involved in S. aureus virulence modulation in both genomes: the two-component iron export system and three autoinducer-like proteins. Moreover, H13 strain harbours a prophage inserted in 1,261,110-1,295,549 (attL-attR), which is absent in H10 strain. The prophage PhiCD38_2 was previously reported for its ability to modulate secretion profile, reinforcing the autoinducer-like hypothesis. In the future, transcriptomics or metaproteomics approaches could be performed to better characterize the H13 strain and possibly identify the underlying mechanism for S. aureus virulence modulation.

Reinfection or relapse? A case study of whole genome sequencing guided genomic characterization of Mycobacterium abscessus chronic infection in a cystic fibrosis patient

A 7-year-old cystic fibrosis patient with increased cough, new pulmonary infiltrate, and declining pulmonary function was diagnosed with clarithromycin resistant Mycobacterium abscessus infection. Treatment was initiated with clofazimine, linezolid and cefoxitin; she responded well to therapy and achieved microbiological clearance after completion of 12-month treatment. One year later, she had re-emergence of worsening symptoms and her sputum culture again grew clarithromycin resistant M. abscessus. Using a laboratory developed whole genome sequencing (WGS) test, the bacterium was determined to be the same strain with the same resistance mechanisms, indicating a relapse. This was deemed a critical element of clinical information, as the isolation of a genetically distinct organism would have indicated a new infection and would have served as evidence that a 12-month regimen was likely sufficient to achieve eradication. The confirmation of a relapse prompted the prolongation of the therapy plan to a goal of 24 months. Reinfection and relapse are great challenges in patients with cystic fibrosis who may acquire new strain of M. abscessus from the environment, may harbor multiple subpopulations of bacteria, or may have persistent infections but intermittent bacteria shedding that could not be eradicated. WGS has emerged as a powerful molecular tool to accurately differentiate re-infection from relapse thus solving this conundrum.

Cross-species transmission potential of chicken anemia virus and avian gyrovirus 2

The Gyrovirus genus consists of nonenveloped, icosahedral viruses with small circular single-stranded DNA genomes. Gyroviruses have been detected in diverse hosts, including humans, chickens, rodents, and cats. Two Gyroviruses were detected in canine serum samples using PCR in this study. The results indicated that four serum samples were positive for CAV (0.28%, 2/700) or AGV2 (0.28%, 2/700). Additionally, recombination analysis showed that AGV2 and CAV might have originated from the recombination of viruses similar to those detected in chickens and humans. We detected a total of 14 mutations in CAV VP1 amino acid sequences and identified new mutations at positions 31, 388, 390, 399, and 421 for the first time. The identification of T390C, C912T, T1230C, and T1297C mutations in AGV2 VP1, R93C mutations in AGV2 VP2, and R58C mutations AGV2 VP3 indicated that the differences might be related to a transboundary movement among hosts, which requires further elucidation. To the best of our knowledge, this study is the first report of an AGV2-infected dog in China, suggesting that the cross-species transmission of viruses with circular single-stranded DNA genomes is a public health concern.

Patient-derived xenograft models capture genomic heterogeneity in endometrial cancer

Background

Endometrial cancer (EC) is a major gynecological cancer with increasing incidence. It comprises four molecular subtypes with differing etiology, prognoses, and responses to chemotherapy. In the future, clinical trials testing new single agents or combination therapies will be targeted to the molecular subtype most likely to respond. As pre-clinical models that faithfully represent the molecular subtypes of EC are urgently needed, we sought to develop and characterize a panel of novel EC patient-derived xenograft (PDX) models.

Methods

Here, we report whole exome or whole genome sequencing of 11 PDX models and their matched primary tumor. Analysis of multiple PDX lineages and passages was performed to study tumor heterogeneity across lineages and/or passages. Based on recent reports of frequent defects in the homologous recombination (HR) pathway in EC, we assessed mutational signatures and HR deficiency scores and correlated these with in vivo responses to the PARP inhibitor (PARPi) talazoparib in six PDXs representing the copy number high/p53-mutant and mismatch-repair deficient molecular subtypes of EC.

Results

PDX models were successfully generated from grade 2/3 tumors, including three uterine carcinosarcomas. The models showed similar histomorphology to the primary tumors and represented all four molecular subtypes of EC, including five mismatch-repair deficient models. The different PDX lineages showed a wide range of inter-tumor and intra-tumor heterogeneity. However, for most PDX models, one arm recapitulated the molecular landscape of the primary tumor without major genomic drift. An in vivo response to talazoparib was detected in four copy number high models. Two models (carcinosarcomas) showed a response consistent with stable disease and two models (one copy number high serous EC and another carcinosarcoma) showed significant tumor growth inhibition, albeit one consistent with progressive disease; however, all lacked the HR deficiency genomic signature.

Conclusions

EC PDX models represent the four molecular subtypes of disease and can capture intra-tumor heterogeneity of the original primary tumor. PDXs of the copy number high molecular subtype showed sensitivity to PARPi; however, deeper and more durable responses will likely require combination of PARPi with other agents.

Supplementary Information

The online version contains supplementary material available at 10.1186/s13073-021-00990-z.

Characterization and evolutionary analysis of Cucumber mosaic virus isolate infecting Salvia sclarea in India

Cucumber mosaic virus (CMV) is one of the most widespread viruses that infects a large number of cultivated crops worldwide and causes severe losses. Besides vegetables and ornamental crops, it is also spreading on medicinal and aromatic plants (MAPs) at an alarming rate. Natural occurrence of green mosaic, leaf curling, necrosis, and distortion were observed on the leaves of Salvia sclarea in the experimental field of CSIR-CIMAP, Lucknow, India. Mechanical transmission assay, morphological features, CP gene, and RNA3-based genomic characterization revealed the association of a new CMV isolate with the present disease. Virus is mechanically sap transmissible to the test plants. Transmission electron microscopy showed the presence of isodiametric particles of ~ 28 nm in diameter. Phylogenetic studies revealed that RNA3 (~ 2.2 kb) belongs to the subgroup IA of CMV. Multiple sequence alignment of amino acid sequences showed that MP possesses two unique changes in the RNA-binding domain and CP was found to be the conserved one with one change only. To the best of our knowledge, this is the first report of a CMV isolate infecting S. sclarea.

Comparison of Genomic Characterization in Upper Tract Urothelial Carcinoma and Urothelial Carcinoma of the Bladder

BACKGROUND

Different genomic characterization in urothelial carcinoma (UC) by site of origin may imply contrasting therapeutic opportunities and pathogenetic mechanisms. The aim of this study was to investigate whether differences between upper tract UC (UTUC) and UC of the bladder (UCB) result from intrinsic biological diversity.

MATERIALS AND METHODS

We prospectively sequenced 118 tumors and matched blood DNA from Chinese patients with UC using next-generation sequencing techniques, including 45 UTUC and 73 UCB. Two hundred twenty-six patients with UTUC and 350 patients with UCB for The Cancer Genome Atlas were acquired from the cbioportal.

RESULTS

There were marked disparities in the mutational landscape for UC according to race and site of origin. Signature 22 for exposure to aristolochic acid was only observed in the UTUC cohort. Conversely, signature 6 for defective DNA mismatch repair only existed in the UCB cohort. Compared with UCB, UTUC had higher clonal and subclonal mutation numbers. TP53, PIK3CA, and FGFR3 mutations may be the driver genes for UTUC, whereas for UCB, the driver gene may be BRCA1. Patients with UTUC had lower PD-L1 than those with UCB. There was no significant difference in the number of DDR mutations, copy number variation counts, and tumor mutational burden between UTUC and UCB.

CONCLUSION

UTUC and UCB exhibit significant differences in the prevalence of genomic landscape and carcinogenesis. Consequently, molecular subtypes differ according to location, and these results may imply the site-specific management of patients with urothelial carcinoma. Mutational signature may be used as a screening tool to assist clinical differential diagnosis between UTUC and UCB.

IMPLICATIONS FOR PRACTICE

This study's findings lay the foundation for a deeper understanding of distinct molecular mechanisms and similar treatment opportunities between upper tract urothelial carcinoma (UTUC) and urothelial carcinoma of the bladder (UCB) and had important implications for the site-specific management of patients with urothelial carcinoma. A comprehensive understanding of the biology of UTUC and UCB is needed to identify new drug targets in order to improve clinical outcomes.

First detection and complete genome analysis of the Lyon IARC polyomavirus in China from samples of diarrheic cats

Lyon IARC polyomavirus (LIPyV), a newly discovered polyomavirus (PyV), was first identified in 2017 in human skin samples in the USA. Later, it was detected in several other countries in samples of human and feline origin. Our aim was to find out if the virus occurs in China. To this end, 100 fecal samples were collected from cats with diarrhea in Guangxi Province during 2016 and 2018 and tested with polymerase chain reaction (PCR). Only 2 samples that originated from two related individuals were found to be positive. Based on the sequence identity of the 240-bp PCR products, the two positive samples supposedly contained identical viruses. Therefore, only one of them, which was designated as LIPyV-GXNN01, was selected for full genome amplification, cloning, sequencing and analysis. LIPyV-GXNN01, which comprises 5,263 nucleotides, has an early region that consists of small T antigen (ST-Ag) and large T antigen (LT-Ag) and a late region coding for the VP1, VP2, and VP3 structural proteins. Moreover, the LIPyV-GXNN01 strain structural proteins share 95.9–99.4%, 97.6–99.2%, and 97.1–99.2% nucleic acid identity with the VP1, VP2, and VP3of other LIPyV reference strains, respectively. A phylogenetic analysis revealed that GXNN01 clustered together with previously reported LIPyV strain. This present study is the first report of LIPyV in China.

Allelic and genotypic frequencies for loci associated with meat quality in Mexican Braunvieh cattle

The objective was to estimate allelic and genotypic frequencies for loci associated with meat quality in a Mexican population of Braunvieh cattle. Information was obtained from 300 animals genotyped with the Genomic Profile Bovine LD chip of 30K and 50K SNPs. After the final edition, including quality control, the data contained information for 12 loci of the CAPN1, CAPN3, CAPN5, CAPN14, DGAT1, DGAT2, TG, ANK1, and MADH3 genes. Allelic and genotypic frequencies and Hardy-Weinberg equilibrium were estimated with the Cervus 3.0.7 software. The studied population markers were in Hardy-Weinberg equilibrium, except for those associated with CAPN1, DGAT1, and MADH3. Frequencies higher than those reported for other breeds were found for genotypes associated with meat softness, higher marbling score, lower quantity of saturated fatty acids, and lower shear force (CAPN1 and DGAT2). There were similarities with frequencies reported for Bos taurus breeds for the CAPN3 and TG genes. For the DGAT1 and ANK1 genes, the frequencies of the desired genotypes were low. A marker for DGAT1 and another for MADH3 were monomorphic. The results of this study are encouraging in terms of the potential of the Braunvieh population studied for breeding programs aiming to increase meat quality. The breed has strengths that could be used either by crossbreeding to generate heterozygous animals or by selection to increase frequencies of valuable alleles.

Isolation and Characterization of Bacteriophages from Inland Saline Aquaculture Environments to Control Vibrio parahaemolyticus Contamination in Shrimp

Among the various bacterial pathogens associated with the aquaculture environment, Vibrio parahaemolyticus the important one from shrimp and human health aspects. Though having been around for several decades, phage-based control of bacterial pathogens (phage therapy) has recently re-emerged as an attractive alternative due to the availability of modern phage characterization tools and the global emergence of antibiotic-resistant bacteria. In the present study, a total of 12 V. parahaemolyticus specific phages were isolated from 264 water samples collected from inland saline shrimp culture farms. During the host range analysis against standard/field isolates of V. parahaemolyticus and other bacterial species, lytic activity was observed against 2.3-45.5% of tested V. parahaemolyticus isolates. No lytic activity was observed against other bacterial species. For genomic characterization, high-quality phage nucleic acid with concentrations ranging from 7.66 to 220 ng/µl was isolated from 9 phages. After digestion treatments with DNase, RNase and S1 nuclease, the nature of phage nucleic acid was determined as ssDNA and dsDNA for 7 and 2 phages respectively. During transmission electron microscopy analysis of phage V5, it was found to have a filamentous shape making it a member of the family Inoviridae. During efficacy study of phage against V. parahaemolyticus in shrimp, 78.1% reduction in bacterial counts was observed within 1 h of phage application. These results indicate the potential of phage therapy for the control of V. parahaemoyticus in shrimp.

Supplementary Information

The online version contains supplementary material available at 10.1007/s12088-021-00934-6.

Genomic characterization and expression profiles of stress-associated proteins (SAPs) in castor bean (Ricinus communis)

Stress-associated proteins (SAPs) are known as response factors to multiple abiotic and biotic stresses in plants. However, the potential physiological and molecular functions of SAPs remain largely unclear. Castor bean (Ricinus communis L.) is one of the most economically valuable non-edible woody oilseed crops, able to be widely cultivated in marginal lands worldwide because of its broad adaptive capacity to soil and climate conditions. Whether SAPs in castor bean plays a key role in adapting diverse soil conditions and stresses remains unknown. In this study, we used the castor bean genome to identify and characterize nine castor bean SAP genes (RcSAP). Structural analysis showed that castor bean SAP gene structures and functional domain types vary greatly, differing in intron number, protein sequence, and functional domain type. Notably, the AN1-C2H2-C2H2 zinc finger domain within RcSAP9 has not been often observed in other plant families. High throughput RNA-seq data showed that castor bean SAP gene profiles varied among different tissues. In addition, castor bean SAP gene expression varied in response to different stresses, including salt, drought, heat, cold and ABA and MeJA, suggesting that the transcriptional regulation of castor bean SAP genes might operate independently of each other, and at least partially independent from ABA and MeJA signal pathways. Cis-element analyses for each castor bean SAP gene showed that no common cis-elements are shared across the nine castor bean SAP genes. Castor bean SAPs were localized to different regions of cells, including the cytoplasm, nucleus, and cytomembrane. This study provides a comprehensive profile of castor bean SAP genes that advances our understanding of their potential physiological and molecular functions in regulating growth and development and their responses to different abiotic stresses.

Genomic characterization and phylogenetic analysis of a new canine picornavirus variant in the mainland of China

A new canine picornavirus (CanPV) variant, designated as dog/SH1901/CHN/2019, was detected in a pool of various canine fecal samples in the mainland of China using a viral metagenomic analysis, and its nearly complete genome sequence was determined and analyzed. Sequence analyses revealed that it had a standard picornavirus genome organization, a type I internal ribosome entry site (IRES) in the 5'UTR. However, dog/SH1901/CHN/2019 has generated a serial of unique aa mutations and 7aa insertion when compared with the closely related CanPVs. Phylogenetic analysis and pairwise sequence comparisons based on the P1, 2C, 3C, and 3D protein sequences showed that dog/SH1901/ CHN/2019 was closely related to CanPV strains 244 F, 325 F and 6D, which clustered into an independent evolutionary clade and distantly related to CanPV strain A128thr of the genus Mischivirus, which indicated the unclassified CanPV strains may belong to a novel species or genus in the family Picornaviridae. This study extends our knowledge on the evolution and genetic diversity of CanPVs.

Genomic characterization of covert mortality nodavirus from farming shrimp: Evidence for a new species within the family Nodaviridae

The prevalence of covert mortality nodavirus (CMNV) has become one of the major threats to the shrimp farming industry in Asia and South America recently. Here, the genomic RNA1 and RNA2 of CMNV were characterized by using transcriptome sequencing and RT-PCR. Our study revealed that RNA1 is 3228 bp in length, and contains two putative Open Reading Frames (ORFs), one encoding the RNA dependent RNA polymerase (RdRp) of length 1043 amino acids and another encoding the protein B2 with a length of 132 amino acids. RNA2 is 1448 bp in length and encodes a capsid protein of 437 amino acids. CMNV shared the highest similarity of 51.78 % for RdRp with the other known nodaviruses. Phylogenetic analyses on the basis of RdRp, B2 and capsid proteins indicated that CMNV might represent a novel viral species in the family Nodaviridae. This study reported the first genome sequence of CMNV and it would be helpful for further studies of CMNV in relation to its evolution, diagnostic technique and control strategy.

Fur Seal Feces-Associated Circular DNA Virus Identified in Pigs in Anhui, China

Fur seal feces-associated circular DNA virus (FSfaCV) is an unclassified circular replication-associated protein (Rep)-encoding single-stranded (CRESS) DNA virus that has been detected in mammals (fur seals and pigs). The biology and epidemiology of the virus remain largely unknown. To investigate the virus diversity among pigs in Anhui Province, China, we pooled 600 nasal samples in 2017 and detected viruses using viral metagenomic methods. From the assembled contigs, 12 showed notably high nucleotide acid sequence similarities to the genome sequences of FSfaCVs. Based on these sequences, a full-length genome sequence of the virus was then obtained using overlapping PCR and sequencing, and the virus was designated as FSfaCV-CHN (GenBank No. MK462122). This virus shared 91.3% and 90.9% genome-wide nucleotide sequence similarities with the New Zealand fur seal strain FSfaCV-as50 and the Japanese pig strain FSfaCV-JPN1, respectively. It also clustered with the two previously identified FSfaCVs in a unique branch in the phylogenetic tree based on the open reading frame 2 (ORF2), Rep-coding gene, and the genome of the reference CRESS DNA viruses. Further epidemiological investigation using samples collected in 2018 showed that the overall positive rate for the virus was 56.4% (111/197) in Anhui Province. This is the first report of FSfaCVs identified in pigs in China, and further epidemiological studies are warranted to evaluate the influence of the virus on pigs.

Genomic characterization of a novel SARS-CoV-2

A new severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) associated with human to human transmission and extreme human sickness has been as of late announced from the city of Wuhan in China. Our objectives were to mutation analysis between recently reported genomes at various times and locations and to characterize the genomic structure of SARS-CoV-2 using bioinformatics programs. Information on the variation of viruses is of considerable medical and biological impacts on the prevention, diagnosis, and therapy of infectious diseases. To understand the genomic structure and variations of the SARS-CoV-2. The study analyzed 95 SARS-CoV-2 complete genome sequences available in GenBank, National MicrobiologyData Center (NMDC) and NGDC Genome Warehouse from December-2019 until 05 of April-2020. The genomic signature analysis demonstrates that a strong association between the time of sample collection, location of sample and accumulation of genetic diversity. We found 116 mutations, the three most common mutations were 8782C>T in ORF1ab gene, 28144T>C in ORF8 gene and 29095C>T in the N gene. The mutations might affect the severity and spread of the SARS-CoV-2. The finding heavily supports an intense requirement for additional prompt, inclusive investigations that combine genomic detail, epidemiological information and graph records of the clinical features of patients with COVID-19.

A divergent canine parvovirus type 2c (CPV-2c) isolate circulating in China

The variability and the intrinsically high mutation rate of canine parvovirus type 2(CPV-2) increased the diversity of CPV-2 in canine populations. Since the first occurrence of CPV-2, three antigenic variants (2a, 2b and 2c) were detected and distributed worldwide. CPV-2c infection has been detected and increasingly reported in China. Here, a CPV-2c strain CPV-SH1516 was isolated and its complete genome sequence was first characterized. Compared with other CPV-2c isolates, CPV-2c isolates from China continued to evolve into divergent CPV-2c variants with specific unique amino acid substitutions under purifying selection. Emergence of CPV-2c isolates from China was driven by the unique gradual point mutations in key sites of VP2 rather than introduction from outside China. Combining sequence comparison with phylogenetic analysis based on the amino acid sequences of VP2, the vast majority of CPV-2c isolates from China formed a monophyletic cluster and CPV-SH1516 was a representative isolate of CPV-2c circulating in China. Overall, our study provides valuable insight into the evolutionary mechanism of CPV-2c.

DNA Microarray-Based Genomic Characterization of the Pathotypes of Escherichia coli O26, O45, O103, O111, and O145 Isolated from Feces of Feedlot Cattle †

Shiga toxin-producing Escherichia coli (STEC) serogroups O26, O45, O103, O111, O121, and O145, referred to as the top six non-O157 serogroups, are responsible for more than 70% of human non-O157 STEC infections in North America. Cattle harbor non-O157 strains in the hindgut and shed them in the feces. The objective of this study was to use the U.S. Food and Drug Administration (FDA) E. coli identification (ECID) DNA microarray to identify the serotype, assess the virulence potential of each, and determine the phylogenetic relationships among five of the six non-O157 E. coli serogroups isolated from feedlot cattle feces. Forty-four strains of STEC, enterohemorrhagic E. coli (EHEC), enteropathogenic E. coli (EPEC), or putative nonpathotype E. coli (NPEC) of cattle origin and five human clinical strains of EHEC were assayed with the FDA-ECID DNA microarray. The cattle strains harbored diverse flagellar genes. The bovine and human strains belonging to serogroups O26, O45, and O103 carried stx₁ only, O111 carried both stx₁ and stx₂, and O145 carried either stx₁ or stx₂. The strains were also positive for various subtypes of intimin and other adhesins (IrgA homologue adhesin, long polar fimbriae, mannose-specific adhesin, and curli). Both human and cattle strains were positive for LEE-encoded type III secretory system genes and non-LEE-encoded effector genes. SplitsTree4, a program used to determine the phylogenetic relationship among the strains, revealed that the strains within each serogroup clustered according to their pathotype. In addition to genes encoding Shiga toxins, bovine non-O157 E. coli strains possessed other major virulence genes, including those for adhesins, type III secretory system proteins, and plasmid-borne virulence genes, similar to human clinical strains. Because virulence factors encoded by these genes are involved in the pathogenesis of various pathotypes of E. coli, the bovine non-O157 strains could cause human illness. The FDA-ECID DNA microarray assay rapidly provided a profile of the virulence genes for assessment of the virulence potential of each strain.

High prevalence of Gammapapillomaviruses (Gamma-PVs) in pre-malignant cutaneous lesions of immunocompetent individuals using a new broad-spectrum primer system, and identification of HPV210, a novel Gamma-PV type

Genus Gammapapillomavirus (Gamma-PV) is the most diverse and largest clade within the Papillomaviridae family. A novel set of degenerate primers targeting the E1 gene was designed and further used in combination with the well-known CUT PCR assay to assess HPV prevalence and genus distribution in a variety of cutaneous samples from 448 immunocompetent individuals. General HPV, Gamma-PV and mixed infections prevalence were significantly higher in actinic keratosis with respect to benign and malignant neoplasms, respectively (p = 0.0047, p = 0.0172, p = 0.00001). Gamma-PVs were significantly more common in actinic keratosis biopsies than Beta- and Alpha-PVs (p = 0.002). The full-length genome sequence of a novel putative Gamma-PV type was amplified by 'hanging droplet' long-range PCR and cloned. The novel virus, designated HPV210, clustered within species Gamma-12. This study provides an additional tool enabling detection of HPV infections in skin and adds new insights about possible early roles of Gamma-PVs in the development of cutaneous malignant lesions.

Complete genome sequencing of Kaisodi virus isolated from ticks in India belonging to Phlebovirus genus, family Phenuiviridae

An unknown virus was repeatedly isolated from hard tick (Haemaphysalis spinigera) during a proactive arbovirus survey in ticks conducted in 1957, in India. The virus remained uncharacterized for a long time. The passages of this virus in different vertebrate and invertebrate cells along with human and monkey-derived cell culture showed no cytopathic effect. It was identified later to be a member of Kaisodi group among Phlebovirus genus in the family Phenuiviridae (Order: Bunyavirales) by serological methods. Due to its genomic diversity, sequencing of this virus was a challenge for a while. In this study, we were able to sequence the complete genome of this virus isolate using next-generation sequencing (NGS) platform. The unknown virus was identified to be Kaisodi virus (KASDV) using NGS analysis. De novo genome assembly derived three genomic segments for the KASDV which encode for RNA-dependent RNA polymerase, glycoprotein precursor, and nucleoprotein. Functional as well as conserved domains for Kaisodi serogroup viruses were predicted and compared to a known representative of the genus Phlebovirus. The phylogenetic tree revealed its closeness to Silverwater virus, of Kaisodi serogroup with nucleotide (69%, 62%, and 61%) and amino acid (52%, 51%, and 62%) identity for L, M, and S segment, respectively. The study demonstrates the presence of a conserved motif (⁷²TRGNK⁷⁶) around the RNA binding motif region in tick-borne phleboviruses. The intergenic region encompassing the S segment of Kaisodi serogroup was GC-rich whereas the other Phlebovirus had AT-rich genome. KASDV has the largest intergenic region and larger loops, suggesting stem-loops formed due to larger loops as a possible factor for instability and cause of transcription termination. This paper also describes the real-time RT-PCR and RT-PCR assays developed and used for the detection of KASDV RNA in ticks from Karnataka, Kerala and Maharashtra State, India. The KASDV positivity observed in the recently collected tick pools indicates that the KASDV, isolated from Karnataka state in 1957, is also circulating in the adjoining Kerala state. On the basis of the current study, it should be possible to develop diagnostic assays which would facilitate an in-depth field survey exploring the veterinary and medical significance of KASDV.

Genomic characterization of colitis-associated colorectal cancer

BACKGROUND

Inflammatory bowel disease (IBD), which includes ulcerative colitis (UC) and Crohn's disease (CD), is a chronic, idiopathic, repeated inflammatory disease. Colorectal cancer (CRC) that develops in patients with IBD is known as colitis-associated colorectal cancer (CAC), but the underlying carcinogenic mechanism remains unclear. Genomic analysis of sporadic CRC has been well described based on next-generation sequencing (NGS) data. Using NGS, we compared all exons of 415 cancer-associated genes in patients in Japan and the USA who had CRC and found similar genomic alteration patterns among the two populations. However, genomic analysis of CAC has not been thoroughly investigated.

MAIN BODY

The molecular pathogenesis of CAC shares many features with sporadic CRC, but there are distinct variations in the time and frequency of some alterations. Gene alterations in CAC are gradually being elucidated using genomic sequencing analyses. Some studies have shown that gene alteration patterns differ between UC and CD. The carcinogenesis of CAC depends on unique environmental, genetic, and immunological factors.

CONCLUSIONS

In this review, we have discussed the differences in genomic alterations between sporadic CRC and CAC. NGS in patients with IBD has the potential to detect early CAC and to suggest therapeutic targets.

Analysis of enterotoxigenic Bacillus cereus strains from dried foods using whole genome sequencing, multi-locus sequence analysis and toxin gene prevalence and distribution using endpoint PCR analysis

Bacillus cereus strains were isolated from dried foods, which included international brands of spices from South East Asia, Mexico and India purchased from several retail stores, samples of powdered infant formula (PIF), medicated fish feed and dietary supplements. The genetic diversity of 64 strains from spices and PIF was determined using a multiplex endpoint PCR assay designed to identify hemolysin BL, nonhemolytic enterotoxin, cytotoxin K, and enterotoxin FM toxin genes. Thirteen different B. cereus toxigenic gene patterns or profiles were identified among the strains. Randomly selected B. cereus strains were sequenced and compared with reference Genomic Groups from National Center Biotechnology Information using bioinformatics tools. A comprehensive multi-loci sequence analysis (MLSA) was designed using alleles from 25 known MLST genes specifically tailored for use with whole genome assemblies. A cohort of representative genomes of strains from a few FDA regulated commodities like dry foods and medicated fish feed was used to demonstrate the utility of the 25-MLSA approach for rapid clustering and identification of Genome Groups. The analysis clustered the strains from medicated fish feed, dry foods, and dietary supplements into phylogenetically-related groups. 25-MLSA also pointed to a greater diversity of B. cereus strains from foods and feed than previously recognized. Our integrated approach of toxin gene PCR, and to our knowledge, whole genome sequencing (WGS) based sequence analysis, may be the first of its kind that demonstrates enterotoxigenic potential and genomic diversity in parallel.

Genomic characterization and evolution of Tacaiuma orthobunyavirus (Peribunyaviridae family) isolated in Brazil

Tacaiuma virus (TCMV) is antigenically characterized as a member of the Anopheles A complex in the Orthobunyavirus genus, Peribunyaviridae family (Bunyavirales order). Clinically, the TCMV infection is characterized by acute febrile illness with myalgia and arthralgia lasting three to five days. However, the genomic and evolutionary aspect of this virus has not been elucidated. In this study, we described the complete coding sequences of three segments of two TCMV strains isolated in Brazil and three complete coding sequences of the small segment of three TCMV strains. All the strains sequenced in this study showed the typical genomic organization of orthobunyaviruses that infect vertebrates, except for the absence of the open reading frame that encodes the well-described non-structural small protein. This study presents the genomic and evolutionary characterization of TCMV strains and would be helpful for diagnostic purposes and epidemiology.

HER2 is not a cancer subtype but rather a pan-cancer event and is highly enriched in AR-driven breast tumors

BACKGROUND

Approximately one in five breast cancers are driven by amplification and overexpression of the human epidermal growth factor receptor 2 (HER2) receptor kinase, and HER2-enriched (HER2E) is one of four major transcriptional subtypes of breast cancer. We set out to understand the genomics of HER2 amplification independent of subtype, and the underlying drivers and biology of HER2E tumors.

METHODS

We investigated published genomic data from 3155 breast tumors and 5391 non-breast tumors.

RESULTS

HER2 amplification is a distinct driver event seen in all breast cancer subtypes, rather than a subtype marker, with major characteristics restricted to amplification and overexpression of HER2 and neighboring genes. The HER2E subtype has a distinctive transcriptional landscape independent of HER2A that reflects androgen receptor signaling as replacement for estrogen receptor (ER)-driven tumorigenesis. HER2 amplification is also an event in 1.8% of non-breast tumors.

CONCLUSIONS

These discoveries reveal therapeutic opportunities for combining anti-HER2 therapy with anti-androgen agents in breast cancer, and highlight the potential for broader therapeutic use of HER2 inhibitors.

Genomic characterization and phylogenetic analysis of Zika virus circulating in the Americas

The rapid spread and potential link with birth defects have made Zika virus (ZIKV) a global public health problem. The virus was discovered 70years ago, yet the knowledge about its genomic structure and the genetic variations associated with current ZIKV explosive epidemics remains not fully understood. In this review, the genome organization, especially conserved terminal structures of ZIKV genome were characterized and compared with other mosquito-borne flaviviruses. It is suggested that major viral proteins of ZIKV share high structural and functional similarity with other known flaviviruses as shown by sequence comparison and prediction of functional motifs in viral proteins. Phylogenetic analysis demonstrated that all ZIKV strains circulating in the America form a unique clade within the Asian lineage. Furthermore, we identified a series of conserved amino acid residues that differentiate the Asian strains including the current circulating American strains from the ancient African strains. Overall, our findings provide an overview of ZIKV genome characterization and evolutionary dynamics in the Americas and point out critical clues for future virological and epidemiological studies.

Molecular characterization of three Zika flaviviruses obtained from sylvatic mosquitoes in the Central African Republic

Zika virus (ZIKV) is an emerging pathogen belonging to the Spondweni serocomplex within the genus Flavivirus. It has been isolated from several mosquito species. Two lineages of ZIKV have been defined by polyprotein homology. Using high-throughput sequencing, we obtained and characterized three complete genomes of ZIKV isolated between 1976 and 1980 in the Central African Republic. The three viruses were isolated from two species of mosquito, Aedes africanus and Ae. opok. Two sequences from Ae. africanus had 99.9% nucleotide sequence identity and 100% amino acid identity, whereas the complete genome obtained from Ae. opok had 98.3% nucleotide identity and 99.4% amino acid identity with the other two genomes. Phylogenetic analysis based on the amino acid sequence of the polyprotein showed that the three ZIKV strains clustered together but diverged from all other ZIKV strains. Our molecular data suggest that a different subtype of West African ZIKV strains circulated in Aedes species in Central Africa.