Genome-wide comparative analysis of digital gene expression tag profiles during maize ear development

Background: Development of the maize (Zea mays L.) female inflorescence (ear) has an important impact on corn yield. However, the molecular mechanisms underlying maize ear development are poorly understood. Results: We profiled and analyzed gene expression of the maize ear at four developmental stages: elongation phase (I), spikelet differentiation phase (II), floret primordium differentiation phase (III), and floret organ differentiation phase (IV). Based on genome-wide profile analysis, we detected differential mRNA of maize genes. Among the ~6,800 differentially expressed genes (DEGs), 3,325 genes were differentially expressed in stage II, 3,765 genes in III, and 1,698 genes in IV, compared to its previous adjacent stages, respectively. Furthermore, some of DEGs were predicted to be potential candidates in maize ear development, such as AGAMOUS (GRMZM2G052890) and ATFP3 (GRMZM2G155281). Meanwhile, some genes were well-known annotated to the mutants during maize inflorescence development such as compact plant2 (ct2), zea AGAMOUS homolog1 (zag1), bearded ear (bde), and silky1 (si1). Some DEGs were predicted targets of microRNAs such as microRNA156. K-means clustering revealed that the DEGs showed 18 major expression patterns. Thirteen transcriptional factors from 10 families were differentially expressed across three comparisons of adjacent stages (II vs. I, III vs. II, IV vs. III). Antisense transcripts were widespread during all four stages, and might play important roles in maize ear development. Finally, we randomly selected 32 DEGs to validate their expression patterns using quantitative reverse-transcription polymerase chain reaction (qRT-PCR). The results were consistent with those from Solexa sequencing. Conclusions: DEGs technique had shown an advantage in detecting candidates, and some transcription factors during maize ear development. RT-PCR data were consistent with our sequencing data and supplied additional information on ear developmental processes. These results provide a molecular foundation for future research on maize ear development

The genome evolution and domestication of tropical fruit mango

Background: Mango is one of the world’s most important tropical fruits. It belongs to the family Anacardiaceae, which includes several other economically important species, notably cashew, sumac and pistachio from other genera. Many species in this family produce family-specific urushiols and related phenols, which can induce contact dermatitis. Results: We generate a chromosome-scale genome assembly of mango, providing a reference genome for the Anacardiaceae family. Our results indicate the occurrence of a recent whole-genome duplication (WGD) event in mango. Duplicated genes preferentially retained include photosynthetic, photorespiration, and lipid metabolic genes that may have provided adaptive advantages to sharp historical decreases in atmospheric carbon dioxide and global temperatures. A notable example of an extended gene family is the chalcone synthase (CHS) family of genes, and particular genes in this family show universally higher expression in peels than in flesh, likely for the biosynthesis of urushiols and related phenols. Genome resequencing reveals two distinct groups of mango varieties, with commercial varieties clustered with India germplasms and demonstrating allelic admixture, and indigenous varieties from Southeast Asia in the second group. Landraces indigenous in China formed distinct clades, and some showed admixture in genomes. Conclusions: Analysis of chromosome-scale mango genome sequences reveals photosynthesis and lipid metabolism are preferentially retained after a recent WGD event, and expansion of CHS genes is likely associated with urushiol biosynthesis in mango. Genome resequencing clarifies two groups of mango varieties, discovers allelic admixture in commercial varieties, and shows distinct genetic background of landraces

The Impact of Genetic Relationship and Linkage Disequilibrium on Genomic Selection

Genomic selection is a promising research area due to its practical application in breeding. In this study, impact of realized genetic relationship and linkage disequilibrium (LD) on marker density and training population size required was investigated and their impact on practical application was further discussed. This study is based on experimental data of two populations derived from the same two founder lines (B73, Mo17). Two populations were genotyped with different marker sets at different density: IBM Syn4 and IBM Syn10. A high-density marker set in Syn10 was imputed into the Syn4 population with low marker density. Seven different prediction scenarios were carried out with a random regression best linear unbiased prediction (RR-BLUP) model. The result showed that the closer the real genetic relationship between training and validation population, the fewer markers were required to reach a good prediction accuracy. Taken the short-term cost for consideration, relationship information is more valuable than LD information. Meanwhile, the result indicated that accuracies based on high LD between QTL and markers were more stable over generations, thus LD information would provide more robust prediction capacity in practical applications

Identification of miRNAs and their target genes in developing maize ears by combined small RNA and degradome sequencing

