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

Dynamic analysis of a fuzzy Bobwhite quail population model under g-division law

This paper is concerned with a kind of Bobwhite quail population modelx n + 1 = A + B x n + x n x n - 1 x n - 2 , n = 0 , 1 , ⋯ , where the parameters and initial values are positive parabolic fuzzy numbers. According to g-division of fuzzy sets and based on the symmetrical parabolic fuzzy numbers, the conditional stability of this model is proved. Besides the existence, boundedness and persistence of its unique positive fuzzy solution. When some fuzzy stability conditions are satisfied, the model evolution exhibits oscillations with return to a fixed fuzzy equilibrium no matter what the initial value is. This phenomena provided a vivid counterexample to Allee effect in density-dependent populations of organisms. As a supplement, two numerical examples with data-table are interspersed to illustrate the effectiveness. Our findings have been verified precise with collected northern bobwhite data in Texas, and will help to form some efficient density estimates for wildlife populations of universal applications.

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.

A high-quality assembly reveals genomic characteristics, phylogenetic status, and causal genes for leucism plumage of Indian peafowl

Background

The dazzling phenotypic characteristics of male Indian peafowl (Pavo cristatus) are attractive both to the female of the species and to humans. However, little is known about the evolution of the phenotype and phylogeny of these birds at the whole-genome level. So far, there are no reports regarding the genetic mechanism of the formation of leucism plumage in this variant of Indian peafowl.

Results

A draft genome of Indian peafowl was assembled, with a genome size of 1.05 Gb (the sequencing depth is 362×), and contig and scaffold N50 were up to 6.2 and 11.4 Mb, respectively. Compared with other birds, Indian peafowl showed changes in terms of metabolism, immunity, and skeletal and feather development, which provided a novel insight into the phenotypic evolution of peafowl, such as the large body size and feather morphologies. Moreover, we determined that the phylogeny of Indian peafowl was more closely linked to turkey than chicken. Specifically, we first identified that PMEL was a potential causal gene leading to the formation of the leucism plumage variant in Indian peafowl.

Conclusions

This study provides an Indian peafowl genome of high quality, as well as a novel understanding of phenotypic evolution and phylogeny of Indian peafowl. These results provide a valuable reference for the study of avian genome evolution. Furthermore, the discovery of the genetic mechanism for the development of leucism plumage is both a breakthrough in the exploration of peafowl plumage and also offers clues and directions for further investigations of the avian plumage coloration and artificial breeding in peafowl.

Current and Forthcoming Approaches for Benchmarking Genetic and Genomic Diversity

The current attrition of biodiversity extends beyond loss of species and unique populations to steady loss of a vast genomic diversity that remains largely undescribed. Yet the accelerating development of new techniques allows us to survey entire genomes ever faster and cheaper, to obtain robust samples from a diversity of sources including degraded DNA and residual DNA in the environment, and to address conservation efforts in new and innovative ways. Here we review recent studies that highlight the importance of carefully considering where to prioritize collection of genetic samples (e.g., organisms in rapidly changing landscapes or along edges of geographic ranges) and what samples to collect and archive (e.g., from individuals of little-known subspecies or populations, even of species not currently considered endangered). Those decisions will provide the sample infrastructure to detect the disappearance of certain genotypes or gene complexes, increases in inbreeding levels, and loss of genomic diversity as environmental conditions change. Obtaining samples from currently endangered, protected, and rare species can be particularly difficult, thus we also focus on studies that use new, non-invasive ways of obtaining genomic samples and analyzing them in these cases where other sampling options are highly constrained. Finally, biological collections archiving such samples face an inherent contradiction: their main goal is to preserve biological material in good shape so it can be used for scientific research for centuries to come, yet the technologies that can make use of such materials are advancing faster than collections can change their standardized practices. Thus, we also discuss current and potential new practices in biological collections that might bolster their usefulness for future biodiversity conservation research.

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.

Dynamic shifts in social network structure and composition within a breeding hybrid population

