Single nucleotide polymorphisms (SNPs) in the open reading frame (ORF) of prion protein gene (PRNP) in Nigerian livestock species

BACKGROUND

Prion diseases, also known as transmissible spongiform encephalopathies (TSEs) remain one of the deleterious disorders, which have affected several animal species. Polymorphism of the prion protein (PRNP) gene majorly determines the susceptibility of animals to TSEs. However, only limited studies have examined the variation in PRNP gene in different Nigerian livestock species. Thus, this study aimed to identify the polymorphism of PRNP gene in Nigerian livestock species (including camel, dog, horse, goat, and sheep). We sequenced the open reading frame (ORF) of 65 camels, 31 village dogs and 12 horses from Nigeria and compared with PRNP sequences of 886 individuals retrieved from public databases.

RESULTS

All the 994 individuals were assigned into 162 haplotypes. The sheep had the highest number of haplotypes (n = 54), and the camel had the lowest (n = 7). Phylogenetic tree further confirmed clustering of Nigerian individuals into their various species. We detected five non-synonymous SNPs of PRNP comprising of G9A, G10A, C11G, G12C, and T669C shared by all Nigerian livestock species and were in Hardy-Weinberg Equilibrium (HWE). The amino acid changes in these five non-synonymous SNP were all "benign" via Polyphen-2 program. Three SNPs G34C, T699C, and C738G occurred only in Nigerian dogs while C16G, G502A, G503A, and C681A in Nigerian horse. In addition, C50T was detected only in goats and sheep.

CONCLUSION

Our study serves as the first to simultaneously investigate the polymorphism of PRNP gene in Nigerian livestock species and provides relevant information that could be adopted in programs targeted at breeding for prion diseases resistance.

Polymorphism of prion protein gene (PRNP) in Nigerian sheep

Polymorphism of the prion protein gene (PRNP) gene determines an animal's susceptibility to scrapie. Three polymorphisms at codons 136, 154, and 171 have been linked to classical scrapie susceptibility, although many variants of PRNP have been reported. However, no study has investigated scrapie susceptibility in Nigerian sheep from the drier agro-climate zones. In this study, we aimed to identify PRNP polymorphism in nucleotide sequences of 126 Nigerian sheep by comparing them with public available studies on scrapie-affected sheep. Further, we deployed Polyphen-2, PROVEAN, and AMYCO analyses to determine the structure changes produced by the non-synonymous SNPs. Nineteen (19) SNPs were found in Nigerian sheep with 14 being non-synonymous. Interestingly, one novel SNP (T718C) was identified. There was a significant difference (P < 0.05) in the allele frequencies of PRNP codon 154 between sheep in Italy and Nigeria. Based on the prediction by Polyphen-2, R154H was probably damaging while H171Q was benign. Contrarily, all SNPs were neutral via PROVEAN analysis while two haplotypes (HYKK and HDKK) had similar amyloid propensity of PRNP with resistance haplotype in Nigerian sheep. Our study provides valuable information that could be possibly adopted in programs targeted at breeding for scrapie resistance in sheep from tropical regions.

Multiple Origins and Genomic Basis of Complex Traits in Sighthounds

Sighthounds, a distinctive group of hounds comprising numerous breeds, have their origins rooted in ancient artificial selection of dogs. In this study, we performed genome sequencing for 123 sighthounds, including one breed from Africa, six breeds from Europe, two breeds from Russia, and four breeds and 12 village dogs from the Middle East. We gathered public genome data of five sighthounds and 98 other dogs as well as 31 gray wolves to pinpoint the origin and genes influencing the morphology of the sighthound genome. Population genomic analysis suggested that sighthounds originated from native dogs independently and were comprehensively admixed among breeds, supporting the multiple origins hypothesis of sighthounds. An additional 67 published ancient wolf genomes were added for gene flow detection. Results showed dramatic admixture of ancient wolves in African sighthounds, even more than with modern wolves. Whole-genome scan analysis identified 17 positively selected genes (PSGs) in the African population, 27 PSGs in the European population, and 54 PSGs in the Middle Eastern population. None of the PSGs overlapped in the three populations. Pooled PSGs of the three populations were significantly enriched in "regulation of release of sequestered calcium ion into cytosol" (gene ontology: 0051279), which is related to blood circulation and heart contraction. In addition, ESR1, JAK2, ADRB1, PRKCE, and CAMK2D were under positive selection in all three selected groups. This suggests that different PSGs in the same pathway contributed to the similar phenotype of sighthounds. We identified an ESR1 mutation (chr1: g.42,177,149 T > C) in the transcription factor (TF) binding site of Stat5a and a JAK2 mutation (chr1: g.93,277,007 T > A) in the TF binding site of Sox5. Functional experiments confirmed that the ESR1 and JAK2 mutation reduced their expression. Our results provide new insights into the domestication history and genomic basis of sighthounds.

Target capture sequencing for the first Nigerian genotype I ASFV genome

African swine fever (ASF) is a contagious viral disease that affects domestic pigs and wild boars, causing significant economic losses globally. After the first Nigerian outbreak in 1997, there have been frequent reports of ASF in pig-producing regions in the country. To facilitate control, it is important to understand the genotype and phylogenetic relationship of ASF viruses (ASFVs). Recent genetic analysis of Nigerian ASFV isolates has revealed the presence of both genotypes I and II; this is based on analysis of a few selected genes. Phylogenetic analysis of ASFV whole genomes highlights virus origins and evolution in greater depth. However, there is currently no information on the ASFV genome from Nigerian isolates. Two ASFV-positive samples were detected during a random survey of 150 Nigerian indigenous pig samples collected in 2016. We assembled near-complete genomes of the two ASFV-positive samples using in-solution hybrid capture sequencing. The genome-wide phylogenetic tree assigned these two genomes into p72 genotype I, particularly close to the virulent Benin 97/1 strain. The two ASFVs share 99.94 and 99.92 % genomic sequence identity to Benin97/1. This provides insight into the origin and relationship of ASFV strains from Nigeria and Italy. The study reports for the first time the determination of near-complete genomes of ASFV using in-solution hybrid capture sequencing, which represents an important advance in understanding the global evolutionary landscape of ASFVs.

