A draft de novo genome assembly for the northern bobwhite (Colinus virginianus) reveals evidence for a rapid decline in effective population size beginning in the Late Pleistocene

Selection of suitable reference genes for normalization of RT-qPCR in three tissues of Northern bobwhite (Colinus virginianus) infected with eyeworm (Oxyspirura petrowi)

BACKGROUND

The Northern bobwhite (Colinus virginianus) is an economically important, and popular game bird in North America. Northern bobwhites have experiencing declines of > 3.5% annually in recent decades due to several factors. The eyeworm Oxyspirura petrowi is a nematode parasite frequently found in the eyes of bobwhites. Although reported frequently in wild bobwhites, there is no research to understand the host-parasite mechanism. Hence, it is important to investigate mechanisms of eyeworm invasion and immune modulation in bobwhite. Cytokine gene expression using RT-PCR is widely used to identify the innate immune response of a host to an infection.

METHODOLOGY

In this study, we evaluated ten reference genes (HMBS, RPL19, RPL32, RPS7, RPS8, TATA, SDHA, YWHAZ, GAPDH, and ACTB) for their stability across three tissues (liver, spleen, and caecal tonsils) of control and O. petrowi infected Northern bobwhites. Primer efficiency and reference genes stability were assessed using GeNorm, NormFinder, and BestKeeper.

RESULTS

Expression of these reference genes with respect to O. petrowi infection in bobwhites showed RPL32 and HMBS were the most stable genes in the liver, HMBS and SDHA were the most stable genes in the spleen, and HMBS and YWHAZ were equally stable reference genes in the caecal tonsils.

CONCLUSION

Based on the geometric mean of all three analyses, our results indicate that the combination of RPL32 and HMBS for the liver, HMBS and SDHA for the spleen, and YWHAZ and HMBS for caecal tonsils might be used as reference genes for normalization in gene expression investigations on Northern bobwhites.

An island 'endemic' born out of hybridization between introduced lineages

Humans have profoundly impacted the distribution of plant and animal species over thousands of years. The most direct example of these effects is human-mediated movement of individuals, either through translocation of individuals within their range or through the introduction of species to new habitats. While human involvement may be suspected in species with obvious range disjunctions, it can be difficult to detect natural versus human-mediated dispersal events for populations at the edge of a species' range, and this uncertainty muddles how we understand the evolutionary history of populations and broad biogeographical patterns. Studies combining genetic data with archaeological, linguistic and historical evidence have confirmed prehistoric examples of human-mediated dispersal; however, it is unclear whether these methods can disentangle recent dispersal events, such as species translocated by European colonizers during the past 500 years. We use genomic DNA from historical museum specimens and historical records to evaluate three hypotheses regarding the timing and origin of Northern Bobwhites (Colinus virginianus) in Cuba, whose status as an endemic or introduced population has long been debated. We discovered that bobwhites from southern Mexico arrived in Cuba between the 12th and 16th centuries, followed by the subsequent introduction of bobwhites from the southeastern USA to Cuba between the 18th and 20th centuries. These dates suggest the introduction of bobwhites to Cuba was human-mediated and concomitant with Spanish colonial shipping routes between Veracruz, Mexico and Havana, Cuba during this period. Our results identify endemic Cuban bobwhites as a genetically distinct population born of hybridization between divergent, introduced lineages.

Molecular phylogenetic and in silico analysis of glyceraldeyde-3-phosphate dehydrogenase (GAPDH) gene from northern bobwhite quail (Colinus virginianus)

