Homoploid hybrid origin of Yucca gloriosa: intersectional hybrid speciation in Yucca (Agavoideae,Asparagaceae)

There is a growing appreciation for the importance of hybrid speciation in angiosperm evolution. Here, we show that Yucca gloriosa (Asparagaceae: Agavoideae) is the product of intersectional hybridization between Y. aloifolia and Y. filamentosa. These species, all named by Carl Linnaeus, exist in sympatry along the southeastern Atlantic coast of the United States. Yucca gloriosa was found to share a chloroplast haplotype with Y. aloifolia in all populations sampled. In contrast, nuclear gene‐based microsatellite markers in Y. gloriosa are shared with both parents. The hybrid origin of Y. gloriosa is supported by multilocus analyses of the nuclear microsatellite markers including principal coordinates analysis (PCO), maximum‐likelihood hybrid index scoring (HINDEX), and Bayesian cluster analysis (STRUCTURE). The putative parental species share only one allele at a single locus, suggesting there is little to no introgressive gene flow occurring between these species and Y. gloriosa. At the same time, diagnostic markers are segregating in Y. gloriosa populations. Lack of variation in the chloroplast of Y. aloifolia, the putative maternal parent, makes it difficult to rule out multiple hybrid origins of Y. gloriosa, but allelic variation at nuclear loci can be explained by a single hybrid origin of Y. gloriosa. Overall, these data provide strong support for the homoploid hybrid origin of Y. gloriosa.     

Comprehensive transcriptome analysis of fluid shear stress altered gene expression in renal epithelial cells

Renal epithelial cells are exposed to mechanical forces due to flow‐induced shear stress within the nephrons. Shear stress is altered in renal diseases caused by tubular dilation, obstruction, and hyperfiltration, which occur to compensate for lost nephrons. Fundamental in regulation of shear stress are primary cilia and other mechano‐sensors, and defects in cilia formation and function have profound effects on development and physiology of kidneys and other organs. We applied RNA sequencing to get a comprehensive overview of fluid‐shear regulated genes and pathways in renal epithelial cells. Functional enrichment‐analysis revealed TGF‐β, MAPK, and Wnt signaling as core signaling pathways up‐regulated by shear. Inhibitors of TGF‐β and MAPK/ERK signaling modulate a wide range of mechanosensitive genes, identifying these pathways as master regulators of shear‐induced gene expression. However, the main down‐regulated pathway, that is, JAK/STAT, is independent of TGF‐β and MAPK/ERK. Other up‐regulated cytokine pathways include FGF, HB‐EGF, PDGF, and CXC. Cellular responses to shear are modified at several levels, indicated by altered expression of genes involved in cell‐matrix, cytoskeleton, and glycocalyx remodeling, as well as glycolysis and cholesterol metabolism. Cilia ablation abolished shear induced expression of a subset of genes, but genes involved in TGF‐β, MAPK, and Wnt signaling were hardly affected, suggesting that other mechano‐sensors play a prominent role in the shear stress response of renal epithelial cells. Modulations in signaling due to variations in fluid shear stress are relevant for renal physiology and pathology, as suggested by elevated gene expression at pathological levels of shear stress compared to physiological shear.     

下一代测序中ChIP-seq数据的处理与分析

将染色质免疫共沉淀技术(ChIP)与下一代高通量测序技术相结合的染色质免疫共沉淀测序(ChIP-seq),已成为功能基因组学、特别是基因表达调控领域研究的关键技术。ChIP-seq实验带来的海量数据向生物信息学研究人员提出了新的挑战。由于此领域数据处理技术的发展大大滞后于实验技术进步,有必要系统地介绍和回顾ChIP-seq数据处理的各个方面,以便更多研究人员进入此领域设计或改进相应的算法。文章结合实例详细介绍了ChIP-seq数据整个流程,并重点讨论了其中的主要问题和关键环节,为这一研究领域的科研人员提供一个快速而深入的认识。     

Phylogenomic analysis clarifies the evolutionary origin of Coffea arabica

Interspecific hybridization events have played a major role in plant speciation, yet the evolutionary origin of hybrid species often remains enigmatic. Here, we inferred the evolutionary origin of the allotetraploid species Coffea arabica, which is widely cultivated for Arabica coffee production. We estimated genetic distances between C. arabica and all species that are known to be closely related to C. arabica using genotyping-by-sequencing (GBS) data. In addition, we reconstructed a time-calibrated multilabeled phylogenetic tree of 24 species to estimate the age of the C. arabica hybridization event. Ancestral states of self-compatibility were also inferred to shed new light on the evolution of self-compatibility in Coffea. Coffea canephora and C. eugenioides were confirmed as the putative progenitor species of C. arabica. These species most likely hybridized between 1.08 million and 543 000 years ago, coinciding with periods of environmental upheaval, which may have induced range shifts of the progenitor species that facilitated the emergence of C. arabica.     

