Genomic resources for a historical collection of cultivated two-row European spring barley genotypes.

Barley genomic resources are increasing rapidly, with the publication of a barley pangenome as one of the latest developments. Two-row spring barley cultivars are intensely studied as they are the source of high-quality grain for malting and distilling. Here we provide data from a European two-row spring barley population containing 209 different genotypes registered for the UK market between 1830 to 2014. The dataset encompasses RNA-sequencing data from six different tissues across a range of barley developmental stages, phenotypic datasets from two consecutive years of field-grown trials in the United Kingdom, Germany and the USA; and whole genome shotgun sequencing from all cultivars, which was used to complement the RNA-sequencing data for variant calling. The outcomes are a filtered SNP marker file, a phenotypic database and a large gene expression dataset providing a comprehensive resource which allows for downstream analyses like genome wide association studies or expression associations.

The Brassica napus boron deficient inflorescence transcriptome resembles a wounding and infection response.

Oilseed rape and other crops of Brassica napus have a high demand for boron (B). Boron deficiencies result in the inhibition of root growth, and eventually premature flower abortion. Understanding the genetic mechanisms underlying flower abortion in B-limiting conditions could provide the basis to enhance B-efficiency and prevent B-deficiency-related yield losses. In this study, we assessed transcriptomic responses to B-deficiency in diverse inflorescence tissues at multiple time points of soil-grown plants that were phenotypically unaffected by B-deficiency until early flowering. Whilst transcript levels of known B transporters were higher in B-deficient samples, these remained remarkably stable as the duration of B-deficiency increased. Meanwhile, GO-term enrichment analysis indicated a growing response resembling that of a pathogen or pest attack, escalating to a huge transcriptome response in shoot heads at mid-flowering. Grouping differentially expressed genes within this tissue into MapMan functional bins indicated enrichment of genes related to wounding, jasmonic acid and WRKY transcription factors. Individual candidate genes for controlling the "flowering-without-seed-setting" phenotype from within MapMan biotic stress bins include those of the metacaspase family, which have been implicated in orchestrating programmed cell death. Overall temporal expression patterns observed here imply a dynamic response to B-deficiency, first increasing expression of B transporters before recruiting various biotic stress-related pathways to coordinate targeted cell death, likely in response to as yet unidentified B-deficiency induced damage-associated molecular patterns (DAMPs). This response indicates new pathways to target and dissect to control B-deficiency-induced flower abortion and to develop more B-efficient crops.

Large-scale genotyping and phenotyping of a worldwide winter wheat genebank for its use in pre-breeding.

Plant genetic resources (PGR) stored at genebanks are humanity's crop diversity savings for the future. Information on PGR contrasted with modern cultivars is key to select PGR parents for pre-breeding. Genotyping-by-sequencing was performed for 7,745 winter wheat PGR samples from the German Federal ex situ genebank at IPK Gatersleben and for 325 modern cultivars. Whole-genome shotgun sequencing was carried out for 446 diverse PGR samples and 322 modern cultivars and lines. In 19 field trials, 7,683 PGR and 232 elite cultivars were characterized for resistance to yellow rust - one of the major threats to wheat worldwide. Yield breeding values of 707 PGR were estimated using hybrid crosses with 36 cultivars - an approach that reduces the lack of agronomic adaptation of PGR and provides better estimates of their contribution to yield breeding. Cross-validations support the interoperability between genomic and phenotypic data. The here presented data are a stepping stone to unlock the functional variation of PGR for European pre-breeding and are the basis for future breeding and research activities.

On the way to plant data commons - a genotyping use case.

Over the last years it has been observed that the progress in data collection in life science has created increasing demand and opportunities for advanced bioinformatics. This includes data management as well as the individual data analysis and often covers the entire data life cycle. A variety of tools have been developed to store, share, or reuse the data produced in the different domains such as genotyping. Especially imputation, as a subfield of genotyping, requires good Research Data Management (RDM) strategies to enable use and re-use of genotypic data. To aim for sustainable software, it is necessary to develop tools and surrounding ecosystems, which are reusable and maintainable. Reusability in the context of streamlined tools can e.g. be achieved by standardizing the input and output of the different tools and adapting to open and broadly used file formats. By using such established file formats, the tools can also be connected with others, improving the overall interoperability of the software. Finally, it is important to build strong communities that maintain the tools by developing and contributing new features and maintenance updates. In this article, concepts for this will be presented for an imputation service.

