The subgenomes show asymmetric expression of alleles in hybrid lineages of Megalobrama amblycephala × Culter alburnus.

Hybridization drives rapid speciation by shaping novel genotypic and phenotypic profiles. Genomic incompatibility and transcriptome shock have been observed in hybrids, although this is rarer in animals than in plants. Using the newly sequenced genomes of the blunt snout bream (Megalobrama amblycephala [BSB]) and the topmouth culter (Culter alburnus [TC]), we focused on the sequence variation and gene expression changes in the reciprocal intergeneric hybrid lineages (F1-F3) of BSB × TC. A genome-wide transcriptional analysis identified 145-974 expressed recombinant genes in the successive generations of hybrid fish, suggesting the rapid emergence of allelic variation following hybridization. Some gradual changes of gene expression with additive and dominance effects and various cis and trans regulations were observed from F1 to F3 in the two hybrid lineages. These asymmetric patterns of gene expression represent the alternative strategies for counteracting deleterious effects of the subgenomes and improving adaptability of novel hybrids. Furthermore, we identified positive selection and additive expression patterns in transforming growth factor, beta 1b (tgfb1b), which may account for the morphological variations of the pharyngeal jaw in the two hybrid lineages. Our current findings provide insights into the evolution of vertebrate genomes immediately following hybridization.

The hard clam genome reveals massive expansion and diversification of inhibitors of apoptosis in Bivalvia.

BACKGROUND: Inhibitors of apoptosis (IAPs) are critical regulators of programmed cell death that are essential for development, oncogenesis, and immune and stress responses. However, available knowledge regarding IAP is largely biased toward humans and model species, while the distribution, function, and evolutionary novelties of this gene family remain poorly understood in many taxa, including Mollusca, the second most speciose phylum of Metazoa. RESULTS: Here, we present a chromosome-level genome assembly of an economically significant bivalve, the hard clam Mercenaria mercenaria, which reveals an unexpected and dramatic expansion of the IAP gene family to 159 members, the largest IAP gene repertoire observed in any metazoan. Comparative genome analysis reveals that this massive expansion is characteristic of bivalves more generally. Reconstruction of the evolutionary history of molluscan IAP genes indicates that most originated in early metazoans and greatly expanded in Bivalvia through both lineage-specific tandem duplication and retroposition, with 37.1% of hard clam IAPs located on a single chromosome. The expanded IAPs have been subjected to frequent domain shuffling, which has in turn shaped their architectural diversity. Further, we observed that extant IAPs exhibit dynamic and orchestrated expression patterns among tissues and in response to different environmental stressors. CONCLUSIONS: Our results suggest that sophisticated regulation of apoptosis enabled by the massive expansion and diversification of IAPs has been crucial for the evolutionary success of hard clam and other molluscan lineages, allowing them to cope with local environmental stresses. This study broadens our understanding of IAP proteins and expression diversity and provides novel resources for studying molluscan biology and IAP function and evolution.

Genomic Mechanisms of Physiological and Morphological Adaptations of Limestone Langurs to Karst Habitats.

Knowledge of the physiological and morphological evolution and adaptation of nonhuman primates is critical to understand hominin origins, physiological ecology, morphological evolution, and applications in biomedicine. Particularly, limestone langurs represent a direct example of adaptations to the challenges of exploiting a high calcium and harsh environment. Here, we report a de novo genome assembly (Tfra_2.0) of a male François's langur (Trachypithecus francoisi) with contig N50 of 16.3 Mb and resequencing data of 23 individuals representing five limestone and four forest langur species. Comparative genomics reveals evidence for functional evolution in genes and gene families related to calcium signaling in the limestone langur genome, probably as an adaptation to naturally occurring high calcium levels present in water and plant resources in karst habitats. The genomic and functional analyses suggest that a single point mutation (Lys1905Arg) in the α1c subunit of the L-type voltage-gated calcium channel Cav1.2 (CACNA1C) attenuates the inward calcium current into the cells in vitro. Population genomic analyses and RNA-sequencing indicate that EDNRB is less expressed in white tail hair follicles of the white-headed langur (T. leucocephalus) compared with the black-colored François's langur and hence might be responsible for species-specific differences in body coloration. Our findings contribute to a new understanding of gene-environment interactions and physiomorphological adaptative mechanisms in ecologically specialized primate taxa.

