Cis-eQTLs in seven duck tissues identify novel candidate genes for growth and carcass traits.

BACKGROUND: Expression quantitative trait loci (eQTL) studies aim to understand the influence of genetic variants on gene expression. The colocalization of eQTL mapping and GWAS strategy could help identify essential candidate genes and causal DNA variants vital to complex traits in human and many farm animals. However, eQTL mapping has not been conducted in ducks. It is desirable to know whether eQTLs within GWAS signals contributed to duck economic traits. RESULTS: In this study, we conducted an eQTL analysis using publicly available RNA sequencing data from 820 samples, focusing on liver, muscle, blood, adipose, ovary, spleen, and lung tissues. We identified 113,374 cis-eQTLs for 12,266 genes, a substantial fraction 39.1% of which were discovered in at least two tissues. The cis-eQTLs of blood were less conserved across tissues, while cis-eQTLs from any tissue exhibit a strong sharing pattern to liver tissue. Colocalization between cis-eQTLs and genome-wide association studies (GWAS) of 50 traits uncovered new associations between gene expression and potential loci influencing growth and carcass traits. SRSF4, GSS, and IGF2BP1 in liver, NDUFC2 in muscle, ELF3 in adipose, and RUNDC1 in blood could serve as the candidate genes for duck growth and carcass traits. CONCLUSIONS: Our findings highlight substantial differences in genetic regulation of gene expression across duck primary tissues, shedding light on potential mechanisms through which candidate genes may impact growth and carcass traits. Furthermore, this availability of eQTL data offers a valuable resource for deciphering further genetic association signals that may arise from ongoing extensive endeavors aimed at enhancing duck production traits.

Genetic fine-mapping reveals single nucleotide polymorphism mutations in the MC1R regulatory region associated with duck melanism.

Birds are among the most colourful terrestrial vertebrates, with various plumage colours and patterns. We conducted a genome-wide association study (GWAS) on an intercross F2 population of Pekin ducks and mallards (n = 722) and identified a 1.57-Mb genetic region (Chr11: 20,176,480-21,750,101 bp) related to duck melanism. Fine mapping by linkage disequilibrium (LD) and FST analysis narrowed the final candidate region to a region of 22,500 bp (Chr11: 20,677,500-20,700,000 bp) including three coding genes, TCF25, MC1R and TUBB3. Combined with transcriptome and qRT-PCR analysis, MC1R was identified as the unique genetic locus responsible for black plumage in ducks, and it was significantly more highly expressed in the feather bulbs of black ducks. We also identified 52G > A (Chr11: 20,696,354G > A) and 376G > A (Chr11: 20,696,678G > A) mutations in the MC1R coding region that have been widely studied in ducks. In addition, structural variations (SVs) were screened by nanopore sequencing, and no significant SV was found to be associated with the duck black plumage trait. However, we identified four novel single nucleotide polymorphisms in the MC1R regulator region (Chr11: 20,678,412G > A, Chr11: 20,679,236G > A, Chr11: 20,692,496 A > G and Chr11: 20,692,791 A > G) that had a strong association with the black plumage phenotype of ducks and combined with potential changes in transcription binding affinities. The luciferase reporter gene assay demonstrated that Chr11: 20,678,412G > A and Chr11: 20,679,236G > A led to significant promoter activity changes. Our research emphasizes the importance of MC1R regulatory region mutation in determining the duck black plumage phenotype, and these results expand our understanding of the genetic mechanism underlying duck plumage colour.

Genome-wide association studies demonstrate the genes associated with perimysial thickness in ducks.

