Ab initio gene finding in Drosophila genomic DNA

The wild sweetpotato (Ipomoea trifida) genome provides insights into storage root development

BACKGROUND

Sweetpotato (Ipomoea batatas (L.) Lam.) is the seventh most important crop in the world and is mainly cultivated for its underground storage root (SR). The genetic studies of this species have been hindered by a lack of high-quality reference sequence due to its complex genome structure. Diploid Ipomoea trifida is the closest relative and putative progenitor of sweetpotato, which is considered a model species for sweetpotato, including genetic, cytological, and physiological analyses.

RESULTS

Here, we generated the chromosome-scale genome sequence of SR-forming diploid I. trifida var. Y22 with high heterozygosity (2.20%). Although the chromosome-based synteny analysis revealed that the I. trifida shared conserved karyotype with Ipomoea nil after the separation, I. trifida had a much smaller genome than I. nil due to more efficient eliminations of LTR-retrotransposons and lack of species-specific amplification bursts of LTR-RTs. A comparison with four non-SR-forming species showed that the evolution of the beta-amylase gene family may be related to SR formation. We further investigated the relationship of the key gene BMY11 (with identity 47.12% to beta-amylase 1) with this important agronomic trait by both gene expression profiling and quantitative trait locus (QTL) mapping. And combining SR morphology and structure, gene expression profiling and qPCR results, we deduced that the products of the activity of BMY11 in splitting starch granules and be recycled to synthesize larger granules, contributing to starch accumulation and SR swelling. Moreover, we found the expression pattern of BMY11, sporamin proteins and the key genes involved in carbohydrate metabolism and stele lignification were similar to that of sweetpotato during the SR development.

CONCLUSIONS

We constructed the high-quality genome reference of the highly heterozygous I. trifida through a combined approach and this genome enables a better resolution of the genomics feature and genome evolutions of this species. Sweetpotato SR development genes can be identified in I. trifida and these genes perform similar functions and patterns, showed that the diploid I. trifida var. Y22 with typical SR could be considered an ideal model for the studies of sweetpotato SR development.

Genome sequence of the non-conventional wine yeast Hanseniaspora guilliermondii UTAD222 unveils relevant traits of this species and of the Hanseniaspora genus in the context of wine fermentation

Hanseanispora species, including H. guilliermondii, are long known to be abundant in wine grape-musts and to play a critical role in vinification by modulating, among other aspects, the wine sensory profile. Despite this, the genetics and physiology of Hanseniaspora species remains poorly understood. The first genomic sequence of a H. guilliermondii strain (UTAD222) and the discussion of its potential significance are presented in this work. Metabolic reconstruction revealed that H. guilliermondii is not equipped with a functional gluconeogenesis or glyoxylate cycle, nor does it harbours key enzymes for glycerol or galactose catabolism or for biosynthesis of biotin and thiamine. Also, no fructose-specific transporter could also be predicted from the analysis of H. guilliermondii genome leaving open the mechanisms underlying the fructophilic character of this yeast. Comparative analysis involving H. guilliermondii, H. uvarum, H. opuntiae and S. cerevisiae revealed 14 H. guilliermondii-specific genes (including five viral proteins and one β-glucosidase). Furthermore, 870 proteins were only found within the Hanseniaspora proteomes including several β-glucosidases and decarboxylases required for catabolism of biogenic amines. The release of H. guilliermondii genomic sequence and the comparative genomics/proteomics analyses performed, is expected to accelerate research focused on Hanseniaspora species and to broaden their application in the wine industry and in other bio-industries in which they could be explored as cell factories.

Identification, characterization and field testing of Brassica napus mutants producing high-oleic oils

Producing healthy, high-oleic oils and eliminating trans-fatty acids from foods are two goals that can be addressed by reducing activity of the oleate desaturase, FAD2, in oilseeds. However, it is essential to understand the consequences of reducing FAD2 activity on the metabolism, cell biology and physiology of oilseed crop plants. Here, we translate knowledge from studies of fad2 mutants in Arabidopsis (Arabidopsis thaliana) to investigate the limits of non-GMO approaches to maximize oleic acid in the seed oil of canola (Brassica napus), a species that expresses three active FAD2 isozymes. A series of hypomorphic and null mutations in the FAD2.A5 isoform were characterized in yeast (Saccharomyes cerevisiae). Then, four of these were combined with null mutations in the other two isozymes, FAD2.C5 and FAD2.C1. The resulting mutant lines contained 71-87% oleic acid in their seed oil, compared with 62% in wild-type controls. All the mutant lines grew well in a greenhouse, but in field experiments we observed a clear demarcation in plant performance. Mutant lines containing less than 80% oleate in the seed oil were indistinguishable from wild-type controls in growth parameters and seed oil content. By contrast, lines with more than 80% oleate in the seed oil had significantly lower seedling establishment and vigor, delayed flowering and reduced plant height at maturity. These lines also had 7-11% reductions in seed oil content. Our results extend understanding of the B. napusFAD2 isozymes and define the practical limit to increasing oil oleate content in this crop species.

Chromosome level comparative analysis of Brassica genomes

We provided a chromosome-length assembly of B. nigra and show the comprehensive chromosome-scale variations among Brassica genomes. Chromosome-level assembly of the Brassica species, which include many important crops, is essential for the agricultural and evolutionary studies. While the present B. nigra chromosomes was connected with genetic map of B. juncea, hindering the comparative analysis of the B chromosomes. Here we present a chromosome-length B. nigra assembly constructed with Hi-C connections and its variations on chromosome level compared with other Brassica species. We produced an assembly of 484 Mb annotated with 51,829 genes, of which 393 Mb were anchored onto 8 chromosomes, taking 81.26% of the assembly. Comparison of the B chromosomes shows high concordance of the two B. nigra assemblies and reveals comprehensive variations of the B chromosomes after polyploidization and gene loss in syntenic regions. Chromosome blocks with variations have lower gene density and higher TE content. Furthermore, we compared the chromosomes of the three major Brassica diploids, which showed that most of the variations between B and A/C had completed before A/C divergence and there are more variations on C chromosomes after their divergence. In summary, our work presents a chromosome-length assembly of B. nigra and comprehensive comparative analysis of the Brassica chromosomes, which provides a useful reference for other studies and comprehensive information of Brassica chromosome evolution.

Chromosome conformation capture resolved near complete genome assembly of broomcorn millet

Broomcorn millet (Panicum miliaceum L.) has strong tolerance to abiotic stresses, and is probably one of the oldest crops, with its earliest cultivation that dated back to ca. ~10,000 years. We report here its genome assembly through a combination of PacBio sequencing, BioNano, and Hi-C (in vivo) mapping. The 18 super scaffolds cover ~95.6% of the estimated genome (~887.8 Mb). There are 63,671 protein-coding genes annotated in this tetraploid genome. About ~86.2% of the syntenic genes in foxtail millet have two homologous copies in broomcorn millet, indicating rare gene loss after tetraploidization in broomcorn millet. Phylogenetic analysis reveals that broomcorn millet and foxtail millet diverged around ~13.1 Million years ago (Mya), while the lineage specific tetraploidization of broomcorn millet may be happened within ~5.91 million years. The genome is not only beneficial for the genome assisted breeding of broomcorn millet, but also an important resource for other Panicum species.

Broomcorn millet is one of the oldest crops cultivated by human that has strong abiotic stress tolerance. To facilitate genome assisted breeding of this and related species, the authors report its genome assembly and conduct comparative genome structure and evolution analyses with foxtail millet.

Metal tolerance protein MTP6 affects mitochondrial iron and manganese homeostasis in cucumber