Rye B chromosomes differently influence the expression of A chromosome-encoded genes depending on the host species.

The B chromosome (B) is a dispensable component of the genome in many species. To evaluate the impact of Bs on the transcriptome of the standard A chromosomes (A), comparative RNA-seq analyses of rye and wheat anthers with and without additional rye Bs were conducted. In both species, 5-6% of the A-derived transcripts across the entire genomes were differentially expressed in the presence of 2Bs. The GO term enrichment analysis revealed that Bs influence A chromosome encoded processes like "gene silencing"; "DNA methylation or demethylation"; "chromatin silencing"; "negative regulation of gene expression, epigenetic"; "post-embryonic development"; and "chromosome organization." 244 B chromosome responsive A-located genes in + 2B rye and + B wheat shared the same biological function. Positively correlated with the number of Bs, 939 and 1391 B-specific transcripts were identified in + 2B and + 4B wheat samples, respectively. 85% of B-transcripts in + 2B were also found in + 4B transcriptomes. 297 B-specific transcripts were identified in + 2B rye, and 27% were common to the B-derived transcripts identified in + B wheat. Bs encode mobile elements and housekeeping genes, but most B-transcripts were without detectable similarity to known genes. Some of these genes are involved in cell division-related functions like Nuf2 and might indicate their importance in maintaining Bs. The transcriptome analysis provides new insights into the complex interrelationship between standard A chromosomes and supernumerary B chromosomes.

Recommendations for the formatting of Variant Call Format (VCF) files to make plant genotyping data FAIR.

In this opinion article, we discuss the formatting of files from (plant) genotyping studies, in particular the formatting of (meta-) data in Variant Call Format (VCF) files. The flexibility of the VCF format specification facilitates its use as a generic interchange format across domains but can lead to inconsistency between files in the presentation of metadata. To enable fully autonomous machine actionable data flow, generic elements need to be further specified. We strongly support the merits of the FAIR principles and see the need to facilitate them also through technical implementation specifications. VCF files are an established standard for the exchange and publication of genotyping data. Other data formats are also used to capture variant call data (for example, the HapMap format and the gVCF format), but none currently have the reach of VCF. In VCF, only the sites of variation are described, whereas in gVCF, all positions are listed, and confidence values are also provided. For the sake of simplicity, we will only discuss VCF and our recommendations for its use. However, the part of the VCF standard relating to metadata (as opposed to the actual variant calls) defines a syntactic format but no vocabulary, unique identifier or recommended content. In practice, often only sparse (if any) descriptive metadata is included. When descriptive metadata is provided, proprietary metadata fields are frequently added that have not been agreed upon within the community which may limit long-term and comprehensive interoperability. To address this, we propose recommendations for supplying and encoding metadata, focusing on use cases from the plant sciences. We expect there to be overlap, but also divergence, with the needs of other domains.

The barley pan-genome reveals the hidden legacy of mutation breeding.

Genetic diversity is key to crop improvement. Owing to pervasive genomic structural variation, a single reference genome assembly cannot capture the full complement of sequence diversity of a crop species (known as the 'pan-genome'1). Multiple high-quality sequence assemblies are an indispensable component of a pan-genome infrastructure. Barley (Hordeum vulgare L.) is an important cereal crop with a long history of cultivation that is adapted to a wide range of agro-climatic conditions2. Here we report the construction of chromosome-scale sequence assemblies for the genotypes of 20 varieties of barley-comprising landraces, cultivars and a wild barley-that were selected as representatives of global barley diversity. We catalogued genomic presence/absence variants and explored the use of structural variants for quantitative genetic analysis through whole-genome shotgun sequencing of 300 gene bank accessions. We discovered abundant large inversion polymorphisms and analysed in detail two inversions that are frequently found in current elite barley germplasm; one is probably the product of mutation breeding and the other is tightly linked to a locus that is involved in the expansion of geographical range. This first-generation barley pan-genome makes previously hidden genetic variation accessible to genetic studies and breeding.