Many recent studies have been focused on prevalence and impact of two helminth parasites, eyeworm Oxyspirura petrowi and caecal worm Aulonocephalus pennula, in the northern bobwhite quail (Colinus virginianus). However, few studies have attempted to examine the effect of these parasites on the bobwhite immune system. This is likely due to the lack of proper reference genes for relative gene expression studies. Glyceraldehyde 3-phosphate dehydrogenase (GAPDH) is a glycolytic enzyme that is often utilized as a reference gene, and in this preliminary study, we evaluated the similarity of bobwhite GAPDH to GAPDH in other avian species to evaluate its potential as a reference gene in bobwhite. GAPDH was identified in the bobwhite full genome sequence and multiple sets of PCR primers were designed to generate overlapping PCR products. These products were then sequenced and then aligned to generate the sequence for the full-length open reading frame (ORF) of bobwhite GAPDH. Utilizing this sequence, phylogenetic analyses and comparative analysis of the exon-intron pattern were conducted that revealed high similarity of GAPDH encoding sequences among bobwhite and other Galliformes. Additionally, This ORF sequence was also used to predict the encoded protein and its three-dimensional structure which like the phylogenetic analyses reveal that bobwhite GAPDH is similar to GAPDH in other Galliformes. Finally, GAPDH qPCR primers were designed, standardized, and tested with bobwhite both uninfected and infected with O. petrowi, and this preliminary test showed no statistical difference in expression of GAPDH between the two groups. These analyses are the first to investigate GAPDH in bobwhite. These efforts in phylogeny, sequence analysis, and protein structure suggest that there is > 97% conservation of GADPH among Galliformes. Furthermore, the results of these in silico tests and the preliminary qPCR indicate that GAPDH is a prospective candidate for use in gene expression analyses in bobwhite.

Accurate Genomic Predictions for Chronic Wasting Disease in U.S. White-Tailed Deer

The geographic expansion of chronic wasting disease (CWD) in U.S. white-tailed deer (Odocoileus virginianus) has been largely unabated by best management practices, diagnostic surveillance, and depopulation of positive herds. Using a custom Affymetrix Axiom single nucleotide polymorphism (SNP) array, we demonstrate that both differential susceptibility to CWD, and natural variation in disease progression, are moderately to highly heritable (h2=0.337±0.079─0.637±0.070) among farmed U.S. white-tailed deer, and that loci other than PRNP are involved. Genome-wide association analyses using 123,987 quality filtered SNPs for a geographically diverse cohort of 807 farmed U.S. white-tailed deer (n = 284 CWD positive; n = 523 CWD non-detect) confirmed the prion gene (PRNP; G96S) as a large-effect risk locus (P-value < 6.3E-11), as evidenced by the estimated proportion of phenotypic variance explained (PVE ≥ 0.05), but also demonstrated that more phenotypic variance was collectively explained by loci other than PRNP. Genomic best linear unbiased prediction (GBLUP; n = 123,987 SNPs) with k-fold cross validation (k = 3; k = 5) and random sampling (n = 50 iterations) for the same cohort of 807 farmed U.S. white-tailed deer produced mean genomic prediction accuracies ≥ 0.81; thereby providing the necessary foundation for exploring a genomically-estimated CWD eradication program.

Conservation Genomics in the Sagebrush Sea: Population Divergence, Demographic History, and Local Adaptation in Sage-Grouse (Centrocercus spp.)

Sage-grouse are two closely related iconic species of the North American West, with historically broad distributions across sagebrush-steppe habitat. Both species are dietary specialists on sagebrush during winter, with presumed adaptations to tolerate the high concentrations of toxic secondary metabolites that function as plant chemical defenses. Marked range contraction and declining population sizes since European settlement have motivated efforts to identify distinct population genetic variation, particularly that which might be associated with local genetic adaptation and dietary specialization of sage-grouse. We assembled a reference genome and performed whole-genome sequencing across sage-grouse from six populations, encompassing both species and including several populations on the periphery of the species ranges. Population genomic analyses reaffirmed genome-wide differentiation between greater and Gunnison sage-grouse, revealed pronounced intraspecific population structure, and highlighted important differentiation of a small isolated population of greater sage-grouse in the northwest of the range. Patterns of genome-wide differentiation were largely consistent with a hypothesized role of genetic drift due to limited gene flow among populations. Inferred ancient population demography suggested persistent declines in effective population sizes that have likely contributed to differentiation within and among species. Several genomic regions with single-nucleotide polymorphisms exhibiting extreme population differentiation were associated with candidate genes linked to metabolism of xenobiotic compounds. In vitro activity of enzymes isolated from sage-grouse livers supported a role for these genes in detoxification of sagebrush, suggesting that the observed interpopulation variation may underlie important local dietary adaptations, warranting close consideration for conservation strategies that link sage-grouse to the chemistry of local sagebrush.

