Introgression of Alectoris chukar genes into a Spanish wild Alectoris rufa population

In order to detect introgression of other Alectoris genus species into wild populations of Spanish Alectoris rufa, we studied a sample of 93 red-legged partridges (supposed to be A. rufa) captured in the island of Majorca. A set of 31 chukar partridges (Alectoris chukar) from Cyprus and 33 red-legged partridges (A. rufa) from one Spanish farm were also studied to provide suitable populations for comparison. Factorial correspondence analysis on microsatellite genotypes supported a clear distinction of birds from Cyprus, whereas partridges from Majorca and the Spanish farm overlapped in a wide area. The existence of A. chukar mitochondrial DNA in 16 individuals from Majorca indicated introgression into their maternal lineage even if their phenotypes were not different from A. rufa. Bayesian inference based on microsatellite analysis indicated a noticeable degree of genetic proximity to A. chukar only for one of these hybrids.     

Parasites, condition, immune responsiveness and carotenoid-based ornamentation in male red-legged partridge Alectoris rufa

Sexual ornaments might indicate better condition, fewer parasites or a greater immune responsiveness. Carotenoid-based ornaments are common sexual signals of birds and often influence mate choice. Skin or beaks pigmented by carotenoids can change colour rapidly, and could be particularly useful as honest indicators of an individual's current condition and/or health. This is because carotenoids must be acquired through diet and/or allocation for ornamental coloration might be to the detriment of self-maintenance needs. Here, we investigated whether the carotenoid-based coloration of eye rings and beak of male red-legged partridges Alectoris rufa predicted condition (mass corrected for size), parasite load (more specifically infection by coccidia, a main avian intestinal parasite) or a greater immune responsiveness (swelling response to a plant lectin, phytohaemagluttinin, or PHA). Redness of beak and eye rings positively correlated with plasma carotenoid levels. Also, males in better condition had fewer coccidia, more circulating carotenoids and a greater swelling response to PHA. Carotenoid-based ornamentation predicted coccidia abundance and immune responsiveness (redder males had fewer coccidia and greater swelling response to PHA), but was only weakly positively related to condition. Thus, the carotenoid pigmentation of beak and eye rings reflected the current health status of individuals. Our results are consistent with the hypothesis that allocation trade-offs (carotenoid use for ornamentation versus parasite defence needs) might ensure reliable carotenoid-based signalling.     

SMRT sequencing only de novo assembly of the sugar beet (Beta vulgaris) chloroplast genome

Background

Third generation sequencing methods, like SMRT (Single Molecule, Real-Time) sequencing developed by Pacific Biosciences, offer much longer read length in comparison to Next Generation Sequencing (NGS) methods. Hence, they are well suited for de novo- or re-sequencing projects. Sequences generated for these purposes will not only contain reads originating from the nuclear genome, but also a significant amount of reads originating from the organelles of the target organism. These reads are usually discarded but they can also be used for an assembly of organellar replicons. The long read length supports resolution of repetitive regions and repeats within the organelles genome which might be problematic when just using short read data. Additionally, SMRT sequencing is less influenced by GC rich areas and by long stretches of the same base.

Results

We describe a workflow for a de novo assembly of the sugar beet (Beta vulgaris ssp. vulgaris) chloroplast genome sequence only based on data originating from a SMRT sequencing dataset targeted on its nuclear genome. We show that the data obtained from such an experiment are sufficient to create a high quality assembly with a higher reliability than assemblies derived from e.g. Illumina reads only. The chloroplast genome is especially challenging for de novo assembling as it contains two large inverted repeat (IR) regions. We also describe some limitations that still apply even though long reads are used for the assembly.

Conclusions

SMRT sequencing reads extracted from a dataset created for nuclear genome (re)sequencing can be used to obtain a high quality de novo assembly of the chloroplast of the sequenced organism. Even with a relatively small overall coverage for the nuclear genome it is possible to collect more than enough reads to generate a high quality assembly that outperforms short read based assemblies. However, even with long reads it is not always possible to clarify the order of elements of a chloroplast genome sequence reliantly which we could demonstrate with Fosmid End Sequences (FES) generated with Sanger technology. Nevertheless, this limitation also applies to short read sequencing data but is reached in this case at a much earlier stage during finishing.

Electronic supplementary material

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

The expression of melanin-based plumage is separately modulated by exogenous oxidative stress and a melanocortin

