Whole-genome Single Nucleotide Polymorphism Research Articles

High-Throughput Genomics Identify Novel FBN1/2 Variants in Severe Neonatal Marfan Syndrome and Congenital Heart Defects.

Fibrillin-1 and fibrillin-2, encoded by FBN1 and FBN2, respectively, play significant roles in elastic fiber assembly, with pathogenic variants causing a diverse group of connective tissue disorders such as Marfan syndrome (MFS) and congenital contractural arachnodactyly (CCD). Different genomic variations may lead to heterogeneous phenotypic features and functional consequences. Recent high-throughput sequencing modalities have allowed detection of novel variants that may guide the care for patients and inform the genetic counseling for their families. We performed clinical phenotyping for two newborn infants with complex congenital heart defects. For genetic investigations, we employed next-generation sequencing strategies including whole-genome Single-Nucleotide Polymorphism (SNP) microarray for infant A with valvular insufficiency, aortic sinus dilatation, hydronephrosis, and dysmorphic features, and Trio whole-exome sequencing (WES) for infant B with dextro-transposition of the great arteries (D-TGA) and both parents. Infant A is a term male with neonatal marfanoid features, left-sided hydronephrosis, and complex congenital heart defects including tricuspid regurgitation, aortic sinus dilatation, patent foramen ovale, patent ductus arteriosus, mitral regurgitation, tricuspid regurgitation, aortic regurgitation, and pulmonary sinus dilatation. He developed severe persistent pulmonary hypertension and worsening acute hypercapnic hypoxemic respiratory failure, and subsequently expired on day of life (DOL) 10 after compassionate extubation. Cytogenomic whole-genome SNP microarray analysis revealed a deletion within the FBN1 gene spanning exons 7-30, which overlapped with the exon deletion hotspot region associated with neonatal Marfan syndrome. Infant B is a term male prenatally diagnosed with isolated D-TGA. He required balloon atrial septostomy on DOL 0 and subsequent atrial switch operation, atrial septal defect repair, and patent ductus arteriosus ligation on DOL 5. Trio-WES revealed compound heterozygous c.518C>T and c.8230T>G variants in the FBN2 gene. Zygosity analysis confirmed each of the variants was inherited from one of the parents who were healthy heterozygous carriers. Since his cardiac repair at birth, he has been growing and developing well without any further hospitalization. Our study highlights novel FBN1/FBN2 variants and signifies the phenotype-genotype association in two infants affected with complex congenital heart defects with and without dysmorphic features. These findings speak to the importance of next-generation high-throughput genomics for novel variant detection and the phenotypic variability associated with FBN1/FBN2 variants, particularly in the neonatal period, which may significantly impact clinical care and family counseling.

Strengthening Bordetella pertussis genomic surveillance by direct sequencing of residual positive specimens.

Whole-genome sequencing (WGS) of microbial pathogens recovered from patients with infectious disease facilitates high-resolution strain characterization and molecular epidemiology. However, increasing reliance on culture-independent methods to diagnose infectious diseases has resulted in few isolates available for WGS. Here, we report a novel culture-independent approach to genome characterization of Bordetella pertussis, the causative agent of pertussis and a paradigm for insufficient genomic surveillance due to limited culture of clinical isolates. Sequencing libraries constructed directly from residual pertussis-positive diagnostic nasopharyngeal specimens were hybridized with biotinylated RNA "baits" targeting B. pertussis fragments within complex mixtures that contained high concentrations of host and microbial background DNA. Recovery of B. pertussis genome sequence data was evaluated with mock and pooled negative clinical specimens spiked with reducing concentrations of either purified DNA or inactivated cells. Targeted enrichment increased the yield of B. pertussis sequencing reads up to 90% while simultaneously decreasing host reads to less than 10%. Filtered sequencing reads provided sufficient genome coverage to perform characterization via whole-genome single nucleotide polymorphisms and whole-genome multilocus sequencing typing. Moreover, these data were concordant with sequenced isolates recovered from the same specimens such that phylogenetic reconstructions from either consistently clustered the same putatively linked cases. The optimized protocol is suitable for nasopharyngeal specimens with diagnostic IS481 Ct < 35 and >10 ng DNA. Routine implementation of these methods could strengthen surveillance and study of pertussis resurgence by capturing additional cases with genomic characterization.

