Whole Genome Sequence Information Research Articles

A Complete Genome Sequence of Podosphaera xanthii Isolate YZU573, the Causal Agent of Powdery Mildew Isolated from Cucumber in China.

Podosphaera xanthii is a well-known obligate biotrophic pathogen that causes powdery mildew (PM) disease on cucurbitaceous plants and is one of the most important limiting factors for cucumber production worldwide. To better understand the avirulence effector proteins in this species that are known to be involved in host-pathogen interaction, the draft genome assembly of P. xanthii isolate YZU573 from cucumber leaves with symptoms of PM was obtained with a hybrid approach, combining nanopore long-read and llumina paired-end sequencing. The final P. xanthii YZU573 genome assembly of 152.7 Mb consists of 58 contigs, with an N50 value of 0.75 Mb and 6491 predicted protein-coding genes. The effector analysis using the whole-genome sequence information revealed a total of 87 putative effector candidates, and 65 of them had their analogs, whereas the remaining 22 were novel ones. The new P. xanthii genome provides valuable resources to better understand plant-microbe interaction in cucumber PM disease.

Whole genome analysis for plant growth promotion profiling of Pantoea agglomerans CPHN2, a non-rhizobial nodule endophyte.

Reduced agricultural production as well as issues like nutrient-depleted soils, eutrophication, and groundwater contamination have drawn attention to the use of endophyte-based bioformulations to restore soil fertility. Pantoea agglomerans CPHN2, a non-rhizobial nodule endophyte isolated from Cicer arietinum, exhibited a variety of plant growth-promoting traits. In this study, we used NextSeq500 technology to analyze whole-genome sequence information of this plant growth-promoting endophytic bacteria. The genome of P. agglomerans CPHN2 has a length of 4,839,532 bp and a G + C content of 55.2%. The whole genome comprises three different genomic fractions, comprising one circular chromosome and two circular plasmids. A comparative analysis between P. agglomerans CPHN2 and 10 genetically similar strains was performed using a bacterial pan-genome pipeline. All the predicted and annotated gene sequences for plant growth promotions (PGPs), such as phosphate solubilization, siderophore synthesis, nitrogen metabolism, and indole-3-acetic acid (IAA) of P. agglomerans CPHN2, were identified. The whole-genome analysis of P. agglomerans CPHN2 provides an insight into the mechanisms underlying PGP by endophytes and its potential applications as a biofertilizer.

Reclassification of Streptococcus ilei as a later heterotypic synonym of Streptococcus koreensis based on whole-genome sequence analysis.

The genus Streptococcus, a member of family Streptococcaceae, is known for its wide range of industrial, clinical and human relevance. Among the species of genus Streptococcus two members, namely Streptococcus koreensis and Streptococcus ilei, were isolated from subgingival dental plaque and human small intestinal fluid, respectively. The 16S rRNA gene sequence similarity of the type strains of these members shows a similarity of 99.87%. In this study, we performed a systematic study to clarify the taxonomic assignment of these two species. Genome similarity assessment based on whole-genome sequence information such as average nucleotide identity using orthoANI and fastANI, digital DNA-DNA hybridization value between S. koreensis and S. ilei were 96.31, 96.60, 86.4 and 97.63, respectively. All these genome similarity values clearly exceeded the species delineation cutoffs. Phylogenetic assessment using 16S rRNA gene and whole-genome information using PhyloPhlAn, which uses around 400 conserved genes across bacterial phyla, provides additional evidence for these members forming a monophyletic clade in the phylogenetic tree. Pan genome analysis suggests a very large core genome (n = 1374) and the presence of no unique gene between the genomes of S. koreensis and S. ilei. Additionally, we found highly syntenic genomes of type strains of these two species. Based on these evidences, we propose S. ilei should be reclassified as a later heterotypic synonym of S. koreensis.

Whole-Genome Sequence and Comparative Analysis of Trichoderma asperellum ND-1 Reveal Its Unique Enzymatic System for Efficient Biomass Degradation

