Lymphocyte pathway analysis using naturally lymphocyte-deficient fish.

Comparative phylogenetic analyses are of potential value to establish the essential components of genetic networks underlying physiological traits. For species that naturally lack particular lymphocyte lineages, we show here that this strategy readily distinguishes trait-specific actors from pleiotropic components of the genetic network governing lymphocyte differentiation. Previously, three of the four members of the DNA polymerase X family have been implicated in the junctional diversification process during the somatic assembly of antigen receptors. Our phylogenetic analysis indicates that the presence of terminal deoxynucleotidyl transferase is strictly associated with the facility of V(D)J recombination, whereas PolL and PolM genes are retained even in species lacking Rag-mediated somatic diversification of antigen receptor genes.

De Novo Genome Assembly of a Plant-Associated Rhodococcus qingshengii Strain (RL1) Isolated from Eruca sativa Mill. and Showing Plant Growth-Promoting Properties.

Rhodococcus qingshengii RL1 was isolated from surface-sterilized leaves of Eruca sativa Mill. and shows plant growth-promoting (PGP) properties. The de novo genome assembly consists of one chromosome with 6,253,838 bp and two plasmids with 144,038 bp and 448,745 bp. Many genes could be identified reflecting its PGP potential.

An aerobic eukaryotic parasite with functional mitochondria that likely lacks a mitochondrial genome.

Dinoflagellates are microbial eukaryotes that have exceptionally large nuclear genomes; however, their organelle genomes are small and fragmented and contain fewer genes than those of other eukaryotes. The genus Amoebophrya (Syndiniales) comprises endoparasites with high genetic diversity that can infect other dinoflagellates, such as those forming harmful algal blooms (e.g., Alexandrium). We sequenced the genome (~100 Mb) of Amoebophrya ceratii to investigate the early evolution of genomic characters in dinoflagellates. The A. ceratii genome encodes almost all essential biosynthetic pathways for self-sustaining cellular metabolism, suggesting a limited dependency on its host. Although dinoflagellates are thought to have descended from a photosynthetic ancestor, A. ceratii appears to have completely lost its plastid and nearly all genes of plastid origin. Functional mitochondria persist in all life stages of A. ceratii, but we found no evidence for the presence of a mitochondrial genome. Instead, all mitochondrial proteins appear to be lost or encoded in the A. ceratii nucleus.

Chromosome-scale comparative sequence analysis unravels molecular mechanisms of genome dynamics between two wheat cultivars.

BACKGROUND: Recent improvements in DNA sequencing and genome scaffolding have paved the way to generate high-quality de novo assemblies of pseudomolecules representing complete chromosomes of wheat and its wild relatives. These assemblies form the basis to compare the dynamics of wheat genomes on a megabase scale. RESULTS: Here, we provide a comparative sequence analysis of the 700-megabase chromosome 2D between two bread wheat genotypes-the old landrace Chinese Spring and the elite Swiss spring wheat line 'CH Campala Lr22a'. Both chromosomes were assembled into megabase-sized scaffolds. There is a high degree of sequence conservation between the two chromosomes. Analysis of large structural variations reveals four large indels of more than 100 kb. Based on the molecular signatures at the breakpoints, unequal crossing over and double-strand break repair were identified as the molecular mechanisms that caused these indels. Three of the large indels affect copy number of NLRs, a gene family involved in plant immunity. Analysis of SNP density reveals four haploblocks of 4, 8, 9 and 48 Mb with a 35-fold increased SNP density compared to the rest of the chromosome. Gene content across the two chromosomes was highly conserved. Ninety-nine percent of the genic sequences were present in both genotypes and the fraction of unique genes ranged from 0.4 to 0.7%. CONCLUSIONS: This comparative analysis of two high-quality chromosome assemblies enabled a comprehensive assessment of large structural variations and gene content. The insight obtained from this analysis will form the basis of future wheat pan-genome studies.