Multiple wheat genomes reveal global variation in modern breeding.

Advances in genomics have expedited the improvement of several agriculturally important crops but similar efforts in wheat (Triticum spp.) have been more challenging. This is largely owing to the size and complexity of the wheat genome1, and the lack of genome-assembly data for multiple wheat lines2,3. Here we generated ten chromosome pseudomolecule and five scaffold assemblies of hexaploid wheat to explore the genomic diversity among wheat lines from global breeding programs. Comparative analysis revealed extensive structural rearrangements, introgressions from wild relatives and differences in gene content resulting from complex breeding histories aimed at improving adaptation to diverse environments, grain yield and quality, and resistance to stresses4,5. We provide examples outlining the utility of these genomes, including a detailed multi-genome-derived nucleotide-binding leucine-rich repeat protein repertoire involved in disease resistance and the characterization of Sm16, a gene associated with insect resistance. These genome assemblies will provide a basis for functional gene discovery and breeding to deliver the next generation of modern wheat cultivars.

Chromosome-scale genome assembly for the duckweed Spirodela intermedia, integrating cytogenetic maps, PacBio and Oxford Nanopore libraries.

Duckweeds are small, free-floating, morphologically highly reduced organisms belonging to the monocot order Alismatales. They display the most rapid growth among flowering plants, vary ~ 14-fold in genome size and comprise five genera. Spirodela is the phylogenetically oldest genus with only two mainly asexually propagating species: S. polyrhiza (2n = 40; 160 Mbp/1C) and S. intermedia (2n = 36; 160 Mbp/1C). This study combined comparative cytogenetics and de novo genome assembly based on PacBio, Illumina and Oxford Nanopore (ON) reads to obtain the first genome reference for S. intermedia and to compare its genomic features with those of the sister species S. polyrhiza. Both species' genomes revealed little more than 20,000 putative protein-coding genes, very low rDNA copy numbers and a low amount of repetitive sequences, mainly Ty3/gypsy retroelements. The detection of a few new small chromosome rearrangements between both Spirodela species refined the karyotype and the chromosomal sequence assignment for S. intermedia.

Tissue-Specific Transcriptome Analysis Reveals Candidate Transcripts Associated with the Process of Programmed B Chromosome Elimination in Aegilops speltoides.

Some eukaryotes exhibit dramatic genome size differences between cells of different organs, resulting from programmed elimination of chromosomes. Here, we present the first transcriptome analysis of programmed chromosome elimination using laser capture microdissection (LCM)-based isolation of the central meristematic region of Aegilops speltoides embryos where B chromosome (B) elimination occurs. The comparative RNA-seq analysis of meristematic cells of embryos with (Bplus) and without Bs (B0) allowed the identification of 14,578 transcript isoforms (35% out of 41,615 analyzed transcript isoforms) that are differentially expressed during the elimination of Bs. A total of 2908 annotated unigenes were found to be up-regulated in Bplus condition. These genes are either associated with the process of B chromosome elimination or with the presence of B chromosomes themselves. GO enrichment analysis categorized up-regulated transcript isoforms into 27 overrepresented terms related to the biological process, nine terms of the molecular function aspect and three terms of the cellular component category. A total of 2726 annotated unigenes were down-regulated in Bplus condition. Based on strict filtering criteria, 341 B-unique transcript isoforms could be identified in central meristematic cells, of which 70 were functionally annotated. Beside others, genes associated with chromosome segregation, kinetochore function and spindle checkpoint activity were retrieved as promising candidates involved in the process of B chromosome elimination.

Depletion of KNL2 Results in Altered Expression of Genes Involved in Regulation of the Cell Cycle, Transcription, and Development in Arabidopsis.