A Highly Contiguous Reference Genome for Northern Bobwhite (Colinus virginianus)

Northern bobwhites (Colinus virginianus) are small quails in the New World Quail family (Odontophoridae) and are one of the most phenotypically diverse avian species. Despite extensive research on bobwhite ecology, genomic studies investigating the evolution of phenotypic diversity in this species are lacking. Here, we present a new, highly contiguous assembly for bobwhites using tissue samples from a vouchered, wild, female bird collected in Louisiana. By performing a de novo assembly and scaffolding the assembly with Dovetail Chicago and HiC libraries and the HiRise pipeline, we produced an 866.8 Mb assembly including 1,512 scaffolds with a scaffold N50 of 66.8 Mb, a scaffold L90 of 17, and a BUSCO completeness score of 90.8%. This new assembly represents approximately 96% of the non-repetitive and 84% of the entire bobwhite genome size, greatly improves scaffold lengths and contiguity compared to an existing draft bobwhite genome, and provides an important tool for future studies of evolutionary and functional genomics in bobwhites.

A High-Quality Draft Genome Assembly of the Black-Necked Crane (Grus nigricollis) Based on Nanopore Sequencing

The black-necked crane (Grus nigricollis) which inhabits high-altitude areas has the largest body size of the world's 15 crane species, and is classified as threatened by the IUCN. To support further studies on population genetics and genomics, we present a high-quality genome assembly based on both Illumina and nanopore sequencing. In total, 54.59 Gb Illumina short reads and 116.5 Gb nanopore long reads were generated. The 1.23 Gb assembled genome has a high contig N50 of 17.89 Mb, and has a longest contig of 87.83 Mb. The completeness (97.7%) of the draft genome was evaluated with single-copy orthologous genes using BUSCO. We identified 17,789 genes and found that 8.11% of the genome is composed of repetitive elements. In total, 84 of the 2,272 one-to-one orthologous genes were under positive selection in the black-necked crane lineage. SNP-based inference indicated two bottlenecks in the recent demographic trajectories of the black-necked crane. The genome information will contribute to future study of crane evolutionary history and provide new insights into the potential adaptation mechanisms of the black-necked crane to its high-altitude habitat.

The Draft Genome of the Endangered Sichuan Partridge (Arborophila rufipectus) with Evolutionary Implications

The Sichuan partridge (Arborophila rufipectus, Phasianidae, Galliformes) is distributed in south-west China, and classified as endangered grade. To examine the evolution and genomic features of Sichuan partridge, we de novo assembled the Sichuan partridge reference genome. The final draft assembly consisted of approximately 1.09 Gb, and had a scaffold N50 of 4.57 Mb. About 1.94 million heterozygous single-nucleotide polymorphisms (SNPs) were detected, 17,519 protein-coding genes were predicted, and 9.29% of the genome was identified as repetitive elements. A total of 56 olfactory receptor (OR) genes were found in Sichuan partridge, and conserved motifs were detected. Comparisons between the Sichuan partridge genome and chicken genome revealed a conserved genome structure, and phylogenetic analysis demonstrated that Arborophila possessed a basal phylogenetic position within Phasianidae. Gene Ontology (GO) enrichment analysis of positively selected genes (PSGs) in Sichuan partridge showed over-represented GO functions related to environmental adaptation, such as energy metabolism and behavior. Pairwise sequentially Markovian coalescent analysis revealed the recent demographic trajectory for the Sichuan partridge. Our data and findings provide valuable genomic resources not only for studying the evolutionary adaptation, but also for facilitating the long-term conservation and genetic diversity for this endangered species.