Members of the cation diffusion facilitator (CDF) family have been identified in all kingdoms of life. They have been divided into three subgroups, namely Zn-CDF, Fe/Zn-CDF, and Mn-CDF, based on their putative specificity to transported metal ions. The plant metal tolerance protein 6 (MTP6) proteins fall into the Fe/Zn-CDF subgroup; however, their function in iron/zinc transport has not yet been confirmed. Here, we characterized the MTP6 protein from cucumber, Cucumis sativus. When expressed in yeast and in protoplasts isolated from Arabidopsis cells, CsMTP6 localized in mitochondria and contributed to the efflux of Fe and Mn from these organelles. Immunolocalization of CsMTP6 in cucumber membranes confirmed this association with mitochondria. Root expression and protein levels of CsMTP6 were significantly up-regulated in conditions of Fe deficiency and excess, but were not affected by Mn availability. These results indicate that MTP6 proteins contribute to the distribution of Fe and Mn between the cytosol and mitochondria of plant cells, and are regulated by Fe to maintain mitochondrial and cytosolic iron homeostasis under varying conditions of Fe availability.

The genomic landscape of molecular responses to natural drought stress in Panicum hallii

Environmental stress is a major driver of ecological community dynamics and agricultural productivity. This is especially true for soil water availability, because drought is the greatest abiotic inhibitor of worldwide crop yields. Here, we test the genetic basis of drought responses in the genetic model for C4 perennial grasses, Panicum hallii, through population genomics, field-scale gene-expression (eQTL) analysis, and comparison of two complete genomes. While gene expression networks are dominated by local cis-regulatory elements, we observe three genomic hotspots of unlinked trans-regulatory loci. These regulatory hubs are four times more drought responsive than the genome-wide average. Additionally, cis- and trans-regulatory networks are more likely to have opposing effects than expected under neutral evolution, supporting a strong influence of compensatory evolution and stabilizing selection. These results implicate trans-regulatory evolution as a driver of drought responses and demonstrate the potential for crop improvement in drought-prone regions through modification of gene regulatory networks.

A New Reference Genome Shows the One-Speed Genome Structure of the Barley Pathogen Ramularia collo-cygni

Ramularia leaf spot has recently emerged as a major threat to barley production world-wide, causing 25% yield loss in many barley growing regions. Here, we provide a new reference genome of the causal agent, the Dothideomycete Ramularia collo-cygni. The assembly of 32 Mb consists of 78 scaffolds. We used RNA-seq to identify 11,622 genes of which 1,303 and 282 are coding for predicted secreted proteins and putative effectors respectively.

The pathogen separated from its nearest sequenced relative, Zymoseptoria tritici ∼27 Ma. We calculated the divergence of the two species on protein level and see remarkably high synonymous and nonsynonymous divergence. Unlike in many other plant pathogens, the comparisons of transposable elements and gene distributions, show a very homogeneous genome for R. collo-cygni. We see no evidence for higher selective pressure on putative effectors or other secreted proteins and repetitive sequences are spread evenly across the scaffolds. These findings could be associated to the predominantly endophytic life-style of the pathogen. We hypothesize that R. collo-cygni only recently became pathogenic and that therefore its genome does not yet show the typical pathogen characteristics. Because of its high scaffold length and improved CDS annotations, our new reference sequence provides a valuable resource for the community for future comparative genomics and population genetics studies.

Phylogenetic and mutational analyses of human LEUTX, a homeobox gene implicated in embryogenesis

Recently, human PAIRED-LIKE homeobox transcription factor (TF) genes were discovered whose expression is limited to the period of embryo genome activation up to the 8-cell stage. One of these TFs is LEUTX, but its importance for human embryogenesis is still subject to debate. We confirmed that human LEUTX acts as a TAATCC-targeting transcriptional activator, like other K50-type PAIRED-LIKE TFs. Phylogenetic comparisons revealed that Leutx proteins are conserved across Placentalia and comprise two conserved domains, the homeodomain, and a Leutx-specific domain containing putative transcriptional activation motifs (9aaTAD). Examination of human genotype resources revealed 116 allelic variants in LEUTX. Twenty-four variants potentially affect function, but they occur only heterozygously at low frequency. One variant affects a DNA-specificity determining residue, mutationally reachable by a one-base transition. In vitro and in silico experiments showed that this LEUTX mutation (alanine to valine at position 54 in the homeodomain) results in a transactivational loss-of-function to a minimal TAATCC-containing promoter and a 36 bp motif enriched in genes involved in embryo genome activation. A compensatory change in residue 47 restores function. The results support the notion that human LEUTX functions as a transcriptional activator important for human embryogenesis.

Genome assembly of the Pink Ipê (Handroanthus impetiginosus, Bignoniaceae), a highly valued, ecologically keystone Neotropical timber forest tree

Background

Handroanthus impetiginosus (Mart. ex DC.) Mattos is a keystone Neotropical hardwood tree widely distributed in seasonally dry tropical forests of South and Mesoamerica. Regarded as the “new mahogany,” it is the second most expensive timber, the most logged species in Brazil, and currently under significant illegal trading pressure. The plant produces large amounts of quinoids, specialized metabolites with documented antitumorous and antibiotic effects. The development of genomic resources is needed to better understand and conserve the diversity of the species, to empower forensic identification of the origin of timber, and to identify genes for important metabolic compounds.

Findings

The genome assembly covers 503.7 Mb (N50 = 81 316 bp), 90.4% of the 557-Mbp genome, with 13 206 scaffolds. A repeat database with 1508 sequences was developed, allowing masking of ∼31% of the assembly. Depth of coverage indicated that consensus determination adequately removed haplotypes assembled separately due to the extensive heterozygosity of the species. Automatic gene prediction provided 31 688 structures and 35 479 messenger RNA transcripts, while external evidence supported a well-curated set of 28 603 high-confidence models (90% of total). Finally, we used the genomic sequence and the comprehensive gene content annotation to identify genes related to the production of specialized metabolites.

Conclusions

This genome assembly is the first well-curated resource for a Neotropical forest tree and the first one for a member of the Bignoniaceae family, opening exceptional opportunities to empower molecular, phytochemical, and breeding studies. This work should inspire the development of similar genomic resources for the largely neglected forest trees of the mega-diverse tropical biomes.

Chromosome-scale assembly of the Monopterus genome

Background

The teleost fish Monopterus albus is emerging as a new model for biological studies due to its natural sex transition and small genome, in addition to its enormous economic and potential medical value. However, no genomic information for the Monopterus is currently available.

Findings

Here, we sequenced and de novo assembled the genome of M. albus and report the de novochromosome assembly by FISH walking assisted by conserved synteny (Cafs). Using Cafs, 328 scaffolds were assembled into 12 chromosomes, which covered genomic sequences of 555 Mb, accounting for 81.3% of the sequences assembled in scaffolds (∼689 Mb). A total of 18 ,660 genes were mapped on the chromosomes and showed a nonrandom distribution along chromosomes.

Conclusions

We report the first reference genome of the Monopterus and provide an efficient Cafs strategy for a de novo chromosome-level assembly of the Monopterus genome, which provides a valuable resource, not only for further studies in genetics, evolution, and development, particularly sex determination, but also for breed improvement of the species.

Heat stress induced alternative splicing in catfish as determined by transcriptome analysis