Centromeres contain specialized nucleosomes at which histone H3 is partially replaced by the centromeric histone H3 variant cenH3 that is required for the assembly, maintenance, and proper function of kinetochores during mitotic and meiotic divisions. Previously, we identified a KINETOCHORE NULL 2 (KNL2) of Arabidopsis thaliana that is involved in the licensing of centromeres for the cenH3 recruitment. We also demonstrated that a knockout mutant for KNL2 shows mitotic and meiotic defects, slower development, reduced growth rate, and fertility. To analyze an effect of KNL2 mutation on global gene transcription of Arabidopsis, we performed RNA-sequencing experiments using seedling and flower bud tissues of knl2 and wild-type plants. The transcriptome data analysis revealed a high number of differentially expressed genes (DEGs) in knl2 plants. The set was enriched in genes involved in the regulation of the cell cycle, transcription, development, and DNA damage repair. In addition to comprehensive information regarding the effects of KNL2 mutation on the global gene expression, physiological changes in plants are also presented, which provides an integrated understanding of the critical role played by KNL2 in plant growth and development.

Complete genome sequence analysis of Archaeoglobus fulgidus strain 7324 (DSM 8774), a hyperthermophilic archaeal sulfate reducer from a North Sea oil field.

Archaeoglobus fulgidus is the type species of genus Archaeoglobus Stetter 1998, a hyperthermophilic sulfate reducing group within the Archaeoglobi class of the euryarchaeota phylum. Members of this genus grow heterotrophically or chemolithoautotrophically with sulfate or thiosulfate as electron acceptors. Except for A. fulgidus strain 7324 and the candidate species "Archaeoglobus lithotrophicus", which both originate from deep oil-fields, the other members of this genus have been recovered from marine hydrothermal systems. Here we describe the features of the A. fulgidus strain 7324 genome as compared to the A. fulgidus VC16 type strain. The 2.3 Mbp genome sequence of strain 7324 shares about 93.5% sequence identity with that of strain VC16T but is about 138 Kbp longer, which is mostly due to two large 'insertions' carrying one extra cdc6 (cell-cycle control protein 6) gene, extra CRISPR elements and mobile genetic elements, a high-GC ncRNA gene (hgcC) and a large number of hypothetical gene functions. A comparison with four other Archaeoglobus spp. genomes identified 1001 core Archaeoglobus genes and more than 2900 pan-genome orthologous genes.

Genome sequence of the pink-pigmented marine bacterium Loktanella hongkongensis type strain (UST950701-009P(T)), a representative of the Roseobacter group.

Loktanella hongkongensis UST950701-009P(T) is a Gram-negative, non-motile and rod-shaped bacterium isolated from a marine biofilm in the subtropical seawater of Hong Kong. When growing as a monospecies biofilm on polystyrene surfaces, this bacterium is able to induce larval settlement and metamorphosis of a ubiquitous polychaete tubeworm Hydroides elegans. The inductive cues are low-molecular weight compounds bound to the exopolymeric matrix of the bacterial cells. In the present study we describe the features of L. hongkongensis strain DSM 17492(T) together with its genome sequence and annotation and novel aspects of its phenotype. The 3,198,444 bp long genome sequence encodes 3104 protein-coding genes and 57 RNA genes. The two unambiguously identified extrachromosomal replicons contain replication modules of the RepB and the Rhodobacteraceae-specific DnaA-like type, respectively.

Comparative genomics of Streptococcus pyogenes M1 isolates differing in virulence and propensity to cause systemic infection in mice.

Streptococcus pyogenes serotype M1 is a frequent cause of severe infections in humans. Some M1 isolates are pathogenic in mice and used in studies on infection pathogenesis. We observed marked differences in murine infections caused by M1 strain SF370, 5448, 5448AP or AP1 which prompted us to sequence the whole genome of isolates 5448 and AP1 for comparative analysis. Strain 5448 is known to acquire inactivating mutations in the CovRS two-component system during mouse infection, producing hypervirulent progeny such as 5448AP. Isolates AP1 and 5448AP, more than 5448, caused disseminating infections that became systemic and lethal. SF370 was not pathogenic. Phages caused gross genetic differences and increased the gene content of AP1 by 8% as compared to 5448 and SF370. Each of six examined M1 genomes contained two CRISPR-Cas systems. Phage insertion destroyed a type II CRISPR-Cas system in AP1 and other strains of serotypes M1, M3, M6 and M24, but not in M1 strains 5448, SF370, MGAS5005, A20 or M1 476. A resulting impaired defence against invading genetic elements could have led to the wealth of phages in AP1. AP1 lacks genetic features of the MGAS5005-like clonal complex including the streptodornase that drives selection for hypervirulent clones with inactivated CovRS system. Still, inactivating mutations in covS were a common genetic feature of AP1 and the MGAS5005-like isolate 5448AP. Abolished expression of the cysteine proteinase SpeB, due to CovRS inactivation could be a common cause for hypervirulence of the two isolates. Moreover, an additional protein H-coding gene and a mutation in the regulator gene rofA distinguished AP1 form other M1 isolates. In conclusion, hypervirulence of S. pyogenes M1 in mice is not limited to the MGAS5005-like genotype.