The pomegranate (Punica granatum L.) genome provides insights into fruit quality and ovule developmental biology.

Pomegranate (Punica granatum L.) has an ancient cultivation history and has become an emerging profitable fruit crop due to its attractive features such as the bright red appearance and the high abundance of medicinally valuable ellagitannin-based compounds in its peel and aril. However, the limited genomic resources have restricted further elucidation of genetics and evolution of these interesting traits. Here, we report a 274-Mb high-quality draft pomegranate genome sequence, which covers approximately 81.5% of the estimated 336-Mb genome, consists of 2177 scaffolds with an N50 size of 1.7 Mb and contains 30 903 genes. Phylogenomic analysis supported that pomegranate belongs to the Lythraceae family rather than the monogeneric Punicaceae family, and comparative analyses showed that pomegranate and Eucalyptus grandis share the paleotetraploidy event. Integrated genomic and transcriptomic analyses provided insights into the molecular mechanisms underlying the biosynthesis of ellagitannin-based compounds, the colour formation in both peels and arils during pomegranate fruit development, and the unique ovule development processes that are characteristic of pomegranate. This genome sequence provides an important resource to expand our understanding of some unique biological processes and to facilitate both comparative biology studies and crop breeding.

Identification, characterization, and utilization of single copy genes in 29 angiosperm genomes.

BACKGROUND: Single copy genes are common across angiosperm genomes. With the sufficiently high quality sequenced genomes, the identification of large-scale single copy genes among multiple species is possible. Although some characteristics have been reported, our study provides novel insights into single copy genes. RESULTS: We identified single copy genes across 29 angiosperm genomes. A significant negative correlation was found between the number of duplicate blocks and the number of single copy genes. We found that a considerable number of single copy genes are located in organelles, showing a preference for binding and catalytic activity. The analysis of effective number of codons (Nc) illustrates that single copy genes have a stronger codon bias than non-single copy genes in eudicots. The relative high expression level of single copy genes was partially confirmed by the RNA-seq data, rather than the Codon Adaptation Index (CAI). Unlike in most other species, a strongly negatively correlation occurs between Nc and GC3 among single copy genes in grass genomes. When compared to all non-single copy genes, single copy genes indicate more conservation (as indicated by Ka and Ks values). But our alternative splicing (AS) results reveal that selective constraints are weaker in single copy genes than in low copy family genes (1-10 in-paralogs) and stronger than high copy family genes (>10 in-paralogs). Using concatenated shared single copy genes, we obtained a well-resolved phylogenetic tree. With the addition of intron sequences, the branch support is improved, but striking incongruences are also evident. Therefore, it is noteworthy that inclusion of intron sequences seems more appropriate for the phylogenetic reconstruction at lower taxonomic levels. CONCLUSIONS: Our analysis provides insight into the evolutionary characteristics of single copy genes across 29 angiosperm genomes. The results suggest that there are key differences in evolutionary constraints between single copy genes and non-single copy genes. And to some extent, these evolutionary constraints show some species-specific differences, especially between eudicots and monocots. Our preliminary evidence also suggests that the concatenated shared single copy genes are well suited for use in resolving phylogenetic relationships.

[Molecular identification of 15 species of Ilex genus based on ITS sequence analysis ].

OBJECTIVE: To establish molecular method for identification of Ilex plants using the internal transcribed spacer(ITS) bar- code. METHODS: DNA was extracted using centrifugal columnar kits. The ITS sequences of all experimental samples were amplified and sequenced using a pair of universal primers. The alignment and the phylogenetic tree were performed in MEGA5. 0. RESULTS: The length of ITS sequence in each species was 759 bp,in which there were 216 sites discrepancy. CONCLUSION: DNA barcoding technology based on ITS sequence can be used to identify plants of Ilex genus.

Comparative genomics reveal shared genomic changes in syngnathid fishes and signatures of genetic convergence with placental mammals.