A Tandem Amino Acid Residue Motif in Guard Cell SLAC1 Anion Channel of Grasses Allows for the Control of Stomatal Aperture by Nitrate.

The latest major group of plants to evolve were the grasses. These became important in the mid-Paleogene about 40 million years ago. During evolution, leaf CO2 uptake and transpirational water loss were optimized by the acquisition of grass-specific stomatal complexes. In contrast to the kidney-shaped guard cells (GCs) typical of the dicots such as Arabidopsis, in the grasses and agronomically important cereals, the GCs are dumbbell shaped and are associated with morphologically distinct subsidiary cells (SCs). We studied the molecular basis of GC action in the major cereal crop barley. Upon feeding ABA to xylem sap of an intact barley leaf, stomata closed in a nitrate-dependent manner. This process was initiated by activation of GC SLAC-type anion channel currents. HvSLAC1 expressed in Xenopus oocytes gave rise to S-type anion currents that increased several-fold upon stimulation with >3 mM nitrate. We identified a tandem amino acid residue motif that within the SLAC1 channels differs fundamentally between monocots and dicots. When the motif of nitrate-insensitive dicot Arabidopsis SLAC1 was replaced by the monocot signature, AtSLAC1 converted into a grass-type like nitrate-sensitive channel. Our work reveals a fundamental difference between monocot and dicot GCs and prompts questions into the selective pressures during evolution that resulted in fundamental changes in the regulation of SLAC1 function.

A chromosome conformation capture ordered sequence of the barley genome.

Cereal grasses of the Triticeae tribe have been the major food source in temperate regions since the dawn of agriculture. Their large genomes are characterized by a high content of repetitive elements and large pericentromeric regions that are virtually devoid of meiotic recombination. Here we present a high-quality reference genome assembly for barley (Hordeum vulgare L.). We use chromosome conformation capture mapping to derive the linear order of sequences across the pericentromeric space and to investigate the spatial organization of chromatin in the nucleus at megabase resolution. The composition of genes and repetitive elements differs between distal and proximal regions. Gene family analyses reveal lineage-specific duplications of genes involved in the transport of nutrients to developing seeds and the mobilization of carbohydrates in grains. We demonstrate the importance of the barley reference sequence for breeding by inspecting the genomic partitioning of sequence variation in modern elite germplasm, highlighting regions vulnerable to genetic erosion.

Construction of a map-based reference genome sequence for barley, Hordeum vulgare L.

Barley (Hordeum vulgare L.) is a cereal grass mainly used as animal fodder and raw material for the malting industry. The map-based reference genome sequence of barley cv. 'Morex' was constructed by the International Barley Genome Sequencing Consortium (IBSC) using hierarchical shotgun sequencing. Here, we report the experimental and computational procedures to (i) sequence and assemble more than 80,000 bacterial artificial chromosome (BAC) clones along the minimum tiling path of a genome-wide physical map, (ii) find and validate overlaps between adjacent BACs, (iii) construct 4,265 non-redundant sequence scaffolds representing clusters of overlapping BACs, and (iv) order and orient these BAC clusters along the seven barley chromosomes using positional information provided by dense genetic maps, an optical map and chromosome conformation capture sequencing (Hi-C). Integrative access to these sequence and mapping resources is provided by the barley genome explorer (BARLEX).

Comparative genomics to explore phylogenetic relationship, cryptic sexual potential and host specificity of Rhynchosporium species on grasses.