Background In plants, microRNAs (miRNAs) are endogenous ~22 nt RNAs that play important regulatory roles in many aspects of plant biology, including metabolism, hormone response, epigenetic control of transposable elements, and stress response. Extensive studies of miRNAs have been performed in model plants such as rice and Arabidopsis thaliana. In maize, most miRNAs and their target genes were analyzed and identified by clearly different treatments, such as response to low nitrate, salt and drought stress. However, little is known about miRNAs involved in maize ear development. The objective of this study is to identify conserved and novel miRNAs and their target genes by combined small RNA and degradome sequencing at four inflorescence developmental stages. Results We used deep-sequencing, miRNA microarray assays and computational methods to identify, profile, and describe conserved and non-conserved miRNAs at four ear developmental stages, which resulted in identification of 22 conserved and 21-maize-specific miRNA families together with their corresponding miRNA*. Comparison of miRNA expression in these developmental stages revealed 18 differentially expressed miRNA families. Finally, a total of 141 genes (251 transcripts) targeted by 102 small RNAs including 98 miRNAs and 4 ta-siRNAs were identified by genomic-scale high-throughput sequencing of miRNA cleaved mRNAs. Moreover, the differentially expressed miRNAs-mediated pathways that regulate the development of ears were discussed. Conclusions This study confirmed 22 conserved miRNA families and discovered 26 novel miRNAs in maize. Moreover, we identified 141 target genes of known and new miRNAs and ta-siRNAs. Of these, 72 genes (117 transcripts) targeted by 62 differentially expressed miRNAs may attribute to the development of maize ears. Identification and characterization of these important classes of regulatory genes in maize may improve our understanding of molecular mechanisms controlling ear development

Comparative genomics reveals adaptive evolution of Asian tapeworm in switching to a new intermediate host

Taenia saginata, Taenia solium and Taenia asiatica (beef, pork and Asian tapeworms, respectively) are parasitic flatworms of major public health and food safety importance. Among them, T. asiatica is a newly recognized species that split from T. saginata via an intermediate host switch ∼1.14 Myr ago. Here we report the 169- and 168-Mb draft genomes of T. saginata and T. asiatica. Comparative analysis reveals that high rates of gene duplications and functional diversifications might have partially driven the divergence between T. asiatica and T. saginata. We observe accelerated evolutionary rates, adaptive evolutions in homeostasis regulation, tegument maintenance and lipid uptakes, and differential/specialized gene family expansions in T. asiatica that may favour its hepatotropism in the new intermediate host. We also identify potential targets for developing diagnostic or intervention tools against human tapeworms. These data provide new insights into the evolution of Taenia parasites, particularly the recent speciation of T. asiatica

The genome evolution and domestication of tropical fruit mango

Background: Mango is one of the world’s most important tropical fruits. It belongs to the family Anacardiaceae, which includes several other economically important species, notably cashew, sumac and pistachio from other genera. Many species in this family produce family-specific urushiols and related phenols, which can induce contact dermatitis. Results: We generate a chromosome-scale genome assembly of mango, providing a reference genome for the Anacardiaceae family. Our results indicate the occurrence of a recent whole-genome duplication (WGD) event in mango. Duplicated genes preferentially retained include photosynthetic, photorespiration, and lipid metabolic genes that may have provided adaptive advantages to sharp historical decreases in atmospheric carbon dioxide and global temperatures. A notable example of an extended gene family is the chalcone synthase (CHS) family of genes, and particular genes in this family show universally higher expression in peels than in flesh, likely for the biosynthesis of urushiols and related phenols. Genome resequencing reveals two distinct groups of mango varieties, with commercial varieties clustered with India germplasms and demonstrating allelic admixture, and indigenous varieties from Southeast Asia in the second group. Landraces indigenous in China formed distinct clades, and some showed admixture in genomes. Conclusions: Analysis of chromosome-scale mango genome sequences reveals photosynthesis and lipid metabolism are preferentially retained after a recent WGD event, and expansion of CHS genes is likely associated with urushiol biosynthesis in mango. Genome resequencing clarifies two groups of mango varieties, discovers allelic admixture in commercial varieties, and shows distinct genetic background of landraces

Leveraging Fecal Bacterial Survey Data to Predict Colorectal Tumors