Mercator: a fast and simple web server for genome scale functional annotation of plant sequence data

Next‐generation technologies generate an overwhelming amount of gene sequence data. Efficient annotation tools are required to make these data amenable to functional genomics analyses. The Mercator pipeline automatically assigns functional terms to protein or nucleotide sequences. It uses the MapMan ‘BIN’ ontology, which is tailored for functional annotation of plant ‘omics’ data. The classification procedure performs parallel sequence searches against reference databases, compiles the results and computes the most likely MapMan BINs for each query. In the current version, the pipeline relies on manually curated reference classifications originating from the three reference organisms (Arabidopsis, Chlamydomonas, rice), various other plant species that have a reviewed SwissProt annotation, and more than 2000 protein domain and family profiles at InterPro, CDD and KOG. Functional annotations predicted by Mercator achieve accuracies above 90% when benchmarked against manual annotation. In addition to mapping files for direct use in the visualization software MapMan, Mercator provides graphical overview charts, detailed annotation information in a convenient web browser interface and a MapMan‐to‐GO translation table to export results as GO terms. Mercator is available free of charge via http://mapman.gabipd.org/web/guest/app/Mercator .     

MicNeSs: genotyping microsatellite loci from a collection of (NGS) reads

Microsatellites are widely used in population genetics to uncover recent evolutionary events. They are typically genotyped using capillary sequencer, which capacity is usually limited to 9, at most 12 loci for each run, and which analysis is a tedious task that is performed by hand. With the rise of next‐generation sequencing (NGS), a much larger number of loci and individuals are available from sequencing: for example, on a single run of a GS Junior, 28 loci from 96 individuals are sequenced with a 30X cover. We have developed an algorithm to automatically and efficiently genotype microsatellites from a collection of reads sorted by individual (e.g. specific PCR amplifications of a locus or a collection of reads that encompass a locus of interest). As the sequencing and the PCR amplification introduce artefactual insertions or deletions, the set of reads from a single microsatellite allele shows several length variants. The algorithm infers, without alignment, the true unknown allele(s) of each individual from the observed distributions of microsatellites length of all individuals. MicNeSs, a python implementation of the algorithm, can be used to genotype any microsatellite locus from any organism and has been tested on 454 pyrosequencing data of several loci from fruit flies (a model species) and red deers (a nonmodel species). Without any parallelization, it automatically genotypes 22 loci from 441 individuals in 11 hours on a standard computer. The comparison of MicNeSs inferences to the standard method shows an excellent agreement, with some differences illustrating the pros and cons of both methods.     

江苏常熟地区汉族人群47个微单倍型的遗传多态性及遗传结构分析

目的 探索江苏常熟地区汉族人群在47个常染色体微单倍型基因座的遗传多态性,评估应用效能及法庭科学参数。方法 采用MHSeqTyper47混合DNA鉴定试剂盒进行基因座复合扩增及文库构建,使用MiSeq FGx测序平台进行测序,测得的数据应用MHTyper数据分析软件进行分析,对获得的样本遗传信息进行评估,结合千人基因组数据（1000 Genomes Project phase 3,1KG）评估群体间遗传分化指数及遗传距离,并计算法庭科学参数。结果 江苏常熟地区汉族人群与1KG中的中国北京群体的遗传分化和遗传距离最小,并得到最接近的有效等位基因数（A_e）,累积随机匹配概率（combined matching probability,CMP）与1KG中东亚参考人群的5个群体均接近,为1.25×10^-36,累积非父排除概率达0.999 999 999 964 1。结论 本研究报告了47个微单倍型基因座在江苏常熟地区汉族人群中的等位基因频率及遗传多态性信息,为47个微单倍型在法医学应用中提供了数据基础。另外,比较了1KG参考人群与江苏常熟地区汉族人群的多态性差异,并揭示了47个微单倍型在江苏常熟地区汉族人群中的遗传结构。总的来说,1KG中的东亚人群参考数据更符合江苏常熟地区汉族人群的遗传特征。     

Comparison of the experimental methods in haplotype sequencing via next generation sequencing

Although the diploid nature has been observed for over 50 years, phasing the diploid is still a laborious task. The speed and throughput of next generation sequencing have largely increased in the past decades. However, the short read-length remains one of the biggest challenges of haplotype analysis. For instance, reads as short as 150 bp span no more than one variant in most cases. Numerous experimental technologies have been developed to overcome this challenge. Distance, complexity and accuracy of the linkages obtained are the main factors to evaluate the efficiency of whole genome haplotyping methods. Here, we review these experimental technologies, evaluating their efficiency in linkages obtaining and system complexity. The technologies are organized into four categories based on its strategy: (i) chromosomes separation, (ii) dilution pools, (iii) crosslinking and proximity ligation, (ix) long-read technologies. Within each category, several subsections are listed to classify each technology. Innovative experimental strategies are expected to have high-quality performance, low cost and be labor-saving, which will be largely desired in the future.