Quantitative analysis of Northern bobwhite (Colinus virginianus) cytokines and TLR expression to eyeworm (Oxyspirura petrowi) and caecal worm (Aulonocephalus pennula) glycoproteins

Helminth parasites have been a popular research topic due to their global prevalence and adverse effects on livestock and game species. The Northern bobwhite (Colinus virginianus), a popular game bird in the USA, is one species subject to helminth infection and has been experiencing a decline of > 4% annually over recent decades. In the Rolling Plains Ecoregion of Texas, the eyeworm (Oxyspirura petrowi) and caecal worm (Aulonocephalus pennula) helminths are found to be highly prevalent in bobwhite. While there have been increasing studies on the prevalence, pathology, and phylogeny of the eyeworm and caecal worm, there is still a need to investigate the bobwhite immune response to infection. This study utilizes previously sequenced bobwhite cytokines and toll-like receptors to develop and optimize qPCR primers and measure gene expression in bobwhite intramuscularly challenged with eyeworm and caecal worm glycoproteins. For the challenge experiments, separate treatments of eyeworm and caecal worm glycoproteins were administered to bobwhite on day 1 and day 21. Measurements of primary and secondary immune responses were taken at day 7 and day 28, respectively. Using the successfully optimized qPCR primers for TLR7, IL1β, IL6, IFNα, IFNγ, IL10, and β-actin, the gene expression analysis from the challenge experiments revealed that there was a measurable immune reaction in bobwhite in response to the intramuscular challenge of eyeworm and caecal worm glycoproteins.

Genome-wide analysis sheds light on the high-altitude adaptation of the buff-throated partridge (Tetraophasis szechenyii)

The buff-throated partridge (Tetraophasis szechenyii) is a hypoxia-tolerant bird living in an extremely inhospitable high-altitude environment, which has high ultraviolet (UV) radiation as well as a low oxygen supply when compared with low-altitude areas. To further understand the molecular genetic mechanisms of the high-altitude adaptation of the buff-throated partridges, we de novo assembled the complete genome of the buff-throated partridge. Comparative genomics revealed that positively selected hypoxia-related genes in the buff-throated partridge were distributed in the HIF-1 signaling pathway (map04066), response to hypoxia (GO:0001666), response to oxygen-containing compound (GO:1901700), ATP binding (GO:0005524), and angiogenesis (GO:0001525). Of these positively selected hypoxia-related genes, one positively selected gene (LONP1) had one buff-throated partridge-specific missense mutation which was classified as deleterious by PolyPhen-2. Moreover, positively selected genes in the buff-throated partridge were enriched in cellular response to DNA damage stimulus (corrected P value: 0.028006) and DNA repair (corrected P value: 0.044549), which was related to the increased exposure of the buff-throated partridge to UV radiation. Compared with other avian genomes, the buff-throated partridge showed expansion in genes associated with steroid hormone receptor activity and contractions in genes related to immune and olfactory perception. Furthermore, comparisons between the buff-throated partridge genome and red junglefowl genome revealed a conserved genome structure and provided strong evidence of the sibling relationship between Tetraophasis and Lophophorus. Our data and analysis contributed to the study of Phasianidae evolutionary history and provided new insights into the potential adaptation mechanisms to the high altitude employed by the buff-throated partridge.

Systematic Identification and Evolution Analysis of Sox Genes in Coturnix japonica Based on Comparative Genomics

Coturnix japonica (Japanese quail) has been extensively used as a model animal for biological studies. The Sox gene family, which was systematically characterized by a high-mobility group (HMG-box) in many animal species, encodes transcription factors that play central roles during multiple developmental processes. However, genome-wide investigations on the Sox gene family in birds are scarce. In the current study, we first performed a genome-wide study to explore the Sox gene family in galliform birds. Based on available genomic sequences retrieved from the NCBI database, we focused on the global identification of the Sox gene family in C. japonica and other species in Galliformes, and the evolutionary relationships of Sox genes. In our result, a total of 35 Sox genes in seven groups were identified in the C. japonica genome. Our results also revealed that dispersed gene duplications contributed the most to the expansion of the Sox gene family in Galliform birds. Evolutionary analyses indicated that Sox genes are an ancient gene family, and strong purifying selections played key roles in the evolution of CjSox genes of C. japonica. More interestingly, we observed that most Sox genes exhibited highly embryo-specific expression in both gonads. Our findings provided new insights into the molecular function and phylogeny of Sox gene family in birds.