The lignocellulosic enzymes of Trichoderma asperellum have been intensely investigated toward efficient conversion of biomass into high-value chemicals/industrial products. However, lack of genome data is a remarkable hurdle for hydrolase systems studies. The secretory enzymes of newly isolated T. asperellum ND-1 during lignocellulose degradation are currently poorly known. Herein, a high-quality genomic sequence of ND-1, obtained by both Illumina HiSeq 2000 sequencing platforms and PacBio single-molecule real-time, has an assembly size of 35.75 Mb comprising 10,541 predicted genes. Secretome analysis showed that 895 proteins were detected, with 211 proteins associated with carbohydrate-active enzymes (CAZymes) responsible for biomass hydrolysis. Additionally, T. asperellum ND-1, T. atroviride IMI 206040, and T. virens Gv-298 shared 801 orthologues that were not identified in T. reesei QM6a, indicating that ND-1 may play critical roles in biological-control. In-depth analysis suggested that, compared with QM6a, the genome of ND-1 encoded a unique enzymatic system, especially hemicellulases and chitinases. Moreover, after comparative analysis of lignocellulase activities of ND-1 and other fungi, we found that ND-1 displayed higher hemicellulases (particularly xylanases) and comparable cellulases activities. Our analysis, combined with the whole-genome sequence information, offers a platform for designing advanced T. asperellum ND-1 strains for industrial utilizations, such as bioenergy production.

Use of gene expression and whole-genome sequence information to improve the accuracy of genomic prediction for carcass traits in Hanwoo cattle

BackgroundIn this study, we assessed the accuracy of genomic prediction for carcass weight (CWT), marbling score (MS), eye muscle area (EMA) and back fat thickness (BFT) in Hanwoo cattle when using genomic best linear unbiased prediction (GBLUP), weighted GBLUP (wGBLUP), and a BayesR model. For these models, we investigated the potential gain from using pre-selected single nucleotide polymorphisms (SNPs) from a genome-wide association study (GWAS) on imputed sequence data and from gene expression information. We used data on 13,717 animals with carcass phenotypes and imputed sequence genotypes that were split in an independent GWAS discovery set of varying size and a remaining set for validation of prediction. Expression data were used from a Hanwoo gene expression experiment based on 45 animals.ResultsUsing a larger number of animals in the reference set increased the accuracy of genomic prediction whereas a larger independent GWAS discovery dataset improved identification of predictive SNPs. Using pre-selected SNPs from GWAS in GBLUP improved accuracy of prediction by 0.02 for EMA and up to 0.05 for BFT, CWT, and MS, compared to a 50 k standard SNP array that gave accuracies of 0.50, 0.47, 0.58, and 0.47, respectively. Accuracy of prediction of BFT and CWT increased when BayesR was applied with the 50 k SNP array (0.02 and 0.03, respectively) and was further improved by combining the 50 k array with the top-SNPs (0.06 and 0.04, respectively). By contrast, using BayesR resulted in limited improvement for EMA and MS. wGBLUP did not improve accuracy but increased prediction bias. Based on the RNA-seq experiment, we identified informative expression quantitative trait loci, which, when used in GBLUP, improved the accuracy of prediction slightly, i.e. between 0.01 and 0.02. SNPs that were located in genes, the expression of which was associated with differences in trait phenotype, did not contribute to a higher prediction accuracy.ConclusionsOur results show that, in Hanwoo beef cattle, when SNPs are pre-selected from GWAS on imputed sequence data, the accuracy of prediction improves only slightly whereas the contribution of SNPs that are selected based on gene expression is not significant. The benefit of statistical models to prioritize selected SNPs for estimating genomic breeding values is trait-specific and depends on the genetic architecture of each trait.

Antibiotic susceptibility of plant-derived lactic acid bacteria conferring health benefits to human.