The thickness of the perimysium has an essential effect on the tenderness of the meat. However, the genetic basis underlying perimysial thickness has not been determined. The objective of this study was to explore the quantitative trait loci (QTL) that influence perimysial thickness in an F2 segregating population generated by Mallard × Pekin duck using the genome-wide association study (GWAS) method. Two QTL identified in chromosomes 27 and 13 displayed significant associations with perimysial thickness traits at the genome-wide level. The strongest association was the QTL located in chromosome 27, and this region had an effect on perimysial thickness and contained a promising candidate gene MAGI3 (Membrane-associated guanylate kinase, WW and PDZ domain containing 3). Meanwhile, association analysis showed that the top SNP within the MAGI3 gene was also associated with intramuscular fat content traits, which showed that perimysial thickness was positively correlated with intramuscular fat content. The second strongest association was the QTL region of chromosome 13. SUCLG2 (Succinate-CoA ligase GDP-forming subunit beta) is proximal to the top SNP and stood out as another candidate gene. Furthermore, the Transposase-Accessible Chromatin using Sequencing result showed that some key transcription factors (MYF5, MYOD1, KLF11) related to muscle development or energy metabolism might bind to the open regions of MAGI3 and SUCLG2. By analyzing the expression of different genotypes of the candidate gene, we speculate that different genotypes of MAGI3 may have an effect on breast muscle development, and then affect the thickness of the perimysium. This study maps two major genes of the duck breast muscle perimysial thickness trait, which helps to characterize muscle development and contributes to the genetic improvement of meat yield and quality in livestock.

Genome-wide association study reveals the genetic determinism of serum biochemical indicators in ducks.

BACKGROUND: The serum is rich in nutrients and plays an essential role in electrolyte and acid-base balance, maintaining cellular homeostasis. In addition, serum parameters have been commonly used as essential biomarkers for clinical diagnosis. However, little is known about the genetic mechanism of the serum parameters in ducks. RESULTS: This study measured 18 serum parameters in 320 samples of the F2 segregating population generated by Mallard × Pekin duck. The phenotypic correlations showed a high correlation between LDH, HBDH, AST, and ALT (0.59-0.99), and higher coefficients were also observed among TP, ALB, HDL-C, and CHO (0.46-0.87). And then, we performed the GWAS to reveal the genetic basis of the 18 serum biochemical parameters in ducks. Fourteen candidate protein-coding genes were identified with enzyme traits (AST, ALP, LDH, HBDH), and 3 protein-coding genes were associated with metabolism and protein-related serum parameters (UA, TG). Moreover, the expression levels of the above candidate protein-coding genes in different stages of breast muscle and different tissues were analyzed. Furthermore, the genes located within the high-LD region (r2 > 0.4 and - log10(P) < 4) neighboring the significant locus also remained. Finally, 86 putative protein-coding genes were used for GO and KEGG enrichment analysis, the enzyme-linked receptor protein signaling pathway and ErbB signaling pathway deserve further focus. CONCLUSIONS: The obtained results can contribute to new insights into blood metabolism and provide new genetic biomarkers for application in duck breeding programs.

Genome-wide association analysis reveal the genetic reasons affect melanin spot accumulation in beak skin of ducks.

BACKGROUND: Skin pigmentation is a broadly appearing phenomenon of most animals and humans in nature. Here we used a bird model to investigate why melanin spot deposits on the skin. RESULTS: Our result showed that growth age and the sunlight might induce melanin deposition in bird beak skin which was determined by genetic factors. GWAS helped us to identify two major loci affecting melanin deposition, located on chromosomes 13 and 25, respectively. The fine mapping works narrowed the candidate regions to 0.98 Mb and 1.0 Mb on chromosomes 13 and 25. The MITF and POU2F3 may be the causative genes and synergistically affect melanin deposition during duck beak skin. Furthermore, our data strongly demonstrated that the pathway of melanin metabolism contributes to melanin deposition on the skin. CONCLUSIONS: We demonstrated that age and sunlight induce melanin deposition in bird beak skin, while heredity is fundamental. The MITF and POU2F3 likely played a synergistic effect on the regulation of melanin synthesis, and their mutations contribute to phenotypic differences in beak melanin deposition among individuals. It is pointed out that melanin deposition in the skin is related to the pathway of melanin metabolism, which provided insights into the molecular regulatory mechanisms and the genetic improvement of the melanin deposition in duck beak.

Severe pantothenic acid deficiency induces alterations in the intestinal mucosal proteome of starter Pekin ducks.

