Xiaoyong Du

Sequencing and automated whole-genome optical mapping of the genome of a domestic goat ( Capra hircus )

We report the ∼2.66-Gb genome sequence of a female Yunnan black goat. The sequence was obtained by combining short-read sequencing data and optical mapping data from a high-throughput whole-genome mapping instrument. The whole-genome mapping data facilitated the assembly of super-scaffolds >5× longer by the N50 metric than scaffolds augmented by fosmid end sequencing (scaffold N50 = 3.06 Mb, super-scaffold N50 = 16.3 Mb). Super-scaffolds are anchored on chromosomes based on conserved synteny with cattle, and the assembly is well supported by two radiation hybrid maps of chromosome 1. We annotate 22,175 protein-coding genes, most of which were recovered in the RNA-seq data of ten tissues. Comparative transcriptomic analysis of the primary and secondary follicles of a cashmere goat reveal 51 genes that are differentially expressed between the two types of hair follicles. This study, whose results will facilitate goat genomics, shows that whole-genome mapping technology can be used for the de novo assembly of large genomes.

Positive selection drives the evolution of endocrine regulatory bone morphogenetic protein system in mammals

The rapid evolution of reproductive proteins might be driven by positive Darwinian selection. The bone morphogenetic protein family is the largest within the transforming growth factor (TGF) superfamily. A little have been known about the molecular evolution of bone morphogenetic proteins exhibiting potential role in mammalian reproduction. In this study we investigated mammalian bone morphogenetic proteins using maximum likelihood approaches of codon substitutions to identify positive Darwinian selection in various species. The proportion of positively selected sites was tested by different likelihood models for individual codon, and M8 were found to be the best model. The percentage of positively elected sites under M8 are 2.20% with ω = 1.089 for BMP2, 1.6% with ω = 1.61 for BMP 4 0.53% for BMP15 with ω = 1.56 and 0.78% for GDF9 with ω = 1.93. The percentage of estimated selection sites under M8 is strong statistical confirmation that divergence of bone morphogenetic proteins is driven by Darwinian selection. For the proteins, model M8 was found significant for all proteins with ω > 1. To further test positive selection on particular amino acids, the evolutionary conservation of amino acid were measured based on phylogenetic linkage among sequences. For exploring the impact of these somatic substitution mutations in the selection region on human cancer, we identified one pathogenic mutation in human BMP4 and one in BMP15, possibly causing prostate cancer and six neutral mutations in BMPs. The comprehensive map of selection results allows the researchers to perform systematic approaches to detect the evolutionary footprints of selection on specific gene in specific species.

Exploring the Genetic Resistance to Gastrointestinal Nematodes Infection in Goat Using RNA-Sequencing

Gastrointestinal nematodes (GINs) are one of the most economically important parasites of small ruminants and a major animal health concern in many regions of the world. However, the molecular mechanisms of the host response to GIN infections in goat are still little known. In this study, two genetically distinct goat populations, one relatively resistant and the other susceptible to GIN infections, were identified in Yichang goat and then four individuals in each group were chosen to compare mRNA expression profiles using RNA-seq. Field experiment showed lower worm burden, delayed and reduced egg production in the relatively resistant group than the susceptible group. The analysis of RNA-seq showed that 2369 genes, 1407 of which were up-regulated and 962 down-regulated, were significantly (p < 0.001) differentially expressed between these two groups. Functional annotation of the 298 genes more highly expressed in the resistant group yielded a total of 46 significant (p < 0.05) functional annotation clusters including 31 genes (9 in innate immunity, 13 in immunity, and 9 in innate immune response) related to immune biosynthetic process as well as transforming growth factor (TGF)-β, mitogen-activated protein kinase (MAPK), and cell adhesion molecules (CAMs) pathways. Our findings provide insights that are immediately relevant for the improvement of host resistance to GIN infections and which will make it possible to know the mechanisms underlying the resistance of goats to GIN infections.

Genome-wide association analysis identifies potential regulatory genes for eumelanin pigmentation in chicken plumage

