Building the sugarcane genome for biotechnology and identifying evolutionary trends

Background

Sugarcane is the source of sugar in all tropical and subtropical countries and is becoming increasingly important for bio-based fuels. However, its large (10 Gb), polyploid, complex genome has hindered genome based breeding efforts. Here we release the largest and most diverse set of sugarcane genome sequences to date, as part of an on-going initiative to provide a sugarcane genomic information resource, with the ultimate goal of producing a gold standard genome.

Results

Three hundred and seventeen chiefly euchromatic BACs were sequenced. A reference set of one thousand four hundred manually-annotated protein-coding genes was generated. A small RNA collection and a RNA-seq library were used to explore expression patterns and the sRNA landscape. In the sucrose and starch metabolism pathway, 16 non-redundant enzyme-encoding genes were identified. One of the sucrose pathway genes, sucrose-6-phosphate phosphohydrolase, is duplicated in sugarcane and sorghum, but not in rice and maize. A diversity analysis of the s6pp duplication region revealed haplotype-structured sequence composition. Examination of hom(e)ologous loci indicate both sequence structural and sRNA landscape variation. A synteny analysis shows that the sugarcane genome has expanded relative to the sorghum genome, largely due to the presence of transposable elements and uncharacterized intergenic and intronic sequences.

Conclusion

This release of sugarcane genomic sequences will advance our understanding of sugarcane genetics and contribute to the development of molecular tools for breeding purposes and gene discovery.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-540) contains supplementary material, which is available to authorized users.     

Sex specific expression and distribution of small RNAs in papaya

Background

Regulatory function of small non-coding RNAs (sRNA) in response to environmental and developmental cues has been established. Additionally, sRNA, also plays an important role in maintaining the heterochromatin and centromere structures of the chromosome. Papaya, a trioecious species with recently evolved sex chromosomes, has emerged as an excellent model system to study sex determination and sex chromosome evolution in plants. However, role of small RNA in papaya sex determination is yet to be explored.

Results

We analyzed the high throughput sRNAs reads in the Illumina libraries prepared from male, female, and hermaphrodite flowers of papaya. Using the sRNA reads, we identified 29 miRNAs that were not previously reported from papaya. Including this and two previous studies, a total of 90 miRNAs has been identified in papaya. We analyzed the expression of these miRNAs in each sex types. A total of 65 miRNAs, including 31 conserved and 34 novel mirNA, were detected in at least one library. Fourteen of the 65 miRNAs were differentially expressed among different sex types. Most of the miRNA expressed higher in male flowers were related to the auxin signaling pathways, whereas the miRNAs expressed higher in female flowers were the potential regulators of the apical meristem identity genes. Aligning the sRNA reads identified the sRNA hotspots adjacent to the gaps of the X and Y chromosomes. The X and Y chromosomes sRNA hotspots has a 7.8 and 4.4 folds higher expression of sRNA, respectively, relative to the chromosome wide average. Approximately 75% of the reads aligned to the X chromosome hotspot was identical to that of the Y chromosome hotspot.

Conclusion

By analyzing the large-scale sRNA sequences from three sex types, we identified the sRNA hotspots flanking the gaps of papaya X, Y, and Y^h chromosome. The sRNAs expression patterns in these regions were reminiscent of the pericentromeric region indicating that the only remaining gap in each of these chromosomes is likely the centromere. We also identified 14 differentially expressed miRNAs in male, female and hermaphrodite flowers of papaya. Our results provide valuable information toward understanding the papaya sex determination.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-20) contains supplementary material, which is available to authorized users.     

The I2 resistance gene homologues in Solanum have complex evolutionary patterns and are targeted by miRNAs

Background

Several resistance traits, including the I2 resistance against tomato fusarium wilt, were mapped to the long arm of chromosome 11 of Solanum. However, the structure and evolution of this locus remain poorly understood.

Results