Targeted enrichment of whole-genome SNPs from highly burned skeletal remains.

Genetic assessment of highly incinerated and/or degraded human skeletal material is a persistent challenge in forensic DNA analysis, including identifying victims of mass disasters. Few studies have investigated the impact of thermal degradation on whole-genome single-nucleotide polymorphism (SNP) quality and quantity using next-generation sequencing (NGS). We present whole-genome SNP data obtained from the bones and teeth of 27 fire victims using two DNA extraction techniques. Extracts were converted to double-stranded DNA libraries then enriched for whole-genome SNPs using unpublished biotinylated RNA baits and sequenced on an Illumina NextSeq 550 platform. Raw reads were processed using the EAGER (Efficient Ancient Genome Reconstruction) pipeline, and the SNPs filtered and called using FreeBayes and GATK (v. 3.8). Mixed-effects modeling of the data suggest that SNP variability and preservation is predominantly determined by skeletal element and burn category, and not by extraction type. Whole-genome SNP data suggest that selecting long bones, hand and foot bones, and teeth subjected to temperatures <350°C are the most likely sources for higher genomic DNA yields. Furthermore, we observed an inverse correlation between the number of captured SNPs and the extent to which samples were burned, as well as a significant decrease in the total number of SNPs measured for samples subjected to temperatures >350°C. Our data complement previous analyses of burned human remains that compare extraction methods for downstream forensic applications and support the idea of adopting a modified Dabney extraction technique when traditional forensic methods fail to produce DNA yields sufficient for genetic identification.

Major impacts of widespread structural variation on sorghum.

Genetic diversity is critical to crop breeding and improvement, and dissection of the genomic variation underlying agronomic traits can both assist breeding and give insight into basic biological mechanisms. Although recent genome analyses in plants reveal many structural variants (SVs), most current studies of crop genetic variation are dominated by single-nucleotide polymorphisms (SNPs). The extent of the impact of SVs on global trait variation, as well as their utility in genome-wide selection, is not yet understood. In this study, we built an SV data set based on whole-genome resequencing of diverse sorghum lines (n = 363), validated the correlation of photoperiod sensitivity and variety type, and identified SV hotspots underlying the divergent evolution of cellulosic and sweet sorghum. In addition, we showed the complementary contribution of SVs for heritability of traits related to sorghum adaptation. Importantly, inclusion of SV polymorphisms in association studies revealed genotype-phenotype associations not observed with SNPs alone. Three-way genome-wide association studies (GWAS) based on whole-genome SNP, SV, and integrated SNP + SV data sets showed substantial associations between SVs and sorghum traits. The addition of SVs to GWAS substantially increased heritability estimates for some traits, indicating their important contribution to functional allelic variation at the genome level. Our discovery of the widespread impacts of SVs on heritable gene expression variation could render a plausible mechanism for their disproportionate impact on phenotypic variation. This study expands our knowledge of SVs and emphasizes the extensive impacts of SVs on sorghum.

Genome analyses reveal the secondary metabolites potentially influence the geographical distribution of Fusarium pseudograminearum populations.

Fusarium crown rot (FCR), caused by Fusarium pseudograminearum, significantly impacts wheat yield and quality in China's Huanghuai region. The rapid F. pseudograminearum epidemic and FCR outbreak within a decade remain unexplained. In this study, two high-quality, chromosome-level genomes of F. pseudograminearum strains producing 3-acetyl-deoxynivalenol (3AcDON) and 15-acetyl-deoxynivalenol (15AcDON) toxins were assembled. Additionally, 38 related strains were resequenced. Genomic differences such as single nucleotide polymorphisms (SNPs), insertions/deletions (indels), and structural variations (SVs) among F. pseudograminearum strains were analyzed. The whole-genome SNP locus based population classification mirrored the toxin chemotype (3AcDON and 15AcDON)-based classification, indicating the presence of genes associated with the trichothecene toxin gene cluster. Further analysis of differential SNP, indel, and SV loci between the 3AcDON and 15AcDON populations revealed a predominant connection to secondary metabolite synthesis genes. Notably, the majority of the secondary metabolite biosynthesis gene cluster (SMGC) loci were located in SNP-dense genomic regions, suggesting high mutability and a possible contribution to F. pseudograminearum population structure and environmental adaptability. This study provides insightful perspectives on the distribution and evolution of F. pseudograminearum, and for forecasting the spread of wheat FCR, thereby aiding in the development of preventive measures and control strategies.

