Dongming Fang

The draft genomes of soft-shell turtle and green sea turtle yield insights into the development and evolution of the turtle-specific body plan

The unique anatomical features of turtles have raised unanswered questions about the origin of their unique body plan. We generated and analyzed draft genomes of the soft-shell turtle (Pelodiscus sinensis) and the green sea turtle (Chelonia mydas); our results indicated the close relationship of the turtles to the bird-crocodilian lineage, from which they split ∼267.9–248.3 million years ago (Upper Permian to Triassic). We also found extensive expansion of olfactory receptor genes in these turtles. Embryonic gene expression analysis identified an hourglass-like divergence of turtle and chicken embryogenesis, with maximal conservation around the vertebrate phylotypic period, rather than at later stages that show the amniote-common pattern. Wnt5a expression was found in the growth zone of the dorsal shell, supporting the possible co-option of limb-associated Wnt signaling in the acquisition of this turtle-specific novelty. Our results suggest that turtle evolution was accompanied by an unexpectedly conservative vertebrate phylotypic period, followed by turtle-specific repatterning of development to yield the novel structure of the shell.The unique anatomical features of turtles have raised unanswered questions about the origin of their unique body plan. We generated and analyzed draft genomes of the soft-shell turtle (Pelodiscus sinensis) and the green sea turtle (Chelonia mydas); our results indicated the close relationship of the turtles to the bird-crocodilian lineage, from which they split ∼267.9-248.3 million years ago (Upper Permian to Triassic). We also found extensive expansion of olfactory receptor genes in these turtles. Embryonic gene expression analysis identified an hourglass-like divergence of turtle and chicken embryogenesis, with maximal conservation around the vertebrate phylotypic period, rather than at later stages that show the amniote-common pattern. Wnt5a expression was found in the growth zone of the dorsal shell, supporting the possible co-option of limb-associated Wnt signaling in the acquisition of this turtle-specific novelty. Our results suggest that turtle evolution was accompanied by an unexpectedly conservative vertebrate phylotypic period, followed by turtle-specific repatterning of development to yield the novel structure of the shell.

Genomes of the rice pest brown planthopper and its endosymbionts reveal complex complementary contributions for host adaptation

BackgroundThe brown planthopper, Nilaparvata lugens, the most destructive pest of rice, is a typical monophagous herbivore that feeds exclusively on rice sap, which migrates over long distances. Outbreaks of it have re-occurred approximately every three years in Asia. It has also been used as a model system for ecological studies and for developing effective pest management. To better understand how a monophagous sap-sucking arthropod herbivore has adapted to its exclusive host selection and to provide insights to improve pest control, we analyzed the genomes of the brown planthopper and its two endosymbionts.ResultsWe describe the 1.14 gigabase planthopper draft genome and the genomes of two microbial endosymbionts that permit the planthopper to forage exclusively on rice fields. Only 40.8% of the 27,571 identified Nilaparvata protein coding genes have detectable shared homology with the proteomes of the other 14 arthropods included in this study, reflecting large-scale gene losses including in evolutionarily conserved gene families and biochemical pathways. These unique genomic features are functionally associated with the animal’s exclusive plant host selection. Genes missing from the insect in conserved biochemical pathways that are essential for its survival on the nutritionally imbalanced sap diet are present in the genomes of its microbial endosymbionts, which have evolved to complement the mutualistic nutritional needs of the host.ConclusionsOur study reveals a series of complex adaptations of the brown planthopper involving a variety of biological processes, that result in its highly destructive impact on the exclusive host rice. All these findings highlight potential directions for effective pest control of the planthopper.

Genomic insights into salt adaptation in a desert poplar

Despite the high economic and ecological importance of forests, our knowledge of the genomic evolution of trees under salt stress remains very limited. Here we report the genome sequence of the desert poplar, Populus euphratica, which exhibits high tolerance to salt stress. Its genome is very similar and collinear to that of the closely related mesophytic congener, P. trichocarpa. However, we find that several gene families likely to be involved in tolerance to salt stress contain significantly more gene copies within the P. euphratica lineage. Furthermore, genes showing evidence of positive selection are significantly enriched in functional categories related to salt stress. Some of these genes, and others within the same categories, are significantly upregulated under salt stress relative to their expression in another salt-sensitive poplar. Our results provide an important background for understanding tree adaptation to salt stress and facilitating the genetic improvement of cultivated poplars for saline soils.

The Sinocyclocheilus cavefish genome provides insights into cave adaptation

BackgroundAn emerging cavefish model, the cyprinid genus Sinocyclocheilus, is endemic to the massive southwestern karst area adjacent to the Qinghai-Tibetan Plateau of China. In order to understand whether orogeny influenced the evolution of these species, and how genomes change under isolation, especially in subterranean habitats, we performed whole-genome sequencing and comparative analyses of three species in this genus, S. grahami, S. rhinocerous and S. anshuiensis. These species are surface-dwelling, semi-cave-dwelling and cave-restricted, respectively.ResultsThe assembled genome sizes of S. grahami, S. rhinocerous and S. anshuiensis are 1.75 Gb, 1.73 Gb and 1.68 Gb, respectively. Divergence time and population history analyses of these species reveal that their speciation and population dynamics are correlated with the different stages of uplifting of the Qinghai-Tibetan Plateau. We carried out comparative analyses of these genomes and found that many genetic changes, such as gene loss (e.g. opsin genes), pseudogenes (e.g. crystallin genes), mutations (e.g. melanogenesis-related genes), deletions (e.g. scale-related genes) and down-regulation (e.g. circadian rhythm pathway genes), are possibly associated with the regressive features (such as eye degeneration, albinism, rudimentary scales and lack of circadian rhythms), and that some gene expansion (e.g. taste-related transcription factor gene) may point to the constructive features (such as enhanced taste buds) which evolved in these cave fishes.ConclusionAs the first report on cavefish genomes among distinct species in Sinocyclocheilus, our work provides not only insights into genetic mechanisms of cave adaptation, but also represents a fundamental resource for a better understanding of cavefish biology.

