Exploring the evolutionary process of alkannin/shikonin O-acyltransferases by a reliable Lithospermum erythrorhizon genome

Increasing genome data are coming out. Genome size estimation plays an essential role in guiding genome assembly. Several months ago, other researchers were the first to publish a draft genome of the red gromwell (i.e. Lithospermum erythrorhizon). However, we considered that the genome size they estimated and assembled was incorrect. This study meticulously estimated the L. erythrorhizon genome size to should be ∼708.74 Mb and further provided a reliable genome version (size ≈ 693.34 Mb; contig_N50 length ≈ 238.08 Kb) to support our objection. Furthermore, according to our genome, we identified a gene family of the alkannin/shikonin O-acyltransferases (i.e. AAT/SAT) that catalysed enantiomer-specific acylations in the alkannin/shikonin biosynthesis (a characteristic metabolic pathway in L. erythrorhizon’s roots) and further explored its evolutionary process. The results indicated that the existing AAT/SAT were not generated from only one round of gene duplication but three rounds; after different rounds of gene duplication, the existing AAT/SAT and their recent ancestors were under positive selection at different amino acid sites. These suggested that a combined power from gene duplication plus positive selection plausibly propelled AAT/SAT’s functional differentiation in evolution.     

Chromosome-level genome assembly and characterization of Sophora Japonica

Sophora japonica is a medium-size deciduous tree belonging to Leguminosae family and famous for its high ecological, economic and medicinal value. Here, we reveal a draft genome of S. japonica, which was ∼511.49 Mb long (contig N50 size of 17.34 Mb) based on Illumina, Nanopore and Hi-C data. We reliably assembled 110 contigs into 14 chromosomes, representing 91.62% of the total genome, with an improved N50 size of 31.32 Mb based on Hi-C data. Further investigation identified 271.76 Mb (53.13%) of repetitive sequences and 31,000 protein-coding genes, of which 30,721 (99.1%) were functionally annotated. Phylogenetic analysis indicates that S. japonica separated from Arabidopsis thaliana and Glycine max ∼107.53 and 61.24 million years ago, respectively. We detected evidence of species-specific and common-legume whole-genome duplication events in S. japonica. We further found that multiple TF families (e.g. BBX and PAL) have expanded in S. japonica, which might have led to its enhanced tolerance to abiotic stress. In addition, S. japonica harbours more genes involved in the lignin and cellulose biosynthesis pathways than the other two species. Finally, population genomic analyses revealed no obvious differentiation among geographical groups and the effective population size continuously declined since 2 Ma. Our genomic data provide a powerful comparative framework to study the adaptation, evolution and active ingredients biosynthesis in S. japonica. More importantly, our high-quality S. japonica genome is important for elucidating the biosynthesis of its main bioactive components, and improving its production and/or processing.     

Hypoxia Adaptations in the Grey Wolf (Canis lupus chanco) from Qinghai-Tibet Plateau

The Tibetan grey wolf (Canis lupus chanco) occupies habitats on the Qinghai-Tibet Plateau, a high altitude (>3000 m) environment where low oxygen tension exerts unique selection pressure on individuals to adapt to hypoxic conditions. To identify genes involved in hypoxia adaptation, we generated complete genome sequences of nine Chinese wolves from high and low altitude populations at an average coverage of 25× coverage. We found that, beginning about 55,000 years ago, the highland Tibetan grey wolf suffered a more substantial population decline than lowland wolves. Positively selected hypoxia-related genes in highland wolves are enriched in the HIF signaling pathway (P = 1.57E-6), ATP binding (P = 5.62E-5), and response to an oxygen-containing compound (P≤5.30E-4). Of these positively selected hypoxia-related genes, three genes (EPAS1, ANGPT1, and RYR2) had at least one specific fixed non-synonymous SNP in highland wolves based on the nine genome data. Our re-sequencing studies on a large panel of individuals showed a frequency difference greater than 58% between highland and lowland wolves for these specific fixed non-synonymous SNPs and a high degree of LD surrounding the three genes, which imply strong selection. Past studies have shown that EPAS1 and ANGPT1 are important in the response to hypoxic stress, and RYR2 is involved in heart function. These three genes also exhibited significant signals of natural selection in high altitude human populations, which suggest similar evolutionary constraints on natural selection in wolves and humans of the Qinghai-Tibet Plateau.     