Syngnathids (seahorses, pipefishes and seadragons) exhibit an array of morphological innovations including loss of pelvic fins, a toothless tubular mouth and male pregnancy. They comprise two subfamilies: Syngnathinae and Nerophinae. Genomes of three Syngnathinae members have been analyzed previously. In this study, we have sequenced the genome of a Nerophinae member, the Manado pipefish (Microphis manadensis), which has a semi-enclosed brood pouch. Comparative genomic analysis revealed that the molecular evolutionary rate of the four syngnathids is higher than that of other teleosts. The loss of all but one P/Q-rich SCPP gene in the syngnathids suggests a role for the lost genes in dentin and enameloid formation in teleosts. Genome-wide comparison identified a set of 118 genes with parallel identical amino acid substitutions in syngnathids and placental mammals. Association of some of these genes with placental and embryonic development in mammals suggests a role for them in syngnathid pregnancy.

Chromosome-level genome assembly of the East Asian common octopus (Octopus sinensis) using PacBio sequencing and Hi-C technology.

The Cephalopoda are a group of highly diverse marine species in the phylum Mollusca, which are distributed worldwide. They have evolved some vertebrate-like biological traits and exhibit complicated behavioural repertoires. Thus, they are interesting species for studying the mechanisms of evolutionary convergence, innovational functional structures and evolutionary adaptation to a highly active, predatory lifestyle in diverse marine environments. Despite the evolutionary placement and biological significance of cephalopods, genomic data on these organisms remain limited. Here, we assembled a chromosome-level genome of a female East Asian common octopus (Octopus sinensis) by combining Pacific Bioscience (PacBio) single-molecule real-time sequencing, Illumina paired-end sequencing and Hi-C technology. An O. sinensis genome of 2.72 Gb was assembled from a total of 245.01 Gb high-quality PacBio sequences. The assembled genome represents 80.2% completeness (BUSCO) with a contig N50 of 490.36 Kb and a scaffold N50 of 105.89 Mb, showing a considerable improvement compared with other sequenced cephalopod genomes. Hi-C scaffolding of the genome resulted in the construction of 30 pseudochromosomes in Cephalopoda, representing 96.41% of the assembled sequences. The genome contained 42.26% repeat sequences and 5,245 noncoding RNAs. A total of 31,676 protein-coding genes were predicted, of which 82.73% were functionally annotated. The comparative genomic analysis identified 17,020 orthologous gene families, including 819 unique gene families and 629 expanded gene families. This genomic information will be an important molecular resource for further investigation of biological function and evolutionary adaptations in octopuses, and facilitate research into their population genetics and comparative evolution.

The Coix Genome Provides Insights into Panicoideae Evolution and Papery Hull Domestication.

Coix is a grass crop domesticated as early as the Neolithic era. It is still widely cultivated for both highly nutritional food and medicinal use. However, the genetic study and breeding of this crop are hindered by the lack of a sequenced genome. Here, we report de novo sequencing and assembly of the 1619-Mb genome of Coix, and annotation of 75.39% repeats and 39 629 protein-coding genes. Comparative genomics analysis showed that Coix is more closely related to sorghum than maize, but intriguingly only Coix and maize had a recent genome duplication event, which was not detected in sorghum. We further constructed a genetic map and mapped several important traits, especially the strength of hull. Selection of papery hull (thin: easy dehulling) from the stony hull (thick: difficult dehulling) in wild progenitors was a key step in Coix domestication. The papery hull makes seed easier to process and germinate. Anatomic and global transcriptome analysis revealed that the papery hull is a result of inhibition of cell division and wall biogenesis. We also successfully demonstrated that seed hull pressure resistance is controlled by two major quantitative trait loci (QTLs), which are associated with hull thickness and color, respectively. The two QTLs were further fine mapped within intervals of 250 kb and 146 kb, respectively. These resources provide a platform for evolutionary studies and will facilitate molecular breeding of this important crop.

Genomic and transcriptomic insights into molecular basis of sexually dimorphic nuptial spines in Leptobrachium leishanense.

