Search

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., Comparison of multiple genome assemblies from wheat reveals extensive diversity that results from the complex breeding history of wheat and provides a basis for further potential improvements to this important food crop.

Background Sequence exchange between homologous chromosomes through crossing over and gene conversion is highly conserved among eukaryotes, contributing to genome stability and genetic diversity. A lack of recombination limits breeding efforts in crops; therefore, increasing recombination rates can reduce linkage drag and generate new genetic combinations. Results We use computational analysis of 13 recombinant inbred mapping populations to assess crossover and gene conversion frequency in the hexaploid genome of wheat (Triticum aestivum). We observe that high-frequency crossover sites are shared between populations and that closely related parents lead to populations with more similar crossover patterns. We demonstrate that gene conversion is more prevalent and covers more of the genome in wheat than in other plants, making it a critical process in the generation of new haplotypes, particularly in centromeric regions where crossovers are rare. We identify quantitative trait loci for altered gene conversion and crossover frequency and confirm functionality for a novel RecQ helicase gene that belongs to an ancient clade that is missing in some plant lineages including Arabidopsis. Conclusions This is the first gene to be demonstrated to be involved in gene conversion in wheat. Harnessing the RecQ helicase has the potential to break linkage drag utilizing widespread gene conversions. Electronic supplementary material The online version of this article (10.1186/s13059-019-1675-6) contains supplementary material, which is available to authorized users.

BackgroundBarley (Hordeum vulgare) is one of the most important crops worldwide and is also considered a research model for the large-genome small grain temperate cereals. Despite genomic resources improving all the time, they are limited for the cv. Golden Promise, the most efficient genotype for genetic transformation.FindingsWe have developed a barley cv. Golden Promise reference assembly integrating Illumina paired-end reads, long mate-pair reads, Dovetail Chicago in vitro proximity ligation libraries and chromosome conformation capture sequencing (Hi-C) libraries into a contiguous reference assembly. The assembled genome of 7 chromosomes and 4.13Gb in size, has a super-scaffold N50 after Chicago libraries of 4.14Mb and contains only 2.2% gaps. Using BUSCO (benchmarking universal single copy orthologous genes) as evaluation the genome assembly contains 95.2% of complete and single copy genes from the plant database.ConclusionsA high-quality Golden Promise reference assembly will be useful and utilised by the whole barley research community but will prove particularly useful for CRISPR-Cas9 experiments.

With the advent of high-throughput DNA sequencing technology, the genomic sequence of many disparate species has led to the relatively new discipline of genomics, the study of genome structure, function and evolution. Much work has been focused on the role of whole genome duplications (WGD) in the architecture of extant vertebrate genomes, particularly those of teleost fishes which underwent a WGD early in the teleost radiation > 230 million years ago (mya). Our past work has focused on the fate of duplicated copies of a multigene family coding for the intracellular lipid-binding protein (iLBP) genes in the teleost fishes. To define the evolutionary processes that determined the fate of duplicated genes and generated the structure of extant fish genomes, however, requires comparative genomic analysis with a fish lineage that diverged before the teleost WGD, such as the spotted gar ( Lepisosteus oculatus ), an ancient, air-breathing, ray-finned fish. Here, we describe the genomic organization, chromosomal location and tissue-specific expression of a subfamily of the iLBP genes that code for fatty acid-binding proteins (Fabps) in spotted gar. Based on this work, we have defined the minimum suite of fabp genes prior to their duplication in the teleost lineages ~ 230–400 mya. Spotted gar, therefore, serves as an appropriate outgroup, or ancestral/ancient fish, that did not undergo the teleost-specific WGD. As such, analyses of the spatio-temporal regulation of spotted gar genes provides a foundation to determine whether the duplicated fabp genes have been retained in teleost genomes owing to either sub- or neofunctionalization.