Population history and genomic signatures for high-altitude adaptation in Tibetan pigs

Background

The Tibetan pig is one of domestic animals indigenous to the Qinghai-Tibet Plateau. Several geographically isolated pig populations are distributed throughout the Plateau. It remained an open question if these populations have experienced different demographic histories and have evolved independent adaptive loci for the harsh environment of the Plateau. To address these questions, we herein investigated ~ 40,000 genetic variants across the pig genome in a broad panel of 678 individuals from 5 Tibetan geographic populations and 34 lowland breeds.

Results

Using a series of population genetic analyses, we show that Tibetan pig populations have marked genetic differentiations. Tibetan pigs appear to be 3 independent populations corresponding to the Tibetan, Gansu and Sichuan & Yunnan locations. Each population is more genetically similar to its geographic neighbors than to any of the other Tibetan populations. By applying a locus-specific branch length test, we identified both population-specific and -shared candidate genes under selection in Tibetan pigs. These genes, such as PLA2G12A, RGCC, C9ORF3, GRIN2B, GRID1 and EPAS1, are involved in high-altitude physiology including angiogenesis, pulmonary hypertension, oxygen intake, defense response and erythropoiesis. A majority of these genes have not been implicated in previous studies of highlanders and high-altitude animals.

Conclusion

Tibetan pig populations have experienced substantial genetic differentiation. Historically, Tibetan pigs likely had admixture with neighboring lowland breeds. During the long history of colonization in the Plateau, Tibetan pigs have developed a complex biological adaptation mechanism that could be different from that of Tibetans and other animals. Different Tibetan pig populations appear to have both distinct and convergent adaptive loci for the harsh environment of the Plateau.

Electronic supplementary material

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

Divergent adaptation to Qinghai-Tibetan Plateau implicated from transciptome study of Gymnocypris dobula and Schizothorax nukiangensis

The Schizothoracine fishes are widely distributed in the Qinghai-Tibetan Plateau (QTP) area and its peripheral regions, which provide a prime example of adaptation in highland aquatic environments. Recent progresses have revealed various genetic adaptations of these fishes by comparing to distantly related lowerland species, however, comparative studies on closely-related species of different altitudes are still lacking. In this study, we sequenced and annotated a primitive Schizothoracine fish Schizothorax nukiangensis Tsao and a highly specialized one Gymnocypris dobula. We performed evolutionary analyses to investigate the candidate genes and signaling pathways involved QTP highland adaptation in both Schizothoracine fishes. Analysis of the 11,007 one-copy orthologs to the primitive cyprinid species, Danio rerio, revealed that both G. dobula and S. nukiangensis showed elevated evolutionary rates. A large number of genes related to hypoxia, including genes involved metabolic processes and cardiovascular system development, exhibited signatures of positive selection in both Schizothoracine fishes, but very few positively selected genes were found overlapping among these Schizothoracines. Our results indicated divergent genetic adaptation to highland environment for aquatic species living in QTP.     

Comparative Genome Analyses Highlight Transposon-Mediated Genome Expansion and the Evolutionary Architecture of 3D Genomic Folding in Cotton

Transposable element (TE) amplification has been recognized as a driving force mediating genome size expansion and evolution, but the consequences for shaping 3D genomic architecture remains largely unknown in plants. Here, we report reference-grade genome assemblies for three species of cotton ranging 3-fold in genome size, namely Gossypium rotundifolium (K₂), G. arboreum (A₂), and G. raimondii (D₅), using Oxford Nanopore Technologies. Comparative genome analyses document the details of lineage-specific TE amplification contributing to the large genome size differences (K₂, 2.44 Gb; A₂, 1.62 Gb; D₅, 750.19 Mb) and indicate relatively conserved gene content and synteny relationships among genomes. We found that approximately 17% of syntenic genes exhibit chromatin status change between active (“A”) and inactive (“B”) compartments, and TE amplification was associated with the increase of the proportion of A compartment in gene regions (∼7,000 genes) in K₂ and A₂ relative to D₅. Only 42% of topologically associating domain (TAD) boundaries were conserved among the three genomes. Our data implicate recent amplification of TEs following the formation of lineage-specific TAD boundaries. This study sheds light on the role of transposon-mediated genome expansion in the evolution of higher-order chromatin structure in plants.     

Enrichment and physiological characterization of a novel comammox Nitrospira indicates ammonium inhibition of complete nitrification