Comparative analysis showed that the structure and evolutionary patterns of the I2 locus vary considerably between potato and tomato. The I2 homologues from different Solanaceae species usually do not have orthologous relationship, due to duplication, deletion and frequent sequence exchanges. At least 154 sequence exchanges were detected among 76 tomato I2 homologues, but sequence exchanges between I2 homologues in potato is less frequent. Previous study showed that I2 homologues in potato were targeted by miR482. However, our data showed that I2 homologues in tomato were targeted by miR6024 rather than miR482. Furthermore, miR6024 triggers phasiRNAs from I2 homologues in tomato. Sequence analysis showed that miR6024 was originated after the divergence of Solanaceae. We hypothesized that miR6024 and miR482 might have facilitated the expansion of the I2 family in Solanaceae species, since they can minimize their potential toxic effects by down-regulating their expression.

Conclusions

The I2 locus represents a most divergent resistance gene cluster in Solanum. Its high divergence was partly due to frequent sequence exchanges between homologues. We propose that the successful expansion of I2 homologues in Solanum was at least partially attributed to miRNA mediated regulation.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-743) contains supplementary material, which is available to authorized users.     

Genomic and transcriptomic analysis of the endophytic fungus Pestalotiopsis fici reveals its lifestyle and high potential for synthesis of natural products

Background

In recent years, the genus Pestalotiopsis is receiving increasing attention, not only because of its economic impact as a plant pathogen but also as a commonly isolated endophyte which is an important source of bioactive natural products. Pestalotiopsis fici Steyaert W106-1/CGMCC3.15140 as an endophyte of tea produces numerous novel secondary metabolites, including chloropupukeananin, a derivative of chlorinated pupukeanane that is first discovered in fungi. Some of them might be important as the drug leads for future pharmaceutics.

Results

Here, we report the genome sequence of the endophytic fungus of tea Pestalotiopsis fici W106-1/CGMCC3.15140. The abundant carbohydrate-active enzymes especially significantly expanding pectinases allow the fungus to utilize the limited intercellular nutrients within the host plants, suggesting adaptation of the fungus to endophytic lifestyle. The P. fici genome encodes a rich set of secondary metabolite synthesis genes, including 27 polyketide synthases (PKSs), 12 non-ribosomal peptide synthases (NRPSs), five dimethylallyl tryptophan synthases, four putative PKS-like enzymes, 15 putative NRPS-like enzymes, 15 terpenoid synthases, seven terpenoid cyclases, seven fatty-acid synthases, and five hybrids of PKS-NRPS. The majority of these core enzymes distributed into 74 secondary metabolite clusters. The putative Diels-Alderase genes have undergone expansion.

Conclusion

The significant expansion of pectinase encoding genes provides essential insight in the life strategy of endophytes, and richness of gene clusters for secondary metabolites reveals high potential of natural products of endophytic fungi.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-014-1190-9) contains supplementary material, which is available to authorized users.     

Identification of Wolbachia-responsive microRNAs in the two-spotted spider mite,Tetranychus urticae

Background

The two-spotted spider mite, Tetranychus urticae, is infected with Wolbachia, which have the ability to manipulate host reproduction and fitness. MicroRNAs (miRNAs) are small non-coding RNAs that are involved in many biological processes such as development, reproduction and host-pathogen interactions. Although miRNA was observed to involve in Wolbachia-host interactions in the other insect systems, its roles have not been fully deciphered in the two-spotted spider mite.

Results

Small RNA libraries of infected and uninfected T. urticae for both sexes (in total four libraries) were constructed. By integrating the mRNA data originated from the same samples, the target genes of the differentially expressed miRNAs were predicted. Then, GO and pathway analyses were performed for the target genes. Comparison of libraries showed that Wolbachia infection significantly regulated 91 miRNAs in females and 20 miRNAs in males, with an overall suppression of miRNAs in Wolbachia-infected libraries. A comparison of the miRNA and mRNA data predicted that the differentially expressed miRNAs negatively regulated 90 mRNAs in females and 9 mRNAs in males. An analysis of target genes showed that Wolbachia-responsive miRNAs regulated genes with function in sphingolipid metabolism, lysosome function, apoptosis and lipid transporting in both sexes, as well as reproduction in females.

Conclusion

Comparisons of the miRNA and mRNA data can help to identify miRNAs and miRNA target genes involving in Wolbachia-host interactions. The molecular targets identified in this study should be useful in further functional studies.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-1122) contains supplementary material, which is available to authorized users.     