Lactic acid bacteria (LAB) confer health benefits to human when administered orally. We have recently isolated several species of LAB strains fromplant sources, such as fruits, vegetables, flowers, and medicinal plants. Since antibiotics used to treat bacterial infection diseases induce the emergence of drug-resistant bacteria in intestinal microflora, it is important to evaluate the susceptibility of LAB strains to antibiotics to ensure the safety and security of processed foods. The aim of the present study is to determine the minimum inhibitory concentration (MIC) of antibiotics against several plant-derived LAB strains. When aminoglycoside antibiotics, such as streptomycin (SM), kanamycin (KM), and gentamicin (GM), were evaluated using LAB susceptibility test medium (LSM), the MIC was higher than when using Mueller-Hinton (MH) medium. Etest, which is an antibiotic susceptibility assay method consisting of a predefined gradient of antibiotic concentrations on a plastic strip, is used to determine the MIC of antibiotics world-wide. In the present study, we demonstrated that Etest was particularly valuable while testing LAB strains. We also show that the low susceptibility of the plant-derived LAB strains against each antibiotic tested is due to intrinsic resistance and not acquired resistance. This finding is based on the whole-genome sequence information reflecting the horizontal spread of the drug-resistance genes in the LAB strains.

Evaluation of the RAPIDEC CARBA NP, Conventional CarbaNP, and the Modified Carbapenem Inactivation Method (mCIM) Tests for Phenotypic Detection of Carbapenemase-Producing Organisms

Decreased susceptibility to carbapenems is being increasingly reported among Gram-negative organisms worldwide. It is of critical importance for patient management and infection control purposes to detect carbapenemase-producing organisms rapidly and reliably. In this study, the performance of the commercially available RAPIDEC CARBA NP test (bioMérieux, Durham, NC) was evaluated against the performance of the conventional CarbaNP and the modified Carbapenem Inactivation Method (mCIM) for the detection of carbapenemase activity in Enterobacteriaceae and Pseudomonas aeruginosa. A total of 73 Enterobacteriaceae and 32 Pseudomonas aeruginosa isolates were tested for carbapenemase activity using the RAPIDEC CARBA NP, conventional CarbaNP, and mCIM methods. Using available whole-genome sequence information, this collection included 72 isolates with and 33 isolates without carbapenemase genes. The isolates harboring a carbapenemase gene encoded 20 class A, 47 class B, and 5 class D–type β-lactamases. All methods were performed using bacteria from a single overnight blood agar plate. The conventional CarbaNP and mCIM were performed according to guidelines from the Clinical and Laboratory Standards Institute. The RAPIDEC CARBA NP was performed per the manufacturer’s instructions. The RAPIDEC CARBA NP detected 65 of 72 (90%) carbapenemase-producing isolates; the seven isolates yielding false-negative results included five Enterobacteriaceae with blaOXA-48-like and two with blaKPC-4 genes. The sensitivity and specificity were 90.2% and 100%, respectively, for the RAPIDEC test; 94.4% and 100% for conventional CarbaNP; and 100% and 81.8% for mCIM. The conventional CarbaNP detected all carbapenemase producers except four OXA-48-like–producing isolates. The mCIM, which accurately detected all carbapenemase producers, yielded false-positive results in six isolates of P aeruginosa. Overall, the RAPIDEC CARBA NP was easy to use and required less hands-on time than the conventional CarbaNP and mCIM methods. However, the RAPIDEC CARBA NP showed decreased sensitivity for detecting low-level activity carbapenemases such as OXA-48-like enzymes and KPC-4.

Power and precision of QTL mapping in simulated multiple porcine F2 crosses using whole-genome sequence information