Melanin-based traits involved in animal communication have been traditionally viewed as occurring under strict genetic control. However, it is generally accepted that both genetic and environmental factors influence melanin production. Medical studies suggest that, among environmental factors influencing melanization, oxidative stress could play a relevant role. On the other hand, genetic control would be exerted by the melanocortin system, and particularly by the alpha-melanocyte-stimulating hormone (α-MSH), which triggers the production of eumelanins (black pigments). To determine how the melanocortin system and an exogenous source of oxidative stress interact in the expression of melanin-based plumage, developing red-legged partridges (Alectoris rufa) were manipulated. Some partridges were injected with α-MSH, while other birds received a pro-oxidant molecule (diquat) in drinking water. Controls and birds receiving both treatments were also studied. Both α-MSH- and diquat-treated individuals presented larger eumelanin-based traits than controls, but α-MSH+diquat-treated birds showed the largest traits, suggesting that oxidative stress and melanocortins promote additive but independent effects. Diquat also induced a decline in the level of a key intracellular antioxidant (glutathione), which is associated with high expression of eumelanin-based signals in other bird species. Some scenarios for the evolution of melanin-based traits in relation to oxidative stress are proposed.     

Mitochondrial DNA and microsatellite markers evidence a different pattern of hybridization in red-legged partridge (Alectoris rufa) populations from NW Italy

Genetic introgression with exotic genomes represents a major conservation concern for the red-legged partridge (Alectoris rufa, Phasianidae). In particular, massive releases of chukar partridges (Alectoris chukar) and/or red-legged × chukar partridge hybrids for hunting purposes have cast serious doubts on the Italian present-day occurrence of wild A. rufa populations not affected by introgressive hybridization. This study investigates the genetic structure of red-legged partridges populations in two ecologically different areas in Northern Italy. Analysis of maternal mitochondrial DNA and biparental microsatellite markers excluded the presence of hybridization in a typical agricultural habitat where hunting and release of reared birds are strictly banned. By contrast, signs of chukar introgression were detected in a perifluvial habitat unusual for the red-legged partridge in Italy. The present study documents the first red-legged partridge population with no genetic evidences of recent chukar introgression presently living in Italy. We recommend that urgent conservation actions are taken to preserve the genetic integrity of this population from the risk of hybridization with farm-reared birds and to support its long-term conservation.     

New single nucleotide polymorphisms in Alectoris identified using chicken genome information allow Alectoris introgression detection

Using the chicken genome, 114 polymorphisms (109 SNPs and 5 INDELs) were identified in the Alectoris genus by polymerase chain reaction-single strand conformation polymorphism. Using these, a panel of SNPs is described, which allows easy detection of introgression of Alectoris chukar in wild Alectoris rufa populations, when used with a primer extension protocol. The selected polymorphisms were genotyped and their allelic frequencies estimated on 98 A. rufa partridges sampled from nonrestocking Spanish areas, and 63 A. chukar partridges from Greek and Spanish farms. Power calculations to determine an optimum subset of markers for a given significance level were performed.     

Contrasted effects of an oxidative challenge and α-melanocyte-stimulating hormone on cellular immune responsiveness: an experiment with red-legged partridges Alectoris rufa

Oxidative stress is increasingly recognized as a key selective force shaping evolutionary trade-offs. One such trade-off involves investing in immunity versus combating oxidative stress. While there is broad evidence that mounting an immune response causes increased oxidative stress, the effect that increased oxidative stress during development has at a later stage on immune responsiveness remains little known. The production of melanin-based coloration in vertebrates is influenced by oxidative stress and by hormones, such as the alpha-melanocyte-stimulating hormone (α-MSH). Oxidative stress could impair immunity, and this might be a cost associated with the production of melanin traits. α-MSH has immunomodulatory effects, with most evidence pointing towards an improvement of immunity (improved pro-inflammatory activity). Here, we investigated the effects of an oxidative challenge (exposure to a pro-oxidant compound, diquat) and of experimentally elevated α-MSH on the cell-mediated immune responses (CMIR) of growing young (1 month old) red-legged partridges Alectoris rufa in captivity. CMIR were assessed in response to primary and secondary challenges with phytohemagglutinin (PHA). We specifically tested whether an oxidative challenge during growth and development had a delayed effect (4 months after exposure) on immunity. We found that the diquat treatment did not affect primary CMIR, but significantly reduced secondary CMIR. Elevated α-MSH increased primary CMIR in males, but not in females. Our experimental results are consistent with a trade-off between investing in activities that generate oxidative stress (e.g., growth, reproduction, production of ornaments) versus investing in immunity, and shed new lights onto the inter-relationships between immunity, oxidative stress and the expression of melanin-based coloration in vertebrates, revealing a novel, delayed physiological cost that can contribute to ensuring honest signaling.     