Hidden diversity within Beauveria and Metarhizium – Comparing morphology, barcoding, multilocus phylogenies and whole-genome sequences

Accurately identifying entomopathogenic fungi is crucially important, but the current approach of analysing four genes might not provide sufficient resolution. In this study, we investigated the different resolution provided by multilocus phylogenies and approaches based on whole genome sequence data. Fungi were isolated from soil samples that were collected from five different vegetation types (dry sclerophyll forest, agricultural grassland, rainforest, suburban parkland, and sugar cane fields) and across four different suburban soil habitat types in southeast Queensland. Three different agricultural pests were used as live baits, cotton stainer bug, diamondback moth, and rust-red flour beetle. Whole-genome sequencing was conducted for all 83 isolates recovered, and the ITS2 region was extracted from the genome assemblies to make initial species identifications with the UNITE database. We also extracted tef1a 3′, tef1a 5’, rpb1 and rpb2 genes from the Metarhizium genomes and the bloc, tef1a and rpb1 genes from Beauveria genomes to construct multilocus phylogenies and obtain species identification. To investigate the genetic relationships across 14 isolates of Beauveria bassiana and (independently) across 43 isolates of Metarhizium based on whole-genome data, we genotyped single nucleotide polymorphism (SNP) markers and conducted principal components analysis on the whole-genome SNP data. The multilocus methods identified isolates to species more precisely than ITS2, except in the one unresolved clade in the Metarhizium phylogeny. The whole-genome approach identified more genetic clusters than the multilocus phylogenies identified species among the isolates, and the morphological results correlated with some of the genetic clusters, so they likely represent distinct species not detected by the other methods. The genetic clusters were not associated with vegetation type or bait insect species. This is the first comparison of the resolution of multilocus phylogenetics with that of whole-genome SNP data for these genera. We suggest how the genetic clusters identified here may be investigated further to determine whether they represent unrecognised species within these groups.

A global view on carbapenem-resistant Acinetobacter baumannii.

To give an update on the molecular epidemiology and global distribution of carbapenemase encoding genes, we subjected 313 carbapenem-resistant Acinetobacter baumannii isolated from 114 study centers in 47 countries in five world regions, Africa, Asia, Europe, Latin America, and North America, to whole genome sequencing. Numbers of isolates investigated were proportional to the population size of the contributing countries. Molecular epidemiology was investigated using seven-loci and core genome multilocus sequence typing, whole-genome single nucleotide polymorphism phylogenies, and the intrinsic blaOXA-51-like variant. Carbapenemase encoding genes were identified by multiplex PCR and ResFinder. Among the total of 313 isolates, 289 (92.3%) were assigned to A. baumannii international clones (IC) IC1-IC8. IC2 predominated with 196 isolates (62.6%) and was spread worldwide, followed by IC5 with 44 isolates (14.1%) mainly confined to Latin America. Six isolates (1.9%) originating from Belgium, Egypt, Italy, and Pakistan represent the novel IC9. Acquired OXA-type carbapenemase genes were found in 300 (96%) isolates with blaOXA-23-like and blaOXA-40-like predominating, which constitutes a significant increase compared to our findings from 2010. Metallo-beta-lactamases were rare with seven isolates (2.2%). The distribution of ICs and carbapenemase determinants can vary widely among different geographical regions. IMPORTANCE Carbapenem-resistant Acinetobacter baumannii are of increasing public health importance, as they are resistant to last-line antibiotics. International clones with well-characterized resistance genes dominate globally; however, locally, other lineages with different properties may be of importance to consider. This study investigated isolates from a broad geographic origin from 114 hospitals in 47 countries and from five world regions ensuring the greatest possible diversity in an organism known for its propensity for clonal epidemic spread and reflecting the current global epidemiology of carbapenem-resistant A. baumannii. In Latin America, a lineage different from other geographic regions circulates, with a different resistance gene profile. This knowledge is important to adjust local infection prevention measures. In a global world with migration and increasing use of antimicrobials, multidrug-resistant bacteria will continue to adapt and challenge our healthcare systems worldwide.