Complete genome sequence and description of Salinispira pacifica gen. nov., sp. nov., a novel spirochaete isolated form a hypersaline microbial mat.

During a study of the anaerobic microbial community of a lithifying hypersaline microbial mat of Lake 21 on the Kiritimati atoll (Kiribati Republic, Central Pacific) strain L21-RPul-D2(T) was isolated. The closest phylogenetic neighbor was Spirochaeta africana Z-7692(T) that shared a 16S rRNA gene sequence identity value of 90% with the novel strain and thus was only distantly related. A comprehensive polyphasic study including determination of the complete genome sequence was initiated to characterize the novel isolate. Cells of strain L21-RPul-D2(T) had a size of 0.2 - 0.25 × 8-9 µm, were helical, motile, stained Gram-negative and produced an orange carotenoid-like pigment. Optimal conditions for growth were 35°C, a salinity of 50 g/l NaCl and a pH around 7.0. Preferred substrates for growth were carbohydrates and a few carboxylic acids. The novel strain had an obligate fermentative metabolism and produced ethanol, acetate, lactate, hydrogen and carbon dioxide during growth on glucose. Strain L21-RPul-D2(T) was aerotolerant, but oxygen did not stimulate growth. Major cellular fatty acids were C14:0, iso-C15:0, C16:0 and C18:0. The major polar lipids were an unidentified aminolipid, phosphatidylglycerol, an unidentified phospholipid and two unidentified glycolipids. Whole-cell hydrolysates contained L-ornithine as diagnostic diamino acid of the cell wall peptidoglycan. The complete genome sequence was determined and annotated. The genome comprised one circular chromosome with a size of 3.78 Mbp that contained 3450 protein-coding genes and 50 RNA genes, including 2 operons of ribosomal RNA genes. The DNA G + C content was determined from the genome sequence as 51.9 mol%. There were no predicted genes encoding cytochromes or enzymes responsible for the biosynthesis of respiratory lipoquinones. Based on significant differences to the uncultured type species of the genus Spirochaeta, S. plicatilis, as well as to any other phylogenetically related cultured species it is suggested to place strain L21-RPul-D2(T) (=DSM 27196(T) = JCM 18663(T)) in a novel species and genus, for which the name Salinispira pacifica gen. nov., sp. nov. is proposed.

Genome sequence of the Roseovarius mucosus type strain (DSM 17069(T)), a bacteriochlorophyll a-containing representative of the marine Roseobacter group isolated from the dinoflagellate Alexandrium ostenfeldii.

Roseovarius mucosus Biebl et al. 2005 is a bacteriochlorophyll a-producing representative of the marine Roseobacter group within the alphaproteobacterial family Rhodobacteraceae, which was isolated from the dinoflagellate Alexandrium ostenfeldii. The marine Roseobacter group was found to be abundant in the ocean and plays an important role for global and biogeochemical processes. Here we describe the features of the R. mucosus strain DFL-24(T) together with its genome sequence and annotation generated from a culture of DSM 17069(T). The 4,247,724 bp containing genome sequence encodes 4,194 protein-coding genes and 57 RNA genes. In addition to the presence of four plasmids, genome analysis revealed the presence of genes associated with host colonization, DMSP utilization, cytotoxins, and quorum sensing that could play a role in the interrelationship of R. mucosus with the dinoflagellate A. ostenfeldii and other marine organisms. Furthermore, the genome encodes genes associated with mixotrophic growth, where both reduced inorganic compounds for lithotrophic growth and a photoheterotrophic lifestyle using light as additional energy source could be used.

