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 (FST), 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

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

Accelerated evolution of constraint elements for hematophagic adaptation in mosquitoes

Comparative genomics is a powerful approach that comprehensively interprets the genome. Herein, we performed whole genome comparative analysis of 16 Diptera genomes, including four mosquitoes and 12 Drosophilae. We found more than 540 000 constraint elements (CEs) in the Diptera genome, with the majority found in the intergenic, coding and intronic regions. Accelerated elements (AEs) identified in mosquitoes were mostly in the protein-coding regions (>93%), which differs from vertebrates in genomic distribution. Some genes functionally enriched in blood digestion, body temperature regulation and insecticide resistance showed rapid evolution not only in the lineage of the recent common ancestor of mosquitoes (RCAM), but also in some mosquito lineages. This may be associated with lineage-specific traits and/or adaptations in comparison with other insects. Our findings revealed that although universally fast evolution acted on biological systems in RCAM, such as hematophagy, same adaptations also appear to have occurred through distinct degrees of evolution in different mosquito species, enabling them to be successful blood feeders in different environments

Complete mitochondrial genome of the Thai Red Junglefowl (Gallus gallus) and phylogenetic analysis

In this study, we sequenced the complete mitochondrial genome (mitogenome) of the Thai Red Junglefowl (RJF; Gallus gallus) using the next-generation sequencing (NGS) platform of the Ion Torrent PGM. Samples were taken from Mae Wang District, Chiang Mai Province, northern Thailand. Our data showed the complete mitogenome to be 16 785 bp in length, composed by 13 protein-coding genes, 22 tRNA genes, two rRNA genes, and one control region. The genome nucleotide composition was 30.3% A, 23.7% T, 32.5% C, and 13.5% G, resulting in a high percentage of A+T (50.4%). Phylogenetic analysis revealed that the mitogenome belonged to haplogroup X, whereas those of all domestic chickens belong to haplogroups A to G. This newly released mitogenome sequence will advance further evolutionary and population genetics study of the RJF and domestic chicken. The availability of the G. gallus mitogenome will also contribute to further conservation genetics research of a unique species, listed as ‘data deficient’ in Thailand

Genetic Diversity Relationship in Azakheli Buffalo Inferred from mtDNA and MC1R Sequences Comparison

Azakheli is relatively smaller riverine breed with a very peculiar characteristics kept under unique traditional husbandry practices in comparison with rest of the Pakistani buffalo breeds; however, milk production is comparable. The present study was conducted to evaluate the genetic diversity of Azakheli breed. A total of sixty-six blood samples were collected for the amplification of mtDNA D-loop region and MC1R gene sequencing analysis. Median-joining network analysis of 191 mtDNA D-loop sequences of Azakheli and eight Indian riverine buffalo breeds clustered into three clades. Ancient Azakheli Region 1 clade was the oldest with the highest mutation steps and was present close to the root of UPGMA phylogenetic tree. There was 5 mutated lines distance between Pakistan buffalo and Indian riverine buffaloes. The populations of neighboring countries did not share any haplotypes with Azakheli buffalo of Pakistan. Possibly, residing for so long in the cold atmosphere and high elevation regions caused the mutation in mtDNA D-loop, though these conditions did not affect the overall performance of Azakheli as milch buffalo breed of Pakistan. MC1R analyses showed high mutations in Azakheli of Albino phenotype and all the black phenotype individuals of Azakheli buffalo share haplotypes with dominant Chinese and Indian black phenotypes buffaloes in MC1R median-joining network, indicating the reason of black coat color is due to MC1R gene. The haplotype diversity and nucleotide diversity was (H. 0.923, Pi: 0.00895) in Azakheli. Current results illustrated Asian ancestry for Azakheli buffalo, and mtDNA and MC1R analyses provided further evidence. Additional genetic analyses and archeological studies may provide further insight into the domestication period and history of Azakheli buffalo breed. The further studies are required on different coat colors with different genes on Azakheli buffalo to understand the phenotype variation

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

Abstract 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

Mitochondrial DNA sequence variation in Iranian native dogs

<p>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 (<i>n</i> = 305) in comparison with published European (<i>n</i> = 443) and Southwest Asian (<i>n</i> = 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 <i>F</i>_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.</p

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

Abstract 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