BACKGROUND: Pantothenic acid deficiency (PAD) results in growth depression and intestinal hypofunction of animals. However, the underlying molecular mechanisms remain to be elucidated. Mucosal proteome might reflect dietary influences on physiological processes. RESULTS: A total of 128 white Pekin ducks of one-day-old were randomly assigned to two groups, fed either a PAD or a pantothenic acid adequate (control, CON) diet. After a 16-day feeding period, two ducks from each replicate were sampled to measure plasma parameters, intestinal morphology, and mucosal proteome. Compared to the CON group, high mortality, growth retardation, fasting hypoglycemia, reduced plasma insulin, and oxidative stress were observed in the PAD group. Furthermore, PAD induced morphological alterations of the small intestine indicated by reduced villus height and villus surface area of duodenum, jejunum, and ileum. The duodenum mucosal proteome of ducks showed that 198 proteins were up-regulated and 223 proteins were down-regulated (> 1.5-fold change) in the PAD group compared to those in the CON group. Selected proteins were confirmed by Western blotting. Pathway analysis of these proteins exhibited the suppression of glycolysis and gluconeogenesis, fatty acid beta oxidation, tricarboxylic acid cycle, oxidative phosphorylation, oxidative stress, and intestinal absorption in the PAD group, indicating impaired energy generation and abnormal intestinal absorption. We also show that nine out of eleven proteins involved in regulation of actin cytoskeleton were up-regulated by PAD, probably indicates reduced intestinal integrity. CONCLUSION: PAD leads to growth depression and intestinal hypofunction of ducks, which are associated with impaired energy generation, abnormal intestinal absorption, and regulation of actin cytoskeleton processes. These findings provide insights into the mechanisms of intestinal hypofunction induced by PAD.

A single nucleotide polymorphism variant located in the cis-regulatory region of the ABCG2 gene is associated with mallard egg colour.

Avian egg coloration is shaped by natural selection, but its genetic basis remains unclear. Here, we used genome-wide association analysis and identity by descent to finely map green egg colour to a 179-kb region of Chr4 based on the resequencing of 352 ducks (Anas platyrhynchos) from a segregating population resulting from the mating of Pekin ducks (white-shelled eggs) and mallards (green-shelled eggs). We further narrowed the candidate region to a 30-kb interval by comparing genome divergence in seven indigenous duck populations. Among the genes located in the finely mapped region, only one transcript of the ABCG2 gene (XM_013093252.2) exhibited higher uterine expression in green-shelled individuals than in white-shelled individuals, as supported by transcriptome data from four populations. ABCG2 has been reported to encode a protein that functions as a membrane transporter for biliverdin. Sanger sequencing of the whole 30-kb candidate region (Chr4: 47.41-47.44 Mb) and a plasmid reporter assay helped to identify a single nucleotide polymorphism (Chr4: 47,418,074 G>A) located in a conserved predicted promoter region whose variation may alter ABCG2 transcription activity. We provide a useful molecular marker for duck breeding and contribute data to the research on ecological evolution based on egg colour patterns among birds.

Genome-wide association studies demonstrate that TASP1 contributes to increased muscle fiber diameter.

Muscle fiber diameter is an economically important trait because it affects meat yield and quality. However, the genetic basis underlying muscle fiber diameter has not been determined. In this study, we collected THREE muscular histological phenotypes in 479 ducks from an F2 segregating population generated by mallard × Pekin duck crosses. We performed genome-wide association studies (GWAS) and identified a quantitative trait locus (QTL) significantly associated with muscle fiber diameter on chromosome 3. Then, we discovered the selection signatures using the fixation index among 40 mallards and 30 Pekin ducks in this QTL region. Furthermore, we characterized the recombination event in this QTL region and identified a 6-kb block located on TASP1 that was significantly associated with muscle fiber diameter. Finally, five SNPs were screened as potential causative mutations within the 6-kb block. In conclusion, we demonstrated that TASP1 contributes to an increase in muscle fiber diameter, which helps to characterize muscle development and contributes to the genetic improvement of meat yield and quality in livestock.

Genome-wide association and selective sweep analyses reveal genetic loci for FCR of egg production traits in ducks.