Advances in genome sequencing and assembly technologies are generating many high quality genome sequences, but assemblies of large, repeat-rich polyploid genomes, such as that of bread wheat, remain fragmented and incomplete. We have generated a new wheat whole-genome shotgun sequence assembly using a combination of optimised data types and an assembly algorithm designed to deal with large and complex genomes. The new assembly represents more than 78% of the genome with a scaffold N50 of 88.8kbp that has a high fidelity to the input data. Our new annotation combines strand-specific Illumina RNAseq and PacBio full-length cDNAs to identify 104,091 high confidence protein-coding genes and 10,156 non-coding RNA genes. We confirmed three known and identified one novel genome rearrangements. Our approach enables the rapid and scalable assembly of wheat genomes, the identification of structural variants, and the definition of complete gene models, all powerful resources for trait analysis and breeding of this key global crop. [Supplemental material is available for this article.]

Primula vulgaris (primrose) exhibits heterostyly: plants produce self-incompatible pin- or thrum-form flowers, with anthers and stigma at reciprocal heights. Darwin concluded that this arrangement promotes insect-mediated cross-pollination; later studies revealed control by a cluster of genes, or supergene, known as the S (Style length) locus. The P. vulgaris S locus is absent from pin plants and hemizygous in thrum plants (thrum-specific); mutation of S locus genes produces self-fertile homostyle flowers with anthers and stigma at equal heights. Here, we present a 411 Mb P. vulgaris genome assembly of a homozygous inbred long homostyle, representing ~87% of the genome. We annotate over 24,000 P. vulgaris genes, and reveal more genes up-regulated in thrum than pin flowers. We show reduced genomic read coverage across the S locus in other Primula species, including P. veris, where we define the conserved structure and expression of the S locus genes in thrum. Further analysis reveals the S locus has elevated repeat content (64%) compared to the wider genome (37%). Our studies suggest conservation of S locus genetic architecture in Primula, and provide a platform for identification and evolutionary analysis of the S locus and downstream targets that regulate heterostyly in diverse heterostylous species.

Summary In Primula vulgaris outcrossing is promoted through reciprocal herkogamy with insect-mediated cross-pollination between pin and thrum form flowers. Development of heteromorphic flowers is coordinated by genes at the S locus. To underpin construction of a genetic map facilitating isolation of these S locus genes, we have characterised Oakleaf, a novel S locus-linked mutant phenotype. We combine phenotypic observation of flower and leaf development, with classical genetic analysis and next-generation sequencing to address the molecular basis of Oakleaf. Oakleaf is a dominant mutation that affects both leaf and flower development; plants produce distinctive lobed leaves, with occasional ectopic meristems on the veins. This phenotype is reminiscent of overexpression of Class I KNOX-homeodomain transcription factors. We describe the structure and expression of all eight P. vulgaris PvKNOX genes in both wild-type and Oakleaf plants, and present comparative transcriptome analysis of leaves and flowers from Oakleaf and wild-type plants. Oakleaf provides a new phenotypic marker for genetic analysis of the Primula S locus. We show that none of the Class I PvKNOX genes are strongly upregulated in Oakleaf leaves and flowers, and identify cohorts of 507 upregulated and 314 downregulated genes in the Oakleaf mutant.

Summary Heteromorphic flower development in Primula is controlled by the S locus. The S locus genes, which control anther position, pistil length and pollen size in pin and thrum flowers, have not yet been characterized. We have integrated S‐linked genes, marker sequences and mutant phenotypes to create a map of the P. vulgaris S locus region that will facilitate the identification of key S locus genes. We have generated, sequenced and annotated BAC sequences spanning the S locus, and identified its chromosomal location.We have employed a combination of classical genetics and three‐point crosses with molecular genetic analysis of recombinants to generate the map. We have characterized this region by Illumina sequencing and bioinformatic analysis, together with chromosome in situ hybridization.We present an integrated genetic and physical map across the P. vulgaris S locus flanked by phenotypic and DNA sequence markers. BAC contigs encompass a 1.5‐Mb genomic region with 1 Mb of sequence containing 82 S‐linked genes anchored to overlapping BACs. The S locus is located close to the centromere of the largest metacentric chromosome pair.These data will facilitate the identification of the genes that orchestrate heterostyly in Primula and enable evolutionary analyses of the S locus.