BackgroundDuring the last two decades, many QTL (quantitative trait locus) mapping experiments in pigs have been conducted using F2 crosses established from two outbred founder breeds. The founder breeds were frequently chosen from the Asian and European type breeds. A combination of next-generation sequencing, SNP (single nucleotide polymorphism) genotyping technology using SNP-chips, and genotype imputation techniques, can be used to infer the sequence information of all F2 individuals in a cost-effective way. The aim of the present simulation study was to analyze the power and precision of genome-wide association studies (GWASs) with whole-genome sequence data in several types of F2 crosses, including pooled crosses.MethodsBased on a common historical population, three breeds representing two European type breeds (EU1 and EU2) and one Asian type breed (AS) were simulated. Two F2 designs of 500 individuals each were simulated. The cross EU1xEU2 (ASxEU2) was simulated using the phylogenetically closely related breeds EU1 and EU2 (or distantly related breeds AS and EU2) as the founder breeds. The simulated genomes comprised ten chromosomes, each with a length of 1 Morgan and whole-genome sequence information. A polygenic trait with a heritability of 0.5, which was affected by approximately 20 QTL per Morgan, was simulated. GWASs were conducted using single marker mixed linear models, either within the crosses or in their pooled datasets. Additionally, the studies were conducted in the breed EU2, which was a founder breed in both simulated crosses.ResultsThe power to map QTL was high (low) in the ASxEU2 (EU1xEU2) cross and was highest when the data of both crosses were analyzed jointly. By contrast, the mapping precision was the highest in the EU1xEU2 cross. Pooling data led to a precision that was in between the precision of the EU1xEU2 cross and the ASxEU2 cross. A higher mapping precision was observed for QTL segregating within a founder breed.ConclusionsThese results suggest that the existing F2 crosses are promising databases for QTL mapping when the founder breeds are closely related or several crosses can be pooled. This conclusion is particularly applicable for QTL that segregate in a founder breed.

Gene Flow between Divergent Cereal- and Grass-Specific Lineages of the Rice Blast Fungus Magnaporthe oryzae.

ABSTRACTDelineating species and epidemic lineages in fungal plant pathogens is critical to our understanding of disease emergence and the structure of fungal biodiversity and also informs international regulatory decisions. Pyricularia oryzae (syn. Magnaporthe oryzae) is a multihost pathogen that infects multiple grasses and cereals, is responsible for the most damaging rice disease (rice blast), and is of growing concern due to the recent introduction of wheat blast to Bangladesh from South America. However, the genetic structure and evolutionary history of M. oryzae, including the possible existence of cryptic phylogenetic species, remain poorly defined. Here, we use whole-genome sequence information for 76 M. oryzae isolates sampled from 12 grass and cereal genera to infer the population structure of M. oryzae and to reassess the species status of wheat-infecting populations of the fungus. Species recognition based on genealogical concordance, using published data or extracting previously used loci from genome assemblies, failed to confirm a prior assignment of wheat blast isolates to a new species (Pyricularia graminis-tritici). Inference of population subdivisions revealed multiple divergent lineages within M. oryzae, each preferentially associated with one host genus, suggesting incipient speciation following host shift or host range expansion. Analyses of gene flow, taking into account the possibility of incomplete lineage sorting, revealed that genetic exchanges have contributed to the makeup of multiple lineages within M. oryzae. These findings provide greater understanding of the ecoevolutionary factors that underlie the diversification of M. oryzae and highlight the practicality of genomic data for epidemiological surveillance in this important multihost pathogen.

Molecular mechanisms of limb regeneration: insights from regenerating legs of the cricket Gryllus bimaculatus.

This review summarizes recent advances in leg regeneration research, focusing on the cricket Gryllus bimaculatus. Recent studies have revealed molecular mechanisms on blastema formation, establishment of positional information, and epigenetic regulation during leg regeneration. Especially, these studies have provided molecular bases in classical conceptual models such as the polar coordinate model, the intercalation model, the boundary model, the steepness model, etc., which were proposed to interpret regeneration processes of the cockroach legs. When a leg is amputated, a blastema is formed through the activation of the Janus-kinase (Jak)/Signal-Transduction-and-Activator-of-Transcription (STAT) pathway. Subsequently, the Hedgehog/Wingless/Decapentaplegic/Epidermal-growth-factor pathways instruct distalization in the blastema, designated as the molecular boundary model. Downstream targets of this pathway are transcription factors Distal-less (Dll) and dachshund (dac), functioning as key regulators of proximodistal pattern formation. Dll and dac specify the distal and proximal regions in the blastema, respectively, through the regulation of tarsal patterning genes. The expression of leg patterning genes during regeneration may be epigenetically controlled by histone H3K27 methylation via Enhancer-of-zeste and Ubiquitously-transcribed-tetratricopeptide-repeat-gene-X-chromosome. For the molecular mechanism of intercalation of the missing structures between the amputated position and the most distal one, Dachsous/Fat (Ds/Ft) steepness model has been proposed, in which the Ds/Ft pathway maintains positional information and determines leg size through dac expression. This model was theoretically verified to interpret the experimental results obtained with cricket legs. Availability of whole-genome sequence information, regeneration-dependent RNA interference, and genome editing technique will have the cricket be an ideal model system to reveal gene functions in leg regeneration.