An improved 2b-RAD approach (I2b-RAD) offering genotyping tested by a rice (Oryza sativa L.) F2 population

Background2b-RAD (type IIB endonucleases restriction-site associated DNA) approach was invented by Wang in 2012 and proven as a simple and flexible method for genome-wide genotyping. However, there is still plenty of room for improvement for the existent 2b-RAD approach. Firstly, it doesn’t include the samples pooling in library preparation as other reduced representation libraries. Secondly, the information of 2b-RAD tags, such as tags numbers and distributions, in most of species are unknown. The purposes of the research are to improve a new 2b-RAD approach which possesses samples pooling, moreover to figure out the characteristic and application potentiality of 2b-RAD tags by bioinformatics analysis.ResultsTwelve adapter1 and an adapter2 were designed. A library approach comprising digestion, ligation, pooling, PCR and size selection were established. For saving costs, we used non-phosphorylated adapters and indexed PCR primers. A F2 population of rice (Oryza sativa .L) was genotyped to validate the new approach. On average, 2000332 high quality reads of each sample were obtained with high evenness. Totally 3598 markers containing 3804 SNPs were discovered and the missing rate was 18.9%. A genetic linkage map of 1385 markers was constructed and 92% of the markers’ orders in the genetic map were in accordance with the orders in chromosomes. Meanwhile, the bioinformatics simulation in 20 species showed that the BsaXI had the most widespread recognition sites, indicating that 2b-RAD tags had a powerful application potentiality for high density genetic map. Using modified adapters with a fix base in 3′end, 2b-RAD was also fit for QTL studies with low costs.ConclusionsAn improved 2b-RAD genotyping approach was established in this research and named as I2b-RAD. The method was a simple, fast, cost-effective and multiplex sequencing library approach. It could be adjusted by selecting different enzymes and adapters to fit for alternative uses including chromosomes assembly, QTL fine mapping and even natural population analysis.

Multiple IgH Isotypes Including IgD, Subclasses of IgM, and IgY Are Expressed in the Common Ancestors of Modern Birds

Although evolutionarily just as ancient as IgM, it has been thought for many years that IgD is not present in birds. Based on the recently sequenced genomes of 48 bird species as well as high-throughput transcriptome sequencing of immune-related tissues, we demonstrate in this work that the ostrich (Struthio camelus) possesses a functional δ gene that encodes a membrane-bound IgD H chain with seven CH domains. Furthermore, δ sequences were clearly identified in many other bird species, demonstrating that the δ gene is widely distributed among birds and is only absent in certain bird species. We also show that the ostrich possesses two μ genes (μ1, μ2) and two υ genes (υ1, υ2), in addition to the δ and α genes. Phylogenetic analyses suggest that subclass diversification of both the μ and υ genes occurred during the early stages of bird evolution, after their divergence from nonavian reptiles. Although the positions of the two υ genes are unknown, physical mapping showed that the remaining genes are organized in the order μ1-δ-α-μ2, with the α gene being inverted relative to the others. Together with previous studies, our data suggest that birds and nonavian reptile species most likely shared a common ancestral IgH gene locus containing a δ gene and an inverted α gene. The δ gene was then evolutionarily lost in selected birds, whereas the α gene lost in selected nonavian reptiles. The data obtained in this study provide significant insights into the understanding of IgH gene evolution in tetrapods.

Genome-wide Target Enrichment-aided Chip Design: a 66 K SNP Chip for Cashmere Goat

Compared with the commercially available single nucleotide polymorphism (SNP) chip based on the Bead Chip technology, the solution hybrid selection (SHS)-based target enrichment SNP chip is not only design-flexible, but also cost-effective for genotype sequencing. In this study, we propose to design an animal SNP chip using the SHS-based target enrichment strategy for the first time. As an update to the international collaboration on goat research, a 66 K SNP chip for cashmere goat was created from the whole-genome sequencing data of 73 individuals. Verification of this 66 K SNP chip with the whole-genome sequencing data of 436 cashmere goats showed that the SNP call rates was between 95.3% and 99.8%. The average sequencing depth for target SNPs were 40X. The capture regions were shown to be 200 bp that flank target SNPs. This chip was further tested in a genome-wide association analysis of cashmere fineness (fiber diameter). Several top hit loci were found marginally associated with signaling pathways involved in hair growth. These results demonstrate that the 66 K SNP chip is a useful tool in the genomic analyses of cashmere goats. The successful chip design shows that the SHS-based target enrichment strategy could be applied to SNP chip design in other species.