Comparison of the Chinese bamboo partridge and red Junglefowl genome sequences highlights the importance of demography in genome evolution

BACKGROUND

Recent large-scale whole genome sequencing efforts in birds have elucidated broad patterns of avian phylogeny and genome evolution. However, despite the great interest in economically important phasianids like Gallus gallus (Red Junglefowl, the progenitor of the chicken), we know little about the genomes of closely related species. Gallus gallus is highly sexually dichromatic and polygynous, but its sister genus, Bambusicola, is smaller, sexually monomorphic, and monogamous with biparental care. We sequenced the genome of Bambusicola thoracicus (Chinese Bamboo Partridge) using a single insert library to test hypotheses about genome evolution in galliforms. Selection acting at the phenotypic level could result in more evidence of positive selection in the Gallus genome than in Bambusicola. However, the historical range size of Bambusicola was likely smaller than Gallus, and demographic effects could lead to higher rates of nonsynonymous substitution in Bambusicola than in Gallus.

RESULTS

We generated a genome assembly suitable for evolutionary analyses. We examined the impact of selection on coding regions by examining shifts in the average nonsynonymous to synonymous rate ratio (dN/dS) and the proportion of sites subject to episodic positive selection. We observed elevated dN/dS in Bambusicola relative to Gallus, which is consistent with our hypothesis that demographic effects may be important drivers of genome evolution in Bambusicola. We also demonstrated that alignment error can greatly inflate estimates of the number of genes that experienced episodic positive selection and heterogeneity in dN/dS. However, overall patterns of molecular evolution were robust to alignment uncertainty. Bambusicola thoracicus has higher estimates of heterozygosity than Gallus gallus, possibly due to migration events over the past 100,000 years.

CONCLUSIONS

Our results emphasized the importance of demographic processes in generating the patterns of variation between Bambusicola and Gallus. We also demonstrated that genome assemblies generated using a single library can provide valuable insights into avian evolutionary history and found that it is important to account for alignment uncertainty in evolutionary inferences from draft genomes.

Annotated Draft Genome Assemblies for the Northern Bobwhite (Colinus virginianus) and the Scaled Quail (Callipepla squamata) Reveal Disparate Estimates of Modern Genome Diversity and Historic Effective Population Size

Northern bobwhite (Colinus virginianus; hereafter bobwhite) and scaled quail (Callipepla squamata) populations have suffered precipitous declines across most of their US ranges. Illumina-based first- (v1.0) and second- (v2.0) generation draft genome assemblies for the scaled quail and the bobwhite produced N50 scaffold sizes of 1.035 and 2.042 Mb, thereby producing a 45-fold improvement in contiguity over the existing bobwhite assembly, and ≥90% of the assembled genomes were captured within 1313 and 8990 scaffolds, respectively. The scaled quail assembly (v1.0 = 1.045 Gb) was ∼20% smaller than the bobwhite (v2.0 = 1.254 Gb), which was supported by kmer-based estimates of genome size. Nevertheless, estimates of GC content (41.72%; 42.66%), genome-wide repetitive content (10.40%; 10.43%), and MAKER-predicted protein coding genes (17,131; 17,165) were similar for the scaled quail (v1.0) and bobwhite (v2.0) assemblies, respectively. BUSCO analyses utilizing 3023 single-copy orthologs revealed a high level of assembly completeness for the scaled quail (v1.0; 84.8%) and the bobwhite (v2.0; 82.5%), as verified by comparison with well-established avian genomes. We also detected 273 putative segmental duplications in the scaled quail genome (v1.0), and 711 in the bobwhite genome (v2.0), including some that were shared among both species. Autosomal variant prediction revealed ∼2.48 and 4.17 heterozygous variants per kilobase within the scaled quail (v1.0) and bobwhite (v2.0) genomes, respectively, and estimates of historic effective population size were uniformly higher for the bobwhite across all time points in a coalescent model. However, large-scale declines were predicted for both species beginning ∼15-20 KYA.