BACKGROUND: As a major economic trait in poultry, egg production efficiency attracts widespread interest in breeding and production. However, limited information is available about the underlying genetic architecture of egg production traits in ducks. In this paper, we analyzed six egg production-related traits in 352 F2 ducks derived from reciprocal crosses between mallard and Pekin ducks. RESULTS: Feed conversation ratio (FCR) was positively correlated with feed intake but negatively correlated with egg-related traits, including egg weight and egg production, both phenotypically and genetically. Estimates of pedigree-based heritability were higher than 0.2 for all traits investigated, except hip-width. Based on whole-genome sequencing data, we conducted genome-wide association studies to identify genomic regions associated with these traits. In total, 11 genomic regions were associated with FCR. No genomic regions were identified as significantly associated with hip-width, total feed intake, average daily feed intake, and total egg production. Analysis of selective sweeps between mallard and Pekin ducks confirmed three of these genomic regions on chromosomes 13, 3 and 6. Within these three regions, variants in candidate genes that were in linkage disequilibrium with the GWAS leader single nucleotide polymorphisms (SNPs) (Chr13:2,196,728, P = 7.05 × 10-14; Chr3:76,991,524, P = 1.06 × 10-12; Chr6:20,356,803, P = 1.14 × 10-10) were detected. Thus, we identified 31 potential candidate genes associated with FCR, among which the strongest candidates are those that are highly expressed in tissues involved in reproduction and nervous system functions of ducks: CNTN4, CRBR, GPR63, KLHL32, FHL5, TRNT1, MANEA, NDUFAF4, and SCD. CONCLUSIONS: For the first time, we report the identification of genomic regions that are associated with FCR in ducks and our results illustrate the genomic changes that occurred during their domestication and are involved in egg production efficiency.

Genetic variations for egg internal quality of ducks revealed by genome-wide association study.

Egg internal quality traits are important traits related to egg production in poultry industry. To better understand the genetic architecture of egg internal quality traits in ducks, we performed genetic parameters estimates and a genome-wide association study (GWAS). The phenotypic values of egg weight, yolk color, albumin height (AH), yolk weight, and Haugh unit (HU) were collected individually from 352 F2 laying ducks produced by reciprocal crosses between mallards and Pekin ducks, and their genotypes were assayed by whole genome re-sequencing. The results showed that the AH and HU traits have a clear coefficient of variance, around 15% for both mallards and Pekin ducks. The pedigree-based genetic parameters estimates rane from 0.26 to 0.71 for all eight egg quality traits, while the highest heritability was 0.71 for egg weight. The GWAS showed that a clear signal was associated with AH and HU traits. The locus zoom analysis and conditional GWAS helped to narrow the candidate region to ~5.8-Mb spanning from 14.7 to 20.5 Mb on Chromosome 5, which harbored 111 candidate genes. MUC6 and LDLRAD3 were finally promised as the major candidate genes affecting albumen composition. Our data revealed the egg internal quality traits for the first time in ducks, which provides a theoretical basis and technological support for improving duck egg internal quality.

Dynamic accumulation of fatty acids in duck (Anas platyrhynchos) breast muscle and its correlations with gene expression.

BACKGROUND: Fatty acid composition contributes greatly to the quality and nutritional value of meat. However, the molecular regulatory mechanisms underlying fatty acid accumulation in poultry have not yet been cleared. The aims of this study were to characterize the dynamics of fatty acid accumulation in duck breast muscle and investigate its correlations with gene expression. RESULTS: Here, we analyzed the fatty acid profile and transcriptome of breast muscle derived from Pekin ducks and mallards at the ages of 2 weeks, 4 weeks, 6 weeks and 8 weeks. Twenty fatty acids were detected in duck breast muscle, with palmitic acid (C16:0, 16.6%~ 21.1%), stearic acid (C18:0, 9.8%~ 17.7%), oleic acid (C18:1n-9, 15.7%~ 33.8%), linoleic acid (C18:2n-6, 10.8%~ 18.9%) and arachidonic acid (C20:4n-6, 11.7%~ 28.9%) as the major fatty acids. Our results showed that fatty acid composition was similar between the two breeds before 6 weeks, but the compositions diverged greatly after this point, mainly due to the stronger capacity for C16:0 and C18:1n-9 deposition in Pekin ducks. By comparing the multistage transcriptomes of Pekin ducks and mallards, we identified 2025 differentially expressed genes (DEGs). Cluster analysis of these DEGs revealed that the genes involved in oxidative phosphorylation, fatty acid degradation and the PPAR signaling pathway were upregulated in mallard at 8 weeks. Moreover, correlation analysis of the DEGs and fatty acid composition traits suggested that the DEGs involved in lipogenesis, lipolysis and fatty acid ß-oxidation may interact to influence the deposition of fatty acids in duck breast muscle. CONCLUSIONS: We reported the temporal progression of fatty acid accumulation and the dynamics of the transcriptome in breast muscle of Pekin ducks and mallards. Our results provide insights into the transcriptome regulation of fatty acid accumulation in duck breast muscle, and will facilitate improvements of fatty acid composition in duck breeding.