The recent discovery of bacteria within the genus Nitrospira capable of complete ammonia oxidation (comammox) demonstrated that the sequential oxidation of ammonia to nitrate via nitrite can also be performed within a single bacterial cell. Although comammox Nitrospira exhibit a wide distribution in natural and engineered ecosystems, information on their physiological properties is scarce due to the limited number of cultured representatives. Additionally, most available genomic information is derived from metagenomic sequencing and high-quality genomes of Nitrospira in general are limited. In this study, we obtained a high (90%) enrichment of a novel comammox species, tentatively named “Candidatus Nitrospira kreftii”, and performed a detailed genomic and physiological characterization. The complete genome of “Ca. N. kreftii” allowed reconstruction of its basic metabolic traits. Similar to Nitrospira inopinata, the enrichment culture exhibited a very high ammonia affinity (K_{m(app)_NH3} ≈ 0.040 ± 0.01 µM), but a higher nitrite affinity (K_{m(app)_NO2}- = 12.5 ± 4.0 µM), indicating an adaptation to highly oligotrophic environments. Furthermore, we observed partial inhibition of ammonia oxidation at ammonium concentrations as low as 25 µM. This inhibition of “Ca. N. kreftii” indicates that differences in ammonium tolerance rather than affinity could potentially be a niche determining factor for different comammox Nitrospira.Subject terms: Bacterial genomics, Environmental microbiology, Bacterial physiology     

A spinach genome assembly with remarkable completeness,and its use for rapid identification of candidate genes for agronomic traits

Spinach (Spinacia oleracea) is grown as a nutritious leafy vegetable worldwide. To accelerate spinach breeding efficiency, a high-quality reference genome sequence with great completeness and continuity is needed as a basic infrastructure. Here, we used long-read and linked-read technologies to construct a de novo spinach genome assembly, designated SOL_r1.1, which was comprised of 287 scaffolds (total size: 935.7 Mb; N₅₀ = 11.3 Mb) with a low proportion of undetermined nucleotides (Ns = 0.34%) and with high gene completeness (BUSCO complete 96.9%). A genome-wide survey of resistance gene analogues identified 695 genes encoding nucleotide-binding site domains, receptor-like protein kinases, receptor-like proteins and transmembrane-coiled coil domains. Based on a high-density double-digest restriction-site associated DNA sequencing-based linkage map, the genome assembly was anchored to six pseudomolecules representing ∼73.5% of the whole genome assembly. In addition, we used SOL_r1.1 to identify quantitative trait loci for bolting timing and fruit/seed shape, which harbour biologically plausible candidate genes, such as homologues of the FLOWERING LOCUS T and EPIDERMAL PATTERNING FACTOR-LIKE genes. The new genome assembly, SOL_r1.1, will serve as a useful resource for identifying loci associated with important agronomic traits and for developing molecular markers for spinach breeding/selection programs.     

Complete telomere-to-telomere de novo assembly of the Plasmodium falciparum genome through long-read (>11?kb), single molecule,real-time sequencing

The application of next-generation sequencing to estimate genetic diversity of Plasmodium falciparum, the most lethal malaria parasite, has proved challenging due to the skewed AT-richness [∼80.6% (A + T)] of its genome and the lack of technology to assemble highly polymorphic subtelomeric regions that contain clonally variant, multigene virulence families (Ex: var and rifin). To address this, we performed amplification-free, single molecule, real-time sequencing of P. falciparum genomic DNA and generated reads of average length 12 kb, with 50% of the reads between 15.5 and 50 kb in length. Next, using the Hierarchical Genome Assembly Process, we assembled the P. falciparum genome de novo and successfully compiled all 14 nuclear chromosomes telomere-to-telomere. We also accurately resolved centromeres [∼90–99% (A + T)] and subtelomeric regions and identified large insertions and duplications that add extra var and rifin genes to the genome, along with smaller structural variants such as homopolymer tract expansions. Overall, we show that amplification-free, long-read sequencing combined with de novo assembly overcomes major challenges inherent to studying the P. falciparum genome. Indeed, this technology may not only identify the polymorphic and repetitive subtelomeric sequences of parasite populations from endemic areas but may also evaluate structural variation linked to virulence, drug resistance and disease transmission.     

Rheumatoid arthritis response to treatment across IgG1 allotype – anti-TNF incompatibility: a case-only study