Linkage mapping and comparative genomics using next-generation RAD sequencing of a non-model organism

Restriction-site associated DNA (RAD) sequencing is a powerful new method for targeted sequencing across the genomes of many individuals. This approach has broad potential for genetic analysis of non-model organisms including genotype-phenotype association mapping, phylogeography, population genetics and scaffolding genome assemblies through linkage mapping. We constructed a RAD library using genomic DNA from a Plutella xylostella (diamondback moth) backcross that segregated for resistance to the insecticide spinosad. Sequencing of 24 individuals was performed on a single Illumina GAIIx lane (51 base paired-end reads). Taking advantage of the lack of crossing over in homologous chromosomes in female Lepidoptera, 3,177 maternally inherited RAD alleles were assigned to the 31 chromosomes, enabling identification of the spinosad resistance and W/Z sex chromosomes. Paired-end reads for each RAD allele were assembled into contigs and compared to the genome of Bombyx mori (n = 28) using BLAST, revealing 28 homologous matches plus 3 expected fusion/breakage events which account for the difference in chromosome number. A genome-wide linkage map (1292 cM) was inferred with 2,878 segregating RAD alleles inherited from the backcross father, producing chromosome and location specific sequenced RAD markers. Here we have used RAD sequencing to construct a genetic linkage map de novo for an organism that has no previous genome data. Comparative analysis of P. xyloxtella linkage groups with B. mori chromosomes shows for the first time, genetic synteny appears common beyond the Macrolepidoptera. RAD sequencing is a powerful system capable of rapidly generating chromosome specific data for non-model organisms.     

Whole Genome Complete Resequencing of Bacillus subtilis Natto by Combining Long Reads with High-Quality Short Reads

De novo microbial genome sequencing reached a turning point with third-generation sequencing (TGS) platforms, and several microbial genomes have been improved by TGS long reads. Bacillus subtilis natto is closely related to the laboratory standard strain B. subtilis Marburg 168, and it has a function in the production of the traditional Japanese fermented food “natto.” The B. subtilis natto BEST195 genome was previously sequenced with short reads, but it included some incomplete regions. We resequenced the BEST195 genome using a PacBio RS sequencer, and we successfully obtained a complete genome sequence from one scaffold without any gaps, and we also applied Illumina MiSeq short reads to enhance quality. Compared with the previous BEST195 draft genome and Marburg 168 genome, we found that incomplete regions in the previous genome sequence were attributed to GC-bias and repetitive sequences, and we also identified some novel genes that are found only in the new genome.     

Human-mediated introgression of exotic chukar (Alectoris chukar, Galliformes) genes from East Asia into native Mediterranean partridges

Mediterranean red-legged (Alectoris rufa) and rock (Alectoris graeca) partridge populations are affected by genetic pollution. The chukar partridge (Alectoris chukar), a species only partly native to Europe, is the most frequently introgressive taxon detected in the genome of hybrid partridges. Both theoretical (evolutionary) and practical (resources management) matters spur to get insight into the geographic origin of the A. chukar hybridizing swarm. The phenotypic A. rufa populations colonizing the easternmost part of the distribution range of this species, the islands of Elba (Italy) and Corsica (France), were investigated. The analysis of both mitochondrial (mtDNA: Cytochrome-b gene plus Control Region: 2,250 characters) and nuclear (Short Tandem Repeats, STR; Random Amplified Polymorphic DNA, RAPD) genomes of 25 wild (Elba) and 20 captive (Corsica) partridges, disclosed spread introgression of chukar origin also in these populations. All mtDNA haplotypes of Elba and Corsica partridges along with those we obtained from other A. rufa (total, n = 111: Italy, Spain, France) and A. graeca (n = 6, Italy), were compared with the mtDNA haplotypes of chukars (n = 205) sampled in 20 countries. It was found that the A. chukar genes detected in red-legged (n = 43) and rock partridges (n = 4) of Spain, France and Italy as well as in either introduced (Italy) or native (Greece, Turkey) chukars (n = 35) were all from East Asia. Hence, a well-defined geographic origin of the exotic chukar genes polluting the genome of native Mediterranean A. rufa and A. graeca (inter-specific level) as well as A. chukar (intra-specific level), was demonstrated.     

Evaluation and Validation of Assembling Corrected PacBio Long Reads for Microbial Genome Completion via Hybrid Approaches