MisMM: An Integrated Pipeline for Misassembly Detection Using Genotyping-by-Sequencing and Its Validation with BAC End Library Sequences and Gene Synteny.

As next-generation sequencing technologies have advanced, enormous amounts of whole-genome sequence information in various species have been released. However, it is still difficult to assemble the whole genome precisely, due to inherent limitations of short-read sequencing technologies. In particular, the complexities of plants are incomparable to those of microorganisms or animals because of whole-genome duplications, repeat insertions, and Numt insertions, etc. In this study, we describe a new method for detecting misassembly sequence regions of Brassica rapa with genotyping-by-sequencing, followed by MadMapper clustering. The misassembly candidate regions were cross-checked with BAC clone paired-ends library sequences that have been mapped to the reference genome. The results were further verified with gene synteny relations between Brassica rapa and Arabidopsis thaliana. We conclude that this method will help detect misassembly regions and be applicable to incompletely assembled reference genomes from a variety of species.

Isolation of Vaccine-Like Poliovirus Strains in Sewage Samples From the United Kingdom.

Environmental surveillance (ES) is a sensitive method for detecting human enterovirus (HEV) circulation, and it is used worldwide to support global polio eradication. We describe a novel ES approach using next-generation sequencing (NGS) to identify HEVs in sewage samples collected in London, United Kingdom, from June 2016 to May 2017. Two different methods were used to process raw sewage specimens: a 2-phase aqueous separation system and size exclusion by filtration and centrifugation. HEVs were isolated using cell cultures and analyzed using NGS. Type 1 and 3 vaccine-like poliovirus (PV) strains were detected in samples collected from September 2016 through January 2017. NGS analysis allowed us to rapidly obtain whole-genome sequences of PV and non-PV HEV strains. As many as 6 virus strains from different HEV serotypes were identified in a single cell culture flask. PV isolates contained only a small number of mutations from vaccine strains commonly seen in early isolates from vaccinees. Our ES setup has high sensitivity for polio and non-PV HEV detection, generating nearly whole-genome sequence information. Such ES systems provide critical information to assist the polio eradication endgame and contribute to the improvement of our understanding of HEV circulation patterns in humans.

Highly accurate sequence imputation enables precise QTL mapping in Brown Swiss cattle

BackgroundWithin the last few years a large amount of genomic information has become available in cattle. Densities of genomic information vary from a few thousand variants up to whole genome sequence information. In order to combine genomic information from different sources and infer genotypes for a common set of variants, genotype imputation is required.ResultsIn this study we evaluated the accuracy of imputation from high density chips to whole genome sequence data in Brown Swiss cattle. Using four popular imputation programs (Beagle, FImpute, Impute2, Minimac) and various compositions of reference panels, the accuracy of the imputed sequence variant genotypes was high and differences between the programs and scenarios were small. We imputed sequence variant genotypes for more than 1600 Brown Swiss bulls and performed genome-wide association studies for milk fat percentage at two stages of lactation. We found one and three quantitative trait loci for early and late lactation fat content, respectively. Known causal variants that were imputed from the sequenced reference panel were among the most significantly associated variants of the genome-wide association study.ConclusionsOur study demonstrates that whole-genome sequence information can be imputed at high accuracy in cattle populations. Using imputed sequence variant genotypes in genome-wide association studies may facilitate causal variant detection.

