Accurate genome relative abundance estimation based on shotgun metagenomic reads

Accurate estimation of microbial community composition based on metagenomic sequencing data is fundamental for subsequent metagenomics analysis. Prevalent estimation methods are mainly based on directly summarizing alignment results or its variants; often result in biased and/or unstable estimates. We have developed a unified probabilistic framework (named GRAMMy) by explicitly modeling read assignment ambiguities, genome size biases and read distributions along the genomes. Maximum likelihood method is employed to compute Genome Relative Abundance of microbial communities using the Mixture Model theory (GRAMMy). GRAMMy has been demonstrated to give estimates that are accurate and robust across both simulated and real read benchmark datasets. We applied GRAMMy to a collection of 34 metagenomic read sets from four metagenomics projects and identified 99 frequent species (minimally 0.5% abundant in at least 50% of the data-sets) in the human gut samples. Our results show substantial improvements over previous studies, such as adjusting the over-estimated abundance for Bacteroides species for human gut samples, by providing a new reference-based strategy for metagenomic sample comparisons. GRAMMy can be used flexibly with many read assignment tools (mapping, alignment or composition-based) even with low-sensitivity mapping results from huge short-read datasets. It will be increasingly useful as an accurate and robust tool for abundance estimation with the growing size of read sets and the expanding database of reference genomes.     

棕果蝠微卫星位点的筛选及其对近缘种的通用性

棕果蝠在我国热带和亚热带地区分布广、种群数量大,对龙眼、荔枝等水果种植业造成巨大危害。近年来,棕果蝠在国内被大范围捕杀,造成其种群数量急剧下降。为了获得用于研究棕果蝠的种群结构和遗传多样性的分子标记,我们利用探针富集法构建了棕果蝠部分微卫星基因组文库,筛选出棕果蝠的4 个高度多态性位点;并结合前期本研究组已经发表的棕果蝠微卫星位点,共检测了10 个棕果蝠特异的微卫星标记在其它4 种果蝠(短耳犬蝠、犬蝠、大长舌果蝠、小长舌果蝠)内的通用性,为不同果蝠种群遗传学研究提供了更多的微卫星分子标记。     

DNA reassociation kinetics of closely related species

The kinetics of reassociation of the DNA of three groups of closely related organisms were examined. The laboratory mouse was compared to an Asiatic mouse, whose chromosome number is the same but whose chromosome organization is different. Chinese hamster (2N=22) was compared to Syrian hamster (2N=44), and Haplopappus gracilis (2N=4) was compared to H. ravenit (2N=8). It was found that the most highly repeated DNA fractions of the three comparative sets of organisms differ in their reaction rates. However, these fractions of the related hamsters, haplopappi, and probably the mice, do not differ in the amount of DNA composing the fractions. The intermediately fast reassociating DNA and the unique DNA do not differ between members of related pairs of organisms. The implication of these results is that a short sequence of DNA may be highly copied in one organism, while in a related organism a longer DNA sequence is repeated a fewer number of times, and the total amount of repeated DNA may be the same in both related organisms.     

Contrasting genetic structuring in the closely related basidiomycetes Trichaptum abietinum and Trichaptum fuscoviolaceum (Hymenochaetales)

Trichaptum abietinum and Trichaptum fuscoviolaceum (Hymenochaetales, Basidiomycota) are closely related saprotrophic fungi, widely distributed on coniferous wood in temperate regions worldwide. Three intersterility groups have previously been detected in T. abietinum, while no prezygotic barriers have been proven within T. fuscoviolaceum. The aim of this study was to reveal the phylogeography and genetic relationship between these two closely related species and to explore whether the previously observed intersterility groups in T. abietinum are reflected in the genetic data. We assembled worldwide fruit body collections of both species (N = 314) and generated DNA sequences from three nuclear (ITS2, LSU, IGS) and one mitochondrial rDNA region (mtLSU). The two species are genetically well separated in all analyses. In correspondence with observations from earlier mating studies, our results revealed that T. fuscoviolaceum is genetically more uniform than T. abietinum. Multiple genetic sub-groups exist in T. abietinum that may correspond to the previously observed intersterility groups. However, there is low consistency across the investigated loci in delimiting the different sub-groups, except for a consistent North American group. As for many other widespread fungi, a complex phylogeographic pattern is found in T. abietinum which may have been formed by geographic, as well as multiple genetic intersterility barriers.     

Construction of a dairy microbial genome catalog opens new perspectives for the metagenomic analysis of dairy fermented products

Background