Spatial and temporal profiles of cytokinin biosynthesis and accumulation in developing caryopses of maize

Background and Aims

Cytokinins are a major group of plant hormones and are associated with various developmental processes. Developing caryopses of maize have high levels of cytokinins, but little is known about their spatial and temporal distribution. The localization and quantification of cytokinins was investigated in maize (Zea mays) caryopsis from 0 to 28 d after pollination together with the expression and localization of isopentenyltransferase ZmIPT1 involved in cytokinin biosynthesis and ZmCNGT, the gene putatively involved in N9-glucosylation.

Methods

Biochemical, cellular and molecular approaches resolved the overall cytokinin profiles, and several gene expression assays were used for two critical genes to assess cytokinin cell-specific biosynthesis and conversion to the biologically inactive form. Cytokinins were immunolocalized for the first time in maize caryopses.

Key Results

During the period 0–28 d after pollination (DAP): (1) large quantities of cytokinins were detected in the maternal pedicel region relative to the filial tissues during the early stages after fertilization; (2) unpollinated ovules did not accumulate cytokinins; (3) the maternal nucellar region showed little or no cytokinin signal; (4) the highest cytokinin concentrations in filial endosperm and embryo were detected at 12 DAP, predominantly zeatin riboside and zeatin-9-glucoside, respectively; and (5) a strong cytokinin immuno-signal was detected in specific cell types in the pedicel, endosperm and embryo.

Conclusions

The cytokinins of developing maize caryopsis may originate from both local syntheses as well as by transport. High levels of fertilization-dependent cytokinins in the pedicel suggest filial control on metabolism in the maternal tissue; they may also trigger developmental programmed cell death in the pedicel.     

Impaired macrophage phagocytosis of bacteria in severe asthma

Background

Bacteria are frequently cultured from sputum samples of severe asthma patients suggesting a defect in bacterial clearance from the airway. We measured the capacity of macrophages from patients with asthma to phagocytose bacteria.

Methods

Phagocytosis of fluorescently-labelled polystyrene beads, Haemophilus influenzae or Staphylococcus aureus by broncholaveolar lavage alveolar macrophages (AM) and by monocyte-derived macrophages (MDM) from non-asthmatics, mild-moderate and severe asthmatic patients was assessed using fluorimetry.

Results

There were no differences in phagocytosis of polystyrene beads by AMs or MDMs from any of the subject groups. There was reduced phagocytosis of Haemophilus influenzae and Staphylococcus aureus in MDMs from patients with severe asthma compared to non-severe asthma (p < 0.05 and p < 0.01, respectively) and healthy subjects (p < 0.01and p < 0.001, respectively). Phagocytosis of Haemophilus influenzae and Staphylococcus aureus by AM was also reduced in severe asthma compared to normal subjects (p < 0.05). Dexamethasone and formoterol did not suppress phagocytosis of bacteria by MDMs from any of the groups.

Conclusions

Persistence of bacteria in the lower airways may result partly from a reduced phagocytic capacity of macrophages for bacteria. This may contribute to increased exacerbations, airway colonization and persistence of inflammation.     

Theoretical evidence for the functional benefit of root cortical aerenchyma in soils with low phosphorus availability

Background and Aims

The formation of root cortical aerenchyma (RCA) reduces root respiration and nutrient content by converting living tissue to air volume. It was hypothesized that RCA increases soil resource acquisition by reducing the metabolic and phosphorus cost of soil exploration.

Methods

To test the quantitative logic of the hypothesis, SimRoot, a functional–structural plant model with emphasis on root architecture and nutrient acquisition, was employed. Sensitivity analyses for the effects of RCA on the initial 40 d of growth of maize (Zea mays) and common bean (Phaseolus vulgaris) were conducted in soils with varying degrees of phosphorus availability. With reference to future climates, the benefit of having RCA in high CO₂ environments was simulated.

Key Results