Mating behaviour and the timing of reproduction can inhibit genetic exchange between closely related species; however, these reproductive barriers are challenging to measure within natural populations. Social network analysis provides promising tools for studying the social context of hybridization, and the exchange of genetic variation, more generally. We test how social networks within a hybrid population of California Callipepla californica and Gambel's quail Callipepla gambelii change over discrete periods of a breeding season. We assess patterns of phenotypic and genotypic assortment, and ask whether altered associations between individuals (association rewiring), or changes to the composition of the population (individual turnover) drive network dynamics. We use genetic data to test whether social associations and relatedness between individuals correlate with patterns of parentage within the hybrid population. To achieve these aims, we combine RFID association data, phenotypic data and genomic measures with social network analyses. We adopt methods from the ecological network literature to quantify shifts in network structure and to partition changes into those due to individual turnover and association rewiring. We integrate genomic data into networks as node-level attributes (ancestry) and edges (relatedness, parentage) to test links between social and parentage networks. We show that rewiring of associations between individuals that persist across network periods, rather than individual turnover, drives the majority of the changes in network structure throughout the breeding season, and that the traits involved in phenotypic/genotypic assortment were highly dynamic over time. Social networks were randomly assorted based on genetic ancestry, suggesting weak behavioural reproductive isolation within this hybrid population. Finally, we show that the strength of associations within the social network, but not levels of genetic relatedness, predicts patterns of parentage. Social networks play an important role in population processes such as the transmission of disease and information, yet there has been less focus on how networks influence the exchange of genetic variation. By integrating analyses of social structure, phenotypic assortment and reproductive outcomes within a hybrid zone, we demonstrate the utility of social networks for analysing links between social context and gene flow within wild populations.

Haplotype-resolved genomes provide insights into structural variation and gene content in Angus and Brahman cattle

Inbred animals were historically chosen for genome analysis to circumvent assembly issues caused by haplotype variation but this resulted in a composite of the two genomes. Here we report a haplotype-aware scaffolding and polishing pipeline which was used to create haplotype-resolved, chromosome-level genome assemblies of Angus (taurine) and Brahman (indicine) cattle subspecies from contigs generated by the trio binning method. These assemblies reveal structural and copy number variants that differentiate the subspecies and that variant detection is sensitive to the specific reference genome chosen. Six genes with immune related functions have additional copies in the indicine compared with taurine lineage and an indicus-specific extra copy of fatty acid desaturase is under positive selection. The haplotyped genomes also enable transcripts to be phased to detect allele-specific expression. This work exemplifies the value of haplotype-resolved genomes to better explore evolutionary and functional variations.

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.

BlobToolKit - Interactive Quality Assessment of Genome Assemblies

Reconstruction of target genomes from sequence data produced by instruments that are agnostic as to the species-of-origin may be confounded by contaminant DNA. Whether introduced during sample processing or through co-extraction alongside the target DNA, if insufficient care is taken during the assembly process, the final assembled genome may be a mixture of data from several species. Such assemblies can confound sequence-based biological inference and, when deposited in public databases, may be included in downstream analyses by users unaware of underlying problems. We present BlobToolKit, a software suite to aid researchers in identifying and isolating non-target data in draft and publicly available genome assemblies. BlobToolKit can be used to process assembly, read and analysis files for fully reproducible interactive exploration in the browser-based Viewer. BlobToolKit can be used during assembly to filter non-target DNA, helping researchers produce assemblies with high biological credibility. We have been running an automated BlobToolKit pipeline on eukaryotic assemblies publicly available in the International Nucleotide Sequence Data Collaboration and are making the results available through a public instance of the Viewer at https://blobtoolkit.genomehubs.org/view. We aim to complete analysis of all publicly available genomes and then maintain currency with the flow of new genomes. We have worked to embed these views into the presentation of genome assemblies at the European Nucleotide Archive, providing an indication of assembly quality alongside the public record with links out to allow full exploration in the Viewer.

The quail genome: insights into social behaviour, seasonal biology and infectious disease response

BACKGROUND

The Japanese quail (Coturnix japonica) is a popular domestic poultry species and an increasingly significant model species in avian developmental, behavioural and disease research.

RESULTS

We have produced a high-quality quail genome sequence, spanning 0.93 Gb assigned to 33 chromosomes. In terms of contiguity, assembly statistics, gene content and chromosomal organisation, the quail genome shows high similarity to the chicken genome. We demonstrate the utility of this genome through three diverse applications. First, we identify selection signatures and candidate genes associated with social behaviour in the quail genome, an important agricultural and domestication trait. Second, we investigate the effects and interaction of photoperiod and temperature on the transcriptome of the quail medial basal hypothalamus, revealing key mechanisms of photoperiodism. Finally, we investigate the response of quail to H5N1 influenza infection. In quail lung, many critical immune genes and pathways were downregulated after H5N1 infection, and this may be key to the susceptibility of quail to H5N1.

CONCLUSIONS

We have produced a high-quality genome of the quail which will facilitate further studies into diverse research questions using the quail as a model avian species.

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.

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.

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.

Computational detection and experimental validation of segmental duplications and associated copy number variations in water buffalo ( Bubalus bubalis )