Identification of New World Quails Susceptible to Infection with Avian Leukosis Virus Subgroup J

The J subgroup of avian leukosis virus (ALV-J) infects domestic chickens, jungle fowl, and turkeys. This virus enters the host cell through a receptor encoded by the tvj locus and identified as Na+/H+ exchanger 1. The resistance to avian leukosis virus subgroup J in a great majority of galliform species has been explained by deletions or substitutions of the critical tryptophan 38 in the first extracellular loop of Na+/H+ exchanger 1. Because there are concerns of transspecies virus transmission, we studied natural polymorphisms and susceptibility/resistance in wild galliforms and found the presence of tryptophan 38 in four species of New World quails. The embryo fibroblasts of New World quails are susceptible to infection with avian leukosis virus subgroup J, and the cloned Na+/H+ exchanger 1 confers susceptibility on the otherwise resistant host. New World quails are also susceptible to new avian leukosis virus subgroup J variants but resistant to subgroups A and B and weakly susceptible to subgroups C and D of avian sarcoma/leukosis virus due to obvious defects of the respective receptors. Our results suggest that the avian leukosis virus subgroup J could be transmitted to New World quails and establish a natural reservoir of circulating virus with a potential for further evolution.

IMPORTANCE

Since its spread in broiler chickens in China and Southeast Asia in 2000, ALV-J remains a major enzootic challenge for the poultry industry. Although the virus diversifies rapidly in the poultry, its spillover and circulation in wild bird species has been prevented by the resistance of most species to ALV-J. It is, nevertheless, important to understand the evolution of the virus and its potential host range in wild birds. Because resistance to avian retroviruses is due particularly to receptor incompatibility, we studied Na+/H+ exchanger 1, the receptor for ALV-J. In New World quails, we found a receptor compatible with virus entry, and we confirmed the susceptibilities of four New World quail species in vitro We propose that a prospective molecular epidemiology study be conducted to identify species with the potential to become reservoirs for ALV-J.

Genome sequencing of Metrosideros polymorpha (Myrtaceae), a dominant species in various habitats in the Hawaiian Islands with remarkable phenotypic variations

Whole genome sequences, which can be provided even for non-model organisms owing to high-throughput sequencers, are valuable in enhancing the understanding of adaptive evolution. Metrosideros polymorpha, a tree species endemic to the Hawaiian Islands, occupies a wide range of ecological habitats and shows remarkable polymorphism in phenotypes among/within populations. The biological functions of genetic variations observed within this species could provide significant insights into the adaptive radiation found in a single species. Here de novo assembled genome sequences of M. polymorpha are presented to reveal basic genomic parameters about this species and to develop our knowledge of ecological divergences. The assembly yielded 304-Mbp genome sequences, half of which were covered by 19 scaffolds with >5 Mbp, and contained 30 K protein-coding genes. Demographic history inferred from the genome-wide heterozygosity indicated that this species experienced a dramatic rise and fall in the effective population size, possibly owing to past geographic or climatic changes in the Hawaiian Islands. This M. polymorpha genome assembly represents a high-quality genome resource useful for future functional analyses of both intra- and interspecies genetic variations or comparative genomics.

Northern Bobwhite (Colinus virginianus) Mitochondrial Population Genomics Reveals Structure, Divergence, and Evidence for Heteroplasmy