Accelerating international trade and climate change make pathogen spread an increasing concern. Hymenoscyphus fraxineus, the causal agent of ash dieback, is a fungal pathogen that has been moving across continents and hosts from Asian to European ash. Most European common ash trees (Fraxinus excelsior) are highly susceptible to H. fraxineus, although a minority (~5%) have partial resistance to dieback. Here, we assemble and annotate a H. fraxineus draft genome, which approaches chromosome scale. Pathogen genetic diversity across Europe and in Japan, reveals a strong bottleneck in Europe, though a signal of adaptive diversity remains in key host interaction genes. We find that the European population was founded by two divergent haploid individuals. Divergence between these haplotypes represents the ancestral polymorphism within a large source population. Subsequent introduction from this source would greatly increase adaptive potential of the pathogen. Thus, further introgression of H. fraxineus into Europe represents a potential threat and Europe-wide biological security measures are needed to manage this disease.

Accelerating international trade and climate change make pathogen spread an increasing concern. Hymenoscyphus fraxineus, the causal agent of ash dieback is one such pathogen, moving across continents and hosts from Asian to European ash. Most European common ash (Fraxinus excelsior) trees are highly susceptible to H. fraxineus although a small minority (~5%) evidently have partial resistance to dieback. We have assembled and annotated a draft of the H. fraxineus genome which approaches chromosome scale. Pathogen genetic diversity across Europe, and in Japan, reveals a tight bottleneck into Europe, though a signal of adaptive diversity remains in key host interaction genes (effectors). We find that the European population was founded by two divergent haploid individuals. Divergence between these haplotypes represents the 'shadow' of a large source population and subsequent introduction would greatly increase adaptive potential and the pathogen's threat. Thus, EU wide biological security measures remain an important part of the strategy to manage this disease.

In the duplication-degeneration-complementation model, duplicated gene-pairs undergo nonfunctionalization (loss from the genome), subfunctionalization (the functions of the ancestral gene are sub-divided between duplicate genes), or neofunctionalization (one of the duplicate genes acquires a new function). These processes occur by loss or gain of regulatory elements in gene promoters. Fatty acid-binding proteins (Fabp) belong to a multigene family composed of orthologous proteins that are highly conserved in sequence and function, but differ in their gene regulation. We previously reported that the zebrafish fabp1a, fabp1b.1, and fabp1b.2 promoters underwent subfunctionalization of PPAR responsiveness. Here, we describe the regulation at the duplicated zebrafish fabp7a/fabp7b, fabp10a/fabp10b and fabp11a/fabp11b gene promoters. Differential control at the duplicated fabp promoters was assessed by DNA sequence analysis, responsiveness to PPAR-isoform specific agonists and NF-κB p50 antagonists in zebrafish liver and intestine explant tissue, and in HEK293A cells transfected with fabp promoter-reporter constructs. Each zebrafish fabp gene displayed unique transcriptional regulation compared to its paralogous duplicate. This work provides a framework to account for the evolutionary trajectories that led to the high retention (57%) of duplicated fabp genes in the zebrafish genome compared to only ~3% of all duplicated genes in the zebrafish genome.

The Sequence Distance Graph (SDG) framework works with genome assembly graphs and raw data from paired, linked and long reads. It includes a simple deBruijn graph module, and can import graphs using the graphical fragment assembly (GFA) format. It also maps raw reads onto graphs, and provides a Python application programming interface (API) to navigate the graph, access the mapped and raw data and perform interactive or scripted analyses. Its complete workspace can be dumped to and loaded from disk, decoupling mapping from analysis and supporting multi-stage pipelines. We present the design and implementation of the framework, and example analyses scaffolding a short read graph with long reads, and navigating paths in a heterozygous graph for a simulated parent-offspring trio dataset. SDG is freely available under the MIT license at https://github.com/bioinfologics/sdg