Sexually dimorphic (SD) traits are important in sexual selection and species survival, yet the molecular basis remains elusive, especially in amphibians where SD traits have evolved repeatedly. We focus on the Leishan moustache toad (Leptobrachium leishanense), in which males develop nuptial spines on their maxillary skin. Here we report a 3.5 Gb genome assembly with a contig N50 of 1.93 Mb. We find a specific expansion of the intermediate filament gene family including numerous keratin genes. Within these genes, a cluster of duplicated hair keratin genes exhibits male-biased and maxillary skin-specific expression, suggesting a role in developing nuptial spines. We identify a module of coexpressed genes significantly associated with spine formation. In addition, we find several hormones likely to be involved in regulating spine development. This study not only presents a high-quality anuran genome but also provides a reference for studying skin-derived SD traits in amphibians.

Whole genomes and transcriptomes reveal adaptation and domestication of pistachio.

BACKGROUND: Pistachio (Pistacia vera), one of the most important commercial nut crops worldwide, is highly adaptable to abiotic stresses and is tolerant to drought and salt stresses. RESULTS: Here, we provide a draft de novo genome of pistachio as well as large-scale genome resequencing. Comparative genomic analyses reveal stress adaptation of pistachio is likely attributable to the expanded cytochrome P450 and chitinase gene families. Particularly, a comparative transcriptomic analysis shows that the jasmonic acid (JA) biosynthetic pathway plays an important role in salt tolerance in pistachio. Moreover, we resequence 93 cultivars and 14 wild P. vera genomes and 35 closely related wild Pistacia genomes, to provide insights into population structure, genetic diversity, and domestication. We find that frequent genetic admixture occurred among the different wild Pistacia species. Comparative population genomic analyses reveal that pistachio was domesticated about 8000 years ago and suggest that key genes for domestication related to tree and seed size experienced artificial selection. CONCLUSIONS: Our study provides insight into genetic underpinning of local adaptation and domestication of pistachio. The Pistacia genome sequences should facilitate future studies to understand the genetic basis of agronomically and environmentally related traits of desert crops.

Hypoxic and Cold Adaptation Insights from the Himalayan Marmot Genome.

The Himalayan marmot (Marmota himalayana) is a hibernating mammal that inhabits the high-elevation regions of the Himalayan mountains. Here we present a draft genome of the Himalayan marmot, with a total assembly length of 2.47 Gb. Phylogenetic analyses showed that the Himalayan marmot diverged from the Mongolian marmot approximately 1.98 million years ago. Transcriptional changes during hibernation included genes responsible for fatty acid metabolism in liver and genes involved in complement and coagulation cascades and stem cell pluripotency pathways in brain. Two selective sweep genes, Slc25a14 and ψAamp, showed apparent genotyping differences between low- and high-altitude populations. As a processed pseudogene, ψAamp may be biologically active to influence the stability of Aamp through competitive microRNA binding. These findings shed light on the molecular and genetic basis underlying adaptation to extreme environments in the Himalayan marmot.

Evolutionary analysis of three gibberellin oxidase genes in rice, Arabidopsis, and soybean.

GAs are plant hormones that play fundamental roles in plant growth and development. GA2ox, GA3ox, and GA20ox are three key enzymes in GA biosynthesis. These enzymes belong to the 2OG-Fe (II) oxygenase superfamily and are independently encoded by different gene families. To date, genome-wide comparative analyses of GA oxidases in plant species have not been thoroughly carried out. In the present work, 61 GA oxidase family genes from rice (Oryza sativa), Arabidopsis, and soybean (Glycine max) were identified and a full study of these genes including phylogenetic tree construction, gene structure, gene family expansion and analysis of functional motifs was performed. Based on phylogeny, most of the GA oxidases were divided into four subgroups that reflected functional classifications. Intron/intron average length of GA oxidase genes in rice analysis revealed that GA oxidase genes in rice experienced substantial evolutionary divergence. Segmental duplication events were mainly found in soybean genome. However, in rice and Arabidopsis, no single expansion pattern exhibited dominance, indicating that GA oxidase genes from these species might have been subjected to a more complex evolutionary mechanism. In addition, special functional motifs were discovered in GA20ox, GA3ox, and GA2ox, which suggested that different functional motifs are associated with differences in protein function. Taken together our results suggest that GA oxidase family genes have undergone divergent evolutionary routes, especially at the monocot-dicot split, with dynamic evolution occurring in Arabidopsis thaliana and soybean.