Rapid evolution of a retro-transposable hotspot of ovine genome underlies the alteration of BMP2 expression and development of fat tails.

BACKGROUND: Sheep have developed the ability to store fat in their tails, which is a unique way of reserving energy to survive a harsh environment. However, the mechanism underlying this adaptive trait remains largely unsolved. RESULTS: In the present study, we provide evidence for the genetic determinants of fat tails, based on whole genome sequences of 89 individual sheep. A genome-wide scan of selective sweep identified several candidate loci including a region at chromosome 13, a haplotype of which underwent rapid evolution and spread through fat-tailed populations in China and the Middle East. Sequence analysis revealed an inter-genic origin of this locus, which later became a hotspot of ruminant-specific retro-transposon named BovB. Additionally, the candidate locus was validated based on a fat- and thin-tailed cross population. The expression of an upstream gene BMP2 was differentially regulated between fat-tailed and thin-tailed individuals in tail adipose and several other tissue types. CONCLUSIONS: Our findings suggest the fixation of fat tails in domestic sheep is caused by a selective sweep near a retro-transposable hotspot at chromosome 13, the diversity of which specifically affects the expression of BMP2. The present study has shed light onto the understanding of fat metabolism.

Functional evolution of phosphatidylethanolamine binding proteins in soybean and Arabidopsis.

Gene duplication provides resources for novel gene functions. Identification of the amino acids responsible for functional conservation and divergence of duplicated genes will strengthen our understanding of their evolutionary course. Here, we conducted a systemic functional investigation of phosphatidylethanolamine binding proteins (PEBPs) in soybean (Glycine max) and Arabidopsis thaliana. Our results demonstrated that after the ancestral duplication, the lineage of the common ancestor of the FLOWERING LOCUS T (FT) and TERMINAL FLOWER1 (TFL1) subfamilies functionally diverged from the MOTHER OF FT AND TFL1 (MFT) subfamily to activate flowering and repress flowering, respectively. They also underwent further specialization after subsequent duplications. Although the functional divergence increased with duplication age, we observed rapid functional divergence for a few pairs of young duplicates in soybean. Association analysis between amino acids and functional variations identified critical amino acid residues that led to functional differences in PEBP members. Using transgenic analysis, we validated a subset of these differences. We report clear experimental evidence for the functional evolution of the PEBPs in the MFT, FT, and TFL1 subfamilies, which predate the origin of angiosperms. Our results highlight the role of amino acid divergence in driving evolutionary novelty after duplication.

Expanding the BLUP alphabet for genomic prediction adaptable to the genetic architectures of complex traits.

Improvement of statistical methods is crucial for realizing the potential of increasingly dense genetic markers. Bayesian methods treat all markers as random effects, exhibit an advantage on dense markers, and offer the flexibility of using different priors. In contrast, genomic best linear unbiased prediction (gBLUP) is superior in computing speed, but only superior in prediction accuracy for extremely complex traits. Currently, the existing variety in the BLUP method is insufficient for adapting to new sequencing technologies and traits with different genetic architectures. In this study, we found two ways to change the kinship derivation in the BLUP method that improve prediction accuracy while maintaining the computational advantage. First, using the settlement under progressively exclusive relationship (SUPER) algorithm, we substituted all available markers with estimated quantitative trait nucleotides (QTNs) to derive kinship. Second, we compressed individuals into groups based on kinship, and then used the groups as random effects instead of individuals. The two methods were named as SUPER BLUP (sBLUP) and compressed BLUP (cBLUP). Analyses on both simulated and real data demonstrated that these two methods offer flexibility for evaluating a variety of traits, covering a broadened realm of genetic architectures. For traits controlled by small numbers of genes, sBLUP outperforms Bayesian LASSO (least absolute shrinkage and selection operator). For traits with low heritability, cBLUP outperforms both gBLUP and Bayesian LASSO methods. We implemented these new BLUP alphabet series methods in an R package, Genome Association and Prediction Integrated Tool (GAPIT), available at http://zzlab.net/GAPIT .