Heat tolerance is increasingly becoming an important trait for aquaculture species with a changing climate. Transcriptional studies on responses to heat stress have been conducted in catfish, one of the most important economic aquaculture species around the world. The molecular mechanisms underlying heat tolerance is still poorly understood, especially at the post-transcriptional level including regulation of alternative splicing. In this study, existing RNA-Seq datasets were utilized to characterize the change of alternative splicing in catfish following heat treatment. Heat-tolerant and -intolerant catfish were differentiated by the time to lost equilibrium after heat stress. With heat stress, alternative splicing was generally increased. In heat-intolerant fish, the thermal stress induced 29.2% increases in alternative splicing events and 25.8% increases in alternatively spliced genes. A total of 282, 189, and 44 differential alternative splicing (DAS) events were identified in control-intolerant, control-tolerant, and intolerant-tolerant comparisons, corresponding to 252, 171, and 42 genes, respectively. Gene ontology analyses showed that genes involved in the molecular function of RNA binding were significantly enriched in DAS gene sets after heat stress in both heat-intolerant and -tolerant catfish compared with the control group. Similar results were also observed in the DAS genes between heat-intolerant and -tolerant catfish, and the biological process of RNA splicing was also enriched in this comparison, indicating the involvement of RNA splicing-related genes underlying heat tolerance. This is the first comprehensive study of alternative splicing in response to heat stress in fish species, providing insights into the molecular mechanisms of responses to the abiotic stress.

Genome Sequence of the Wine Yeast Saccharomycodes ludwigii UTAD17

This work describes, for the first time, the genome sequence of a Saccharomycodes ludwigii strain. Although usually seen as a wine spoilage yeast, S. ludwigii has been of interest for the production of fermented beverages because it harbors several interesting properties, including the production of beneficial aroma compounds.

This work describes, for the first time, the genome sequence of a Saccharomycodes ludwigii strain. Although usually seen as a wine spoilage yeast, S. ludwigii has been of interest for the production of fermented beverages because it harbors several interesting properties, including the production of beneficial aroma compounds.

Chromosome-level genome assembly of the spotted sea bass, Lateolabrax maculatus

Background

The spotted sea bass (Lateolabrax maculatus) is a valuable commercial fish that is widely cultured in China. While analyses using molecular markers and population genetics have been conducted, genomic resources are lacking.

Findings

Here, we report a chromosome-scale assembly of the spotted sea bass genome by high-depth genome sequencing, assembly, and annotation. The genome scale was 0.67 Gb with contig and scaffold N50 length of 31 Kb and 1,040 Kb, respectively. Hi-C scaffolding of the genome resulted in 24 pseudochromosomes containing 77.68% of the total assembled sequences. A total of 132.38 Mb repeat sequences were detected, accounting for 20.73% of the assembled genome. A total of 22, 015 protein-coding genes were predicted, of which 96.52% were homologous to proteins in various databases. In addition, we constructed a phylogenetic tree using 1,586 single-copy gene families and identified 125 unique gene families in the spotted sea bass genome.

Conclusions

We assembled a spotted sea bass genome that will be a valuable genomic resource to understanding the biology of the spotted sea bass and will also lead to the development of molecular breeding techniques to generate spotted sea bass with better economic traits.

Analysis of the Trichuris suis excretory/secretory proteins as a function of life cycle stage and their immunomodulatory properties

Parasitic worms have a remarkable ability to modulate host immune responses through several mechanisms including excreted/secreted proteins (ESP), yet the exact nature of these proteins and their targets often remains elusive. Here, we performed mass spectrometry analyses of ESP (TsESP) from larval and adult stages of the pig whipworm Trichuris suis (Ts) and identified ~350 proteins. Transcriptomic analyses revealed large subsets of differentially expressed genes in the various life cycle stages of the parasite. Exposure of bone marrow-derived macrophages and dendritic cells to TsESP markedly diminished secretion of the pro-inflammatory cytokines TNFα and IL-12p70. Conversely, TsESP exposure strongly induced release of the anti-inflammatory cytokine IL-10, and also induced high levels of nitric oxide (NO) and upregulated arginase activity in macrophages. Interestingly, TsESP failed to directly induce CD4+ CD25+ FoxP3+ regulatory T cells (Treg cells), while OVA-pulsed TsESP-treated dendritic cells suppressed antigen-specific OT-II CD4+ T cell proliferation. Fractionation of TsESP identified a subset of proteins that promoted anti-inflammatory functions, an activity that was recapitulated using recombinant T. suis triosephosphate isomerase (TPI) and nucleoside diphosphate kinase (NDK). Our study helps illuminate the intricate balance that is characteristic of parasite-host interactions at the immunological interface, and further establishes the principle that specific parasite-derived proteins can modulate immune cell functions.

The Aquilegia genome provides insight into adaptive radiation and reveals an extraordinarily polymorphic chromosome with a unique history

The columbine genus Aquilegia is a classic example of an adaptive radiation, involving a wide variety of pollinators and habitats. Here we present the genome assembly of A. coerulea 'Goldsmith', complemented by high-coverage sequencing data from 10 wild species covering the world-wide distribution. Our analyses reveal extensive allele sharing among species and demonstrate that introgression and selection played a role in the Aquilegia radiation. We also present the remarkable discovery that the evolutionary history of an entire chromosome differs from that of the rest of the genome - a phenomenon that we do not fully understand, but which highlights the need to consider chromosomes in an evolutionary context.

Modeling one thousand intron length distributions with fitild

Motivation

Intron length distribution (ILD) is a specific feature of a genome that exhibits extensive species-specific variation. Whereas ILD contributes to up to 30% of the total information content for intron recognition in some species, rendering it an important component of computational gene prediction, very few studies have been conducted to quantitatively characterize ILDs of various species.

Results

We developed a set of computer programs (fitild, compild, etc.) to build statistical models of ILDs and compare them with one another. Each ILD of more than 1000 genomes was fitted with fitild to a statistical model consisting of one, two, or three components of Frechet distributions. Several measures of distances between ILDs were calculated by compild. A theoretical model was presented to better understand the origin of the observed shape of an ILD.

Availability and implementation

The C++ source codes are available at https://github.com/ogotoh/fitild.git/.

Supplementary information

Supplementary data are available at Bioinformatics online.

Engineering Kluyveromyces marxianus as a Robust Synthetic Biology Platform Host

The yeast Kluyveromyces marxianus grows at high temperatures and on a wide range of carbon sources, making it a promising host for industrial biotechnology to produce renewable chemicals from plant biomass feedstocks. However, major genetic engineering limitations have kept this yeast from replacing the commonly used yeast Saccharomyces cerevisiae in industrial applications. Here, we describe genetic tools for genome editing and breeding K. marxianus strains, which we use to create a new thermotolerant strain with promising fatty acid production. These results open the door to using K. marxianus as a versatile synthetic biology platform organism for industrial applications.

Throughout history, the yeast Saccharomyces cerevisiae has played a central role in human society due to its use in food production and more recently as a major industrial and model microorganism, because of the many genetic and genomic tools available to probe its biology. However, S. cerevisiae has proven difficult to engineer to expand the carbon sources it can utilize, the products it can make, and the harsh conditions it can tolerate in industrial applications. Other yeasts that could solve many of these problems remain difficult to manipulate genetically. Here, we engineered the thermotolerant yeast Kluyveromyces marxianus to create a new synthetic biology platform. Using CRISPR-Cas9 (clustered regularly interspaced short palindromic repeats with Cas9)-mediated genome editing, we show that wild isolates of K. marxianus can be made heterothallic for sexual crossing. By breeding two of these mating-type engineered K. marxianus strains, we combined three complex traits—thermotolerance, lipid production, and facile transformation with exogenous DNA—into a single host. The ability to cross K. marxianus strains with relative ease, together with CRISPR-Cas9 genome editing, should enable engineering of K. marxianus isolates with promising lipid production at temperatures far exceeding those of other fungi under development for industrial applications. These results establish K. marxianus as a synthetic biology platform comparable to S. cerevisiae, with naturally more robust traits that hold potential for the industrial production of renewable chemicals.