Microbial communities of traditional cheeses are complex and insufficiently characterized. The origin, safety and functional role in cheese making of these microbial communities are still not well understood. Metagenomic analysis of these communities by high throughput shotgun sequencing is a promising approach to characterize their genomic and functional profiles. Such analyses, however, critically depend on the availability of appropriate reference genome databases against which the sequencing reads can be aligned.

Results

We built a reference genome catalog suitable for short read metagenomic analysis using a low-cost sequencing strategy. We selected 142 bacteria isolated from dairy products belonging to 137 different species and 67 genera, and succeeded to reconstruct the draft genome of 117 of them at a standard or high quality level, including isolates from the genera Kluyvera, Luteococcus and Marinilactibacillus, still missing from public database. To demonstrate the potential of this catalog, we analysed the microbial composition of the surface of two smear cheeses and one blue-veined cheese, and showed that a significant part of the microbiota of these traditional cheeses was composed of microorganisms newly sequenced in our study.

Conclusions

Our study provides data, which combined with publicly available genome references, represents the most expansive catalog to date of cheese-associated bacteria. Using this extended dairy catalog, we revealed the presence in traditional cheese of dominant microorganisms not deliberately inoculated, mainly Gram-negative genera such as Pseudoalteromonas haloplanktis or Psychrobacter immobilis, that may contribute to the characteristics of cheese produced through traditional methods.

Electronic supplementary material

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

Discrimination between closely related Triticeae species using genomic DNA as a probe

Summary Labelled total genomic DNA was used as a probe in combination with blocking DNA to discriminate between taxonomically closely related species in the genera Hordeum and Secale. Discrimination was possible both by Southern hybridization to size-fractionated restriction enzyme digests of genomic DNA and by in situ hybridization to chromosome preparations. To distinguish between two species (e.g. H. vulgare and H. bulbosum), genomic DNA from one species was used as the labelled probe, while unlabelled DNA from the other species was applied at a much higher concentration as a block. The blocking DNA presumably hybridized to sequences in common between the block and the labelled probe, and between the block and DNA sequences on the membrane or chromosomes in situ. If so, mainly species-specific sequences would remain as sites for probe hybridization. These species-specific sequences are dispersed and represent a substantial proportion of the genome (unlike many cloned, species-specific sequences). Consequently, rapid nonradioactive methods detected probe hybridization sites satisfactorily. The method was able to confirm the parentage of hybrid plants. It has potentially wide application in plant breeding for the detection of alien DNA transfer, and it can be easily adapted to many species.     

促旋酶(gyrase) B亚单位基因gyrB在鉴别细菌近缘种中的应用

单拷贝的gyrB是普遍存在于细菌中编码促旋酶B亚单位的基因,该基因进化速率快,每100万年的平均碱基替换率为0.7%～0.8%.研究证明,该区段基因能对假单胞菌、芽孢杆菌、弧菌、肠杆菌、分枝杆菌、气单胞菌、乳酸菌等不同属或科内的近缘种进行区分鉴定,还可通过设计种特异性引物进行定量PCR或者限制性片段分析,同时也能结合变性梯度凝胶电泳技术(DGGE)追踪微生物的动态.gyrB基因弥补了非蛋白编码基因16S rDNA或者基因间隔序列(ITS)DNA无法区分近缘种的缺陷,给近缘种的鉴别或者其群体指纹图谱的分析带来了新的希望,为当前国内外用于近缘种研究的热点,是细菌系统发育分析上非常值得重视的新靶标.     

Limits to the estimation of species richness: The use of relative abundance distributions

The present study demonstrates the possibility of estimating species numbers of animal or plant communities from samples using relative abundance distributions. We use log‐abundance–species‐rank order plots and derive two new estimators that are based on log‐series and lognormal distributions. At small to moderate sample sizes these estimators appear to be more precise than previous parametric and nonparametric estimators. We test our estimators using samples from 171 published medium‐sized to large animal and plant communities taken from the literature. By this we show that our new estimators define also limits of precision.     

Nuclear simple sequence repeat markers are superior to DNA barcodes for identification of closely related Rhododendron species on the same mountain