BACKGROUND: The Rhynchosporium species complex consists of hemibiotrophic fungal pathogens specialized to different sweet grass species including the cereal crops barley and rye. A sexual stage has not been described, but several lines of evidence suggest the occurrence of sexual reproduction. Therefore, a comparative genomics approach was carried out to disclose the evolutionary relationship of the species and to identify genes demonstrating the potential for a sexual cycle. Furthermore, due to the evolutionary very young age of the five species currently known, this genus appears to be well-suited to address the question at the molecular level of how pathogenic fungi adapt to their hosts. RESULTS: The genomes of the different Rhynchosporium species were sequenced, assembled and annotated using ab initio gene predictors trained on several fungal genomes as well as on Rhynchosporium expressed sequence tags. Structures of the rDNA regions and genome-wide single nucleotide polymorphisms provided a hypothesis for intra-genus evolution. Homology screening detected core meiotic genes along with most genes crucial for sexual recombination in ascomycete fungi. In addition, a large number of cell wall-degrading enzymes that is characteristic for hemibiotrophic and necrotrophic fungi infecting monocotyledonous hosts were found. Furthermore, the Rhynchosporium genomes carry a repertoire of genes coding for polyketide synthases and non-ribosomal peptide synthetases. Several of these genes are missing from the genome of the closest sequenced relative, the poplar pathogen Marssonina brunnea, and are possibly involved in adaptation to the grass hosts. Most importantly, six species-specific genes coding for protein effectors were identified in R. commune. Their deletion yielded mutants that grew more vigorously in planta than the wild type. CONCLUSION: Both cryptic sexuality and secondary metabolites may have contributed to host adaptation. Most importantly, however, the growth-retarding activity of the species-specific effectors suggests that host adaptation of R. commune aims at extending the biotrophic stage at the expense of the necrotrophic stage of pathogenesis. Like other apoplastic fungi Rhynchosporium colonizes the intercellular matrix of host leaves relatively slowly without causing symptoms, reminiscent of the development of endophytic fungi. Rhynchosporium may therefore become an object for studying the mutualism-parasitism transition.

Genetic Factors of the Disease Course after Sepsis: A Genome-Wide Study for 28Day Mortality.

Sepsis is the dysregulated host response to an infection which leads to life-threatening organ dysfunction that varies by host genomic factors. We conducted a genome-wide association study (GWAS) in 740 adult septic patients and focused on 28day mortality as outcome. Variants with suggestive evidence for an association (p≤10-5) were validated in two additional GWA studies (n=3470) and gene coding regions related to the variants were assessed in an independent exome sequencing study (n=74). In the discovery GWAS, we identified 243 autosomal variants which clustered in 14 loci (p≤10-5). The best association signal (rs117983287; p=8.16×10-8) was observed for a missense variant located at chromosome 9q21.2 in the VPS13A gene. VPS13A was further supported by additional GWAS (p=0.03) and sequencing data (p=0.04). Furthermore, CRISPLD2 (p=5.99×10-6) and a region on chromosome 13q21.33 (p=3.34×10-7) were supported by both our data and external biological evidence. We found 14 loci with suggestive evidence for an association with 28day mortality and found supportive, converging evidence for three of them in independent data sets. Elucidating the underlying biological mechanisms of VPS13A, CRISPLD2, and the chromosome 13 locus should be a focus of future research activities.

Genetic Factors of the Disease Course After Sepsis: Rare Deleterious Variants Are Predictive.

Sepsis is a life-threatening organ dysfunction caused by dysregulated host response to infection. For its clinical course, host genetic factors are important and rare genomic variants are suspected to contribute. We sequenced the exomes of 59 Greek and 15 German patients with bacterial sepsis divided into two groups with extremely different disease courses. Variant analysis was focusing on rare deleterious single nucleotide variants (SNVs). We identified significant differences in the number of rare deleterious SNVs per patient between the ethnic groups. Classification experiments based on the data of the Greek patients allowed discrimination between the disease courses with estimated sensitivity and specificity>75%. By application of the trained model to the German patients we observed comparable discriminatory properties despite lower population-specific rare SNV load. Furthermore, rare SNVs in genes of cell signaling and innate immunity related pathways were identified as classifiers discriminating between the sepsis courses. Sepsis patients with favorable disease course after sepsis, even in the case of unfavorable preconditions, seem to be affected more often by rare deleterious SNVs in cell signaling and innate immunity related pathways, suggesting a protective role of impairments in these processes against a poor disease course.

The multicellularity genes of dictyostelid social amoebas.