The maize W22 genome provides a foundation for functional genomics and transposon biology

The maize W22 inbred has served as a platform for maize genetics since the mid twentieth century. To streamline maize genome analyses, we have sequenced and de novo assembled a W22 reference genome using short-read sequencing technologies. We show that significant structural heterogeneity exists in comparison to the B73 reference genome at multiple scales, from transposon composition and copy number variation to single-nucleotide polymorphisms. The generation of this reference genome enables accurate placement of thousands of Mutator (Mu) and Dissociation (Ds) transposable element insertions for reverse and forward genetics studies. Annotation of the genome has been achieved using RNA-seq analysis, differential nuclease sensitivity profiling and bisulfite sequencing to map open reading frames, open chromatin sites and DNA methylation profiles, respectively. Collectively, the resources developed here integrate W22 as a community reference genome for functional genomics and provide a foundation for the maize pan-genome.

Evolutionary and expression analyses reveal a pattern of ancient duplications and functional specializations in the diversification of the Downstream of Kinase (DOK) genes

Downstream of Kinase (DOK) proteins represent a multigenic family of adaptors that includes negative regulators of immune cell signaling. Using phylogenetics and intron/exon structure data, we show here that the seven human DOK genes (DOK1 to DOK7) form three highly divergent groups that emerged before the protostome-deuterostome split: DOK1/2/3, DOK4/5/6, and DOK7. For two of these three groups (DOK1/2/3 and DOK4/5/6), further gene duplications occurred in vertebrates and so while chordates only have three DOK genes, vertebrates have seven DOK genes over the three groups. From our expression analysis in humans, we show that each group of DOK genes has a distinct pattern of expression. The DOK1/2/3 group is immune specific, yet each of the three genes in the group has a distinct pattern of expression in immune cells. This immune specificity could thus be ancestral, with the DOK1/2/3 gene also being immune-related in protostomes. The DOK4/5/6 and DOK7 groups represent genes that are much less expressed in immune system than the DOK1/2/3 group. Interestingly, we identify a novel tyrosine based motif that is specific to the vertebrate DOK4/5/6 sequences. The evolution of the DOK genes is thus marked by a pattern of ancient duplications and functional specializations.

Aedes aegypti antiviral adaptive response against DENV-2

Priming is the conceptual term defining memory phenomenon in innate immune response. Numerous examples of enhanced secondary immune response have been described in diverse taxa of invertebrates; which naturally lacks memory response. In mosquitoes, a previous non-lethal challenge with some specific pathogens modify their immune response against the same microorganism; developing an improved antimicrobial reaction. In this work, we explore the ability of Aedes aegypti to mount a higher antiviral response upon a second oral DENV challenge. When previously challenged with inactive virus, we observed that the posterior infection showed a diminished number of DENV infectious particles in midguts and carcasses. In challenged tissues, we detected higher de novo midgut DNA synthesis than control group, as determined by DNA incorporation of 5-bromo-2-deoxyuridine. We demonstrated that inactive DENV particle are capable to induce DNA synthesis levels comparable to infective DENV. We considered the Drosophila melanogaster hindsight and Delta-Notch mosquitoes orthologues as potential de novo DNA synthesis pathway components (as observed in fly oocyte development and midgut tissue renewal). We showed that Aedes aegypti hindsight transcript relative expression levels were higher than control during DENV infection and inactive DENV particle alimentation. Also, Aedes aegypti second challenge with active DENV induced higher hindsight, Delta and Notch transcriptions in the primed mosquitoes (compared with the primary infection levels). Considering that the mosquito de novo DNA synthesis is concomitant to viral particle reduction, this finding opens a new perspective on the mechanisms underlying the vector antiviral immune response and the effector molecules involved.

De Novo DNA Synthesis in Aedes aegypti Midgut Cells as a Complementary Strategy to Limit Dengue Viral Replication

Aedes aegypti is the main vector of Dengue Virus, carrying the virus during the whole mosquito life post-infection. Few mosquito fitness costs have been associated to the virus infection, thereby allowing for a swift dissemination. In order to diminish the mosquito population, public health agency use persistent chemicals with environmental impact for disease control. Most countries barely use biological controls, if at all. With the purpose of developing novel Dengue control strategies, a detailed understanding of the unexplored virus-vector interactions is urgently needed. Damage induced (through tissue injury or bacterial invasion) DNA duplication (endoreplication) has been described in insects during epithelial cells renewal. Here, we delved into the mosquito midgut tissue ability to synthesize DNA de novo; postulating that Dengue virus infection could trigger a protective endoreplication mechanism in some mosquito cells. We hypothesized that the Aedes aegypti orthologue of the Drosophila melanogaster hindsight gene (not previously annotated in Aedes aegypti transcriptome/genome) is part of the Delta-Notch pathway. The activation of this transcriptional cascade leads to genomic DNA endoreplication. The amplification of the genomic copies of specific genes ultimately limits the viral spreading during infection. Conversely, inhibiting DNA synthesis capacity, hence endoreplication, leads to a higher viral replication.

Draft genome of the Peruvian scallop Argopecten purpuratus

Background

The Peruvian scallop, Argopecten purpuratus, is mainly cultured in southern Chile and Peru was introduced into China in the last century. Unlike other Argopecten scallops, the Peruvian scallop normally has a long life span of up to 7 to 10 years. Therefore, researchers have been using it to develop hybrid vigor. Here, we performed whole genome sequencing, assembly, and gene annotation of the Peruvian scallop, with an important aim to develop genomic resources for genetic breeding in scallops.

Findings

A total of 463.19-Gb raw DNA reads were sequenced. A draft genome assembly of 724.78 Mb was generated (accounting for 81.87% of the estimated genome size of 885.29 Mb), with a contig N50 size of 80.11 kb and a scaffold N50 size of 1.02 Mb. Repeat sequences were calculated to reach 33.74% of the whole genome, and 26,256 protein-coding genes and 3,057 noncoding RNAs were predicted from the assembly.

Conclusions

We generated a high-quality draft genome assembly of the Peruvian scallop, which will provide a solid resource for further genetic breeding and for the analysis of the evolutionary history of this economically important scallop.

Genome-wide sequencing of longan (Dimocarpus longan Lour.) provides insights into molecular basis of its polyphenol-rich characteristics

Longan (Dimocarpus longan Lour.), an important subtropical fruit in the family Sapindaceae, is grown in more than 10 countries. Longan is an edible drupe fruit and a source of traditional medicine with polyphenol-rich traits. Tree size, alternate bearing, and witches' broom disease still pose serious problems. To gain insights into the genomic basis of longan traits, a draft genome sequence was assembled. The draft genome (about 471.88 Mb) of a Chinese longan cultivar, "Honghezi," was estimated to contain 31 007 genes and 261.88 Mb of repetitive sequences. No recent whole-genome-wide duplication event was detected in the genome. Whole-genome resequencing and analysis of 13 cultivated D. longan accessions revealed the extent of genetic diversity. Comparative transcriptome studies combined with genome-wide analysis revealed polyphenol-rich and pathogen resistance characteristics. Genes involved in secondary metabolism, especially those from significantly expanded (DHS, SDH, F3΄H, ANR, and UFGT) and contracted (PAL, CHS, and F3΄5΄H) gene families with tissue-specific expression, may be important contributors to the high accumulation levels of polyphenolic compounds observed in longan fruit. The high number of genes encoding nucleotide-binding site leucine-rich repeat (NBS-LRR) and leucine-rich repeat receptor-like kinase proteins, as well as the recent expansion and contraction of the NBS-LRR family, suggested a genomic basis for resistance to insects, fungus, and bacteria in this fruit tree. These data provide insights into the evolution and diversity of the longan genome. The comparative genomic and transcriptome analyses provided information about longan-specific traits, particularly genes involved in its polyphenol-rich and pathogen resistance characteristics.