Despite the ever-increasing output of next-generation sequencing data along with developing assemblers, dozens to hundreds of gaps still exist in de novo microbial assemblies due to uneven coverage and large genomic repeats. Third-generation single-molecule, real-time (SMRT) sequencing technology avoids amplification artifacts and generates kilobase-long reads with the potential to complete microbial genome assembly. However, due to the low accuracy (~85%) of third-generation sequences, a considerable amount of long reads (>50X) are required for self-correction and for subsequent de novo assembly. Recently-developed hybrid approaches, using next-generation sequencing data and as few as 5X long reads, have been proposed to improve the completeness of microbial assembly. In this study we have evaluated the contemporary hybrid approaches and demonstrated that assembling corrected long reads (by runCA) produced the best assembly compared to long-read scaffolding (e.g., AHA, Cerulean and SSPACE-LongRead) and gap-filling (SPAdes). For generating corrected long reads, we further examined long-read correction tools, such as ECTools, LSC, LoRDEC, PBcR pipeline and proovread. We have demonstrated that three microbial genomes including Escherichia coli K12 MG1655, Meiothermus ruber DSM1279 and Pdeobacter heparinus DSM2366 were successfully hybrid assembled by runCA into near-perfect assemblies using ECTools-corrected long reads. In addition, we developed a tool, Patch, which implements corrected long reads and pre-assembled contigs as inputs, to enhance microbial genome assemblies. With the additional 20X long reads, short reads of S. cerevisiae W303 were hybrid assembled into 115 contigs using the verified strategy, ECTools + runCA. Patch was subsequently applied to upgrade the assembly to a 35-contig draft genome. Our evaluation of the hybrid approaches shows that assembling the ECTools-corrected long reads via runCA generates near complete microbial genomes, suggesting that genome assembly could benefit from re-analyzing the available hybrid datasets that were not assembled in an optimal fashion.     

QC-Chain: Fast and Holistic Quality Control Method for Next-Generation Sequencing Data

Next-generation sequencing (NGS) technologies have been widely used in life sciences. However, several kinds of sequencing artifacts, including low-quality reads and contaminating reads, were found to be quite common in raw sequencing data, which compromise downstream analysis. Therefore, quality control (QC) is essential for raw NGS data. However, although a few NGS data quality control tools are publicly available, there are two limitations: First, the processing speed could not cope with the rapid increase of large data volume. Second, with respect to removing the contaminating reads, none of them could identify contaminating sources de novo, and they rely heavily on prior information of the contaminating species, which is usually not available in advance. Here we report QC-Chain, a fast, accurate and holistic NGS data quality-control method. The tool synergeticly comprised of user-friendly tools for (1) quality assessment and trimming of raw reads using Parallel-QC, a fast read processing tool; (2) identification, quantification and filtration of unknown contamination to get high-quality clean reads. It was optimized based on parallel computation, so the processing speed is significantly higher than other QC methods. Experiments on simulated and real NGS data have shown that reads with low sequencing quality could be identified and filtered. Possible contaminating sources could be identified and quantified de novo, accurately and quickly. Comparison between raw reads and processed reads also showed that subsequent analyses (genome assembly, gene prediction, gene annotation, etc.) results based on processed reads improved significantly in completeness and accuracy. As regard to processing speed, QC-Chain achieves 7–8 time speed-up based on parallel computation as compared to traditional methods. Therefore, QC-Chain is a fast and useful quality control tool for read quality process and de novo contamination filtration of NGS reads, which could significantly facilitate downstream analysis. QC-Chain is publicly available at: http://www.computationalbioenergy.org/qc-chain.html.     

Identification of Optimum Sequencing Depth Especially for De Novo Genome Assembly of Small Genomes Using Next Generation Sequencing Data