Scrapie-associated polymorphisms of the prion protein gene (PRNP) in Nigerian native goats

Scrapie is a fatal prion protein disease stiffly associated with single nucleotide polymorphism (SNPs) of the prion protein gene (PRNP). The prevalence of this deadly disease has been reported in small ruminants, including goats. The Nigerian goats are hardy, trypano-tolerant, and contribute to the protein intake of the increasing population. Although scrapie has been reported in Nigerian goats, there is no study on the polymorphism of the PRNP gene. Herein, we evaluated the genetic and allele distributions of PRNP polymorphism in 132 Nigerian goats and compared them with publicly available studies on scrapie-affected goats. We utilized Polyphen-2, PROVEAN and AMYCO programs to examine structural variations produced by the non-synonymous SNPs. Our study revealed 29 SNPs in Nigerian goats, of which 14 were non-synonymous, and 23 were novel. There were significant differences (P < 0.001) in the allele frequencies of PRNP codons 139, 146, 154 and 193 in Nigerian goats compared with scrapie-affected goats, except for Northern Italian goats at codon 154. Based on the prediction by Polyphen-2, R139S and N146S were 'benign', R154H was 'probably damaging', and T193I was 'possibly damaging'. In contrast, PROVEAN predicted 'neutral' for all non-synonymous SNPs, while AMYCO showed a similar amyloid propensity of PRNP for resistant haplotype and two haplotypes of Nigerian goats. Our study is the first to investigate the polymorphism of scrapie-related genes in Nigerian goats.

Long amplicon HiFi sequencing for mitochondrial DNA genomes