Improving eukaryotic genome annotation using single molecule mRNA sequencing

BACKGROUND

The advantages of Pacific Biosciences (PacBio) single-molecule real-time (SMRT) technology include long reads, low systematic bias, and high consensus read accuracy. Here we use these attributes to improve on the genome annotation of the parasitic hookworm Ancylostoma ceylanicum using PacBio RNA-Seq.

RESULTS

We sequenced 192,888 circular consensus sequences (CCS) derived from cDNAs generated using the CloneTech SMARTer system. These SMARTer-SMRT libraries were normalized and size-selected providing a robust population of expressed structural genes for subsequent genome annotation. We demonstrate PacBio mRNA sequences based genome annotation improvement, compared to genome annotation using conventional sequencing-by-synthesis alone, by identifying 1609 (9.2%) new genes, extended the length of 3965 (26.7%) genes and increased the total genomic exon length by 1.9 Mb (12.4%). Non-coding sequence representation (primarily from UTRs based on dT reverse transcription priming) was particularly improved, increasing in total length by fifteen-fold, by increasing both the length and number of UTR exons. In addition, the UTR data provided by these CCS allowed for the identification of a novel SL2 splice leader sequence for A. ceylanicum and an increase in the number and proportion of functionally annotated genes. RNA-seq data also confirmed some of the newly annotated genes and gene features.

CONCLUSION

Overall, PacBio data has supported a significant improvement in gene annotation in this genome, and is an appealing alternative or complementary technique for genome annotation to the other transcript sequencing technologies.

Whole genome duplications have provided teleosts with many roads to peptide loaded MHC class I molecules

Background

In sharks, chickens, rats, frogs, medaka and zebrafish there is haplotypic variation in MHC class I and closely linked genes involved in antigen processing, peptide translocation and peptide loading. At least in chicken, such MHCIa haplotypes of MHCIa, TAP2 and Tapasin are shown to influence the repertoire of pathogen epitopes being presented to CD8+ T-cells with subsequent effect on cell-mediated immune responses.

Results

Examining MHCI haplotype variation in Atlantic salmon using transcriptome and genome resources we found little evidence for polymorphism in antigen processing genes closely linked to the classical MHCIa genes. Looking at other genes involved in MHCI assembly and antigen processing we found retention of functional gene duplicates originating from the second vertebrate genome duplication event providing cyprinids, salmonids, and neoteleosts with the potential of several different peptide-loading complexes. One of these gene duplications has also been retained in the tetrapod lineage with orthologs in frogs, birds and opossum.

Conclusion

We postulate that the unique salmonid whole genome duplication (SGD) is responsible for eliminating haplotypic content in the paralog MHCIa regions possibly due to frequent recombination and reorganization events at early stages after the SGD. In return, multiple rounds of whole genome duplications has provided Atlantic salmon, other teleosts and even lower vertebrates with alternative peptide loading complexes. How this affects antigen presentation remains to be established.

Electronic supplementary material

The online version of this article (10.1186/s12862-018-1138-9) contains supplementary material, which is available to authorized users.

Genome sequencing and comparative genomics reveal the potential pathogenic mechanism of Cercospora sojina Hara on soybean

Frogeye leaf spot, caused by Cercospora sojina Hara, is a common disease of soybean in most soybean-growing countries of the world. In this study, we report a high-quality genome sequence of C. sojina by Single Molecule Real-Time sequencing method. The 40.8-Mb genome encodes 11,655 predicated genes, and 8,474 genes are revealed by RNA sequencing. Cercospora sojina genome contains large numbers of gene clusters that are involved in synthesis of secondary metabolites, including mycotoxins and pigments. However, much less carbohydrate-binding module protein encoding genes are identified in C. sojina genome, when compared with other phytopathogenic fungi. Bioinformatics analysis reveals that C. sojina harbours about 752 secreted proteins, and 233 of them are effectors. During early infection, the genes for metabolite biosynthesis and effectors are significantly enriched, suggesting that they may play essential roles in pathogenicity. We further identify 13 effectors that can inhibit BAX-induced cell death. Taken together, our results provide insights into the infection mechanisms of C. sojina on soybean.

Pan-genome analysis highlights the extent of genomic variation in cultivated and wild rice

The rich genetic diversity in Oryza sativa and Oryza rufipogon serves as the main sources in rice breeding. Large-scale resequencing has been undertaken to discover allelic variants in rice, but much of the information for genetic variation is often lost by direct mapping of short sequence reads onto the O. sativa japonica Nipponbare reference genome. Here we constructed a pan-genome dataset of the O. sativa-O. rufipogon species complex through deep sequencing and de novo assembly of 66 divergent accessions. Intergenomic comparisons identified 23 million sequence variants in the rice genome. This catalog of sequence variations includes many known quantitative trait nucleotides and will be helpful in pinpointing new causal variants that underlie complex traits. In particular, we systemically investigated the whole set of coding genes using this pan-genome data, which revealed extensive presence and absence of variation among rice accessions. This pan-genome resource will further promote evolutionary and functional studies in rice.

The Medicago truncatula Genome: Genomic Data Availability

Medicago truncatula emerged in 1990 as a model for legumes, comprising the third largest land plant family. Most legumes form symbiotic nitrogen-fixing root nodules with compatible soil bacteria and thus are important contributors to the global nitrogen cycle and sustainable agriculture. Legumes and legume products are important sources for human and animal protein as well as for edible and industrial oils. In the years since M. truncatula was chosen as a legume model, many genetic, genomic, and molecular resources have become available, including reference quality genome sequences for two widely used genotypes. Accessibility of genomic data is important for many different types of studies with M. truncatula as well as for research involving crop and forage legumes. In this chapter, we discuss strategies to obtain archived M. truncatula genomic data originally deposited into custom databases that are no longer maintained but are now accessible in general databases. We also review key current genomic databases that are specific to M. truncatula as well as those that contain M. truncatula data in addition to data from other plants.

Manual Gene Curation and Functional Annotation

No genome sequencing project is complete without structural and functional annotation. Gene models and functional predictions for these models can be obtained relatively easily using computational methods, but they are prone to errors. We describe herein the steps we use to manually curate gene models and functionally annotate them. Our approach is to examine each gene model carefully, and improve its structure if necessary, using a comprehensive set of experimental and computational data as evidence. Then, functional predictions are assigned to the gene models based on conserved protein domains and sequence similarities. We use stringent sequence similarity cutoffs and reviewed sequence-database records as external sources for our annotations. By methodically choosing which evidence to use for each annotation, we minimize the risk of adopting and assigning false predictions to the gene models.

Contribution of NtZIP1-Like to the Regulation of Zn Homeostasis