Next Generation Sequencing (NGS) is a disruptive technology that has found widespread acceptance in the life sciences research community. The high throughput and low cost of sequencing has encouraged researchers to undertake ambitious genomic projects, especially in de novo genome sequencing. Currently, NGS systems generate sequence data as short reads and de novo genome assembly using these short reads is computationally very intensive. Due to lower cost of sequencing and higher throughput, NGS systems now provide the ability to sequence genomes at high depth. However, currently no report is available highlighting the impact of high sequence depth on genome assembly using real data sets and multiple assembly algorithms. Recently, some studies have evaluated the impact of sequence coverage, error rate and average read length on genome assembly using multiple assembly algorithms, however, these evaluations were performed using simulated datasets. One limitation of using simulated datasets is that variables such as error rates, read length and coverage which are known to impact genome assembly are carefully controlled. Hence, this study was undertaken to identify the minimum depth of sequencing required for de novo assembly for different sized genomes using graph based assembly algorithms and real datasets. Illumina reads for E.coli (4.6 MB) S.kudriavzevii (11.18 MB) and C.elegans (100 MB) were assembled using SOAPdenovo, Velvet, ABySS, Meraculous and IDBA-UD. Our analysis shows that 50X is the optimum read depth for assembling these genomes using all assemblers except Meraculous which requires 100X read depth. Moreover, our analysis shows that de novo assembly from 50X read data requires only 6–40 GB RAM depending on the genome size and assembly algorithm used. We believe that this information can be extremely valuable for researchers in designing experiments and multiplexing which will enable optimum utilization of sequencing as well as analysis resources.     

Rapid microsatellite marker development in the endangered neotropical freshwater turtle Podocnemis lewyana (Testudines: Podocnemididae) using 454 sequencing

Using high throughput sequencing we obtained a large number of microsatellites from Podocnemis lewyana, an endemic turtle from northwestern South America. We used 454 Genome Sequence FLX platform of sheared genomic DNA from randomly sampling approximately 17% of the haploid genome. We identified 86,501 reads (8.1% of all reads) that contained our definition of microsatellite loci. AC and TC were the most abundant motifs in the P. lewyana genome. TGC and AAAC were most abundant tri and tetra-nucleotide motifs respectively. 72.7% of microsatellite reads had flanking sequence regions suitable for primer design and PCR amplification. We validated the identified potentially amplifiable loci (PAL) and tested for polymorphism by selecting 15 loci corresponding to tetranucleotides. Twelve loci showed polymorphism in eight individuals. These findings demonstrates that microsatellite detection using next-generation sequencing is an efficient way of getting a lot of loci for listed taxa and in turn will have a large impact on future genetic studies aiming to understand and implement conservation plans for this highly threatened freshwater turtle.     

The De Novo Assembly of Mitochondrial Genomes of the Extinct Passenger Pigeon (Ectopistes migratorius) with Next Generation Sequencing

The information from ancient DNA (aDNA) provides an unparalleled opportunity to infer phylogenetic relationships and population history of extinct species and to investigate genetic evolution directly. However, the degraded and fragmented nature of aDNA has posed technical challenges for studies based on conventional PCR amplification. In this study, we present an approach based on next generation sequencing to efficiently sequence the complete mitochondrial genome (mitogenome) of two extinct passenger pigeons (Ectopistes migratorius) using de novo assembly of massive short (90 bp), paired-end or single-end reads. Although varying levels of human contamination and low levels of postmortem nucleotide lesion were observed, they did not impact sequencing accuracy. Our results demonstrated that the de novo assembly of shotgun sequence reads could be a potent approach to sequence mitogenomes, and offered an efficient way to infer evolutionary history of extinct species.     

The First Complete Chloroplast Genome Sequences in Actinidiaceae: Genome Structure and Comparative Analysis

Actinidia chinensis is an important economic plant belonging to the basal lineage of the asterids. Availability of a complete Actinidia chloroplast genome sequence is crucial to understanding phylogenetic relationships among major lineages of angiosperms and facilitates kiwifruit genetic improvement. We report here the complete nucleotide sequences of the chloroplast genomes for Actinidia chinensis and A. chinensis var deliciosa obtained through de novo assembly of Illumina paired-end reads produced by total DNA sequencing. The total genome size ranges from 155,446 to 157,557 bp, with an inverted repeat (IR) of 24,013 to 24,391 bp, a large single copy region (LSC) of 87,984 to 88,337 bp and a small single copy region (SSC) of 20,332 to 20,336 bp. The genome encodes 113 different genes, including 79 unique protein-coding genes, 30 tRNA genes and 4 ribosomal RNA genes, with 16 duplicated in the inverted repeats, and a tRNA gene (trnfM-CAU) duplicated once in the LSC region. Comparisons of IR boundaries among four asterid species showed that IR/LSC borders were extended into the 5’ portion of the psbA gene and IR contraction occurred in Actinidia. The clap gene has been lost from the chloroplast genome in Actinidia, and may have been transferred to the nucleus during chloroplast evolution. Twenty-seven polymorphic simple sequence repeat (SSR) loci were identified in the Actinidia chloroplast genome. Maximum parsimony analyses of a 72-gene, 16 taxa angiosperm dataset strongly support the placement of Actinidiaceae in Ericales within the basal asterids.