The draft genome of the pear (Pyrus bretschneideri) using a combination of BAC-by-BAC and next-generation sequencing is reported. A 512.0-Mb sequence corresponding to 97.1% of the estimated genome size of this highly heterozygous species is assembled with 1943 coverage. High-density genetic maps comprising 2005 SNP markers anchored 75.5% of the sequence to all 17 chromosomes. The pear genome encodes 42,812 protein-coding genes, and of these,~28.5% encode multiple isoforms. Repetitive sequences of 271.9 Mb in length, accounting for 53.1% of the pear genome, are identified. Simulation of eudicots to the ancestor of Rosaceae has reconstructed nine ancestral chromosomes. Pear and apple diverged from each other~5.4-21.5 million years ago, and a recent whole-genome duplication (WGD) event must have occurred 30-45 MYA prior to their divergence, but following divergence from strawberry. When compared with the apple genome sequence, size differences between the apple and pear genomes are confirmed mainly due to the presence of repetitive sequences predominantly contributed by transposable elements (TEs), while genic regions are similar in both species. Genes critical for self-incompatibility, lignified stone cells (a unique feature of pear fruit), sorbitol metabolism, and volatile compounds of fruit have also been identified. Multiple candidate SFB genes appear as tandem repeats in the S-locus region of pear; while lignin synthesis-related gene family expansion and highly expressed gene families of HCT, C39H, and CCOMT contribute to high accumulation of both G-lignin and S-lignin. Moreover, alpha-linolenic acid metabolism is a key pathway for aroma in pear fruit.

BackgroundPear (Pyrus) is a globally grown fruit, with thousands of cultivars in five domesticated species and dozens of wild species. However, little is known about the evolutionary history of these pear species and what has contributed to the distinct phenotypic traits between Asian pears and European pears.ResultsWe report the genome resequencing of 113 pear accessions from worldwide collections, representing both cultivated and wild pear species. Based on 18,302,883 identified SNPs, we conduct phylogenetics, population structure, gene flow, and selective sweep analyses. Furthermore, we propose a model for the divergence, dissemination, and independent domestication of Asian and European pears in which pear, after originating in southwest China and then being disseminated throughout central Asia, has eventually spread to western Asia, and then on to Europe. We find evidence for rapid evolution and balancing selection for S-RNase genes that have contributed to the maintenance of self-incompatibility, thus promoting outcrossing and accounting for pear genome diversity across the Eurasian continent. In addition, separate selective sweep signatures between Asian pears and European pears, combined with co-localized QTLs and differentially expressed genes, underline distinct phenotypic fruit traits, including flesh texture, sugar, acidity, aroma, and stone cells.ConclusionsThis study provides further clarification of the evolutionary history of pear along with independent domestication of Asian and European pears. Furthermore, it provides substantive and valuable genomic resources that will significantly advance pear improvement and molecular breeding efforts.Electronic supplementary materialThe online version of this article (10.1186/s13059-018-1452-y) contains supplementary material, which is available to authorized users.

Universal oligonucleotide hybridization probes targeting the small-subunit rRNA are commonly used to quantify total microbial representation in environmental samples. Universal probes also serve to normalize results obtained with probes targeting specific phylogenetic groups of microorganisms. In this study, six universal probes were evaluated for stability of probe-target duplexes by using rRNA from nine organisms representing the three domains of Bacteria, Archaea, and Eucarya. Domain-specific variations in dissociation temperatures were observed for all probes. This could lead to a significant bias when these probes are used to quantify microbial populations in environmental samples. We suggest lowering the posthybridization wash stringency for two of the universal probes (S-*-Univ-1390-a-A-18 and S-*-Univ-1392-a-A-15) examined. These two probes were evaluated with traditional and modified hybridization conditions to characterize defined mixtures of rRNAs extracted from pure cultures and rRNA samples obtained from anaerobic digester samples. Probe S-*-Univ-1390-a-A-18 provided excellent estimations of domain-level community composition of these samples and is recommended for future use in microbial ecology studies.

The microbial community structure of twenty-one single-phase and one two-phase full-scale anaerobic sewage sludge digesters was evaluated using oligonucleotide probes complementary to conserved tracts of the 16S rRNAs of phylogenetically defined groups of methanogens and sulfate-reducing bacteria. These probe results were interpreted in combination with results from traditional chemical analyses and metabolic activity assays. It was determined that methanogens in "healthy" mesophilic, single-phase sewage sludge digesters accounted for approximately 8-12% of the total community and that Methanosarcinales and Methanomicrobiales constituted the majority of the total methanogen population. Methanobacteriales and Methanococcales played a relatively minor role in the digesters. Phylogenetic groups of mesophilic, Gram-negative sulfate-reducing bacteria were consistently present at significant levels: Desulfovibrio and Desulfobulbus spp. were the dominant sulfate-reducing populations, Desulfobacter and Desulfobacterium spp. were present at lower levels, and Desulfosarcina, Desulfococcus, and Desulfobotulus spp. were absent. Sulfate reduction by one or more of these populations played a significant role in all digesters evaluated in this study. In addition, sulfate-reducing bacteria played a role in favoring methanogenesis by providing their substrates. The analysis of the two-phase digester indicated that true phase separation was not accomplished: significant levels of active methanogens were present in the first phase. It was determined that the dominant populations in the second phase were different from those in the single-phase digesters.