Tobacco has frequently been suggested as a candidate plant species for use in phytoremediation of metal contaminated soil but knowledge on the regulation of its metal-homeostasis is still in the infancy. To identify new tobacco metal transport genes that are involved in Zn homeostasis a bioinformatics study using the tobacco genome information together with expression analysis was performed. Ten new tobacco metal transport genes from the ZIP, NRAMP, MTP, and MRP/ABCC families were identified with expression levels in leaves that were modified by exposure to Zn excess. Following exposure to high Zn there was upregulation of NtZIP11-like, NtNRAMP3, three isoforms of NtMTP2, three MRP/ABCC genes (NtMRP5-like, NtMRP10-like, and NtMRP14 like) and downregulation of NtZIP1-like and NtZIP4. This suggests their involvement in several processes governing the response to Zn-related stress and in the efficiency of Zn accumulation (uptake, sequestration, and redistribution). Further detailed analysis of NtZIP1-like provided evidence that it is localized at the plasma membrane and is involved in Zn but not Fe and Cd transport. NtZIP1-like is expressed in the roots and shoots, and is regulated developmentally and in a tissue-specific manner. It is highly upregulated by Zn deficiency in the leaves and the root basal region but not in the root apical zone (region of maturation and absorption containing root hairs). Thus NtZIP1-like is unlikely to be responsible for Zn uptake by the root apical region but rather in the uptake by root cells within the already mature basal zone. It is downregulated by Zn excess suggesting it is involved in a mechanism to protect the root and leaf cells from accumulating excess Zn.

Fungal Genome Annotation

The term "genome annotation" includes identification of protein-coding and noncoding sequences (e.g., repeats, rDNA, and ncRNA) in genome assemblies and attaching functional information (metadata) to these annotated features. Here, we describe the basic outline of fungal nuclear and mitochondrial genome annotation as performed at the US Department of Energy Joint Genome Institute (JGI).

Bacterial endosymbionts influence host sexuality and reveal reproductive genes of early divergent fungi

Many heritable mutualisms, in which beneficial symbionts are transmitted vertically between host generations, originate as antagonisms with parasite dispersal constrained by the host. Only after the parasite gains control over its transmission is the symbiosis expected to transition from antagonism to mutualism. Here, we explore this prediction in the mutualism between the fungus Rhizopus microsporus (Rm, Mucoromycotina) and a beta-proteobacterium Burkholderia, which controls host asexual reproduction. We show that reproductive addiction of Rm to endobacteria extends to mating, and is mediated by the symbiont gaining transcriptional control of the fungal ras2 gene, which encodes a GTPase central to fungal reproductive development. We also discover candidate G-protein-coupled receptors for the perception of trisporic acids, mating pheromones unique to Mucoromycotina. Our results demonstrate that regulating host asexual proliferation and modifying its sexual reproduction are sufficient for the symbiont's control of its own transmission, needed for antagonism-to-mutualism transition in heritable symbioses. These properties establish the Rm-Burkholderia symbiosis as a powerful system for identifying reproductive genes in Mucoromycotina.

A modified GC-specific MAKER gene annotation method reveals improved and novel gene predictions of high and low GC content in Oryza sativa

BACKGROUND

Accurate structural annotation depends on well-trained gene prediction programs. Training data for gene prediction programs are often chosen randomly from a subset of high-quality genes that ideally represent the variation found within a genome. One aspect of gene variation is GC content, which differs across species and is bimodal in grass genomes. When gene prediction programs are trained on a subset of grass genes with random GC content, they are effectively being trained on two classes of genes at once, and this can be expected to result in poor results when genes are predicted in new genome sequences.

RESULTS

We find that gene prediction programs trained on grass genes with random GC content do not completely predict all grass genes with extreme GC content. We show that gene prediction programs that are trained with grass genes with high or low GC content can make both better and unique gene predictions compared to gene prediction programs that are trained on genes with random GC content. By separately training gene prediction programs with genes from multiple GC ranges and using the programs within the MAKER genome annotation pipeline, we were able to improve the annotation of the Oryza sativa genome compared to using the standard MAKER annotation protocol. Gene structure was improved in over 13% of genes, and 651 novel genes were predicted by the GC-specific MAKER protocol.

CONCLUSIONS

We present a new GC-specific MAKER annotation protocol to predict new and improved gene models and assess the biological significance of this method in Oryza sativa. We expect that this protocol will also be beneficial for gene prediction in any organism with bimodal or other unusual gene GC content.

Nucleotide patterns aiding in prediction of eukaryotic promoters

Computational analysis of promoters is hindered by the complexity of their architecture. In less studied genomes with complex organization, false positive promoter predictions are common. Accurate identification of transcription start sites and core promoter regions remains an unsolved problem. In this paper, we present a comprehensive analysis of genomic features associated with promoters and show that probabilistic integrative algorithms-driven models allow accurate classification of DNA sequence into “promoters” and “non-promoters” even in absence of the full-length cDNA sequences. These models may be built upon the maps of the distributions of sequence polymorphisms, RNA sequencing reads on genomic DNA, methylated nucleotides, transcription factor binding sites, as well as relative frequencies of nucleotides and their combinations. Positional clustering of binding sites shows that the cells of Oryza sativa utilize three distinct classes of transcription factors: those that bind preferentially to the [-500,0] region (188 “promoter-specific” transcription factors), those that bind preferentially to the [0,500] region (282 “5′ UTR-specific” TFs), and 207 of the “promiscuous” transcription factors with little or no location preference with respect to TSS. For the most informative motifs, their positional preferences are conserved between dicots and monocots.

Architecture and Distribution of Introns in Core Genes of Four Fusarium Species

Removal of introns from transcribed RNA represents a crucial step during the production of mRNA in eukaryotes. Available whole-genome sequences and expressed sequence tags (ESTs) have increased our knowledge of this process and revealed various commonalities among eukaryotes. However, certain aspects of intron structure and diversity are taxon-specific, which can complicate the accuracy of in silico gene prediction methods. Using core genes, we evaluated the distribution and architecture of Fusarium circinatum spliceosomal introns, and linked these characteristics to the accuracy of the predicted gene models of the genome of this fungus. We also evaluated intron distribution and architecture in F. verticillioides, F. oxysporum, and F. graminearum, and made comparisons with F. circinatum. Results indicated that F. circinatum and the three other Fusarium species have canonical 5′ and 3′ splice sites, but with subtle differences that are apparently not shared with those of other fungal genera. The polypyrimidine tract of Fusarium introns was also found to be highly divergent among species and genes. Furthermore, the conserved adenosine nucleoside required during the first step of splicing is contained within unique branch site motifs in certain Fusarium introns. Data generated here show that introns of F. circinatum, as well as F. verticillioides, F. oxysporum, and F. graminearum, are characterized by a number of unique features such as the CTHAH and ACCAT motifs of the branch site. Incorporation of such information into genome annotation software will undoubtedly improve the accuracy of gene prediction methods used for Fusarium species and related fungi.

Describing Genomic and Epigenomic Traits Underpinning Emerging Fungal Pathogens

An unprecedented number of pathogenic fungi are emerging and causing disease in animals and plants, putting the resilience of wild and managed ecosystems in jeopardy. While the past decades have seen an increase in the number of pathogenic fungi, they have also seen the birth of new big data technologies and analytical approaches to tackle these emerging pathogens. We review how the linked fields of genomics and epigenomics are transforming our ability to address the challenge of emerging fungal pathogens. We explore the methodologies and bioinformatic toolkits that currently exist to rapidly analyze the genomes of unknown fungi, then discuss how these data can be used to address key questions that shed light on their epidemiology. We show how genomic approaches are leading a revolution into our understanding of emerging fungal diseases and speculate on future approaches that will transform our ability to tackle this increasingly important class of emerging pathogens.

Comparative genomics of geographically distant Fusarium fujikuroi isolates revealed two distinct pathotypes correlating with secondary metabolite profiles