Selection response and estimation of the genetic parameters for multidimensional measured breast meat yield related traits in a long-term breeding Pekin duck line.

OBJECTIVE: This study was conducted to estimate the genetic parameters and breeding values of breast meat related traits of Pekin ducks. Selection response was also determined by using ultrasound breast muscle thickness (BMT) measurements in combination with bosom breadth (BB) and keel length (KL) values. METHODS: The traits analyzed were breast meat weight (BMW), body weight (BW), breast meat percentage (BMP) and the three parameters of breast meat (BB, KL, and BMT). These measurements were derived from studying 15,781 Pekin ducks selected from 10 generations based on breast meat weight. Genetic parameters and breeding value were estimated for the analysis of the breeding process. RESULTS: Estimated heritability of BMW and BMP were moderate (0.23 and 0.16, respectively), and heritability of BW was high (0.48). Other traits such as BB, KL, and BMT indicated moderate heritability ranging between 0.11 and 0.28. Significant phenotypic correlations of BMW with BW and BMP were discovered (p<0.05), and genetic correlations of BMW with BW and BMP were positive and high (0.83 and 0.66, respectively). It was noted that BMW had positive correlations with all the other traits. Generational average estimated breeding values of all traits increased substantially over the course of selection, which demonstrated that the ducks responded efficiently to increased breast meat yield after 10 generations of breeding. CONCLUSION: The results indicated that duck BMW had the potential to be increased through genetic selection with positive effects on BW and BMP. The ultrasound BMT, in combination with the measurement of BB and KL, is shown to be essential and effective in the process of high breast meat yield duck breeding.

Global investigation of the co-evolution of MIRNA genes and microRNA targets during soybean domestication.

Although the selection of coding genes during plant domestication has been well studied, the evolution of MIRNA genes (MIRs) and the interaction between microRNAs (miRNAs) and their targets in this process are poorly understood. Here, we present a genome-wide survey of the selection of MIRs and miRNA targets during soybean domestication and improvement. Our results suggest that, overall, MIRs have higher evolutionary rates than miRNA targets. Nonetheless, they do demonstrate certain similar evolutionary patterns during soybean domestication: MIRs and miRNA targets with high expression and duplication status, and with greater numbers of partners, exhibit lower nucleotide divergence than their counterparts without these characteristics, suggesting that expression level, duplication status, and miRNA-target interaction are essential for evolution of MIRs and miRNA targets. Further investigation revealed that miRNA-target pairs that are subjected to strong purifying selection have greater similarities than those that exhibited genetic diversity. Moreover, mediated by domestication and improvement, the similarities of a large number of miRNA-target pairs in cultivated soybean populations were increased compared to those in wild soybeans, whereas a small number of miRNA-target pairs exhibited decreased similarity, which may be associated with the adoption of particular domestication traits. Taken together, our results shed light on the co-evolution of MIRs and miRNA targets during soybean domestication.

Global dissection of alternative splicing in paleopolyploid soybean.

Alternative splicing (AS) is common in higher eukaryotes and plays an important role in gene posttranscriptional regulation. It has been suggested that AS varies dramatically among species, tissues, and duplicated gene families of different sizes. However, the genomic forces that govern AS variation remain poorly understood. Here, through genome-wide identification of AS events in the soybean (Glycine max) genome using high-throughput RNA sequencing of 28 samples from different developmental stages, we found that more than 63% of multiexonic genes underwent AS. More AS events occurred in the younger developmental stages than in the older developmental stages for the same type of tissue, and the four main AS types, exon skipping, intron retention, alternative donor sites, and alternative acceptor sites, exhibited different characteristics. Global computational analysis demonstrated that the variations of AS frequency and AS types were significantly correlated with the changes of gene features and gene transcriptional level. Further investigation suggested that the decrease of AS within the genome-wide duplicated genes were due to the diminution of intron length, exon number, and transcriptional level. Altogether, our study revealed that a large number of genes were alternatively spliced in the soybean genome and that variations in gene structure and transcriptional level may play important roles in regulating AS.