A total of 355 simple sequence repeat (SSR) markers were developed, based on expressed sequence tag (EST) and bacterial artificial chromosome (BAC)-end sequence databases, and successfully used to construct an SSR-based genetic linkage map of the apple. The consensus linkage map spanned 1143 cM, with an average density of 2.5 cM per marker. Newly developed SSR markers along with 279 SSR markers previously published by the HiDRAS project were further used to integrate physical and genetic maps of the apple using a PCR-based BAC library screening approach. A total of 470 contigs were unambiguously anchored onto all 17 linkage groups of the apple genome, and 158 contigs contained two or more molecular markers. The genetically mapped contigs spanned ∼421 Mb in cumulative physical length, representing 60.0% of the genome. The sizes of anchored contigs ranged from 97 kb to 4.0 Mb, with an average of 995 kb. The average physical length of anchored contigs on each linkage group was ∼24.8 Mb, ranging from 17.0 Mb to 37.73 Mb. Using BAC DNA as templates, PCR screening of the BAC library amplified fragments of highly homologous sequences from homoeologous chromosomes. Upon integrating physical and genetic maps of the apple, the presence of not only homoeologous chromosome pairs, but also of multiple locus markers mapped to adjacent sites on the same chromosome was detected. These findings demonstrated the presence of both genome-wide and segmental duplications in the apple genome and provided further insights into the complex polyploid ancestral origin of the apple.

Three genes encoding flavonoid 3'-hydroxylase (F3'H) in apple (Malus x domestica), designated MdF3'HI, MdF3'HIIa, and MdF3'HIIb, have been identified. MdF3'HIIa and MdF3'HIIb are almost identical in amino acid sequences, and they are allelic, whereas MdF3'HI has 91% nucleotide sequence identity in the coding region to both MdF3'HIIa and MdF3'HIIb. MdF3'HI and MdF3'HII genes are mapped onto linkage groups 14 and 6, respectively, of the apple genome. Throughout the development of apple fruit, transcriptional levels of MdF3'H genes along with other anthocyanin biosynthesis genes are higher in the red-skinned cv Red Delicious than that in the yellow-skinned cv Golden Delicious. Moreover, patterns of MdF3'H gene expression correspond to accumulation patterns of flavonoids in apple fruit. These findings suggest that MdF3'H genes are coordinately expressed with other genes in the anthocyanin biosynthetic pathway in apple. The functionality of these apple F3'H genes has been demonstrated via their ectopic expression in both the Arabidopsis (Arabidopsis thaliana) transparent testa7-1 (tt7) mutant and tobacco (Nicotiana tabacum). When grown under nitrogen-deficient conditions, transgenic Arabidopsis tt7 seedlings expressing apple F3'H regained red color pigmentation and significantly accumulated both 4'-hydrylated pelargonidin and 3',4'-hydrylated cyanidin. When compared with wild-type plants, flowers of transgenic tobacco lines overexpressing apple F3'H genes exhibited enhanced red color pigmentation. This suggests that the F3'H enzyme may coordinately interact with other flavonoid enzymes in the anthocyanin biosynthesis pathway.

Microbial population dynamics were investigated during start-up and during periods of overload conditions in anaerobic co-digesters treating municipal solid waste and sewage sludge. Changes in community structure were monitored using ribosomal RNA-based oligonucleotide probe hybridization to measure the abundance of syntrophic propionate-oxidizing bacteria (SPOB), saturated fatty acid-beta-oxidizing syntrophs (SFAS), and methanogens. These changes were linked to traditional performance parameters such as biogas production and volatile fatty acid (VFA) concentrations. Digesters with high levels of Archaea started up successfully. Methanosaeta concilii was the dominant aceticlastic methanogen in these systems. In contrast, digesters that experienced a difficult start-up period had lower levels of Archaea with proportionally more abundant Methanosarcina spp. Syntrophic propionate-oxidizing bacteria and saturated fatty acid-beta-oxidizing syntrophs were present at low levels in all digesters, and SPOB appeared to play a role in stabilizing propionate levels during start-up of one digester. Digesters with a history of poor performance tolerated a severe organic overload event better than digesters that had previously performed well. It is hypothesized that higher levels of SPOB and SFAS and their methanogenic partners in previously unstable digesters are responsible for this behavior.

The sugar transporter (ST) plays an important role in plant growth, development and fruit quality. In this study, a total of 75 ST genes were identified in the pear (Pyrus bretschneideri Rehd) genome based on systematic analysis. Furthermore, all ST genes identified were grouped into eight subfamilies according to conserved domains and phylogenetic analysis. Analysis of cis-regulatory element sequences of all ST genes identified the MYBCOREATCYCB1 promoter in sucrose transporter (SUT) and monosaccharide transporter (MST) genes of pear, while in grape it is exclusively found in SUT subfamily members, indicating divergent transcriptional regulation in different species. Gene duplication event analysis indicated that whole-genome duplication (WGD) and segmental duplication play key roles in ST gene amplification, followed by tandem duplication. Estimation of positive selection at codon sites of ST paralog pairs indicated that all plastidic glucose translocator (pGlcT) subfamily members have evolved under positive selection. In addition, the evolutionary history of ST gene duplications indicated that the ST genes have experienced significant expansion in the whole ST gene family after the second WGD, especially after apple and pear divergence. According to the global RNA sequencing results of pear fruit development, gene expression profiling showed the expression of 53 STs. Combined with quantitative real-time PCR (qRT-PCR) analysis, two polyol/monosaccharide transporter (PLT) and three tonoplast monosaccharide transporter (tMT) members were identified as candidate genes, which may play important roles in sugar accumulation during pear fruit development and ripening. Identification of highly expressed STs in fruit is important for finding novel genes contributing to enhanced levels of sugar content in pear fruit.