Genomic prediction based on selective linkage disequilibrium pruning of low-coverage whole-genome sequence variants in a pure Duroc population

BackgroundAlthough the accumulation of whole-genome sequencing (WGS) data has accelerated the identification of mutations underlying complex traits, its impact on the accuracy of genomic predictions is limited. Reliable genotyping data and pre-selected beneficial loci can be used to improve prediction accuracy. Previously, we reported a low-coverage sequencing genotyping method that yielded 11.3 million highly accurate single-nucleotide polymorphisms (SNPs) in pigs. Here, we introduce a method termed selective linkage disequilibrium pruning (SLDP), which refines the set of SNPs that show a large gain during prediction of complex traits using whole-genome SNP data.ResultsWe used the SLDP method to identify and select markers among millions of SNPs based on genome-wide association study (GWAS) prior information. We evaluated the performance of SLDP with respect to three real traits and six simulated traits with varying genetic architectures using two representative models (genomic best linear unbiased prediction and BayesR) on samples from 3579 Duroc boars. SLDP was determined by testing 180 combinations of two core parameters (GWAS P-value thresholds and linkage disequilibrium r2). The parameters for each trait were optimized in the training population by five fold cross-validation and then tested in the validation population. Similar to previous GWAS prior-based methods, the performance of SLDP was mainly affected by the genetic architecture of the traits analyzed. Specifically, SLDP performed better for traits controlled by major quantitative trait loci (QTL) or a small number of quantitative trait nucleotides (QTN). Compared with two commercial SNP chips, genotyping-by-sequencing data, and an unselected whole-genome SNP panel, the SLDP strategy led to significant improvements in prediction accuracy, which ranged from 0.84 to 3.22% for real traits controlled by major or moderate QTL and from 1.23 to 11.47% for simulated traits controlled by a small number of QTN.ConclusionsThe SLDP marker selection method can be incorporated into mainstream prediction models to yield accuracy improvements for traits with a relatively simple genetic architecture, however, it has no significant advantage for traits not controlled by major QTL. The main factors that affect its performance are the genetic architecture of traits and the reliability of GWAS prior information. Our findings can facilitate the application of WGS-based genomic selection.

A Genome-Wide Association Study Reveals Region Associated with Seed Protein Content in Cowpea.