Plumage color in chicken is determined by the proportion of eumelanin and pheomelanin pigmentation. As the main ingredient in plumage melanin, eumelanin plays a key role in the dark black, brown and grey coloration. However, very few studies have been performed to identify the related genes and mutations on a genome-wide scale. Herein, a resource family consisting of one backcross population and two F2 cross populations between a black roster and Yukou Brown I parent stockbreed was constructed for identification of genes related to eumelanin pigmentation. Chickens with eumelanin in their plumage were classified as the case group, and the rest were considered the control group. A genome-wide association study of this phenotype and genotypes using Affymetrix 600K HD SNP arrays in this F2 family revealed 13 significantly associated SNPs and in 10 separate genes on chromosomes 1, 2, 3 and 5. Based on previous studies in model species, we inferred that genes, including NUAK family kinase 1 (NUAK1) and sonic hedgehog (SHH), may play roles in the development of neural crest cells or melanoblasts during the embryonic period, which may also affect the eumelanin pigmentation. Our results facilitate the understanding of the genetic basis of eumelanin pigmentation in chicken plumage.

Discovering All Transcriptome Single-Nucleotide Polymorphisms and Scanning for Selection Signatures in Ducks (Anas platyrhynchos)

The duck is one of the most economically important waterfowl as a source of meat, eggs, and feathers. Characterizing the genetic variation in duck species is an important step toward linking genes or genomic regions with phenotypes. Human-driven selection during duck domestication and subsequent breed formation has likely left detectable signatures in duck genome. In this study, we employed a panel of >1.4 million single-nucleotide polymorphisms (SNPs) identified from the RNA sequencing (RNA-seq) data of 15 duck individuals. The density of the resulting SNPs is significantly positively correlated with the density of genes across the duck genome, which demonstrates that the usage of the RNA-seq data allowed us to enrich variant functional categories, such as coding exons, untranslated regions (UTRs), introns, and downstream/upstream. We performed a complete scan of selection signatures in the ducks using the composite likelihood ratio (CLR) and found 76 candidate regions of selection, many of which harbor genes related to phenotypes relevant to the function of the digestive system and fat metabolism, including TCF7L2, EIF2AK3, ELOVL2, and fatty acid-binding protein family. This study illustrates the potential of population genetic approaches for identifying genomic regions affecting domestication-related phenotypes and further helps to increase the known genetic information about this economically important animal.

Transcriptome sequencing reveals key potential long non-coding RNAs related to duration of fertility trait in the uterovaginal junction of egg-laying hens

Duration of fertility, (DF) is an important functional trait in poultry production and lncRNAs have emerged as important regulators of various process including fertility. In this study we applied a genome-guided strategy to reconstruct the uterovaginal junction (UVJ) transcriptome of 14 egg-laying birds with long- and short-DF (n = 7); and sought to uncover key lncRNAs related to duration of fertility traits by RNA-sequencing technology. Examination of RNA-seq data revealed a total of 9977 lncRNAs including 2576 novel lncRNAs. Differential expression (DE) analysis of lncRNA identified 223 lncRNAs differentially expressed between the two groups. DE-lncRNA target genes prediction uncovered over 200 lncRNA target genes and functional enrichment tests predict a potential function of DE-lncRNAs. Gene ontology classification and pathway analysis revealed 8 DE-lncRNAs, with the majority of their target genes enriched in biological functions such as reproductive structure development, developmental process involved in reproduction, response to cytokine, carbohydrate binding, chromatin organization, and immune pathways. Differential expression of lncRNAs and target genes were confirmed by qPCR. Together, these results significantly expand the utility of the UVJ transcriptome and our analysis identification of key lncRNAs and their target genes regulating DF will form the baseline for understanding the molecular functions of lncRNAs regulating DF.

De novo assembly of mitochondrial genomes provides insights into genetic diversity and molecular evolution in wild boars and domestic pigs

Up to date, the scarcity of publicly available complete mitochondrial sequences for European wild pigs hampers deeper understanding about the genetic changes following domestication. Here, we have assembled 26 de novo mtDNA sequences of European wild boars from next generation sequencing (NGS) data and downloaded 174 complete mtDNA sequences to assess the genetic relationship, nucleotide diversity, and selection. The Bayesian consensus tree reveals the clear divergence between the European and Asian clade and a very small portion (10 out of 200 samples) of maternal introgression. The overall nucleotides diversities of the mtDNA sequences have been reduced following domestication. Interestingly, the selection efficiencies in both European and Asian domestic pigs are reduced, probably caused by changes in both selection constraints and maternal population size following domestication. This study suggests that de novo assembled mitogenomes can be a great boon to uncover the genetic turnover following domestication. Further investigation is warranted to include more samples from the ever-increasing amounts of NGS data to help us to better understand the process of domestication.

Detecting the Population Structure and Scanning for Signatures of Selection in Horses (Equus caballus) From Whole-Genome Sequencing Data