Accurate species delimitation of sampled biological material is critical for a range of studies. Although the DNA barcodes developed in recent years are useful for identifying numerous well differentiated species that have not experienced frequent gene flow, they fail to delimit recently diverged species, especially those with extensive introgressions. Here we use five Rhododendron species growing together on the same mountain as a model system to compare the species delimitation effectiveness of the DNA barcodes (internal transcribed spacer, matK, psbA‐trnH, and rbcL) previously proposed versus 15 pairs of microsatellite markers. Using these markers, we genotyped 129 individuals, which were members of five species according to morphological identification. We identified five simple sequence repeat genetic clusters (independently evolving lineages) corresponding to the morphological identification. However, we found that numerous individuals contained cryptic hybrid introgressions from the other species. The four DNA barcodes could not delimit three out of four closely related species that showed clear morphological differentiation and cryptic introgressions. Even after excluding all cryptic hybrids, two closely related species could not be successfully identified. The low discrimination ability of the DNA barcodes for closely related Rhododendron species could result from two, not mutually exclusive factors: introgressive hybridization and incomplete lineage sorting. Our results highlight the importance of simple sequence repeat markers in delimiting closely related species and identifying cryptic introgressions in the absence of morphological changes.     

A new statistical approach for assessing similarity of species composition with incidence and abundance data

The classic Jaccard and Sørensen indices of compositional similarity (and other indices that depend upon the same variables) are notoriously sensitive to sample size, especially for assemblages with numerous rare species. Further, because these indices are based solely on presence–absence data, accurate estimators for them are unattainable. We provide a probabilistic derivation for the classic, incidence‐based forms of these indices and extend this approach to formulate new Jaccard‐type or Sørensen‐type indices based on species abundance data. We then propose estimators for these indices that include the effect of unseen shared species, based on either (replicated) incidence‐ or abundance‐based sample data. In sampling simulations, these new estimators prove to be considerably less biased than classic indices when a substantial proportion of species are missing from samples. Based on species‐rich empirical datasets, we show how incorporating the effect of unseen shared species not only increases accuracy but also can change the interpretation of results.     

Determining the relative roles of species replacement and species richness differences in generating beta‐diversity patterns

Aim To determine the relative contribution of species replacement and species richness differences to the emergence of beta‐diversity patterns. Innovation A novel method that disentangles all compositional differences (β_cc, overall beta diversity) in its two components, species replacement (β_‐3) and species richness differences (β_rich) is proposed. The performance of the method was studied with ternary plots, which allow visualization of the influence of the relative proportions of shared and unique species of two sites over each metric. The method was also tested in different hypothetical gradients and with real datasets. The novel method was compared with a previous proposal based on the partitioning of overall compositional differences (β_sor) in replacement (β_sim) and nestedness (β_nes). The linear response of β_cc contrasts with the curvilinear response of β_sor to linear gradients of dissimilarity. When two sites did not share any species, β_sim was always 1 and β_‐3 only reached 1 when the number of exclusive species of both sites was equal. β_‐3 remained constant along gradients of richness differences with constant replacement, while β_sim decreased. β_rich had a linear response to a linear gradient of richness differences with constant species replacement, whereas β_nes exhibited a hump‐shaped response. Moreover, β_sim > β_nes when clearly almost all species of one site were lost, whereas β_‐3 < β_rich in the same circumstances. Main conclusions The behaviour of the partition of β_cc into β_‐3 and β_rich is consistent with the variation of replacement and richness differences. The partitioning of β_sor into β_sim and β_nes overestimates the replacement component and underestimates richness differences. The novel methodology allows the discrimination of different causes of beta‐diversity patterns along latitudinal, biogeographic or ecological gradients, by estimating correctly the relative contributions of replacement and richness differences.     

Ambiguous species boundaries: Hybridization and morphological variation in two closely related Rubus species along altitudinal gradients

Although hybridization frequently occurs among plant species, hybrid zones of divergent lineages formed at species boundaries are less common and may not be apparent in later generations of hybrids with more parental‐like phenotypes, as a consequence of backcrossing. To determine the effects of dispersal and selection on species boundaries, we compared clines in leaf traits and molecular hybrid index along two hybrid zones on Yakushima Island, Japan, in which a temperate (Rubus palmatus) and subtropical (Rubus grayanus) species of wild raspberry are found. Leaf sinus depth in the two hybrid zones had narrower clines at 600 m a.s.l. than the molecular hybrid index and common garden tests confirmed that some leaf traits, including leaf sinus depth that is a major trait used in species identification, are genetically divergent between these closely related species. The sharp transition in leaf phenotypic traits compared to molecular markers indicated divergent selection pressure on the hybrid zone structure. We suggest that species boundaries based on neutral molecular data may differ from those based on observed morphological traits.     

Novel microsatellite markers for Dalechampia scandens (Euphorbiaceae) and closely related taxa: application to studying a species complex