Cowpea (Vigna unguiculata L. Walp., 2n = 2x = 22) is a protein-rich crop that complements staple cereals for humans and serves as fodder for livestock. It is widely grown in Africa and other developing countries as the primary source of protein in the diet; therefore, it is necessary to identify the protein-related loci to improve cowpea breeding. In the current study, we conducted a genome-wide association study (GWAS) on 161 cowpea accessions (151 USDA germplasm plus 10 Arkansas breeding lines) with a wide range of seed protein contents (21.8~28.9%) with 110,155 high-quality whole-genome single-nucleotide polymorphisms (SNPs) to identify markers associated with protein content, then performed genomic prediction (GP) for future breeding. A total of seven significant SNP markers were identified using five GWAS models (single-marker regression (SMR), the general linear model (GLM), Mixed Linear Model (MLM), Fixed and Random Model Circulating Probability Unification (FarmCPU), and Bayesian-information and Linkage-disequilibrium Iteratively Nested Keyway (BLINK), which are located at the same locus on chromosome 8 for seed protein content. This locus was associated with the gene Vigun08g039200, which was annotated as the protein of the thioredoxin superfamily, playing a critical function for protein content increase and nutritional quality improvement. In this study, a genomic prediction (GP) approach was employed to assess the accuracy of predicting seed protein content in cowpea. The GP was conducted using cross-prediction with five models, namely ridge regression best linear unbiased prediction (rrBLUP), Bayesian ridge regression (BRR), Bayesian A (BA), Bayesian B (BB), and Bayesian least absolute shrinkage and selection operator (BL), applied to seven random whole genome marker sets with different densities (10 k, 5 k, 2 k, 1 k, 500, 200, and 7), as well as significant markers identified through GWAS. The accuracies of the GP varied between 42.9% and 52.1% across the seven SNPs considered, depending on the model used. These findings not only have the potential to expedite the breeding cycle through early prediction of individual performance prior to phenotyping, but also offer practical implications for cowpea breeding programs striving to enhance seed protein content and nutritional quality.

Implementation of routine genomic surveillance provided insights into a locally acquired outbreak caused by a rare clade of Salmonella enterica serovar Enteritidis in Queensland, Australia.

Salmonellosis is a significant public health problem globally. In Australia, Salmonella enterica serovar Enteritidis is one of the main causes of salmonellosis. This study reports how the implementation of routine genetic surveillance of isolates from human S. Enteritidis cases enabled identification of the likely source of an outbreak that occurred in a remote town in Far North Queensland, Australia. This study included patient, food and water samples collected during an outbreak investigation. S. Enteritidis of the novel sequence type 5438 was isolated from all seven patient samples and one bore water sample but not any of the food samples. Both whole-genome single nucleotide polymorphism (SNP) and core-genome multilocus sequence typing analysis revealed that S. Enteritidis isolated from outbreak-related patient samples and the bore water isolates clustered together with fewer than five SNP differences and ten allelic differences. This genetic relatedness informed the outbreak response team around public health interventions and no further cases were identified post-treatment of the bore water. This disease cluster was identified through the routine sequencing of S. Enteritidis performed by the state public health laboratory in an actionable time frame. Additionally, genomic surveillance captured a case with unknown epidemiological links to the affected community, ruled out a simultaneous outbreak in an adjacent state as the source and provided evidence for the likely source preventing further transmission. Therefore, this report provides compelling support for the implementation of whole-genome sequencing based genotyping methods in public health microbiology laboratories for better outbreak detection and management.

A clinically applicable molecular classification of oncocytic cell thyroid nodules.

Whole-chromosome instability with near-whole genome haploidization (GH) and subsequent endoreduplication is considered a main genomic driver in the tumorigenesis of oncocytic cell thyroid neoplasms (OCN). These copy number alterations (CNA) occur less frequently in oncocytic thyroid adenoma (OA) than in oncocytic carcinoma (OCA), suggesting a continuous process. The current study described the CNA patterns in a cohort of 30 benign and malignant OCN, observed using a next generation sequencing (NGS) panel that assesses genome-wide loss of heterozygosity (LOH) and chromosomal imbalances using 1,500 single nucleotide polymorphisms (SNP) across all autosomes and the X chromosome in DNA derived from cytological and histological samples. Observed CNA patterns were verified using multiparameter DNA flow cytometry with or without whole-genome SNP array analysis and lesser-allele intensity-ratio (LAIR) analysis. On CNA-LOH analysis using the NGS panel, GH type CNA were observed 4 of 11 (36%) OA and 14 of 16 OCA (88%). Endoreduplication was suspected in 8 of 16 (50%) OCA, all with more extensive GH type CNA (p<0.001). Reciprocal chromosomal imbalance type CNA, characterized by (imbalanced) chromosomal copy number gains and associated with benign disease, were observed in 6 of 11 (55%) OA and one equivocal case of OCA. CNA patterns were different between the histopathological subgroups (p<0.001). By applying the structured interpretation and considerations provided by the current study, CNA-LOH analysis using an NGS panel that is feasible for daily practice may be of great added value to the widespread application of molecular diagnostics in the diagnosis and risk stratification of OCN.

Genome Investigation of Brucella melitensis Strains Isolated From Brucellosis Patients in Kuwait

Background: Human brucellosis is present on all inhabited continents with high prevalence in many areas of the world, including Kuwait and the Middle East. To implement proper control measures, the identification and characterization of Brucella species and genotypes are required through a reliable and rapid subtyping method. In previous studies, whole-genome sequencing (WGS) has shown its potential as an epidemiological typing tool. Using WGS data, this study aimed to identify the species, phage sequences, putative antibiotic resistance genes, virulence factors, and genotypes of Brucella melitensis strains isolated from patients in Kuwait and other countries. Methods: Five B. melitensis isolates of Kuwaiti origin and 31 other isolates of B. melitensis originating from 28 countries were analyzed using whole genome-based approaches for genotypic identification and typing. In-silico techniques were used to identify the sequences for phages, antibiotic resistance genes, virulence factors, and genotypes using multilocus sequence typing (MLST) and whole-genome single-nucleotide polymorphism (wgSNP). Results: The analysis of WGS data demonstrated that all five Kuwaiti isolates belonged to the non-vaccine strains of B. melitensis. Furthermore, the data represented the presence of two phage sequences, two antibiotic resistance genes, and 51 virulence factors in Kuwaiti isolates. Eventually, the genotypes of all isolates were identified based on MLST and wgSNP analysis, and wgSNP analysis suggested the possible areas/countries of origin of Kuwaiti isolates. Conclusion: WGS data can be used to characterize Brucella isolates, and molecular techniques can be applied in-silico to rapidly identify and classify Brucella into species and genotypes and trace the possible origin of the isolates.

Whole genome resequencing reveal patterns of genetic variation within Colletotrichum acutatum species complex from rubber trees in China

The Colletotrichum acutatum species complex possesses a diverse number of important traits, such as a wide host range and host preference, different modes of reproduction, and different strategies of host infection. Research using comparative genomics has attempted to find correlations between these traits. Here, we used multi-locus techniques and gene genealogical concordance analysis to investigate the phylogenetic relationships and taxonomic status of the Colletotrichum acutatum species complex using field isolates obtained from rubber trees. The results revealed that the dominant species was C. australisinense, followed by C. bannaense, while strain YNJH17109 was identified as C. laticiphilum. The taxonomic status of strains YNLC510 and YNLC511 was undetermined. Using whole-genome single nucleotide polymorphism data to analyze population structure, 18 strains of C. australisinense were subsequently divided into four populations, one of which was derived from an admixture of two populations. In addition, the strains LD1687, GD1628, and YNLC516, did not belong to any populations, and were considered to be admixtures of two or more populations. A split decomposition network analysis also provided evidence for genetic recombination within Colletotrichum acutatum species complex from rubber trees in China. Overall, a weak phylogeographic sub-structure was observed. Analysis also revealed significant differences in morphological characters and levels of virulence between populations.

Association of periodontal disease with migraine: A large‐scale community‐based cross‐sectional study

Periodontal disease may drive a systemic inflammatory response that triggers migraine; however, the association between periodontal disease and migraine has rarely been investigated in a community-based setting. This cross-sectional study included 66,109 participants aged 30 to 70 years from Taiwan Biobank (TWB). A structured questionnaire was administered to participants, who were also subjected to whole-genome single nucleotide polymorphism genotyping using the customized Axiom-TWB array. To identify subjects with periodontal disease and migraine, the computerized linkage of data obtained from TWB and the National Health Insurance Research Database was performed. Participants were evaluated for their genetic predisposition to migraine using a polygenic risk score. We examined and estimated the magnitude of associations between periodontal disease and migraine. In this study, 4618 (4618/66,109; 7%) participants with migraine and 61,491 (61,491/66,109; 83%) participants without migraine were included. Participants with migraine exhibited a higher prevalence of periodontal disease than participants without migraine (4324/4618; 94% vs. 56,036/61,491; 91%). A significant positive association was observed between periodontal disease and migraine, with an adjusted odds ratio (ORadj ) of 1.40 (95% confidence interval [CI]=1.24-1.59; p <0.001). The association remained consistent even after excluding participants with other comorbidities (ORadj =1.34; 95% CI=1.16-1.55; p <0.001). Moreover, the positive association between periodontal disease and migraine remained significant across the subgroups of age, sex, other comorbidities, and classified polygenic risk scores of migraine, with the ORadj ranging from 1.26 to 1.78. A significant positive association was observed between periodontal disease and migraine. Future studies need to explore the biological mechanisms of how periodontal disease might affect migraine.

Brucella abortus in Kazakhstan, population structure and comparison with worldwide genetic diversity.

Brucella abortus is the main causative agent of brucellosis in cattle, leading to severe economic consequences in agriculture and affecting public health. The zoonotic nature of the infection increases the need to control the spread and dynamics of outbreaks in animals with the incorporation of high resolution genotyping techniques. Based on such methods, B. abortus is currently divided into three clades, A, B, and C. The latter includes subclades C1 and C2. This study presents the results of whole-genome sequencing of 49 B. abortus strains isolated in Kazakhstan between 1947 and 2015 and of 36 B. abortus strains of various geographic origins isolated from 1940 to 2004. In silico Multiple Locus Sequence Typing (MLST) allowed to assign strains from Kazakhstan to subclades C1 and to a much lower extend C2. Whole-genome Single-Nucleotide Polymorphism (wgSNP) analysis of the 46 strains of subclade C1 with strains of worldwide origins showed clustering with strains from neighboring countries, mostly North Caucasia, Western Russia, but also Siberia, China, and Mongolia. One of the three Kazakhstan strains assigned to subclade C2 matched the B. abortus S19 vaccine strain used in cattle, the other two were genetically close to the 104 M vaccine strain. Bayesian phylodynamic analysis dated the introduction of B. abortus subclade C1 into Kazakhstan to the 19th and early 20th centuries. We discuss this observation in view of the history of population migrations from Russia to the Kazakhstan steppes.

EGFR and MMP-9 are associated with neointimal hyperplasia in systemic-to-pulmonary shunts in children with complex cyanotic heart disease

Systemic-to-pulmonary shunt malfunction contributes to morbidity in children with complex congenital heart disease after palliative procedure. Neointimal hyperplasia might play a role in the pathogenesis increasing risk for shunt obstruction. The aim was to evaluate the role of epidermal growth factor receptor (EGFR) and matrix-metalloproteinase 9 (MMP-9) in the formation of neointimal within shunts. Immunohistochemistry was performed with anti-EGFR and anti-MMP-9 on shunts removed at follow-up palliative or corrective procedure. Whole-genome single-nucleotide polymorphisms genotyping was performed on DNA extracted from patients´ blood samples and allele frequencies were compared between the group of patients with shunts displaying severe stenosis (≥ 40% of lumen) and the remaining group. Immunohistochemistry detected EGFR and MMP-9 in 24 of 31 shunts, located mainly in the luminal area. Cross-sectional area of EGFR and MMP-9 measured in median 0.19 mm2 (IQR 0.1–0.3 mm2) and 0.04 mm2 (IQR 0.03–0.09 mm2), respectively, and correlated positively with the area of neointimal measured on histology (r = 0.729, p < 0.001 and r = 0.0479, p = 0.018, respectively). There was a trend of inverse correlation between the dose of acetylsalicylic acid and the degree of EGFR, but not MMP-9, expression within neointima. Certain alleles in epidermal growth factor (EGF) and tissue inhibitor of metalloproteinases 1 (TIMP-1) were associated with increased stenosis and neointimal hyperplasia within shunts. EGFR and MMP-9 contribute to neointimal proliferation in SP shunts of children with complex cyanotic heart disease. SP shunts from patients carrying certain risk alleles in the genes encoding for EGF and TIMP-1 displayed increased neointima.

Whole-genome sequencing reveals genomic characterization of Listeria monocytogenes from food in China.

Listeria monocytogenes is a foodborne bacterium that could persist in food and food processing environments for a long time. Understanding the population structure and genomic characterization of foodborne L. monocytogenes is essential for the prevention and control of listeriosis. A total of 322 foodborne L. monocytogenes isolates from 13 geographical locations and four food sources in China between 2000 and 2018 were selected for whole-genome sequencing. In silico subtyping divided the 322 isolates into five serogroups, 35 sequence types (STs), 26 clonal complexes (CCs) and four lineages. Serogroup IIa was the most prevalent serogroup and ST9 was the most prevalent ST of foodborne L. monocytogenes strains isolated in China. The in-depth phylogenetic analysis on CC9 revealed that ST122 clone might be original from ST9 clone. Furthermore, 23 potentially relevant clusters were identified by pair-wised whole-genome single nucleotide polymorphism analysis, indicating that persistent- and/or cross-contamination had occurred in markets in China. ST8 and ST121 were the second and third top STs of L. monocytogenes in China, which had heterogeneity with that of L. monocytogenes isolates from other countries. The antibiotic resistance genes aacA4, tetM, tetS, dfrG carried by different mobile elements were found in L. monocytogenes strains. One lineage II strain carrying Listeria Pathogenicity Island 3 was first reported. In addition, a novel type of premature stop codon in inlA gene was identified in this study. These findings revealed the genomic characteristics and evolutionary relationship of foodborne L. monocytogenes in China on a scale larger than previous studies, which further confirmed that whole-genome sequencing analysis would be a helpful tool for routine surveillance and source-tracing investigation.

Comparative whole-genome analysis of diphtheria toxin-producing Corynebacterium ulcerans isolates from sheltered cats in Osaka, Japan.

Diphtheria toxin-producing Corynebacterium ulcerans is a zoonotic pathogen that causes human diphtheria-like symptoms. After performing whole-genome analysis of the five isolates from sheltered cats in Osaka, Japan, we compared them with genome sequences of 25 strains of C. ulcerans from a public database. The five isolates from cats harbored 14 genes encoding possible virulence factors in diphtheria-toxin-producing C. ulcerans. These isolates also had diphtheria toxin gene-encoding prophage in their chromosome, although differences were found in other prophages possession. Whole-genome single-nucleotide polymorphism analysis showed that cats' isolates belonged to ST337 branch, as were strains from Japanese human patients, with 41 or more single-nucleotide polymorphisms variations. High-resolution single-nucleotide polymorphism analysis of C. ulcerans was sufficient to distinguish cats' isolates clearly as not different by conventional genotyping methods.

Genetic Diversity and Population Structure of Fusarium oxysporum f. sp. conglutinans Race 1 in Northern China Samples.

Fusarium oxysporum f. sp. conglutinans (FOC), the causal agent of cabbage fusarium wilt, is a serious threat to cabbage production in northern China, and most Chinese FOC isolates were identified as FOC race 1 (FOC1). To better understand the genetic diversity of FOC1 in northern China, we collected FOC isolates from five provinces in northern China and identified them as FOC1 through pathogenicity and race test. To evaluate the genome-level diversity of FOC1, we performed a genome assembly for a FOC1 isolate (FoYQ-1) collected from Yanqing, Beijing, where cabbage fusarium wilt was first reported in China. Using resequencing data of FOC1 isolates, we conducted a genome-wide SNP (single nucleotide polymorphism) analysis to investigate the genetic diversity and population structure of FOC1 isolates in northern China. Our study indicated that Chinese FOC1 can be grouped into four populations and revealed that the genetic diversity of FOC1 were closely associated with geographical locations. Our study further suggests that genetic differentiation occurred when FOC1 spread to the northwest provinces from Beijing Province in China. The FOC1 genetic diversity based on whole-genome SNPs could deepen our understanding of FOC1 variation and provide clues for the control of cabbage fusarium wilt in China.

Comparative Analysis of Molecular RFLP and SNP Markers in Assessing and Understanding the Genetic Diversity of Various Chicken Breeds.

Monitoring the genetic diversity of small populations is important with respect to conserving rare and valuable chicken breeds, as well as discovery and innovation in germplasm research and application. Restriction fragment length polymorphisms (RFLPs), the molecular markers that underlie multilocus DNA fingerprinting (MLDF), have historically been employed for this purpose, but over the past two decades, there has been an irreversible shift toward high-throughput single-nucleotide polymorphisms (SNPs). In this study, we conducted a comparative analysis of archived MLDF results and new data from whole-genome SNP genotyping (SNPg) among 18 divergently selected breeds representing a large sample of the world gene pool. As a result, we obtained data that fit the general concept of the phylogenetic distribution of the studied breeds and compared them with RFLP and SNP markers. RFLPs were found to be useful markers for retrospective assessment of changes in the genetic architecture and variability underlying the phenotypic variation in chicken populations, especially when samples from previous generations used for MLDF are unavailable for SNPg. These results can facilitate further research necessary to assess the possibility of extrapolating previous MLDF results to study the long-term dynamics of genetic diversity in various small chicken germplasm populations over time. In general, the whole-genome characterization of populations and breeds by multiple SNP loci will further form the basis for the development and implementation of genomic selection with the aim of effective use of the genetic potential of the domestic gene pool in the poultry industry.