The evolution of multicellularity enabled specialization of cells, but required novel signalling mechanisms for regulating cell differentiation. Early multicellular organisms are mostly extinct and the origins of these mechanisms are unknown. Here using comparative genome and transcriptome analysis across eight uni- and multicellular amoebozoan genomes, we find that 80% of proteins essential for the development of multicellular Dictyostelia are already present in their unicellular relatives. This set is enriched in cytosolic and nuclear proteins, and protein kinases. The remaining 20%, unique to Dictyostelia, mostly consists of extracellularly exposed and secreted proteins, with roles in sensing and recognition, while several genes for synthesis of signals that induce cell-type specialization were acquired by lateral gene transfer. Across Dictyostelia, changes in gene expression correspond more strongly with phenotypic innovation than changes in protein functional domains. We conclude that the transition to multicellularity required novel signals and sensors rather than novel signal processing mechanisms.

Longitudinal RNA-Seq Analysis of Vertebrate Aging Identifies Mitochondrial Complex I as a Small-Molecule-Sensitive Modifier of Lifespan.

Mutations and genetic variability affect gene expression and lifespan, but the impact of variations in gene expression within individuals on their aging-related mortality is poorly understood. We performed a longitudinal study in the short-lived killifish, Nothobranchius furzeri, and correlated quantitative variations in gene expression during early adult life with lifespan. Shorter- and longer-lived individuals differ in their gene expression before the onset of aging-related mortality; differences in gene expression are more pronounced early in life. We identified mitochondrial respiratory chain complex I as a hub in a module of genes whose expression is negatively correlated with lifespan. Accordingly, partial pharmacological inhibition of complex I by the small molecule rotenone reversed aging-related regulation of gene expression and extended lifespan in N. furzeri by 15%. These results support the use of N. furzeri as a vertebrate model for identifying the protein targets, pharmacological modulators, and individual-to-individual variability associated with aging.

Multiplex sequencing of bacterial artificial chromosomes for assembling complex plant genomes.

Hierarchical shotgun sequencing remains the method of choice for assembling high-quality reference sequences of complex plant genomes. The efficient exploitation of current high-throughput technologies and powerful computational facilities for large-insert clone sequencing necessitates the sequencing and assembly of a large number of clones in parallel. We developed a multiplexed pipeline for shotgun sequencing and assembling individual bacterial artificial chromosomes (BACs) using the Illumina sequencing platform. We illustrate our approach by sequencing 668 barley BACs (Hordeum vulgare L.) in a single Illumina HiSeq 2000 lane. Using a newly designed parallelized computational pipeline, we obtained sequence assemblies of individual BACs that consist, on average, of eight sequence scaffolds and represent >98% of the genomic inserts. Our BAC assemblies are clearly superior to a whole-genome shotgun assembly regarding contiguity, completeness and the representation of the gene space. Our methods may be employed to rapidly obtain high-quality assemblies of a large number of clones to assemble map-based reference sequences of plant and animal species with complex genomes by sequencing along a minimum tiling path.

Insights into Sex Chromosome Evolution and Aging from the Genome of a Short-Lived Fish.