Fusarium fujikuroi causes bakanae ("foolish seedling") disease of rice which is characterized by hyper-elongation of seedlings resulting from production of gibberellic acids (GAs) by the fungus. This plant pathogen is also known for production of harmful mycotoxins, such as fusarins, fusaric acid, apicidin F and beauvericin. Recently, we generated the first de novo genome sequence of F. fujikuroi strain IMI 58289 combined with extensive transcriptional, epigenetic, proteomic and chemical product analyses. GA production was shown to provide a selective advantage during infection of the preferred host plant rice. Here, we provide genome sequences of eight additional F. fujikuroi isolates from distant geographic regions. The isolates differ in the size of chromosomes, most likely due to variability of subtelomeric regions, the type of asexual spores (microconidia and/or macroconidia), and the number and expression of secondary metabolite gene clusters. Whilst most of the isolates caused the typical bakanae symptoms, one isolate, B14, caused stunting and early withering of infected seedlings. In contrast to the other isolates, B14 produced no GAs but high amounts of fumonisins during infection on rice. Furthermore, it differed from the other isolates by the presence of three additional polyketide synthase (PKS) genes (PKS40, PKS43, PKS51) and the absence of the F. fujikuroi-specific apicidin F (NRPS31) gene cluster. Analysis of additional field isolates confirmed the strong correlation between the pathotype (bakanae or stunting/withering), and the ability to produce either GAs or fumonisins. Deletion of the fumonisin and fusaric acid-specific PKS genes in B14 reduced the stunting/withering symptoms, whereas deletion of the PKS51 gene resulted in elevated symptom development. Phylogenetic analyses revealed two subclades of F. fujikuroi strains according to their pathotype and secondary metabolite profiles.

Computational deciphering of biotic stress associated genes in tomato (Solanum lycopersicum)

Tomato (Solanum lycopersicum) is one of the major vegetable plant and a model system for fruit development. Its global importance is due to its lycopene pigment which has anti-oxidative and anti-cancerous properties. Though > 1.5 M biotic stress associated ESTs of tomato are available but cumulative analysis to predict genes is warranted. Availability of whole genome de novo assembly can advantageously be used to map them over different chromosome. Further, available 0.14 M catalogued markers can be used to introgress specific desirable genes in varietal improvement program. We report here 57 novel genes associated with biotic stress of tomato along with 50 genes having physical location over different chromosomes. We also report 52 cis-regulating elements and 69 putative miRNAs which are involved in regulation of these biotic stresses associated genes. These putative candidate genes associated with biotic stress can be used in molecular breeding in the endeavor of tomato productivity along with its sustainable germplasm management.

•

ESts related to biotic stress were collected and assembled into contigs.

•

Total 57 novel genes were computationally mined from the assembled contigs.

•

Among the predicted novel gene, 50 genes were mapped on tomato genome.

•

52 cis-regulating elements and 69 putative miRNAs were predicted for the novel genes.

An optimized approach for annotation of large eukaryotic genomic sequences using genetic algorithm

BACKGROUND

Detection of important functional and/or structural elements and identification of their positions in a large eukaryotic genomic sequence are an active research area. Gene is an important functional and structural unit of DNA. The computation of gene prediction is, therefore, very essential for detailed genome annotation.

RESULTS

In this paper, we propose a new gene prediction technique based on Genetic Algorithm (GA) to determine the optimal positions of exons of a gene in a chromosome or genome. The correct identification of the coding and non-coding regions is difficult and computationally demanding. The proposed genetic-based method, named Gene Prediction with Genetic Algorithm (GPGA), reduces this problem by searching only one exon at a time instead of all exons along with its introns. This representation carries a significant advantage in that it breaks the entire gene-finding problem into a number of smaller sub-problems, thereby reducing the computational complexity. We tested the performance of the GPGA with existing benchmark datasets and compared the results with well-known and relevant techniques. The comparison shows the better or comparable performance of the proposed method. We also used GPGA for annotating the human chromosome 21 (HS21) using cross-species comparisons with the mouse orthologs.

CONCLUSION

It was noted that the GPGA predicted true genes with better accuracy than other well-known approaches.

Echinochloa crus-galli genome analysis provides insight into its adaptation and invasiveness as a weed

Barnyardgrass (Echinochloa crus-galli) is a pernicious weed in agricultural fields worldwide. The molecular mechanisms underlying its success in the absence of human intervention are presently unknown. Here we report a draft genome sequence of the hexaploid species E. crus-galli, i.e., a 1.27 Gb assembly representing 90.7% of the predicted genome size. An extremely large repertoire of genes encoding cytochrome P450 monooxygenases and glutathione S-transferases associated with detoxification are found. Two gene clusters involved in the biosynthesis of an allelochemical 2,4-dihydroxy-7-methoxy-1,4-benzoxazin-3-one (DIMBOA) and a phytoalexin momilactone A are found in the E. crus-galli genome, respectively. The allelochemical DIMBOA gene cluster is activated in response to co-cultivation with rice, while the phytoalexin momilactone A gene cluster specifically to infection by pathogenic Pyricularia oryzae. Our results provide a new understanding of the molecular mechanisms underlying the extreme adaptation of the weed.

Insights into Land Plant Evolution Garnered from the Marchantia polymorpha Genome

The evolution of land flora transformed the terrestrial environment. Land plants evolved from an ancestral charophycean alga from which they inherited developmental, biochemical, and cell biological attributes. Additional biochemical and physiological adaptations to land, and a life cycle with an alternation between multicellular haploid and diploid generations that facilitated efficient dispersal of desiccation tolerant spores, evolved in the ancestral land plant. We analyzed the genome of the liverwort Marchantia polymorpha, a member of a basal land plant lineage. Relative to charophycean algae, land plant genomes are characterized by genes encoding novel biochemical pathways, new phytohormone signaling pathways (notably auxin), expanded repertoires of signaling pathways, and increased diversity in some transcription factor families. Compared with other sequenced land plants, M. polymorpha exhibits low genetic redundancy in most regulatory pathways, with this portion of its genome resembling that predicted for the ancestral land plant. PAPERCLIP.

Molecular cloning and expression analysis of duplicated polyphenol oxidase genes reveal their functional differentiations in sorghum

Polyphenol oxidase (PPO) is believed to play a role in plant growth, reproduction, and resistance to pathogens and pests. PPO causes browning of grains in cereals. In this study, genetic mapping of sorghum grain for phenol color reaction (PHR) was performed using a recombinant inbred line population. Only one locus was detected between SSR markers SM06072 and Xtxp176 on chromosome 6. Two linked orthologous genes (Sb06PPO1 and Sb06PPO2) within the mapped region were discovered and cloned. Transformation experiments using Nipponbare (a PHR negative rice cultivar) showed that Sb06PPO1 from LTR108 and two Sb06PPO2 alleles from both varieties could complement Nipponbare, whereas Sb06PPO1 from 654 could not. Subsequent quantitative real-time PCR (qPCR) experiments showed that Sb06PPO1 and Sb06PPO2 functioned diversely, Sb06PPO1 was mainly expressed in young panicles before flowering. Sb06PPO2 was strongly expressed in flowering panicles, especially in hulls and branches at filling stage. Moreover, the expression of Sb06PPO1 was found to be significantly up-regulated by exogenous ABA and salt, whereas Sb06PPO2 was not changed significantly, further demonstrating functional differentiation between the two genes.

Genome sequencing of the staple food crop white Guinea yam enables the development of a molecular marker for sex determination

Background