Darwin's studies on heterostyly in Primula described two floral morphs, pin and thrum, with reciprocal anther and stigma heights that promote insect-mediated cross-pollination. This key innovation evolved independently in several angiosperm families. Subsequent studies on heterostyly in Primula contributed to the foundation of modern genetic theory and the neo-Darwinian synthesis. The established genetic model for Primula heterostyly involves a diallelic S locus comprising several genes, with rare recombination events that result in self-fertile homostyle flowers with anthers and stigma at the same height. Here we reveal the S locus supergene as a tightly linked cluster of thrum-specific genes that are absent in pins. We show that thrums are hemizygous not heterozygous for the S locus, which suggests that homostyles do not arise by recombination between S locus haplotypes as previously proposed. Duplication of a floral homeotic gene 51.7 million years (Myr) ago, followed by its neofunctionalization, created the current S locus assemblage which led to floral heteromorphy in Primula. Our findings provide new insights into the structure, function and evolution of this archetypal supergene.

Increasing organismal complexity during the evolution of life has been attributed to the duplication of genes and entire genomes. More recently, theoretical models have been proposed that postulate the fate of duplicated genes, among them the duplication-degeneration-complementation (DDC) model. In the DDC model, the common fate of a duplicated gene is lost from the genome owing to nonfunctionalization. Duplicated genes are retained in the genome either by subfunctionalization, where the functions of the ancestral gene are sub-divided between the sister duplicate genes, or by neofunctionalization, where one of the duplicate genes acquires a new function. Both processes occur either by loss or gain of regulatory elements in the promoters of duplicated genes. Here, we review the genomic organization, evolution, and transcriptional regulation of the multigene family of intracellular lipid-binding protein (iLBP) genes from teleost fishes. Teleost fishes possess many copies of iLBP genes owing to a whole genome duplication (WGD) early in the teleost fish radiation. Moreover, the retention of duplicated iLBP genes is substantially higher than the retention of all other genes duplicated in the teleost genome. The fatty acid-binding protein genes, a subfamily of the iLBP multigene family in zebrafish, are differentially regulated by peroxisome proliferator-activated receptor (PPAR) isoforms, which may account for the retention of iLBP genes in the zebrafish genome by the process of subfunctionalization of cis-acting regulatory elements in iLBP gene promoters.

MotivationDe novo assembly of whole genome shotgun (WGS) next-generation sequencing (NGS) data benefits from high-quality input with high coverage. However, in practice, determining the quality and quantity of useful reads quickly and in a reference-free manner is not trivial. Gaining a better understanding of the WGS data, and how that data is utilised by assemblers, provides useful insights that can inform the assembly process and result in better assemblies.ResultsWe present the K-mer Analysis Toolkit (KAT): a multi-purpose software toolkit for reference-free quality control (QC) of WGS reads and de novo genome assemblies, primarily via their k-mer frequencies and GC composition. KAT enables users to assess levels of errors, bias and contamination at various stages of the assembly process. In this paper we highlight KAT’s ability to provide valuable insights into assembly composition and quality of genome assemblies through pairwise comparison of k-mers present in both input reads and the assemblies.AvailabilityKAT is available under the GPLv3 license at: https://github.com/TGAC/KAT.Contactbernardo.clavijo@earlham.ac.ukSupplementary InformationSupplementary Information (SI) is available at Bioinformatics online. In addition, the software documentation is available online at: http://kat.readthedocs.io/en/latest/.