Animal domestication gives rise to gradual changes at the genomic level through selection in populations. Selective sweeps have been traced in the genomes of many animal species, including humans, cattle, and dogs. However, little is known regarding positional candidate genes and genomic regions that exhibit signatures of selection in domestic horses. In addition, an understanding of the genetic processes underlying horse domestication, especially the origin of Chinese native populations, is still lacking. In our study, we generated whole genome sequences from 4 Chinese native horses and combined them with 48 publicly available full genome sequences, from which 15 341 213 high-quality unique single-nucleotide polymorphism variants were identified. Kazakh and Lichuan horses are 2 typical Asian native breeds that were formed in Kazakh or Northwest China and South China, respectively. We detected 1390 loss-of-function (LoF) variants in protein-coding genes, and gene ontology (GO) enrichment analysis revealed that some LoF-affected genes were overrepresented in GO terms related to the immune response. Bayesian clustering, distance analysis, and principal component analysis demonstrated that the population structure of these breeds largely reflected weak geographic patterns. Kazakh and Lichuan horses were assigned to the same lineage with other Asian native breeds, in agreement with previous studies on the genetic origin of Chinese domestic horses. We applied the composite likelihood ratio method to scan for genomic regions showing signals of recent selection in the horse genome. A total of 1052 genomic windows of 10 kB, corresponding to 933 distinct core regions, significantly exceeded neutral simulations. The GO enrichment analysis revealed that the genes under selective sweeps were overrepresented with GO terms, including “negative regulation of canonical Wnt signaling pathway,” “muscle contraction,” and “axon guidance.” Frequent exercise training in domestic horses may have resulted in changes in the expression of genes related to metabolism, muscle structure, and the nervous system.

Characterization of the Promoter Regions of Two Sheep Keratin-Associated Protein Genes for Hair Cortex-Specific Expression.

The keratin-associated proteins (KAPs) are the structural proteins of hair fibers and are thought to play an important role in determining the physical properties of hair fibers. These proteins are activated in a striking sequential and spatial pattern in the keratinocytes of hair fibers. Thus, it is important to elucidate the mechanism that underlies the specific transcriptional activity of these genes. In this study, sheep KRTAP 3–3 and KRTAP11-1 genes were found to be highly expressed in wool follicles in a tissue-specific manner. Subsequently, the promoter regions of the two genes that contained the 5′ flanking/5′ untranslated regions and the coding regions were cloned. Using an in vivo transgenic approach, we found that the promoter regions from the two genes exhibited transcriptional activity in hair fibers. A much stronger and more uniformly expressed green fluorescent signal was observed in the KRTAP11-1-ZsGreen1 transgenic mice. In situ hybridization revealed the symmetrical expression of sheep KRTAP11-1 in the entire wool cortex. Consistently, immunohistochemical analysis demonstrated that the pattern of ZsGreen1 expression in the hair cortex of transgenic mice matches that of the endogenous KRTAP11-1 gene, indicating that the cloned promoter region contains elements that are sufficient to govern the wool cortex-specific transcription of KRTAP11-1. Furthermore, regulatory regions in the 5′ upstream sequence of the sheep KRTAP11-1 gene that may regulate the observed hair keratinocyte specificity were identified using in vivo reporter assays.

ROLE OF CYTOKINES AND CORRELATED CELLS IN HOST IMMUNITY AGAINST NEMATODE IN RUMINANTS

Nematode contagion is a core impediment to the profitable dairy production for livestock farms. Gastrointestinal (GI) parasitism causes weight loss and low milk production, along with high mortality in sheep and goat. The current intensive anthelmintics therapies to control parasitic burden resulted in the appearance of drug resistant parasitic strains. Due to high demand of unadulterated animal products free from drugs motivated alternative strategies for improvement; including breeding plans for parasite control in ruminants. The increase of protective resistance to nematode infections, host genome show particular expression that is frequently confused by mechanisms concurrently essential to control multiple nematodes species as well as protozoan ecto-parasites, viral and microbial pathogens. The involved molecular mechanisms under these developments correspond to crucial steps toward improvement of efficient new parasite control strategies. Knowledge of various immunity methods of host and regulation of development of parasite, physiology, and virulence is able to identify the objective of parasite control. This review recapitulate current evolution and restrictions of optimistic regulatory biological pathways and genetic networks that concern with susceptibility and host resistance to infection of GI nematode in ruminants.