We developed novel microsatellite markers for D alechampia scandens L. (Euphorbiaceae). The target plants belong to a distinct, but undescribed, species in the D . scandens species complex, characterized by small resin‐producing glands. In total, 110 alleles over 36 novel markers were identified across 39 individuals from three populations. The number of alleles varied from one to seven, with an average of 3.06 ± 0.26 alleles per locus. The developed markers, along with previously developed ones for a large‐glanded D . scandens species, were tested for amplification in 11 additional species of the genus D alechampia. Four markers did not produce any detectable allele in 37 individuals from two populations of the large‐glanded species. Average expected heterozygosity across all small‐ and large‐glanded specific loci was 0.36 and 0.15, for the small and large glanded populations, respectively. Cross‐species amplification showed that 89% of all markers were successfully amplified in at least one of the 11 other D alechampia species. These microsatellite markers may be useful for detecting undescribed species in the D . scandens species complex, and can be used for comparative analyses of genetic structure, mating system and phylogeography of other D alechampia species.     

涝渍胁迫下5种菊花近缘种属植物生理特性

采用土培模拟涝害的方法,对5种菊花近缘种属植物的根系活力、叶绿素含量和光合生理特性以及涝渍胁迫后的形态表现进行研究和观察,比较了其耐涝能力.结果表明:涝渍胁迫下,5种菊花近缘种属植物根系活力急剧下降,但紫花野菊在胁迫初期表现出强的根系活力;5种植物叶绿素含量总体呈现先升后降的趋势;涝渍胁迫15d时纪伊潮菊、泡黄金菊的光合强度明显降低,CO2同化作用下降,并出现负值,叶片已丧失光合作用能力,其次是大岛野路菊和那贺川野菊,而紫花野菊维持了一定水平的净光合速率.结合叶片外观形态观察发现,紫花野菊最为耐涝,泡黄金菊最不耐涝,而大岛野路菊、纪伊潮菊和那贺川野菊的耐涝性介于两者之间.     

不同地理种群美洲斑潜蝇及近缘种的rDNA-ITS1序列分析和比较

【目的】通过对美洲斑潜蝇Liriomyza sativae Blanchard 不同地理种群及近缘种间的核糖体DNA第一内转录间隔区（rDNA-ITS1）进行比较,分析美洲斑潜蝇不同地理种群间的遗传分化情况,并为美洲斑潜蝇与近缘种间提供分子鉴别标记。【方法】用PCR产物直接测序法及克隆测序法对我国美洲斑潜蝇8个地理种群的rDNA-ITS1序列进行测序,并调用GenBank中3个近缘种的rDNA-ITS序列,运用软件MEGA3.1对美洲斑潜蝇不同地理种群及近缘种间的rDNA-ITS1序列进行分析。【结果】美洲斑潜蝇8个地理种群间的分化程度较低,只有8个变异位点,遗传距离都在0.02以下,但4个近缘种间的碱基差异显著,遗传距离为0.149～0.390,有126个变异位点,12个美洲斑潜蝇特异性识别位点。【结论】虽然基于rDNA-ITS1序列所显示的美洲斑潜蝇各地理种群之间的遗传分化很小,但是其分化趋势与地理分布基本相吻合;得到的12个特异性识别位点不仅可以作为美洲斑潜蝇与其近缘种间鉴别的分子标记,而且可为今后设计鉴别性PCR引物提供重要的参考依据。     

转基因水稻外源基因向近缘种群漂流和渐渗的研究进展

水稻是我国最重要的粮食作物之一,我国有8亿以上的人口以稻米作为主食。但在水稻生产中,由于病、虫、草害及不良气候等逆境因子的影响,严重制约了水稻的高产、稳产。转基因生物技术的迅速发展,为水稻抗性育种提供了新途径。自20世纪80年代以来,我国全方位地开展了转基因水稻的研发,目前已经培育出大量的抗病、抗虫、抗除草剂和抗逆的转基因水稻品种,这将为提高我国水稻的生产力和确保粮食安全做出重要的贡献。但转基因水稻的基因漂流及其可能带来的生物安全问题备受关注。已有报道证明,外源转基因可以通过异交向非转基因品种和野生近缘种漂流。在不同的试验条件下,抗除草剂基因有0．05％-0．53％逃逸的可能,其向不育系的最大漂移频率可达4．518％。抗虫基因向相邻非转基因水稻的平均漂移频率最高为0．875％。因此,本文对水稻与其近缘野生种的杂交情况,转基因水稻外源基因向非转基因品种、野生近缘种以及野生非近缘种的漂流和渐渗及其潜在的生态环境风险等方面进行了简要分析,并对转基因水稻的发展进行了展望,以期为转基因水稻的安全应用提供参考。