Root and tuber crops are a major food source in tropical Africa. Among these crops are several species in the monocotyledonous genus Dioscorea collectively known as yam, a staple tuber crop that contributes enormously to the subsistence and socio-cultural lives of millions of people, principally in West and Central Africa. Yam cultivation is constrained by several factors, and yam can be considered a neglected “orphan” crop that would benefit from crop improvement efforts. However, the lack of genetic and genomic tools has impeded the improvement of this staple crop.

Results

To accelerate marker-assisted breeding of yam, we performed genome analysis of white Guinea yam (Dioscorea rotundata) and assembled a 594-Mb genome, 76.4% of which was distributed among 21 linkage groups. In total, we predicted 26,198 genes. Phylogenetic analyses with 2381 conserved genes revealed that Dioscorea is a unique lineage of monocotyledons distinct from the Poales (rice), Arecales (palm), and Zingiberales (banana). The entire Dioscorea genus is characterized by the occurrence of separate male and female plants (dioecy), a feature that has limited efficient yam breeding. To infer the genetics of sex determination, we performed whole-genome resequencing of bulked segregants (quantitative trait locus sequencing [QTL-seq]) in F1 progeny segregating for male and female plants and identified a genomic region associated with female heterogametic (male = ZZ, female = ZW) sex determination. We further delineated the W locus and used it to develop a molecular marker for sex identification of Guinea yam plants at the seedling stage.

Conclusions

Guinea yam belongs to a unique and highly differentiated clade of monocotyledons. The genome analyses and sex-linked marker development performed in this study should greatly accelerate marker-assisted breeding of Guinea yam. In addition, our QTL-seq approach can be utilized in genetic studies of other outcrossing crops and organisms with highly heterozygous genomes. Genomic analysis of orphan crops such as yam promotes efforts to improve food security and the sustainability of tropical agriculture.

Electronic supplementary material

The online version of this article (doi:10.1186/s12915-017-0419-x) contains supplementary material, which is available to authorized users.

Pearl millet genome sequence provides a resource to improve agronomic traits in arid environments

Draft genome, 994 re-sequenced lines and GWAS for yield-traits provide a resource of genetics and genomics tools for pearl millet researchers and breeders.

Pearl millet [Cenchrus americanus (L.) Morrone] is a staple food for more than 90 million farmers in arid and semi-arid regions of sub-Saharan Africa, India and South Asia. We report the ∼1.79 Gb draft whole genome sequence of reference genotype Tift 23D₂B₁-P1-P5, which contains an estimated 38,579 genes. We highlight the substantial enrichment for wax biosynthesis genes, which may contribute to heat and drought tolerance in this crop. We resequenced and analyzed 994 pearl millet lines, enabling insights into population structure, genetic diversity and domestication. We use these resequencing data to establish marker trait associations for genomic selection, to define heterotic pools, and to predict hybrid performance. We believe that these resources should empower researchers and breeders to improve this important staple crop.

Evidence-based gene models for structural and functional annotations of the oil palm genome

Background

Oil palm is an important source of edible oil. The importance of the crop, as well as its long breeding cycle (10-12 years) has led to the sequencing of its genome in 2013 to pave the way for genomics-guided breeding. Nevertheless, the first set of gene predictions, although useful, had many fragmented genes. Classification and characterization of genes associated with traits of interest, such as those for fatty acid biosynthesis and disease resistance, were also limited. Lipid-, especially fatty acid (FA)-related genes are of particular interest for the oil palm as they specify oil yields and quality. This paper presents the characterization of the oil palm genome using different gene prediction methods and comparative genomics analysis, identification of FA biosynthesis and disease resistance genes, and the development of an annotation database and bioinformatics tools.

Results

Using two independent gene-prediction pipelines, Fgenesh++ and Seqping, 26,059 oil palm genes with transcriptome and RefSeq support were identified from the oil palm genome. These coding regions of the genome have a characteristic broad distribution of GC₃ (fraction of cytosine and guanine in the third position of a codon) with over half the GC₃-rich genes (GC₃ ≥ 0.75286) being intronless. In comparison, only one-seventh of the oil palm genes identified are intronless. Using comparative genomics analysis, characterization of conserved domains and active sites, and expression analysis, 42 key genes involved in FA biosynthesis in oil palm were identified. For three of them, namely EgFABF, EgFABH and EgFAD3, segmental duplication events were detected. Our analysis also identified 210 candidate resistance genes in six classes, grouped by their protein domain structures.

Conclusions

We present an accurate and comprehensive annotation of the oil palm genome, focusing on analysis of important categories of genes (GC₃-rich and intronless), as well as those associated with important functions, such as FA biosynthesis and disease resistance. The study demonstrated the advantages of having an integrated approach to gene prediction and developed a computational framework for combining multiple genome annotations. These results, available in the oil palm annotation database (http://palmxplore.mpob.gov.my), will provide important resources for studies on the genomes of oil palm and related crops.

Reviewers

This article was reviewed by Alexander Kel, Igor Rogozin, and Vladimir A. Kuznetsov.

Electronic supplementary material

The online version of this article (doi:10.1186/s13062-017-0191-4) contains supplementary material, which is available to authorized users.

The influence of the mating type on virulence of Mucor irregularis

Mucor irregularis is an emerging fungal pathogen that cause cutaneous infection and could cause death. However, little is known about its mechanism of pathogenesis. There is evidence suggesting virulence vary with mating types in fungi, including the Mucorales. Here, we characterized the mating type locus of M. irregularis and the mating type ratio of 17 clinical isolates in China. Genomic data indicated M. irregularis is heterothallic having two mating types – bearing either SexP or SexM allele. Also, we employed a mice model to study the inflammation and pathological effects of different mating types. The comparison of the inflammatory response, cytokine profiles and Th-1, Th-2 and Th-17 cells numbers in each mating type treated mice showed that the severity and disease progress were enhanced in (+) mating type treated mice. One (+/0) mutant strain, with multiple mutations at the mating locus, had defects in sexual mating ability but appeared to be more virulent than the (−) mating type. Although (+) mating type appeared to be more virulent, most of our clinical isolates presented belonged to (−) mating type. Our findings support the involvement of MAT genes in sexual fertility, and the influence of mating type on the severity of cutaneous infection.

Subgenome Dominance in an Interspecific Hybrid, Synthetic Allopolyploid, and a 140-Year-Old Naturally Established Neo-Allopolyploid Monkeyflower

Recent studies have shown that one of the parental subgenomes in ancient polyploids is generally more dominant, having retained more genes and being more highly expressed, a phenomenon termed subgenome dominance. The genomic features that determine how quickly and which subgenome dominates within a newly formed polyploid remain poorly understood. To investigate the rate of emergence of subgenome dominance, we examined gene expression, gene methylation, and transposable element (TE) methylation in a natural, <140-year-old allopolyploid (Mimulus peregrinus), a resynthesized interspecies triploid hybrid (M. robertsii), a resynthesized allopolyploid (M. peregrinus), and progenitor species (M. guttatus and M. luteus). We show that subgenome expression dominance occurs instantly following the hybridization of divergent genomes and significantly increases over generations. Additionally, CHH methylation levels are reduced in regions near genes and within TEs in the first-generation hybrid, intermediate in the resynthesized allopolyploid, and are repatterned differently between the dominant and recessive subgenomes in the natural allopolyploid. Subgenome differences in levels of TE methylation mirror the increase in expression bias observed over the generations following hybridization. These findings provide important insights into genomic and epigenomic shock that occurs following hybridization and polyploid events and may also contribute to uncovering the mechanistic basis of heterosis and subgenome dominance.