Duplicated sequences are an important source of gene evolution and structural variation within mammalian genomes. Using a read depth approach based on next-generation sequencing, we performed a genome-wide analysis of segmental duplications (SDs) and associated copy number variations (CNVs) in the water buffalo (Bubalus bubalis). By aligning short reads of Olimpia (the reference water buffalo) to the UMD3.1 cattle genome, we identified 1,038 segmental duplications comprising 44.6 Mb (equivalent to ~1.73% of the cattle genome) of the autosomal and X chromosomal sequence in the buffalo genome. We experimentally validated 70.3% (71/101) of these duplications using fluorescent in situ hybridization. We also detected a total of 1,344 CNV regions across 14 additional water buffaloes, amounting to 59.8 Mb of variable sequence or the equivalent of 2.2% of the cattle genome. The CNV regions overlap 1,245 genes that are significantly enriched for specific biological functions including immune response, oxygen transport, sensory system and signal transduction. Additionally, we performed array Comparative Genomic Hybridization (aCGH) experiments using the 14 water buffaloes as test samples and Olimpia as the reference. Using a linear regression model, a high Pearson correlation (r = 0.781) was observed between the log₂ ratios between copy number estimates and the log₂ ratios of aCGH probes. We further designed Quantitative PCR assays to confirm CNV regions within or near annotated genes and found 74.2% agreement with our CNV predictions. These results confirm sub-chromosome-scale structural rearrangements present in the cattle and water buffalo. The information on genome variation that will be of value for evolutionary and phenotypic studies, and may be useful for selective breeding of both species.

Comparative genomics and transcriptomics of Chrysolophus provide insights into the evolution of complex plumage coloration

Background

As one of the most recognizable characteristics in birds, plumage color has a high impact on understanding the evolution and mechanisms of coloration. Feather and skin are ideal tissues to explore the genomics and complexity of color patterns in vertebrates. Two species of the genus Chrysolophus, golden pheasant (Chrysolophus pictus) and Lady Amherst's pheasant (Chrysolophus amherstiae), exhibit brilliant colors in their plumage, but with extreme phenotypic differences, making these two species great models to investigate plumage coloration mechanisms in birds.

Results

We sequenced and assembled a genome of golden pheasant with high coverage and annotated 15,552 protein-coding genes. The genome of Lady Amherst's pheasant is sequenced with low coverage. Based on the feather pigment identification, a series of genomic and transcriptomic comparisons were conducted to investigate the complex features of plumage coloration. By identifying the lineage-specific sequence variations in Chrysolophus and golden pheasant against different backgrounds, we found that four melanogenesis biosynthesis genes and some lipid-related genes might be candidate genomic factors for the evolution of melanin and carotenoid pigmentation, respectively. In addition, a study among 47 birds showed some candidate genes related to carotenoid coloration in a broad range of birds. The transcriptome data further reveal important regulators of the two colorations, particularly one splicing transcript of the microphthalmia-associated transcription factor gene for pheomelanin synthesis.

Conclusions

Analysis of the golden pheasant and its sister pheasant genomes, as well as comparison with other avian genomes, are helpful to reveal the underlying regulation of their plumage coloration. The present study provides important genomic information and insights for further studies of avian plumage evolution and diversity.

Genomic structural differences between cattle and River Buffalo identified through comparative genomic and transcriptomic analysis

Water buffalo (Bubalus bubalis L.) is an important livestock species worldwide. Like many other livestock species, water buffalo lacks high quality and continuous reference genome assembly, required for fine-scale comparative genomics studies. In this work, we present a dataset, which characterizes genomic differences between water buffalo genome and the extensively studied cattle (Bos taurus Taurus) reference genome. This data set is obtained after alignment of 14 river buffalo whole genome sequencing datasets to the cattle reference. This data set consisted of 13,444 deletion CNV regions, and 11,050 merged mobile element insertion (MEI) events within the upstream regions of annotated cattle genes. Gene expression data from cattle and buffalo were also presented for genes impacted by these regions. Public assessment of this dataset will allow for further analyses and functional annotation of genes that are potentially associated with phenotypic difference between cattle and water buffalo.

Comparative sequence alignment reveals River Buffalo genomic structural differences compared with cattle

This study sought to characterize differences in gene content, regulation and structure between taurine cattle and river buffalo (one subspecies of domestic water buffalo) using the extensively annotated UMD3.1 cattle reference genome as a basis for comparisons. We identified 127 deletion CNV regions in river buffalo representing 5 annotated cattle genes. We also characterized 583 merged mobile element insertion (MEI) events within the upstream regions of annotated cattle genes. Transcriptome analysis in various tissue types on river buffalo confirmed the absence of four cattle genes. Four genes which may be related to phenotypic differences in meat quality and color, had upstream MEI predictions and were found to have significantly elevated expression in river buffalo compared with cattle. Our comparative alignment approach and gene expression analyses suggested a functional role for many genomic structural variations, which may contribute to the unique phenotypes of river buffalo.