Retrotransposon-based insertion polymorphism markers in mango

Retrotransposons are major components of eukaryotic genomes and are present in high copy numbers. We developed retrotransposon-based insertion polymorphism (RBIP) markers based on long terminal repeat (LTR) sequences and flanking genome regions by using shotgun genome sequence data of mango (Mangifera indica L.). Three novel LTR sequences were identified based on two LTR retrotransposon structural features; a 5′ LTR located upstream of the primer binding site and a 3′ LTR showing high sequence similarity to the 5′ LTR. Starting with 377 unique sequences containing both 3′ LTR and downstream genome region sequences, we developed 82 RBIP markers that were applied to DNA fingerprinting of 16 mango accession. Five RBIP markers were enough to distinguish all 16 accessions. Our result showed that LTR identification from shotgun genome sequences was effective for development of retrotransposon-based DNA markers without whole-genome sequence information. We discuss application of the developed RBIP markers for identification of genetic diversity and construction of a genetic linkage map.

Using Whole-Genome Sequence Information to Foster Conservation Efforts for the European Dark Honey Bee, Apis mellifera mellifera

Pollination is a key ecosystem service for agricultural systems and Western honey bees, Apis mellifera, are the most important managed pollinators. Major losses of managed honey bee colonies reinforced the need to take advantage of locally adapted subspecies and ecotypes to buffer populations against various stressors. However, introductions of non-native honey bees from distant lineages are likely to undermine respective conservation efforts unless reliable and cost effective tools can be used to identify hybridization. The purpose of this study is to characterize current population structure and genetic diversity, and to assess the degree of admixture between native and introduced honey bees. Moreover, we aim to select a reduced number of genetic markers to improve conservation management strategies. We take advantage of recent developments in next-generation sequencing and network-based clustering to investigate conservation efforts for the native European Dark honey bee, A. m. mellifera, which is threatened by introgression in most of its range. We collected whole-genome sequence information from haploid drones of A. m. mellifera, A. m. carnica and Buckfast sampled throughout Switzerland (N=81), as well as from four Swiss A. m. mellifera conservation areas (N=39) and from one conservatory in the French Alps (N = 31). Population structure analyses based upon 3.375 M genome-wide SNPs discerned samples by subspecies and geographic origin (Switzerland or France). Ancestry inference indicated admixed individuals in all of the protected areas, calling for improved management efforts. After testing different subsets of ancestry informative SNPs using three different selection strategies (FST, PCA-based or at random), as few as 50 SNPs are found to be sufficient to differentiate native from introduced honey bees. Therefore our data suggests that a low-density SNP panel can be a precise and cost-effective tool to support conservation management efforts for managed pollinators.

Xylella taiwanensis sp. nov., causing pear leaf scorch disease.

A Gram-stain-negative, nutritionally fastidious bacterium (PLS229T) causing pear leaf scorch was identified in Taiwan and previously grouped into Xylella fastidiosa. Yet, significant variations between PLS229T and Xylellafastidiosa were noted. In this study, PLS229T was evaluated phenotypically and genotypically against representative strains of Xylellafastidiosa, including strains of the currently known subspecies of Xylellafastidiosa, Xylella fastidiosa subsp. multiplex and 'Xylella fastidiosasubsp.pauca'. Because of the difficulty of in vitro culture characterization, emphases were made to utilize the available whole-genome sequence information. The average nucleotide identity (ANI) values, an alternative for DNA-DNA hybridization relatedness, between PLS229T and Xylellafastidiosa were 83.4-83.9 %, significantly lower than the bacterial species threshold of 95 %. In contrast, sequence similarity of 16S rRNA genes was greater than 98 %, higher than the 97 % threshold to justify if two bacterial strains belong to different species. The uniqueness of PLS229T was also evident by observing only about 87 % similarity in the sequence of the 16S-23S internal transcribed spacer (ITS) between PLS229T and strains of Xylellafastidiosa, discovering significant single nucleotide polymorphisms at 18 randomly selected housekeeping gene loci, observing a distinct fatty acid profile for PLS229T compared with Xylellafastidiosa, and PLS229T having different observable phenotypes, such as different susceptibility to antibiotics. A phylogenetic tree derived from 16S rRNA gene sequences showed a distinct PLS229T phyletic lineage positioning it between Xylellafastidiosa and members of the genus Xanthomonas. On the basis of these data, a novel species, Xylella taiwanensis sp. nov. is proposed. The type strain is PLS229T (=BCRC 80915T=JCM 31187T).