Herein, we evaluated the concordance of population inferences and conclusions resulting from the analysis of short mitochondrial fragments (i.e., partial or complete D-Loop nucleotide sequences) versus complete mitogenome sequences for 53 bobwhites representing six ecoregions across TX and OK (USA). Median joining (MJ) haplotype networks demonstrated that analyses performed using small mitochondrial fragments were insufficient for estimating the true (i.e., complete) mitogenome haplotype structure, corresponding levels of divergence, and maternal population history of our samples. Notably, discordant demographic inferences were observed when mismatch distributions of partial (i.e., partial D-Loop) versus complete mitogenome sequences were compared, with the reduction in mitochondrial genomic information content observed to encourage spurious inferences in our samples. A probabilistic approach to variant prediction for the complete bobwhite mitogenomes revealed 344 segregating sites corresponding to 347 total mutations, including 49 putative nonsynonymous single nucleotide variants (SNVs) distributed across 12 protein coding genes. Evidence of gross heteroplasmy was observed for 13 bobwhites, with 10 of the 13 heteroplasmies involving one moderate to high frequency SNV. Haplotype network and phylogenetic analyses for the complete bobwhite mitogenome sequences revealed two divergent maternal lineages (d_XY = 0.00731; F_ST = 0.849; P < 0.05), thereby supporting the potential for two putative subspecies. However, the diverged lineage (n = 103 variants) almost exclusively involved bobwhites geographically classified as Colinus virginianus texanus, which is discordant with the expectations of previous geographic subspecies designations. Tests of adaptive evolution for functional divergence (MKT), frequency distribution tests (D, F_S) and phylogenetic analyses (RAxML) provide no evidence for positive selection or hybridization with the sympatric scaled quail (Callipepla squamata) as being explanatory factors for the two bobwhite maternal lineages observed. Instead, our analyses support the supposition that two diverged maternal lineages have survived from pre-expansion to post-expansion population(s), with the segregation of some slightly deleterious nonsynonymous mutations.

Avoiding Missing Data Biases in Phylogenomic Inference: An Empirical Study in the Landfowl (Aves: Galliformes)

Production of massive DNA sequence data sets is transforming phylogenetic inference, but best practices for analyzing such data sets are not well established. One uncertainty is robustness to missing data, particularly in coalescent frameworks. To understand the effects of increasing matrix size and loci at the cost of increasing missing data, we produced a 90 taxon, 2.2 megabase, 4,800 locus sequence matrix of landfowl using target capture of ultraconserved elements. We then compared phylogenies estimated with concatenated maximum likelihood, quartet-based methods executed on concatenated matrices and gene tree reconciliation methods, across five thresholds of missing data. Results of maximum likelihood and quartet analyses were similar, well resolved, and demonstrated increasing support with increasing matrix size and sparseness. Conversely, gene tree reconciliation produced unexpected relationships when we included all informative loci, with certain taxa placed toward the root compared with other approaches. Inspection of these taxa identified a prevalence of short average contigs, which potentially biased gene tree inference and caused erroneous results in gene tree reconciliation. This suggests that the more problematic missing data in gene tree-based analyses are partial sequences rather than entire missing sequences from locus alignments. Limiting gene tree reconciliation to the most informative loci solved this problem, producing well-supported topologies congruent with concatenation and quartet methods. Collectively, our analyses provide a well-resolved phylogeny of landfowl, including strong support for previously problematic relationships such as those among junglefowl (Gallus), and clarify the position of two enigmatic galliform genera (Lerwa, Melanoperdix) not sampled in previous molecular phylogenetic studies.

Lineage-specific loss of FGF17 within the avian orders Galliformes and Passeriformes