The model shows that RCA may increase the growth of plants faced with suboptimal phosphorus availability up to 70 % for maize and 14 % for bean after 40 d of growth. Maximum increases were obtained at low phosphorus availability (3 µm). Remobilization of phosphorus from dying cells had a larger effect on plant growth than reduced root respiration. The benefit of both these functions was additive and increased over time. Larger benefits may be expected for mature plants. Sensitivity analysis for light-use efficiency showed that the benefit of having RCA is relatively stable, suggesting that elevated CO₂ in future climates will not significantly effect the benefits of having RCA.

Conclusions

The results support the hypothesis that RCA is an adaptive trait for phosphorus acquisition by remobilizing phosphorus from the root cortex and reducing the metabolic costs of soil exploration. The benefit of having RCA in low-phosphorus soils is larger for maize than for bean, as maize is more sensitive to low phosphorus availability while it has a more ‘expensive’ root system. Genetic variation in RCA may be useful for breeding phosphorus-efficient crop cultivars, which is important for improving global food security.     

The diversity of small non-coding RNAs in the diatom Phaeodactylum tricornutum

Background

Marine diatoms constitute a major component of eukaryotic phytoplankton and stand at the crossroads of several evolutionary lineages. These microalgae possess peculiar genomic features and novel combinations of genes acquired from bacterial, animal and plant ancestors. Furthermore, they display both DNA methylation and gene silencing activities. Yet, the biogenesis and regulatory function of small RNAs (sRNAs) remain ill defined in diatoms.

Results

Here we report the first comprehensive characterization of the sRNA landscape and its correlation with genomic and epigenomic information in Phaeodactylum tricornutum. The majority of sRNAs is 25 to 30 nt-long and maps to repetitive and silenced Transposable Elements marked by DNA methylation. A subset of this population also targets DNA methylated protein-coding genes, suggesting that gene body methylation might be sRNA-driven in diatoms. Remarkably, 25-30 nt sRNAs display a well-defined and unprecedented 180 nt-long periodic distribution at several highly methylated regions that awaits characterization. While canonical miRNAs are not detectable, other 21-25 nt sRNAs of unknown origin are highly expressed. Besides, non-coding RNAs with well-described function, namely tRNAs and U2 snRNA, constitute a major source of 21-25 nt sRNAs and likely play important roles under stressful environmental conditions.

Conclusions

P. tricornutum has evolved diversified sRNA pathways, likely implicated in the regulation of largely still uncharacterized genetic and epigenetic processes. These results uncover an unexpected complexity of diatom sRNA population and previously unappreciated features, providing new insights into the diversification of sRNA-based processes in eukaryotes.

Electronic supplementary material

The online version of this article (doi:10.1186/1471-2164-15-698) contains supplementary material, which is available to authorized users.     

Evolutionary patchwork of an insecticidal toxin shared between plant-associated pseudomonads and the insect pathogens Photorhabdus and Xenorhabdus

Background

Root-colonizing fluorescent pseudomonads are known for their excellent abilities to protect plants against soil-borne fungal pathogens. Some of these bacteria produce an insecticidal toxin (Fit) suggesting that they may exploit insect hosts as a secondary niche. However, the ecological relevance of insect toxicity and the mechanisms driving the evolution of toxin production remain puzzling.

Results

Screening a large collection of plant-associated pseudomonads for insecticidal activity and presence of the Fit toxin revealed that Fit is highly indicative of insecticidal activity and predicts that Pseudomonas protegens and P. chlororaphis are exclusive Fit producers. A comparative evolutionary analysis of Fit toxin-producing Pseudomonas including the insect-pathogenic bacteria Photorhabdus and Xenorhadus, which produce the Fit related Mcf toxin, showed that fit genes are part of a dynamic genomic region with substantial presence/absence polymorphism and local variation in GC base composition. The patchy distribution and phylogenetic incongruence of fit genes indicate that the Fit cluster evolved via horizontal transfer, followed by functional integration of vertically transmitted genes, generating a unique Pseudomonas-specific insect toxin cluster.

Conclusions

Our findings suggest that multiple independent evolutionary events led to formation of at least three versions of the Mcf/Fit toxin highlighting the dynamic nature of insect toxin evolution.

Electronic supplementary material

The online version of this article (doi:10.1186/s12864-015-1763-2) contains supplementary material, which is available to authorized users.