Virulence Gene Profiles and Clonal Relationships of Escherichia coli O26:H11 Isolates from Feedlot Cattle as Determined by Whole-Genome Sequencing

Escherichia coli O26 is the second most important enterohemorrhagic E. coli (EHEC) serogroup worldwide. Serogroup O26 strains are categorized mainly into two groups: enteropathogenic (EPEC) O26, carrying a locus of enterocyte effacement (LEE) and mostly causing mild diarrhea, and Shiga-toxigenic (STEC) O26, which carries the Shiga toxin (STX) gene (stx), responsible for more severe outcomes. stx-negative O26 strains can be further split into two groups. One O26 group differs significantly from O26 EHEC, while the other O26 EHEC-like group shows all the characteristics of EHEC O26 except production of STX. In order to determine the different populations of O26 E. coli present in U.S. cattle, we sequenced 42 O26:H11 strains isolated from feedlot cattle and compared them to 37 O26:H11 genomes available in GenBank. Phylogenetic analysis by whole-genome multilocus sequence typing (wgMLST) showed that O26:H11/H(-) strains in U.S. cattle were highly diverse. Most strains were sequence type 29 (ST29). By wgMLST, two clear lineages could be distinguished among cattle strains. Lineage 1 consisted of O26:H11 EHEC-like strains (ST29) (4 strains) and O26:H11 EHEC strains (ST21) (2 strains), and lineage 2 (36 strains) consisted of O26:H11 EPEC strains (ST29). Overall, our analysis showed U.S. cattle carried pathogenic (ST21; stx1 (+) ehxA(+) toxB(+)) and also potentially pathogenic (ST29; ehxA(+) toxB(+)) O26:H11 E. coli strains. Furthermore, in silico analysis showed that 70% of the cattle strains carried at least one antimicrobial resistance gene. Our results showed that whole-genome sequence analysis is a robust and valid approach to identify and genetically characterize E. coli O26:H11, which is of importance for food safety and public health. Escherichia coli O26 is the second most important type of enterohemorrhagic E. coli (EHEC) worldwide. Serogroup O26 strains are categorized into two groups: enteropathogenic (EPEC) carrying LEE, causing mild diarrhea, and Shiga toxigenic (STEC) carrying the stx gene, responsible for more severe outcomes. However, there are currently problems in distinguishing one group from the other. Furthermore, several O26 stx-negative strains are consistently misidentified as either EHEC-like or EPEC. The use of whole-genome sequence (WGS) analysis of O26 strains from cattle in the United States (i) allowed identification of O26 strains present in U.S. cattle, (ii) determined O26 strain diversity, (iii) solved the misidentification problem, and (iv) screened for the presence of antimicrobial resistance and virulence genes in the strains. This study provided a framework showing how to easily and rapidly use WGS information to identify and genetically characterize E. coli O26:H11, which is important for food safety and public health.

A Burden of Rare Variants Associated with Extremes of Gene Expression in Human Peripheral Blood