The research progress of farm animal genomics based on sequencing technologies.

Various farm animal breeds have been domesticated and bred for thousands years, and they provide adequate animal-derived proteins to meet the human nutrition requirement. Although quantitative genetics was applied in animal breeding, which launched a technological revolution in the past century, a number of complex traits remain difficult to be selected based on pedigree derived breeding, due to complicated animal genetics and development mechanisms. Farm animal's genetic potential hasn't yet to be fully exploited. The concept and technology from the Human Genome Project have greatly promoted farm animal genomic researches. It is possible to fine map the causal variations at the whole genome level and then exploit their biological functions, thus providing the theoretical basis for molecular designed breeding. In this review, we summarize the genomics research progress of main farm animals during the past decade, including pigs, cattle, yaks, goats, sheep, chickens, ducks and geese. We focus on the reference genome sequencing and follow-up population-level genomic studies based on high throughput resequencing technologies, and meanwhile envision the future work of farm animal genomics.

Mapping the four-horned locus and testing the polled locus in three Chinese sheep breeds.

Four-horned sheep are an ideal animal model for illuminating the genetic basis of horn development. The objective of this study was to locate the genetic region responsible for the four-horned phenotype and to verify a previously reported polled locus in three Chinese breeds. A genome-wide association study (GWAS) was performed using 34 two-horned and 32 four-horned sheep from three Chinese indigenous breeds: Altay, Mongolian and Sishui Fur sheep. The top two significant single nucleotide polymorphisms (SNPs) associated with the four-horned phenotype were both located in a region spanning positions 132.6 to 132.7 Mb on sheep chromosome 2. Similar locations for the four-horned trait were previously identified in Jacob, Navajo-Churro, Damara and Sishui Fur sheep, suggesting a common genetic component underlying the four-horned phenotype. The two identified SNPs were both downstream of the metaxin 2 (MTX2) gene and the HOXD gene cluster. For the top SNP-OAR2:g.132619300G>A-the strong associations of the AA and AG genotypes with the four-horned phenotype and the GG genotype with the two-horned phenotype indicated the dominant inheritance of the four-horned trait. No significant SNPs for the polled phenotype were identified in the GWAS analysis, and a PCR analysis for the detection of the 1.8-kb insertion associated with polled sheep in other breeds failed to verify the association with polledness in the three Chinese breeds. This study supports the hypothesis that two different loci are responsible for horn existence and number. This study contributes to the understanding of the molecular regulation of horn development and enriches the knowledge of qualitative traits in domestic animals.

Concerted evolution of D1 and D2 to regulate chlorophyll degradation in soybean.

Polyploidy is a common phenomenon, particularly in plants. The soybean (Glycine max [L.] Merr.) genome has undergone two whole genome duplication (WGD) events. The conservation and divergence of duplicated gene pairs are major contributors to genome evolution. D1 and D2 are two unlinked, paralogous nuclear genes, whose double-recessive mutant (d1d1d2d2) results in chlorophyll retention, called 'stay-green'. Through molecular cloning and functional analyses, we demonstrated that D1 and D2 are homologs of the STAY-GREEN (SGR) genes from other plant species and were duplicated as a result of the most recent WGD in soybean. Transcriptional analysis showed that both D1 and D2 were more highly expressed in older tissues, and chlorophyll degradation and programmed cell death-related genes were suppressed in a d1d2 double mutant, this situation indicated that these genes are probably involved in the early stages of tissue senescence. Investigation of genes that flank D1 and D2 revealed that evolution within collinear duplicated blocks may affect the conservation of individual gene pairs within the blocks. Moreover, we found that a long terminal repeat retrotransposon, GmD2IN, resulted in the d2 mutation. Further analysis of this retrotransposon family showed that insertion in or near the coding regions can affect gene expression or splicing patterns, and may be an important force to promote the divergence of duplicated gene pairs.