The killifish Nothobranchius furzeri is the shortest-lived vertebrate that can be bred in the laboratory. Its rapid growth, early sexual maturation, fast aging, and arrested embryonic development (diapause) make it an attractive model organism in biomedical research. Here, we report a draft sequence of its genome that allowed us to uncover an intra-species Y chromosome polymorphism representing-in real time-different stages of sex chromosome formation that display features of early mammalian XY evolution "in action." Our data suggest that gdf6Y, encoding a TGF-ß family growth factor, is the master sex-determining gene in N. furzeri. Moreover, we observed genomic clustering of aging-related genes, identified genes under positive selection, and revealed significant similarities of gene expression profiles between diapause and aging, particularly for genes controlling cell cycle and translation. The annotated genome sequence is provided as an online resource (http://www.nothobranchius.info/NFINgb).

Comprehensive insights in the Mycobacterium avium subsp. paratuberculosis genome using new WGS data of sheep strain JIII-386 from Germany.

Mycobacterium avium (M. a.) subsp. paratuberculosis (MAP) - the etiologic agent of Johne's disease - affects cattle, sheep and other ruminants worldwide. To decipher phenotypic differences among sheep and cattle strains (belonging to MAP-S [Type-I/III] respectively MAP-C [Type-II]) comparative genome analysis needs data from diverse isolates originating from different geographic regions of the world. The current study presents the so far best assembled genome of a MAP-S-strain: sheep isolate JIII-386 from Germany. One newly sequenced cattle isolate (JII-1961, Germany), four published MAP strains of MAP-C and MAP-S from U.S. and Australia and M. a. subsp. hominissuis (MAH) strain 104 were used for assembly improvement and comparisons. All genomes were annotated by BacProt and results compared with NCBI annotation. Corresponding protein-coding sequences (CDSs) were detected, but also CDSs that were exclusively determined either by NCBI or BacProt. A new Shine-Dalgarno sequence motif (5'AGCTGG3') was extracted. Novel CDSs including PE-PGRS family protein genes and about 80 non-coding RNAs exhibiting high sequence conservation are presented. Previously found genetic differences between MAP-types are partially revised. Four out of ten assumed MAP-S-specific large sequence polymorphism regions (LSPSs) are still present in MAP-C strains; new LSPSs were identified. Independently of the regional origin of the strains, the number of individual CDSs and single nucleotide variants confirm the strong similarity of MAP-C strains and show higher diversity among MAP-S strains. This study gives ambiguous results regarding the hypothesis that MAP-S is the evolutionary intermediate between MAH and MAP-C, but it clearly shows a higher similarity of MAP to MAH than to M. intracellulare.

Comparative analysis of mitochondrial genomes from closely related Rhynchosporium species reveals extensive intron invasion.

We sequenced and annotated the complete mitochondrial (mt) genomes of four closely related Rhynchosporium species that diverged ∼14,000-35,000years ago. During this time frame, three of the mt genomes expanded significantly due to an invasion of introns into three genes (cox1, cox2, and nad5). The enlarged mt genomes contained ∼40% introns compared to 8.1% in uninvaded relatives. Many intron gains were accompanied by co-conversion of flanking exonic regions. The comparative analysis revealed a highly variable set of non-intronic, free-standing ORFs of unknown function (uORFs). This is consistent with a rapidly evolving accessory compartment in the mt genome of these closely related species. Only one free-standing uORF was shared among all mt genomes analyzed. This uORF had a mutation rate similar to the core mt protein-encoding genes, suggesting conservation of function among the species. The nucleotide composition of the core protein-encoding genes significantly differed from those of introns and uORFs. The mt mutation rate was 77 times higher than the nuclear mutation rate, indicating that the phylogeny inferred from mt genes may better resolve the phylogenetic relationships among closely related Rhynchosporium species than phylogenies inferred from nuclear genes.

A sequence-ready physical map of barley anchored genetically by two million single-nucleotide polymorphisms.

Barley (Hordeum vulgare) is an important cereal crop and a model species for Triticeae genomics. To lay the foundation for hierarchical map-based sequencing, a genome-wide physical map of its large and complex 5.1 billion-bp genome was constructed by high-information content fingerprinting of almost 600,000 bacterial artificial chromosomes representing 14-fold haploid genome coverage. The resultant physical map comprises 9,265 contigs with a cumulative size of 4.9 Gb representing 96% of the physical length of the barley genome. The reliability of the map was verified through extensive genetic marker information and the analysis of topological networks of clone overlaps. A minimum tiling path of 66,772 minimally overlapping clones was defined that will serve as a template for hierarchical clone-by-clone map-based shotgun sequencing. We integrated whole-genome shotgun sequence data from the individuals of two mapping populations with published bacterial artificial chromosome survey sequence information to genetically anchor the physical map. This novel approach in combination with the comprehensive whole-genome shotgun sequence data sets allowed us to independently validate and improve a previously reported physical and genetic framework. The resources developed in this study will underpin fine-mapping and cloning of agronomically important genes and the assembly of a draft genome sequence.

GenColors-based comparative genome databases for small eukaryotic genomes.

Many sequence data repositories can give a quick and easily accessible overview on genomes and their annotations. Less widespread is the possibility to compare related genomes with each other in a common database environment. We have previously described the GenColors database system (http://gencolors.fli-leibniz.de) and its applications to a number of bacterial genomes such as Borrelia, Legionella, Leptospira and Treponema. This system has an emphasis on genome comparison. It combines data from related genomes and provides the user with an extensive set of visualization and analysis tools. Eukaryote genomes are normally larger than prokaryote genomes and thus pose additional challenges for such a system. We have, therefore, adapted GenColors to also handle larger datasets of small eukaryotic genomes and to display eukaryotic gene structures. Further recent developments include whole genome views, genome list options and, for bacterial genome browsers, the display of horizontal gene transfer predictions. Two new GenColors-based databases for two fungal species (http://fgb.fli-leibniz.de) and for four social amoebas (http://sacgb.fli-leibniz.de) were set up. Both new resources open up a single entry point for related genomes for the amoebozoa and fungal research communities and other interested users. Comparative genomics approaches are greatly facilitated by these resources.

Unusual combinatorial involvement of poly-A/T tracts in organizing genes and chromatin in Dictyostelium.

Dictyostelium discoideum is an amoebozoa that exists in both a free-living unicellular and a multicellular form. It is situated in a deep branch in the evolutionary tree and is particularly noteworthy in having a very A/T-rich genome. Dictyostelium provides an ideal system to examine the extreme to which nucleotide bias may be employed in organizing promoters, genes, and nucleosomes across a genome. We find that Dictyostelium genes are demarcated precisely at their 5' ends by poly-T tracts and precisely at their 3' ends by poly-A tracts. These tracts are also associated with nucleosome-free regions and are embedded with precisely positioned TATA boxes. Homo- and heteropolymeric tracts of A and T demarcate nucleosome border regions. Together, these findings reveal the presence of a variety of functionally distinct polymeric A/T elements. Strikingly, Dictyostelium chromatin may be organized in di-nucleosome units but is otherwise organized as in animals. This includes a +1 nucleosome in a position that predicts the presence of a paused RNA polymerase II. Indeed, we find a strong phylogenetic relationship between the presence of the NELF pausing factor and positioning of the +1 nucleosome. Pausing and +1 nucleosome positioning may have coevolved in animals.

Sequencing of BAC pools by different next generation sequencing platforms and strategies.

BACKGROUND: Next generation sequencing of BACs is a viable option for deciphering the sequence of even large and highly repetitive genomes. In order to optimize this strategy, we examined the influence of read length on the quality of Roche/454 sequence assemblies, to what extent Illumina/Solexa mate pairs (MPs) improve the assemblies by scaffolding and whether barcoding of BACs is dispensable. RESULTS: Sequencing four BACs with both FLX and Titanium technologies revealed similar sequencing accuracy, but showed that the longer Titanium reads produce considerably less misassemblies and gaps. The 454 assemblies of 96 barcoded BACs were improved by scaffolding 79% of the total contig length with MPs from a non-barcoded library.Assembly of the unmasked 454 sequences without separation by barcodes revealed chimeric contig formation to be a major problem, encompassing 47% of the total contig length. Masking the sequences reduced this fraction to 24%. CONCLUSION: Optimal BAC pool sequencing should be based on the longest available reads, with barcoding essential for a comprehensive assessment of both repetitive and non-repetitive sequence information. When interest is restricted to non-repetitive regions and repeats are masked prior to assembly, barcoding is non-essential. In any case, the assemblies can be improved considerably by scaffolding with non-barcoded BAC pool MPs.