The genomic and developmental complexity of vertebrates is commonly attributed to two rounds of whole genome duplications which occurred at the base of the vertebrate radiation. These duplications led to the rise of several, multi-gene families of developmental proteins like the fibroblast growth factors (FGFs); a signaling protein family which functions at various stages of embryonic development. One of the major FGF assemblages arising from these duplications is the FGF8 subfamily, which includes FGF8, FGF17, and FGF18 in tetrapods. While FGF8 and FGF18 are found in all tetrapods and are critical for embryonic survival, genomic analyses suggest putative loss of FGF17 in various lineages ranging from frogs and fish, to the chicken. This study utilizes 27 avian genomes in conjunction with molecular analyses of chicken embryos to confirm the loss of FGF17 in chicken as a true, biological occurrence. FGF17 is also missing in the turkey, black grouse, Japanese quail and northern bobwhite genomes. These species, along with chicken, form a monophyletic clade in the order Galliformes. Four additional species, members of the clade Passeroidea, within the order Passeriformes, are also missing FGF17. Additionally, analysis of intact FGF17 in other avian lineages reveals that it is still under strong purifying selection, despite being seemingly dispensable. Thus, FGF17 likely represents a molecular spandrel arising from a genome duplication event and due to its high connectivity with FGF8/FGF18, and potential for interference with their function, is retained under strong purifying selection, despite itself not having a strong selective advantage.

Genomic sequencing and microsatellite marker development for Boswellia papyrifera, an economically important but threatened tree native to dry tropical forests

Microsatellite (or simple sequence repeat, SSR) markers are highly informative DNA markers often used in conservation genetic research. Next-generation sequencing enables efficient development of large numbers of SSR markers at lower costs. Boswellia papyrifera is an economically important tree species used for frankincense production, an aromatic resinous gum exudate from bark. It grows in dry tropical forests in Africa and is threatened by a lack of rejuvenation. To help guide conservation efforts for this endangered species, we conducted an analysis of its genomic DNA sequences using Illumina paired-end sequencing. The genome size was estimated at 705 Mb per haploid genome. The reads contained one microsatellite repeat per 5.7 kb. Based on a subset of these repeats, we developed 46 polymorphic SSR markers that amplified 2-12 alleles in 10 genotypes. This set included 30 trinucleotide repeat markers, four tetranucleotide repeat markers, six pentanucleotide markers and six hexanucleotide repeat markers. Several markers were cross-transferable to Boswellia pirrotae and B. popoviana. In addition, retrotransposons were identified, the reads were assembled and several contigs were identified with similarity to genes of the terpene and terpenoid backbone synthesis pathways, which form the major constituents of the bark resin.

The Genome 10K Project: a way forward

The Genome 10K Project was established in 2009 by a consortium of biologists and genome scientists determined to facilitate the sequencing and analysis of the complete genomes of 10,000 vertebrate species. Since then the number of selected and initiated species has risen from ∼26 to 277 sequenced or ongoing with funding, an approximately tenfold increase in five years. Here we summarize the advances and commitments that have occurred by mid-2014 and outline the achievements and present challenges of reaching the 10,000-species goal. We summarize the status of known vertebrate genome projects, recommend standards for pronouncing a genome as sequenced or completed, and provide our present and future vision of the landscape of Genome 10K. The endeavor is ambitious, bold, expensive, and uncertain, but together the Genome 10K Consortium of Scientists and the worldwide genomics community are moving toward their goal of delivering to the coming generation the gift of genome empowerment for many vertebrate species.

A field guide to whole-genome sequencing, assembly and annotation

Genome sequencing projects were long confined to biomedical model organisms and required the concerted effort of large consortia. Rapid progress in high-throughput sequencing technology and the simultaneous development of bioinformatic tools have democratized the field. It is now within reach for individual research groups in the eco-evolutionary and conservation community to generate de novo draft genome sequences for any organism of choice. Because of the cost and considerable effort involved in such an endeavour, the important first step is to thoroughly consider whether a genome sequence is necessary for addressing the biological question at hand. Once this decision is taken, a genome project requires careful planning with respect to the organism involved and the intended quality of the genome draft. Here, we briefly review the state of the art within this field and provide a step-by-step introduction to the workflow involved in genome sequencing, assembly and annotation with particular reference to large and complex genomes. This tutorial is targeted at scientists with a background in conservation genetics, but more generally, provides useful practical guidance for researchers engaging in whole-genome sequencing projects.