Background In the duplication-degeneration-complementation (DDC) model, a duplicated gene has three possible fates: it may lose functionality through the accumulation of mutations (nonfunctionalization), acquire a new function (neofunctionalization), or each duplicate gene may retain a subset of functions of the ancestral gene (subfunctionalization). The role that promoter evolution plays in retention of duplicated genes in eukaryotic genomes is not well understood. Fatty acid-binding proteins (Fabp) belong to a multigene family that are highly conserved in sequence and function, but differ in their gene regulation, suggesting selective pressure is exerted via regulatory elements in the promoter. Results In this study, we describe the PPAR regulation of zebrafish fabp1a, fabp1b.1, and fabp1b.2 promoters and compare them to the PPAR regulation of the spotted gar fabp1 promoter, representative of the ancestral fabp1 gene. Evolution of the fabp1 promoter was inferred by sequence analysis, and differential PPAR-agonist activation of fabp1 promoter activity in zebrafish liver and intestine explant cells, and in HEK293A cells transiently transfected with wild-type and mutated fabp1promoter-reporter gene constructs. The promoter activity of spotted gar fabp1, representative of the ancestral fabp1, was induced by both PPARα- and PPARγ-specific agonists, but displayed a biphasic response to PPARα activation. Zebrafish fabp1a was PPARα-selective, fabp1b.1 was PPARγ-selective, and fabp1b.2 was not regulated by PPAR. Conclusions The zebrafish fabp1 promoters underwent two successive rounds of subfunctionalization with respect to PPAR regulation leading to retention of three zebrafish fabp1 genes with stimuli-specific regulation. Using a pharmacological approach, we demonstrated here the divergent regulation of the zebrafish fabp1a, fabp1b.1, and fabp1b.2 with regard to subfunctionalization of PPAR regulation following two rounds of gene duplication. Electronic supplementary material The online version of this article (doi:10.1186/s12862-016-0717-x) contains supplementary material, which is available to authorized users.

Gene duplication is thought to facilitate increasing complexity in the evolution of life. The fate of most duplicated genes is nonfunctionalization: functional decay resulting from the accumulation of mutations. According to the duplication–degeneration–complementation (DDC) model, duplicated genes are retained by subfunctionalization, where the functions of the ancestral gene are sub-divided between duplicate genes, or by neofunctionalization, where one of the duplicates acquires a new function. Here, we report the differential regulation of the zebrafish tandemly duplicated fatty acid-binding protein genes, fabp1b.1 and fabp1b.2, by peroxisome proliferator-activated receptors (PPAR). fabp1b.1 mRNA levels were induced in tissue explants of liver, but not intestine, by PPAR agonists. fabp1b.1 promoter activity was induced to a greater extent by rosiglitazone (PPARγ-selective agonist) compared to WY 14,643 (PPARα-selective agonist) in HEK293A cells. Mutation of a peroxisome proliferator response element (PPRE) at –1232 bp in the fabp1b.1 promoter reduced PPAR-dependent activation. fabp1b.2 promoter activity was not affected by PPAR agonists. Differential regulation of the duplicated fabp1b promoters may be the result of PPRE loss in fabp1b.2 during a meiotic crossing-over event. Retention of PPAR inducibility in fabp1b.1 and not fabp1b.2 suggests unique regulation and function of the fabp1b duplicates.

Aim: To review and compare methods for indirect comparison of aflibercept and ranibizumab in patients with diabetic macular edema. Methods: Post-stratification, inverse probability weighting based on simulated data, weight optimization, and regression model techniques were used to compare pooled individual patient-level data from the RESTORE and RESPOND (ranibizumab 0.5 mg as needed after 3 initial monthly doses) studies with summary-level data from the VIVID and VISTA (aflibercept 2.0 mg every 8 weeks after 5 initial monthly doses, 2q8) studies. The impact of adjusting for up to two baseline characteristics was assessed. Results: All methods provided similar results. After adjustment for baseline best-corrected visual acuity and central retinal thickness, no statistically significant difference in average gain in baseline best-corrected visual acuity from baseline at month 12 was found between ranibizumab 0.5 mg and aflibercept 2q8. Conclusions: Weight optimization and regression methods are useful options to adjust for more than one baseline characteristic.