IntroductionWe have hypothesized that incompatibility between the G1m genotype of the patient and the G1m1 and G1m17 allotypes carried by infliximab (INX) and adalimumab (ADM) could decrease the efficacy of these anti-tumor necrosis factor (anti-TNF) antibodies in the treatment of rheumatoid arthritis (RA).MethodsThe G1m genotypes were analyzed in three collections of patients with RA totaling 1037 subjects. The first, used for discovery, comprised 215 Spanish patients. The second and third were successively used for replication. They included 429 British and Greek patients and 393 Spanish and British patients, respectively. Two outcomes were considered: change in the Disease Activity Score in 28 joint (ΔDAS28) and the European League Against Rheumatism (EULAR) response criteria.ResultsAn association between less response to INX and incompatibility of the G1m1,17 allotype was found in the discovery collection at 6 months of treatment (P = 0.03). This association was confirmed in the replications (P = 0.02 and 0.08, respectively) leading to a global association (P = 0.001) that involved a mean difference in ΔDAS28 of 0.4 units between compatible and incompatible patients (2.3 ± 1.5 in compatible patients vs. 1.9 ± 1.5 in incompatible patients) and an increase in responders and decrease in non-responders according to the EULAR criteria (P = 0.03). A similar association was suggested for patients treated with ADM in the discovery collection, but it was not supported by replication.ConclusionsOur results suggest that G1m1,17 allotypes are associated with response to INX and could aid improved therapeutic targeting in RA.

Electronic supplementary material

The online version of this article (doi:10.1186/s13075-015-0571-z) contains supplementary material, which is available to authorized users.     

Preliminary investigation on otolith microchemistry of naked carp (Gymnocypris przewalskii) in Lake Qinghai, China

The paper represents the first otolith study on trace element (Sr, Mg) concentrations of naked carp (Gymnocypris przewalskii), an endemic and endangered migratory species in high-altitude Lake Qinghai and its tributaries of China, using electron probe microanalysis (EPMA). Results from life-history transect analyses and X-ray intensity maps suggested that the general reference criterion for discrimination of the salinity habitats of freshwater, brackish and seawater by the otolith Sr:Ca ratio of diadromous fish cannot be used as a reliable indicator for reconstructing the migration history of naked carp between different habitats. On the contrary, Mg concentrations and Mg:Ca ratios were found inconsistent in all otolith samples of G. przewalskii. The otolith regions with higher Mg concentrations (greenish) might coincide with the period that yolk provides nutrition to G. przewalskii. The lower Mg:Ca ratios in most regions of otolith showed that Mg in G. przewalskii might be regulated to that low level during its lifespan in Lake Qinghai. This suggests that the distinct EPMA approach can provide new visualization insight into the ontogenetic and migration patterns of G. przewalskii.     

Rare instances of haploid inducer DNA in potato dihaploids and ploidy-dependent genome instability

In cultivated tetraploid potato (Solanum tuberosum), reduction to diploidy (dihaploidy) allows for hybridization to diploids and introgression breeding and may facilitate the production of inbreds. Pollination with haploid inducers (HIs) yields maternal dihaploids, as well as triploid and tetraploid hybrids. Dihaploids may result from parthenogenesis, entailing the development of embryos from unfertilized eggs, or genome elimination, entailing missegregation and the loss of paternal chromosomes. A sign of genome elimination is the occasional persistence of HI DNA in some dihaploids. We characterized the genomes of 919 putative dihaploids and 134 hybrids produced by pollinating tetraploid clones with three HIs: IVP35, IVP101, and PL-4. Whole-chromosome or segmental aneuploidy was observed in 76 dihaploids, with karyotypes ranging from 2n = 2x − 1 = 23 to 2n = 2x + 3 = 27. Of the additional chromosomes in 74 aneuploids, 66 were from the non-inducer parent and 8 from the inducer parent. Overall, we detected full or partial chromosomes from the HI parent in 0.87% of the dihaploids, irrespective of parental genotypes. Chromosomal breaks commonly affected the paternal genome in the dihaploid and tetraploid progeny, but not in the triploid progeny, correlating instability to sperm ploidy and to haploid induction. The residual HI DNA discovered in the progeny is consistent with genome elimination as the mechanism of haploid induction.

A large potato progeny population produced by crossing tetraploid cultivated clones to diploid Phureja lines displays rare instances of haploid inducer chromosomes, which are frequently damaged.