Long-read sequencing technology is a powerful approach with application in various genetic and genomic research. Herein, we developed the pipeline for long amplicon high-fidelity (HiFi) sequencing and then applied it for sequencing mitochondrial DNA (mtDNA) genomes from pools of 79 Tibetan Mastiffs. We amplified the mtDNA genome with long-range PCR using two pairs of primers. Two rounds of circular consensus sequencing (CCS) were conducted and their accuracy was evaluated. The results indicate that the second round of CCS can improve the accuracy of HiFi reads. In addition, the analysis of 79 high-quality mtDNA genomes shows the Tibetan Mastiffs from outside of the Tibetan Plateau experienced hybridization with other dogs. The high quality reads generator (HQGR) software is provided to facilitate data analyses, which is publicly accessible on GitHub (https://github.com/Caizf-script/HQGR). Our long amplicon HiFi sequencing pipeline can also be applied in various target enrichment strategies for small genomes and candidate genes.

Population structure and evolutionary history of the greater cane rat (Thryonomys swinderianus) from the Guinean Forests of West Africa

Grasscutter (Thryonomys swinderianus) is a large-body old world rodent found in sub-Saharan Africa. The body size and the unique taste of the meat of this major crop pest have made it a target of intense hunting and a potential consideration as a micro-livestock. However, there is insufficient knowledge on the genetic diversity of its populations across African Guinean forests. Herein, we investigated the genetic diversity, population structures and evolutionary history of seven Nigerian wild grasscutter populations together with individuals from Cameroon, Republic of Benin, and Ghana, using five mitochondrial fragments, including D-loop and cytochrome b (CYTB). D-loop haplotype diversity ranged from 0.571 (± 0.149) in Republic of Benin to 0.921 (± 0.013) in Ghana. Within Nigeria, the haplotype diversity ranged from 0.659 (± 0.059) in Cross River to 0.837 (± 0.075) in Ondo subpopulation. The fixation index (F_ST), haplotype frequency distribution and analysis of molecular variance revealed varying levels of population structures across populations. No significant signature of population contraction was detected in the grasscutter populations. Evolutionary analyses of CYTB suggests that South African population might have diverged from other populations about 6.1 (2.6-10.18, 95% CI) MYA. Taken together, this study reveals the population status and evolutionary history of grasscutter populations in the region.

African Suid Genomes Provide Insights into the Local Adaptation to Diverse African Environments

African wild suids consist of several endemic species that represent ancient members of the family Suidae and have colonized diverse habitats on the African continent. However, limited genomic resources for African wild suids hinder our understanding of their evolution and genetic diversity. In this study, we assembled high-quality genomes of a common warthog (Phacochoerus africanus), a red river hog (Potamochoerus porcus), as well as an East Asian Diannan small-ear pig (Sus scrofa). Phylogenetic analysis showed that common warthog and red river hog diverged from their common ancestor around the Miocene/Pliocene boundary, putatively predating their entry into Africa. We detected species-specific selective signals associated with sensory perception and interferon signaling pathways in common warthog and red river hog, respectively, which contributed to their local adaptation to savannah and tropical rainforest environments, respectively. The structural variation and evolving signals in genes involved in T-cell immunity, viral infection, and lymphoid development were identified in their ancestral lineage. Our results provide new insights into the evolutionary histories and divergent genetic adaptations of African suids.

Chromosome-level genome assembly of the Muscovy duck provides insight into fatty liver susceptibility

The Muscovy duck (Cairina moschata) is an economically important poultry species, which is susceptible to fatty liver. Thus, the Muscovy duck may serve as an excellent candidate animal model of non-alcoholic fatty liver disease. However, the mechanisms underlying fatty liver development in this species are poorly understood. In this study, we report a chromosome-level genome assembly of the Muscovy duck, with a contig N50 of 11.8 Mb and scaffold N50 of 83.16 Mb. The susceptibility of Muscovy duck to fatty liver was mainly attributed to weak lipid catabolism capabilities (fatty acid β-oxidation and lipolysis). Furthermore, conserved noncoding elements (CNEs) showing accelerated evolution contributed to fatty liver formation by down-regulating the expression of genes involved in hepatic lipid catabolism. We propose that the susceptibility of Muscovy duck to fatty liver is an evolutionary by-product. In conclusion, this study revealed the potential mechanisms underlying the susceptibility of Muscovy duck to fatty liver.

Population Genetic Structure of Three Cichlids in ‎Ilorin, North-Central Nigeria

Cichlids are among the economically important which serve as a source of food for people ‎around the world. A deep understanding of the population structure and genetic diversity of ‎cichlids are vital for initiation of conservation policies and sustainable aquaculture. There is ‎paucity of information on the patterns of genetic variations among and within cichlids in ‎North-central Nigeria. This study, therefore, investigated population genetic structure of ‎Coptodon zillii, Oreochromis niloticus, and Hemichromis fasciatus collected from different ‎freshwater bodies in North-central Nigeria. Genomic DNA was extracted, and five highly ‎polymorphic RAPD primers were used for RAPD-PCR amplification and genotyping of the ‎fish. Genetic polymorphism within and between the three tilapia species were examined. ‎Percentages of polymorphism loci, pairwise population matrix, analysis of molecular variance ‎‎(AMOVA), and genetic distances of cichlid populations were determined using standard ‎methods, and dendrograms were constructed using an un-weighted pair group method of ‎arithmetic mean (UPGMA). Overall, percentages of estimated molecular variance within and ‎among C. zillii, H. fasciatus and O. niloticus populations were 5% and 95%; 4% and 96% and ‎‎13% and 87%, respectively. Our results suggest that the three cichlids have close evolutionary ‎relationship and there were no distinct genetic differences on the basis on sampling locations. ‎C. zillii and H. fasciatus are more genetically closer than O. niloticus. This study concludes that ‎RAPD is useful in studying the population genetic structure of cichlids. This study therefore ‎recommends conservation of genetic pool of cichlid species through proper maintenance and ‎restoration of polluted habitat to guarantee sustainable fishery production. However, markers ‎such as microsatellite DNA can be assayed in further studies for better results‎.

Genetic diversity and population structure of muscovy duck (Cairina moschata) from Nigeria

The domestic Muscovy duck (Cairina moschata) provide unique genetic resources patterned by both tropical environmental conditions and human activities, the evaluation of their genetic diversity and population structure will shade light on the mechanism of their remarkable adaptive capacities. We therefore analyzed the variation in mtDNA cytochrome b and nuclear DNA CYP2U1 sequences of 378 Nigerian Muscovy ducks (comprising of 287 de novo and 91 downloaded) plus 80 published sequences of Muscovy ducks from India. The results showed high haplotype diversity (0.800 ± 0.023) among Nigerian Muscovy duck populations with 91 distinct haplotypes for the nuclear DNA CYP2U1 gene but low (0.266 ± 0.033) for cytochrome b with 31 haplotypes. The median-joining networks of both markers grouped Nigerian Muscovy ducks into two; the first group consisting of only Nigerian Muscovy duck populations, and the second group Nigerian with Indian populations. Neutrality test results indicated that Nigerian populations experienced recent population expansion and/or genetic hitchhiking. A geographic signal was absent in line with previously studied poultry species in Nigeria. The most prominent haplotype dominated across all regions in Nigeria, which may be due to extensive genetic intermixing except for the Indian population (F _ST = 0.02550, P = 0.01075). This indicated low genetic differentiation between and within Nigerian Muscovy duck as revealed by the suitability of the nuclear DNA CYP2U1 gene.

Genome-wide investigations reveal the population structure and selection signatures of Nigerian cattle adaptation in the sub-Saharan tropics

Background

Cattle are considered to be the most desirable livestock by small scale farmers. In Africa, although comprehensive genomic studies have been carried out on cattle, the genetic variations in indigenous cattle from Nigeria have not been fully explored. In this study, genome-wide analysis based on genotyping-by-sequencing (GBS) of 193 Nigerian cattle was used to reveal new insights on the history of West African cattle and their adaptation to the tropical African environment, particularly in sub-Saharan region. 

Results

The GBS data were evaluated against whole-genome sequencing (WGS) data and high rate of variant concordance between the two platforms was evident with high correlated genetic distance matrices genotyped by both methods suggestive of the reliability of GBS applicability in population genetics. The genetic structure of Nigerian cattle was observed to be homogenous and unique from other African cattle populations. Selection analysis for the genomic regions harboring imprints of adaptation revealed genes associated with immune responses, growth and reproduction, efficiency of feeds utilization, and heat tolerance. Our findings depict potential convergent adaptation between African cattle, dogs and humans with adaptive genes SPRY2 and ITGB1BP1 possibly involved in common physiological activities.

Conclusion

The study presents unique genetic patterns of Nigerian cattle which provide new insights on the history of cattle in West Africa based on their population structure and the possibility of parallel adaptation between African cattle, dogs and humans in Africa which require further investigations.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08512-w.

Single-cell RNA-sequencing Reveals Thoracolumbar Vertebra Heterogeneity and Rib-genesis in Pigs

Thoracolumbar vertebra (TLV) and rib primordium (RP) development is a common evolutionary feature across vertebrates, although whole-organism analysis of TLV and RP gene expression dynamics has been lacking. Here we investigated the single-cell transcriptomic landscape of thoracic vertebra (TV), lumbar vertebra (LV), and RP cells from a pig embryo at 27 days post-fertilization (dpf) and identified six cell types with distinct gene-expression signatures. In-depth dissection of the gene-expression dynamics and RNA velocity revealed a coupled process of osteogenesis and angiogenesis during TLV and rib development. Further analysis of cell-type-specific and strand-specific expression uncovered the extremely high levels of Hoxa10 3'-untranslated region (UTR) sequence specific to osteoblast of LV cells, which may function as anti-Hoxa10-antisense by counteracting the Hoxa10-antisense effect to determine TLV transition. Thus, this work provides a valuable resource for understanding embryonic osteogenesis and angiogenesis underlying vertebrate TLV and RP development at the cell-type-specific resolution, which serves as a comprehensive view on the transcriptional profile of animal embryo development.

Genomes reveal selective sweeps in kiang and donkey for high-altitude adaptation

Over the last several hundred years, donkeys have adapted to high-altitude conditions on the Tibetan Plateau. Interestingly, the kiang, a closely related equid species, also inhabits this region. Previous reports have demonstrated the importance of specific genes and adaptive introgression in divergent lineages for adaptation to hypoxic conditions on the Tibetan Plateau. Here, we assessed whether donkeys and kiangs adapted to the Tibetan Plateau via the same or different biological pathways and whether adaptive introgression has occurred. We assembled a de novo genome from a kiang individual and analyzed the genomes of five kiangs and 93 donkeys (including 24 from the Tibetan Plateau). Our analyses suggested the existence of a strong hard selective sweep at the EPAS1 locus in kiangs. In Tibetan donkeys, however, another gene, i.e., EGLN1, was likely involved in their adaptation to high altitude. In addition, admixture analysis found no evidence for interspecific gene flow between kiangs and Tibetan donkeys. Our findings indicate that despite the short evolutionary time scale since the arrival of donkeys on the Tibetan Plateau, as well as the existence of a closely related species already adapted to hypoxia, Tibetan donkeys did not acquire adaptation via admixture but instead evolved adaptations via a different biological pathway.

Genomic Analyses of Unveil Helmeted Guinea Fowl (Numida meleagris) Domestication in West Africa

Domestication of the helmeted guinea fowl (HGF; Numida meleagris) in Africa remains elusive. Here we report a high-quality de novo genome assembly for domestic HGF generated by long- and short-reads sequencing together with optical and chromatin interaction mapping. Using this assembly as the reference, we performed population genomic analyses for newly sequenced whole-genomes for 129 birds from Africa, Asia, and Europe, including domestic animals (n = 89), wild progenitors (n = 34), and their closely related wild species (n = 6). Our results reveal domestication of HGF in West Africa around 1,300-5,500 years ago. Scanning for selective signals characterized the functional genes in behavior and locomotion changes involved in domestication of HGF. The pleiotropy and linkage in genes affecting plumage color and fertility were revealed in the recent breeding of Italian domestic HGF. In addition to presenting a missing piece to the jigsaw puzzle of domestication in poultry, our study provides valuable genetic resources for researchers and breeders to improve production in this species.

Population Genomics Reveals Incipient Speciation, Introgression, and Adaptation in the African Mona Monkey (Cercopithecus mona)

Guenons (tribe Cercopithecini) are the most widely distributed nonhuman primate in the tropical forest belt of Africa and show considerable phenotypic, taxonomic, and ecological diversity. However, genomic information for most species within this group is still lacking. Here, we present a high-quality de novo genome (total 2.90 Gb, contig N50 equal to 22.7 Mb) of the mona monkey (Cercopithecus mona), together with genome resequencing data of 13 individuals sampled across Nigeria. Our results showed differentiation between populations from East and West of the Niger River ∼84 ka and potential ancient introgression in the East population from other mona group species. The PTPRK, FRAS1, BNC2, and EDN3 genes related to pigmentation displayed signals of introgression in the East population. Genomic scans suggest that immunity genes such as AKT3 and IL13 (possibly involved in simian immunodeficiency virus defense), and G6PD, a gene involved in malaria resistance, are under positive natural selection. Our study gives insights into differentiation, natural selection, and introgression in guenons.

Genetic variation of Nigerian cattle inferred from maternal and paternal genetic markers

The African cattle provide unique genetic resources shaped up by both diverse tropical environmental conditions and human activities, the assessment of their genetic diversity will shade light on the mechanism of their remarkable adaptive capacities. We therefore analyzed the genetic diversity of cattle samples from Nigeria using both maternal and paternal DNA markers. Nigerian cattle can be assigned to 80 haplotypes based on the mitochondrial DNA (mtDNA) D-loop sequences and haplotype diversity was 0.985 + 0.005. The network showed two major matrilineal clustering: the dominant cluster constituting the Nigerian cattle together with other African cattle while the other clustered Eurasian cattle. Paternal analysis indicates only zebu haplogroup in Nigerian cattle with high genetic diversity 1.000 ± 0.016 compared to other cattle. There was no signal of maternal genetic structure in Nigerian cattle population, which may suggest an extensive genetic intermixing within the country. The absence of Bos indicus maternal signal in Nigerian cattle is attributable to vulnerability bottleneck of mtDNA lineages and concordance with the view of male zebu genetic introgression in African cattle. Our study shades light on the current genetic diversity in Nigerian cattle and population history in West Africa.

Patterns of ichthyofaunal diversity and distribution across Jebba Hydro-Electric Power (HEP) dam, Jebba, north-central Nigeria

The ichthyofauna diversity of the Jebba Hydroelectric Power (HEP) Dam, Jebba, North-central Nigeria was studied. Fishes were sampled for 24 months using gill net, hook and line, and cast net. Individuals were identified using morphological and molecular (mitochondrial Cytochrome c Oxidase subunit I) data. A total of 9605 freshwater fishes were recorded during the sampling period. The use of an integrative taxonomic approach enabled the identification of 83 species belonging to 42 genera. Additionally, the study recorded three unidentified species - Ctenopoma sp, Malapterurus sp., and Protopterus sp. Analyses showed that individuals belonging to families Cichlidae and Mochokidae dominated the dam. The diversity analyses revealed relatively high fish diversity during the rainy season at the downstream section of Jebba HEP dam compared to the upstream section. The study, therefore, showed the presence of a diverse fish community comprising high species richness and diversity across the Jebba HEP dam. Finally, we recommend proper biodiversity monitoring and assessment of freshwater fish diversity across Nigeria. In addition, the use of an integrated taxonomic approach is recommended for appropriate species' identification and studies of freshwater fishes from Nigeria.

Whole-genome sequencing reveals origin and evolution of influenza A(H1N1)pdm09 viruses in Lincang, China, from 2014 to 2018

The continuous variation of the seasonal influenza viruses, particularly A(H1N1)pdm09, persistently threatens human life and health around the world. In local areas of southwest china, the large time-scale genomic research on A(H1N1)pdm09 is still insufficient. Here, we sequenced 45 whole-genome sequences of influenza A(H1N1)pdm09 viruses in Lincang, China, from 2014 to 2018, by next-generation sequencing technology to characterize molecular mechanisms of their origin and evolution. Our phylogenetic analyses suggest that the A(H1N1)pdm09 strains circulating in Lincang belong to clade 6B and the subclade 6B.1A predominates in 2018. Further, the strains in 2018 possess elevated evolutionary rate as compared to strains in other years. Several newly emerged mutations for HA (hemagglutinin) in 2018 are revealed (i.e., S183P and R221K). Intriguingly, the substitution R221K falls into the RBS (receptor binding site) of HA protein, which could affect antigenic properties of influenza A(H1N1)pdm09 viruses, and another substitution S183P near to RBS with a high covering frequency (11/14 strains) in 2018 is exactly located at the epitope B. Notably, the NA (neuraminidase) protein harbors a new mutation I23T, potentially involved in N-glycosylation. Based on the background with a higher evolutionary rate in 2018 strains, we deeply evaluate the potential vaccine efficacy against Lincang strains and discover a substantive decline of the vaccine efficacy in 2018. Our analyses reaffirm that the real-time molecular surveillance and timely updated vaccine strains for prevention and control of influenza A(H1N1)pdm09 are crucial in the future.

Evolution and transition of expression trajectory during human brain development

BACKGROUND

The remarkable abilities of the human brain are distinctive features that set us apart from other animals. However, our understanding of how the brain has changed in the human lineage remains incomplete, but is essential for understanding cognition, behavior, and brain disorders in humans. Here, we compared the expression trajectory in brain development between humans and rhesus macaques (Macaca mulatta) to explore their divergent transcriptome profiles.

RESULTS

Results showed that brain development could be divided into two stages, with a demarcation date in a range between 25 and 26 postconception weeks (PCW) for humans and 17-23PCWfor rhesus macaques, rather than birth time that have been widely used as a uniform demarcation time of neurodevelopment across species. Dynamic network biomarker (DNB) analysis revealed that the two demarcation dates were transition phases during brain development, after which the brain transcriptome profiles underwent critical transitions characterized by highly fluctuating DNB molecules. We also found that changes between early and later brain developmental stages (as defined by the demarcation points) were substantially greater in the human brain than in the macaque brain. To explore the molecular mechanism underlying prolonged timing during early human brain development, we carried out expression heterochrony tests. Results demonstrated that compared to macaques, more heterochronic genes exhibited neoteny during early human brain development, consistent with the delayed demarcation time in the human lineage, and proving that neoteny in human brain development could be traced to the prenatal period. We further constructed transcriptional networks to explore the profile of early human brain development and identified the hub gene RBFOX1 as playing an important role in regulating early brain development. We also found RBFOX1 evolved rapidly in its non-coding regions, indicating that this gene played an important role in human brain evolution. Our findings provide evidence that RBFOX1 is a likely key hub gene in early human brain development and evolution.

CONCLUSIONS

By comparing gene expression profiles between humans and macaques, we found divergent expression trajectories between the two species, which deepens our understanding of the evolution of the human brain.

DNA Barcoding Silver Butter Catfish (Schilbe intermedius) Reveals Patterns of Mitochondrial Genetic Diversity Across African River Systems

The silver butter catfish (Schilbe intermedius) is widely distributed across African river systems. To date, information on its mitochondrial genetic diversity, population structure, and historical demography are not well-established. Herein, we combined newly generated mitochondrial cytochrome c oxidase (COI) subunit I gene sequences with previously published COI sequences in the global databases to reconstruct its phylogeography, population genetic structure, and historical demography. Results from the mtDNA phylogeography and species delimitation tests (Cluster algorithm – Species Identifier, Automatic Barcode Gap Discovery and Poison Tree Process model) revealed that S. intermedius comprises at least seven geographically defined matrilines. Although the overall haplotype diversity of S. intermedius was high (h = 0.90), results showed that East (Kenya) and West (Nigeria) African populations had low levels of haplotype diversity (h = ~0.40). In addition, population genetic polymorphism and historical demographics showed that S. intermedius populations in both East and West Africa underwent severe contractions as a result of biogeographic influences. The patterns of genetic diversity and population structure were consistent with adaptive responses to historical biogeographic factors and contemporary environmental variations across African river systems. This is suggestive of the influence of historical biogeographic factors and climatic conditions on population divergence of S. intermedius across African river systems. Given our discovery of previously underappreciated diversity within S. intermedius, we recommend that this species be considered for increased conservation and management.

Whole genome resequencing reveals an association of ABCC4 variants with preaxial polydactyly in pigs

Background

Polydactyly is one of the most common congenital limb dysplasia in many animal species. Although preaxial polydactyly (PPD) has been comprehensively studied in humans as a common abnormality, the genetic variations in other animal species have not been fully understood. Herein, we focused on the pig, as an even-toed ungulate mammal model with its unique advantages in medical and genetic researches, two PPD families consisting of four affected and 20 normal individuals were sequenced.

Results

Our results showed that the PPD in the sampled pigs were not related to previously reported variants. A strong association was identified at ABCC4 and it encodes a transmembrane protein involved in ciliogenesis. We found that the affected and normal individuals were highly differentiated at ABCC4, and all the PPD individuals shared long haplotype stretches as compared with the unaffected individuals. A highly differentiated missense mutation (I85T) in ABCC4 was observed at a residue from a transmembrane domain highly conserved among a variety of organisms.

Conclusions

This study reports ABCC4 as a new candidate gene and identifies a missense mutation for PPD in pigs. Our results illustrate a putative role of ciliogenesis process in PPD, coinciding with an earlier observation of ciliogenesis abnormality resulting in pseudo-thumb development in pandas. These results expand our knowledge on the genetic variations underlying PPD in animals.

Potential dual expansion of domesticated donkeys revealed by worldwide analysis on mitochondrial sequences

Molecular studies on donkey mitochondrial sequences have clearly defined two distinct maternal lineages involved in domestication. However, domestication histories of these two lineages remain enigmatic. We therefore compared several population characteristics between these two lineages based on global sampling, which included 171 sequences obtained in this study (including Middle Asian, East Asian, and African samples) plus 536 published sequences (including European, Asian, and African samples). The two lineages were clearly separated from each other based on whole mitochondrial genomes and partial non-coding displacement loop (D-loop) sequences, respectively. The Clade I lineage experienced an increase in population size more than 8 000 years ago and shows a complex haplotype network. In contrast, the population size of the Clade II lineage has remained relatively constant, with a simpler haplotype network. Although the distribution of the two lineages was almost equal across the Eurasian mainland, they still presented discernible but complex geographic bias in most parts of Africa, which are known as their domestication sites. Donkeys from sub-Saharan Africa tended to descend from the Clade I lineage, whereas the Clade II lineage was dominant along the East and North coasts of Africa. Furthermore, the migration routes inferred from diversity decay suggested different expansion across China between the two lineages. Altogether, these differences indicated non-simultaneous domestication of the two lineages, which was possibly influenced by the response of pastoralists to the desertification of the Sahara and by the social expansion and trade of ancient humans in Northeast Africa, respectively.

Population Genomics Analysis Revealed Origin and High-altitude Adaptation of Tibetan Pigs

Tibetan pig is native to the Qinghai-Tibet Plateau and has adapted to the high-altitude environmental condition such as hypoxia. However, its origin and genetic mechanisms underlying high-altitude adaptation still remain controversial and enigmatic. Herein, we analyze 229 genomes of wild and domestic pigs from Eurasia, including 63 Tibetan pigs, and detect 49.6 million high-quality variants. Phylogenomic and structure analyses show that Tibetan pigs have a close relationship with low-land domestic pigs in China, implying a common domestication origin. Positively selected genes in Tibetan pigs involved in high-altitude physiology, such as hypoxia, cardiovascular systems, UV damage, DNA repair. Three of loci with strong signals of selection are associated with EPAS1, CYP4F2, and THSD7A genes, related to hypoxia and circulation. We validated four non-coding mutations nearby EPAS1 and CYP4F2 showing reduced transcriptional activity in Tibetan pigs. A high-frequency missense mutation is found in THSD7A (Lys561Arg) in Tibetan pigs. The selective sweeps in Tibetan pigs was found in association with selection against non-coding variants, indicating an important role of regulatory mutations in Tibetan pig evolution. This study is important in understanding the evolution of Tibetan pigs and advancing our knowledge on animal adaptation to high-altitude environments.

Genetic variation and cryptic lineage diversity of the Nigerian red-headed rock agama Agama agama associate with eco-geographic zones

Nigeria is an Afrotropical region with considerable ecological heterogeneity and levels of biotic endemism. Among its vertebrate fauna, reptiles have broad distributions, thus, they constitute a compelling system for assessing the impact of ecological variation and geographic isolation on species diversification. The red-headed rock agama, Agama agama, lives in a wide range of habitats and, thus, it may show genetic structuring and diversification. Herein, we tested the hypothesis that ecology affects its genetic structure and population divergence. Bayesian inference phylogenetic analysis of a mitochondrial DNA (mtDNA) gene recovered four well-supported matrilines with strong evidence of genetic structuring consistent with eco-geographic regions. Genetic differences among populations based on the mtDNA also correlated with geographic distance. The ecological niche model for the matrilines had a good fit and robust performance. Population divergence along the environmental axes was associated with climatic conditions, and temperature ranked highest among all environmental variables for forest specialists, while precipitation ranked highest for the forest/derived savanna, and savanna specialists. Our results cannot reject the hypothesis that niche conservatism promotes geographic isolation of the western populations of Nigerian A. agama. Thus, ecological gradients and geographic isolation impact the genetic structure and population divergence of the lizards. This species might be facing threats due to recent habitat fragmentation, especially in western Nigeria. Conservation actions appear necessary.

DNA barcoding of economically important freshwater fish species from north-central Nigeria uncovers cryptic diversity

This study examines the utility of morphology and DNA barcoding in species identification of freshwater fishes from north-central Nigeria. We compared molecular data (mitochondrial cytochrome c oxidase subunit I (COI) sequences) of 136 de novo samples from 53 morphologically identified species alongside others in GenBank and BOLD databases. Using DNA sequence similarity-based (≥97% cutoff) identification technique, 50 (94.30%) and 24 (45.30%) species were identified to species level using GenBank and BOLD databases, respectively. Furthermore, we identified cases of taxonomic problems in 26 (49.00%) morphologically identified species. There were also four (7.10%) cases of mismatch in DNA barcoding in which our query sequence in GenBank and BOLD showed a sequence match with different species names. Using DNA barcode reference data, we also identified four unknown fish samples collected from fishermen to species level. Our Neighbor-joining (NJ) tree analysis recovers several intraspecific species clusters with strong bootstrap support (≥95%). Analysis uncovers two well-supported lineages within Schilbe intermedius. The Bayesian phylogenetic analyses of Nigerian S. intermedius with others from GenBank recover four lineages. Evidence of genetic structuring is consistent with geographic regions of sub-Saharan Africa. Thus, cryptic lineage diversity may illustrate species' adaptive responses to local environmental conditions. Finally, our study underscores the importance of incorporating morphology and DNA barcoding in species identification. Although developing a complete DNA barcode reference library for Nigerian ichthyofauna will facilitate species identification and diversity studies, taxonomic revisions of DNA sequences submitted in databases alongside voucher specimens are necessary for a reliable taxonomic and diversity inventory.

Genomic and transcriptomic analyses reveal selection of genes for puberty in Bama Xiang pigs

The Bama Xiang pig (BMX) is a famous early-maturing Chinese indigenous breed with a two-end black coat. To uncover the genetic basis of the BMX phenotype, we conducted comparative genomic analyses between BMX and East Asian wild boars and Laiwu pigs, respectively. Genes under positive selection were enriched in pathways associated with gonadal hormone and melanin synthesis, consistent with the phenotypic changes observed during development in BMX pigs. We also performed differentially expressed gene analysis based on RNA-seq data from pituitary tissues of BMX and Large White pigs. The CTTNBP2NL, FRS2, KANK4, and KATNAL1 genes were under selection and exhibited expressional changes in the pituitary tissue, which may affect BMX pig puberty. Our study demonstrated the positive selection of early maturity in the development of BMX pigs and advances our knowledge on the role of regulatory elements in puberty evolution in pigs.

PigVar: a database of pig variations and positive selection signatures

Pigs are excellent large-animal models for medical research and a promising organ donor source for transplant patients. Next-generation sequencing technology has yielded a dramatic increase in the volume of genomic data for pigs. However, the limited amount of variation data provided by dbSNP, and non-congruent criteria used for calling variation, present considerable hindrances to the utility of this data. We used a uniform pipeline, based on GATK, to identify non-redundant, high-quality, whole-genome SNPs from 280 pigs and 6 outgroup species. A total of 64.6 million SNPs were identified in 280 pigs and 36.8 million in the outgroups. We then used LUMPY to identify a total of 7 236 813 structural variations (SVs) in 211 pigs. Positively selected loci were identified through five statistical tests of different evolutionary attributes of the SNPs. Combining the non-redundant variations and the evolutionary selective scores, we built the first pig-specific variation database, PigVar (http://www.ibiomedical.net/pigvar/), which is a web-based open-access resource. PigVar collects parameters of the variations including summary lists of the locations of the variations within protein-coding and long intergenic non-coding RNA (lincRNA) genes, whether the SNPs are synonymous or non-synonymous, their ancestral and derived states, geographic sampling locations, as well as breed information. The PigVar database will be kept operational and updated to facilitate medical research using the pig as model and agricultural research including pig breeding.

Database URL: http://www.ibiomedical.net/pigvar/

Complete mitochondrial genome of Sri Lankan Junglefowl (Gallus lafayetti) and phylogenetic study

The complete mitochondrial genomes of two Sri Lankan junglefowl (Gallus lafayetti: CJF) individuals were sequenced by using next-generation sequencing technique. Samples were collected from Rathnapura and Pelmadulla areas in Sri Lanka. The complete mitochondrial DNA is 16,839 bp in length, with a typical mitogenome structure composed of a non-coding control region, 22 tRNA, two rRNA, and 13 protein-coding genes. Overall base composition is 30% A, 23.9% T, 32.3% C, and 13.6% G indicating high content of 54.0% A + T for both individuals. Phylogenetic analysis reveals that CJF samples cluster with the clade of the green junglefowl (Gallus varius) and red junglefowl (Gallus gallus) than to grey junglefowl (Gallus sonerattii: GyJF). This result can be subsequently used to provide essential information for junglefowl evolution.

Capture Hybridization of Long-Range DNA Fragments for High-Throughput Sequencing

Capture hybridization coupled with high-throughput sequencing (HTS) has become one of the most popular approaches to address some scientific problems not only for fundamental evolution but also for ecology and human disease in recent years. However, the technical problem of limited probe capture ability affects its widespread application. Here, we propose to capture hybridize long-range DNA fragments for HTS (termed LR-LCH). We provide a case of three amphibian samples to examine LR-LCH with 2 kb libraries and comparison of standard capture hybridization with 480 bp libraries. Capture sensitivity increased from an average 13.57% of standard capture hybridization to an average 19.80% of LR-LCH; capture efficiency also increased from an average 72.56% of standard capture hybridization to an average 97.71% of LR-LCH. These indicate that longer fragments in the library generally contain both relatively variable regions and relatively conservative regions. The divergent parts of target DNA are enriched along with conservative parts of DNA sequence that effectively captured during hybridization. We present a protocol that allows users to overcome the low capture sensitivity problem for high divergent regions.

Analysis of the genetic variation in mitochondrial DNA, Y-chromosome sequences, and MC1R sheds light on the ancestry of Nigerian indigenous pigs

Background

The history of pig populations in Africa remains controversial due to insufficient evidence from archaeological and genetic data. Previously, a Western ancestry for West African pigs was reported based on loci that are involved in the determination of coat color. We investigated the genetic diversity of Nigerian indigenous pigs (NIP) by simultaneously analyzing variation in mitochondrial DNA (mtDNA), Y-chromosome sequence and the melanocortin receptor 1 (MC1R) gene.

Results

Median-joining network analysis of mtDNA D-loop sequences from 201 NIP and previously characterized loci clustered NIP with populations from the West (Europe/North Africa) and East/Southeast Asia. Analysis of partial sequences of the Y-chromosome in 57 Nigerian boars clustered NIP into lineage HY1. Finally, analysis of MC1R in 90 NIP resulted in seven haplotypes, among which the European wild boar haplotype was carried by one individual and the European dominant black by most of the other individuals (93%). The five remaining unique haplotypes differed by a single synonymous substitution from European wild type, European dominant black and Asian dominant black haplotypes.

Conclusions

Our results demonstrate a European and East/Southeast Asian ancestry for NIP. Analyses of MC1R provide further evidence. Additional genetic analyses and archaeological studies may provide further insights into the history of African pig breeds. Our findings provide a valuable resource for future studies on whole-genome analyses of African pigs.

Electronic supplementary material

The online version of this article (doi:10.1186/s12711-017-0326-1) contains supplementary material, which is available to authorized users.

Mitochondrial DNA sequence variation in Iranian native dogs

The dog mtDNA diversity picture from wide geographical sampling but from a small number of individuals per region or breed, displayed little geographical correlation and high degree of haplotype sharing between very distant breeds. For a clear picture, we extensively surveyed Iranian native dogs (n = 305) in comparison with published European (n = 443) and Southwest Asian (n = 195) dogs. Twelve haplotypes related to haplogroups A, B and C were shared by Iranian, European, Southwest Asian and East Asian dogs. In Iran, haplotype and nucleotide diversities were highest in east, southeast and northwest populations while western population had the least. Sarabi and Saluki dog populations can be assigned into haplogroups A, B, C and D; Qahderijani and Kurdi to haplogroups A, B and C, Torkaman to haplogroups A, B and D while Sangsari and Fendo into haplogroups A and B, respectively. Evaluation of population differentiation using pairwise F_ST generally revealed no clear population structure in most Iranian dog populations. The genetic signal of a recent demographic expansion was detected in East and Southeast populations. Further, in accordance with previous studies on dog-wolf hybridization for haplogroup d2 origin, the highest number of d2 haplotypes in Iranian dog as compared to other areas of Mediterranean basin suggests Iran as the probable center of its origin. Historical evidence showed that Silk Road linked Iran to countries in South East Asia and other parts of the world, which might have probably influenced effective gene flow within Iran and these regions. The medium nucleotide diversity observed in Iranian dog calls for utilization of appropriate management techniques in increasing effective population size.

Genetic Variation of Goat Interferon Regulatory Factor 3 Gene and Its Implication in Goat Evolution

The immune systems are fundamentally vital for evolution and survival of species; as such, selection patterns in innate immune loci are of special interest in molecular evolutionary research. The interferon regulatory factor (IRF) gene family control many different aspects of the innate and adaptive immune responses in vertebrates. Among these, IRF3 is known to take active part in very many biological processes. We assembled and evaluated 1356 base pairs of the IRF3 gene coding region in domesticated goats from Africa (Nigeria, Ethiopia and South Africa) and Asia (Iran and China) and the wild goat (Capra aegagrus). Five segregating sites with θ value of 0.0009 for this gene demonstrated a low diversity across the goats’ populations. Fu and Li tests were significantly positive but Tajima’s D test was significantly negative, suggesting its deviation from neutrality. Neighbor joining tree of IRF3 gene in domesticated goats, wild goat and sheep showed that all domesticated goats have a closer relationship than with the wild goat and sheep. Maximum likelihood tree of the gene showed that different domesticated goats share a common ancestor and suggest single origin. Four unique haplotypes were observed across all the sequences, of which, one was particularly common to African goats (MOCH-K14-0425, Poitou and WAD). In assessing the evolution mode of the gene, we found that the codon model d_N/d_S ratio for all goats was greater than one. Phylogenetic Analysis by Maximum Likelihood (PAML) gave a ω₀ (d_N/d_S) value of 0.067 with LnL value of -6900.3 for the first Model (M1) while ω2 = 1.667 in model M2 with LnL value of -6900.3 with positive selection inferred in 3 codon sites. Mechanistic empirical combination (MEC) model for evaluating adaptive selection pressure on particular codons also confirmed adaptive selection pressure in three codons (207, 358 and 408) in IRF3 gene. Positive diversifying selection inferred with recent evolutionary changes in domesticated goat IRF3 led us to conclude that the gene evolution may have been influenced by domestication processes in goats.