Colorectal cancer (CRC) ranks second in cancer-associated mortality and third in the incidence worldwide. Most of CRC follow adenoma-carcinoma sequence, and have more than 90% chance of survival if diagnosed at early stage. But the recommended screening by colonoscopy is invasive, expensive, and poorly adhered to. Recently, several studies reported that the fecal bacteria might provide non-invasive biomarkers for CRC and precancerous tumors. Therefore, we collected and uniformly re-analyzed these published fecal 16S rDNA sequencing datasets to verify the association and identify biomarkers to classify and predict colorectal tumors by random forest method. A total of 1674 samples (330 CRC, 357 advanced adenoma, 141 adenoma, and 846 control) from 7 studies were analyzed in this study. By random effects model and fixed effects model, we observed significant differences in alpha-diversity and beta-diversity between individuals with CRC and the normal colon, but not between adenoma and the normal. We identified various bacterial genera with significant odds ratios for colorectal tumors at different stages. Through building random forest model with 10-fold cross-validation as well as new test datasets, we classified individuals with CRC, advanced adenoma, adenoma and normal colon. All approaches obtained comparable performance at entire OTU level, entire genus level, and the common genus level as measured using AUC. When combined all samples, the AUC of random forest model based on 12 common genera reached 0.846 for CRC, although the predication performed poorly for advance adenoma and adenoma

Evolview v2: an online visualization and management tool for customized and annotated phylogenetic trees

Evolview is an online visualization and management tool for customized and annotated phylogenetic trees. It allows users to visualize phylogenetic trees in various formats, customize the trees through built-in functions and user-supplied datasets and export the customization results to publication-ready figures. Its 'dataset system' contains not only the data to be visualized on the tree, but also 'modifiers' that control various aspects of the graphical annotation. Evolview is a single-page application (like Gmail); its carefully designed interface allows users to upload, visualize, manipulate and manage trees and datasets all in a single webpage. Developments since the last public release include a modern dataset editor with keyword highlighting functionality, seven newly added types of annotation datasets, collaboration support that allows users to share their trees and datasets and various improvements of the web interface and performance. In addition, we included eleven new 'Demo' trees to demonstrate the basic functionalities of Evolview, and five new 'Showcase' trees inspired by publications to showcase the power of Evolview in producing publication-ready figures. Evolview is freely available at: http://www.evolgenius.info/evolview/

Chromosome level high-density integrated genetic maps improve the Pyrus bretschneideri ‘DangshanSuli’ v1.0 genome

Abstract Background Chromosomal level reference genomes provide a crucial foundation for genomics research such as genome-wide association studies (GWAS) and whole genome selection. The chromosomal-level sequences of both the European (Pyrus communis) and Chinese (P. bretschneideri) pear genomes have not been published in public databases so far. Results To anchor the scaffolds of P. bretschneideri ‘DangshanSuli’ (DS) v1.0 genome into pseudo-chromosomes, two genetic maps (MH and YM maps) were constructed using half sibling populations of Chinese pear crosses, ‘Mantianhong’ (MTH) × ‘Hongxiangsu’ (HXS) and ‘Yuluxiang’ (YLX) × MTH, from 345 and 162 seedlings, respectively, which were prepared for SNP discovery using genotyping-by-sequencing (GBS) technology. The MH and YM maps, each with 17 linkage groups (LGs), were constructed from 2606 and 2489 SNP markers and spanned 1847 and 1668 cM, respectively, with average marker intervals of 0.7. The two maps were further merged with a previously published genetic map (BD) based on the cross ‘Bayuehong’ (BYH) × ‘Dangshansuli’ (DS) to build a new integrated MH-YM-BD map. By using 7757 markers located on the integrated MH-YM-BD map, 898 scaffolds (400.57 Mb) of the DS v1.0 assembly were successfully anchored into 17 pseudo-chromosomes, accounting for 78.8% of the assembled genome size. About 88.31% of them (793 scaffolds) were directionally anchored with two or more markers on the pseudo-chromosomes. Furthermore, the errors in each pseudo-chromosome (especially 1, 5, 7 and 11) were manually corrected and pseudo-chromosomes 1, 5 and 7 were extended by adding 19, 12 and 14 scaffolds respectively in the newly constructed DS v1.1 genome. Synteny analyses revealed that the DS v1.1 genome had high collinearity with the apple genome, and the homologous fragments between pseudo-chromosomes were similar to those found in previous studies. Moreover, the red-skin trait of Asian pear was mapped to an identical locus as identified previously. Conclusions The accuracy of DS v1.1 genome was improved by using larger mapping populations and merged genetic map. With more than 400 MB anchored to 17 pseudo-chromosomes, the new DS v1.1 genome provides a critical tool that is essential for studies of pear genetics, genomics and molecular breeding