Genome Sequence of Gammaproteobacterial Pseudohaliea rubra Type Strain DSM 19751, Isolated from Coastal Seawater of the Mediterranean Sea.

Pseudohaliea rubra strain DSM 19751(T) is an aerobic marine gammaproteobacterium that was isolated from surface coastal seawater of the Mediterranean Sea. Here, we present its genome sequence and annotation. Genome analysis revealed the presence of genes involved in the synthesis of bacteriochlorophyll-a and the reserve compound glycogen.

Complete genome sequence of the bacteriochlorophyll a-containing Roseibacterium elongatum type strain (DSM 19469(T)), a representative of the Roseobacter group isolated from Australian coast sand.

Roseibacterium elongatum Suzuki et al. 2006 is a pink-pigmented and bacteriochlorophyll a-producing representative of the Roseobacter group within the alphaproteobacterial family Rhodobacteraceae. Representatives of the marine 'Roseobacter group' were found to be abundant in the ocean and play an important role in global and biogeochemical processes. In the present study we describe the features of R. elongatum strain OCh 323(T) together with its genome sequence and annotation. The 3,555,102 bp long genome consists of one circular chromosome with no extrachromosomal elements and is one of the smallest known Roseobacter genomes. It contains 3,540 protein-coding genes and 59 RNA genes. Genome analysis revealed the presence of a photosynthetic gene cluster, which putatively enables a photoheterotrophic lifestyle. Gene sequences associated with quorum sensing, motility, surface attachment, and thiosulfate and carbon monoxide oxidation could be detected. The genome was sequenced as part of the activities of the Transregional Collaborative Research Centre 51 (TRR51) funded by the German Research Foundation (DFG).

Genome sequence of the Wenxinia marina type strain (DSM 24838(T)), a representative of the Roseobacter group isolated from oilfield sediments.

Wenxinia marina Ying et al. 2007 is the type species of the genus Wenxinia, a representative of the Roseobacter group within the alphaproteobacterial family Rhodobacteraceae, isolated from oilfield sediments of the South China Sea. This family was shown to harbor the most abundant bacteria especially from coastal and polar waters, but was also found in microbial mats, sediments and attached to different kind of surfaces. Here we describe the features of W. marina strain HY34(T) together with the genome sequence and annotation of strain DSM 24838(T) and novel aspects of its phenotype. The 4,181,754 bp containing genome sequence encodes 4,047 protein-coding genes and 59 RNA genes. The genome of W. marina DSM 24838(T) was sequenced as part of the activities of the Genomic Encyclopedia of Type Strains, Phase I: the one thousand microbial genomes (KMG) project funded by the DoE and the Transregional Collaborative Research Centre 51 (TRR51) funded by the German Research Foundation (DFG).

Genome sequence of the pink to light reddish-pigmented Rubellimicrobium mesophilum type strain (DSM 19309(T)), a representative of the Roseobacter group isolated from soil, and emended description of the species.

Rubellimicrobium mesophilum Dastager et al. 2008 is a mesophilic and light reddish-pigmented representative of the Roseobacter group within the alphaproteobacterial family Rhodobacteraceae. Representatives of the Roseobacter group play an important role in the marine biogeochemical cycles and were found in a broad variety of marine environments associated with algal blooms, different kinds of sediments, and surfaces of invertebrates and vertebrates. Roseobacters were shown to be widely distributed, especially within the total bacterial community found in coastal waters, as well as in mixed water layers of the open ocean. Here we describe the features of R. mesophilum strain MSL-20(T) together with its genome sequence and annotation generated from a culture of DSM 19309(T). The 4,927,676 bp genome sequence consists of one chromosome and probably one extrachromosomal element. It contains 5,082 protein-coding genes and 56 RNA genes. As previously reported, the G+C content is significantly different from the actual genome sequence-based G+C content and as the type strain tests positively for oxidase, the species description is emended accordingly. The genome was sequenced as part of the activities of the Transregional Collaborative Research Centre 51 (TRR51) funded by the German Research Foundation (DFG).