In order to evaluate whether rare regulatory variants in the vicinity of promoters are likely to impact gene expression, we conducted a novel burden test for enrichment of rare variants at the extremes of expression. After sequencing 2-kb promoter regions of 472 genes in 410 healthy adults, we performed a quadratic regression of rare variant count on bins of peripheral blood transcript abundance from microarrays, summing over ranks of all genes. After adjusting for common eQTLs and the major axes of gene expression covariance, a highly significant excess of variants with minor allele frequency less than 0.05 at both high and low extremes across individuals was observed. Further enrichment was seen in sites annotated as potentially regulatory by RegulomeDB, but a deficit of effects was associated with known metabolic disease genes. The main result replicates in an independent sample of 75 individuals with RNA-seq and whole-genome sequence information. Three of four predicted large-effect sites were validated by CRISPR/Cas9 knockdown in K562 cells, but simulations indicate that effect sizes need not be unusually large to produce the observed burden. Unusually divergent low-frequency promoter haplotypes were observed at 31 loci, at least 9 of which appear to be derived from Neandertal admixture, but these were not associated with divergent gene expression in blood. The overall burden test results are consistent with rare and private regulatory variants driving high or low transcription at specific loci, potentially contributing to disease.

Fine genetic mapping of the white immature fruit color gene w to a 33.0-kb region in cucumber (Cucumis sativus L.).

The white immature fruit color gene w was rapidly mapped to a 33.0-kb region to identify a valuable candidate gene that encodes peroxidase. The skin color of immature fruit is a crucial external trait of cucumbers, and white skin is shared by limited numbers of commercial cultivars. Herein, one BC1 population and two F2 segregating populations were constructed using four inbred parental lines (WD3×B-2-2 and Q30×Q24) to investigate the inheritance patterns and chromosomal locations of immature fruit color genes in cucumbers. Consequently, a single recessive gene, w, was identified that controls white immature fruit color. A total of 526 markers, which were derived from published genetic maps, two reference cucumber genomes ("9930" and GY14), and two parents (Q30 and Q24) for which whole-genome sequence information is available, were used to map the target gene w to a 33.0-kb region flanked by two SNP-based markers, ASPCR39262 and ASPCR39229, which are physically located at 39262450 and 39229482 of chromosome 3 ("9930" draft genome assembly), respectively. Gene prediction indicated that four potential genes were located in the target region. One gene that encodes peroxidase is likely to be a valuable candidate gene because quantitative real-time PCR revealed an eightfold difference in its transcriptional level, and several amino acid variations were found when the deduced amino acid sequence was aligned. A co-segregating marker was used synergistically to test its ability to predict the skin colors of 83 dark green/white germplasms, and the validity of its utility in marker-assisted selection was confirmed. Fine mapping of this locus will assist in cloning the gene and in marker-assisted breeding to develop dark green/white cucumber cultivars.

Genome-wide macrosynteny among Fusarium species in the Gibberella fujikuroi complex revealed by amplified fragment length polymorphisms.

The Gibberella fujikuroi complex includes many Fusarium species that cause significant losses in yield and quality of agricultural and forestry crops. Due to their economic importance, whole-genome sequence information has rapidly become available for species including Fusarium circinatum, Fusarium fujikuroi and Fusarium verticillioides, each of which represent one of the three main clades known in this complex. However, no previous studies have explored the genomic commonalities and differences among these fungi. In this study, a previously completed genetic linkage map for an interspecific cross between Fusarium temperatum and F. circinatum, together with genomic sequence data, was utilized to consider the level of synteny between the three Fusarium genomes. Regions that are homologous amongst the Fusarium genomes examined were identified using in silico and pyrosequenced amplified fragment length polymorphism (AFLP) fragment analyses. Homology was determined using BLAST analysis of the sequences, with 777 homologous regions aligned to F. fujikuroi and F. verticillioides. This also made it possible to assign the linkage groups from the interspecific cross to their corresponding chromosomes in F. verticillioides and F. fujikuroi, as well as to assign two previously unmapped supercontigs of F. verticillioides to probable chromosomal locations. We further found evidence of a reciprocal translocation between the distal ends of chromosome 8 and 11, which apparently originated before the divergence of F. circinatum and F. temperatum. Overall, a remarkable level of macrosynteny was observed among the three Fusarium genomes, when comparing AFLP fragments. This study not only demonstrates how in silico AFLPs can aid in the integration of a genetic linkage map to the physical genome, but it also highlights the benefits of using this tool to study genomic synteny and architecture.