Small intestinal epithelium homeostasis involves four principal cell types: enterocytes, goblet, enteroendocrine and Paneth cells. Epidermal growth factor (EGF) has been shown to affect enterocyte differentiation. This study determined the effect of dietary EGF on goblet, enteroendocrine and Paneth cell differentiation in piglet small intestine and potential mechanisms. Forty-two weaned piglets were used in a 2 × 3 factorial design; the major factors were time post-weaning (days 7 and 14) and dietary treatment (0, 200 or 400 µg/kg EGF supplementation). The numbers of goblet and enteroendocrine cells were generally greater with the increase in time post-weaning. Moreover, the supplementation of 200 µg/kg EGF increased (P < 0.01) the number of goblet and enteroendocrine cells in villus and crypt of the piglet small intestine as compared with the control. Dietary supplementation with 200 µg/kg EGF enhanced (P < 0.05) abundances of differentiation-related genes atonal homologue 1, mucin 2 and intestinal trefoil factor 3 messenger RNA (mRNA) as compared with the control. Piglets fed 200 or 400 µg/kg EGF diet had increased (P < 0.05) abundances of growth factor-independent 1, SAM pointed domain containing ETS transcription factor and pancreatic and duodenal homeobox 1 mRNA, but decreased the abundance (P < 0.01) of E74 like ETS transcription factor 3 mRNA as compared with the control. Animals receiving 400 µg/kg EGF diets had enhanced (P < 0.05) abundances of neurogenin3 and SRY-box containing gene 9 mRNA as compared with the control. The mRNA abundance and protein expression of lysozyme, a marker of Paneth cell, were also increased (P < 0.05) in those animals. As compared with the control, dietary supplementation with 200 µg/kg EGF increased the abundance of EGF receptor mRNA and the ratio of non-phospho(p)-β-catenin/β-catenin (P < 0.05) in villus epithelial cells at days 7 and 14. This ratio in crypt epithelial cells was higher (P < 0.05) on the both 200 and 400 µg/kg EGF groups during the same period. Our results demonstrated that dietary EGF stimulated goblet, enteroendocrine and Paneth cell differentiation in piglets during the post-weaning period, partly through EGFR and Wnt/β-catenin signalling.

Grape seed procyanidins (GSPs), widely known for their beneficial health properties, fail to bring about the expected improvement in piglets’ growth performance. The effects of dietary supplementation with GSPs on nutrient utilisation may be a critical influencing factor. Hence, the purpose of this study was to investigate the effects of dietary supplementation with GSPs on nutrient utilisation and gut function in weaned piglets. One hundred and twenty crossbred piglets were allocated randomly to four treatment groups, with three replicate pens per treatment and 10 piglets per pen. Each group was given one of the four dietary treatments: the basal diet (control group) or the basal diet with the addition of 50-, 100- or 150-mg/kg GSPs. The trial lasted 28 days. Faeces were collected from d 12 to 14 and from d 26 to 28 for measuring the coefficient of total tract apparent digestibility (CTTAD) of the nutrients. Blood samples were collected on d 14 and 28 for detecting the blood biochemical parameters. Two piglets per pen were slaughtered to collect the pancreas and intestinal digesta for evaluating the digestive enzyme activity and the coefficient of ileal apparent digestibility (CIAD) of the nutrients. On d 14 and 28, supplementation with 150-mg/kg GSPs significantly decreased the CTTAD of DM and CP in piglets. On d 14, GSPs supplementation at a concentration of 150 mg/kg led to a remarkable decrease in the CIAD of CP and gross energy (GE). On d 28, GSPs supplementation at a dose of 150 mg/kg generated a marked decline in the CIAD of DM, GE, CP and ether extract. Grape seed procyanidins supplementation at concentrations of 100 or 150 mg/kg inhibited the activities of lipase and amylase. In contrast, the jejunum mucosa maltase and sucrase activities increased due to the inclusion of GSPs at a concentration of 100 mg/kg in the piglet diet. Compared with the levels of the control group, the serum glucose and total protein levels were enhanced significantly by supplementation with GSPs at 100 mg/kg and reduced dramatically at 150 mg/kg. The serum diamine oxidase activity and endotoxin levels were decreased by GSPs supplementation in piglet diets. In conclusion, higher concentrations of GSPs in weaned piglet diets attenuated nutrient digestion and inhibited digestive enzyme activity; however, suitable concentrations of GSPs could promote brush-border enzyme activity, enhance serum glucose and total protein concentrations and decrease epithelial permeability.

The COllaborative project of Development of Anthropometrical measures in Twins (CODATwins) project is a large international collaborative effort to analyze individual-level phenotype data from twins in multiple cohorts from different environments. The main objective is to study factors that modify genetic and environmental variation of height, body mass index (BMI, kg/m2) and size at birth, and additionally to address other research questions such as long-term consequences of birth size. The project started in 2013 and is open to all twin projects in the world having height and weight measures on twins with information on zygosity. Thus far, 54 twin projects from 24 countries have provided individual-level data. The CODATwins database includes 489,981 twin individuals (228,635 complete twin pairs). Since many twin cohorts have collected longitudinal data, there is a total of 1,049,785 height and weight observations. For many cohorts, we also have information on birth weight and length, own smoking behavior and own or parental education. We found that the heritability estimates of height and BMI systematically changed from infancy to old age. Remarkably, only minor differences in the heritability estimates were found across cultural–geographic regions, measurement time and birth cohort for height and BMI. In addition to genetic epidemiological studies, we looked at associations of height and BMI with education, birth weight and smoking status. Within-family analyses examined differences within same-sex and opposite-sex dizygotic twins in birth size and later development. The CODATwins project demonstrates the feasibility and value of international collaboration to address gene-by-exposure interactions that require large sample sizes and address the effects of different exposures across time, geographical regions and socioeconomic status.

We apply two methods to estimate the 21-cm bispectrum from data taken within the Epoch of Reionisation (EoR) project of the Murchison Widefield Array (MWA). Using data acquired with the Phase II compact array allows a direct bispectrum estimate to be undertaken on the multiple redundantly spaced triangles of antenna tiles, as well as an estimate based on data gridded to the uv-plane. The direct and gridded bispectrum estimators are applied to 21 h of high-band (167–197 MHz; z = 6.2–7.5) data from the 2016 and 2017 observing seasons. Analytic predictions for the bispectrum bias and variance for point-source foregrounds are derived. We compare the output of these approaches, the foreground contribution to the signal, and future prospects for measuring the bispectra with redundant and non-redundant arrays. We find that some triangle configurations yield bispectrum estimates that are consistent with the expected noise level after 10 h, while equilateral configurations are strongly foreground-dominated. Careful choice of triangle configurations may be made to reduce foreground bias that hinders power spectrum estimators, and the 21-cm bispectrum may be accessible in less time than the 21-cm power spectrum for some wave modes, with detections in hundreds of hours.

Potential planting area for tuber mustard was simulated using the Maxent model under current and future conditions based on 591 coordinates and 22 environmental layers. Model accuracy was excellent, with area under the receiving operator curve values of 0.967 and 0.958 for model training and testing, respectively. Dominant factors were mean diurnal range, mean temperature of the coldest quarter, annual mean temperature and minimum temperature of the coldest month, with thresholds of 6.5–7.5, 5.5–9, 16–19 and 2.0–6.5 °C, respectively. Under current conditions, suitable habitat areas (2.16% of total land in China) were concentrated mainly in Central, Southwest and East China, which can be defined as three occurrence and diffusion centres. In the 2050s and 2070s, suitable habitat areas are predicted to change to 3.72 and 3.92%, and 3.60 and 3.73% under scenarios RCP4.5 and RCP6.0, respectively, indicating that suitable habitat areas will increase slightly. However, future distribution of tuber mustard was predicted to differ among provinces or cities, i.e. predicted suitable habitat areas in Sichuan Province increased up to the 2050s but remained relatively unchanged between the 2050s and 2070s; in Chongqing city they first increased and then decreased; in Hunan, Anhui, Jiangsu, Zhejiang and Fujian Provinces they increased continuously; and in Guizhou, Hubei, Jiangxi Provinces and Shanghai city they first decreased, and then increased. The results from the current study provide useful information for management decisions of tuber mustard.

Laser interaction with an ultra-thin pre-structured target is investigated with the help of both two-dimensional and three-dimensional particle-in-cell simulations. With the existence of a periodic structure on the target surface, the laser seems to penetrate through the target at its fundamental frequency even if the plasma density of the target is much higher than the laser’s relativistically critical density. The particle-in-cell simulations show that the transmitted laser energy behind the pre-structured target is increased by about two orders of magnitude compared to that behind the flat target. Theoretical analyses show that the transmitted energy behind the pre-structured target is actually re-emitted by electron ‘islands’ formed by the surface plasma waves on the target surfaces. In other words, the radiation with the fundamental frequency is actually ‘surface emission’ on the target rear surface. Besides the intensity of the component with the fundamental frequency, the intensity of the high-order harmonics behind the pre-structured target is also much enhanced compared to that behind the flat target. The enhancement of the high-order harmonics is also related to the surface plasma waves generated on the target surfaces.

Sodium and chloride are the key factors maintaining normal osmotic pressure (OSM) and volume of the extracellular fluid, and influencing the acid–base balance of body fluids. The experiment was conducted to investigate the effects of dietary Na+ and Cl− level on growth performance, excreta moisture, blood biochemical parameters, intestinal Na+–glucose transporter 1 (SGLT1) messenger RNA (mRNA), and Na+–H+ exchanger 2 (NHE2) mRNA, and to estimate the optimal dietary sodium and chlorine level for yellow-feathered chickens from 22 to 42days. A total of 900 22-day-old Lingnan yellow-feathered male chickens were randomly allotted to five treatments, each of which included six replicates of 30 chickens per floor pen. The basal control diet was based on corn and soybean meal (without added NaCl and NaHCO3). Treatments 2 to 5 consisted of the basal diet supplemented with equal weights of Na+ and Cl−, constituting 0.1%, 0.2%, 0.3% and 0.4% of the diets. Supplemental dietary Na+ and Cl− improved the growth performance (P<0.05). Average daily gain (ADG) showed a quadratic broken-line regression to increasing dietary Na+ and Cl− (R2=0.979, P<0.001), and reached a plateau at 0.1%. Supplemental Na+ and Cl− increased (P<0.05) serum Na+ and OSM in serum and showed a quadratic broken-line regression (R2=0.997, P=0.004) at 0.11%. However, supplemental Na+ and Cl− decreased (P<0.05) serum levels of K+, glucose (GLU) and triglyceride. Higher levels of Na+and Cl− decreased duodenal NHE2 transcripts (P<0.05), but had no effect on ileal SGLT1 transcripts. The activity of Na+ /K+-ATPase in the duodenum decreased (P<0.05) with higher levels of dietary Na+ and Cl−. In conclusion, the optimal dietary Na+ and Cl− requirements for yellow-feathered chickens in the grower phase, from 22 to 42 days of age, to optimize ADG, serum Na+, OSM, K+ and GLU were 0.10%, 0.11%, 0.11%,0.17% and 0.16%, respectively, by regression analysis.

A multichannel calorimeter system is designed and constructed which is capable of delivering single-shot and broad-band spectral measurement of terahertz (THz) radiation generated in intense laser–plasma interactions. The generation mechanism of backward THz radiation (BTR) is studied by using the multichannel calorimeter system in an intense picosecond laser–solid interaction experiment. The dependence of the BTR energy and spectrum on laser energy, target thickness and pre-plasma scale length is obtained. These results indicate that coherent transition radiation is responsible for the low-frequency component ( ${<}$ 1 THz) of BTR. It is also observed that a large-scale pre-plasma primarily enhances the high-frequency component ( ${>}$ 3 THz) of BTR.

A stable reference gene is a key prerequisite for accurate assessment of gene expression. At present, the real-time reverse transcriptase quantitative polymerase chain reaction has been widely used in the analysis of gene expression in a variety of organisms. Neoseiulus barkeri Hughes (Acari: Phytoseiidae) is a major predator of mites on many important economically crops. Until now, however, there are no reports evaluating the stability of reference genes in this species. In view of this, we used GeNorm, NormFinder, BestKeeper, and RefFinder software tools to evaluate the expression stability of 11 candidate reference genes in developmental stages and under various abiotic stresses. According to our results, β-ACT and Hsp40 were the top two stable reference genes in developmental stages. The Hsp60 and Hsp90 were the most stable reference genes in various acaricides stress. For alterations in temperature, Hsp40 and α-TUB were the most suitable reference genes. About UV stress, EF1α and α-TUB were the best choice, and for the different prey stress, β-ACT and α-TUB were best suited. In normal conditions, the β-ACT and α-TUB were the two of the highest stable reference genes to respond to all kinds of stresses. The current study provided a valuable foundation for the further analysis of gene expression in N. barkeri.

Cultivated pastures in southern China are being used to improve forage productivity and animal performance, but studies on grazing behaviour of goats in these cultivated pastures are still rare. In the current study, the grazing behaviour of Yunling black goats under low (5 goats/ha) and high (15 goats/ha) stocking rates (SRs) was evaluated. Data showed that the proportion of time goats spent on activities was: eating (0.59–0.87), ruminating (0.05–0.35), walking (0.03–0.06) and resting (0.01–0.03). Compared with low SR, goats spent more time eating and walking, and less time ruminating and resting under high SR. Goats had similar diet preferences under both SR and preferred to eat grasses (ryegrass and cocksfoot) more than a legume (white clover). The distribution of eating time on each forage species was more uniform under high v. low SR. Bites/step, bite weight and daily intake were greater under low than high SR. Results suggest that the SR affects grazing behaviour of goats on cultivated pasture, and identifying an optimal SR is critical for increasing bite weight and intake.

Using an age-structured process-based simulation model for diamondback moth (DBM), we model the population dynamics of this major Brassica pest using the cropping practices and climate of Guangdong, China. The model simulates two interacting sub-populations (demes), each representing a short season crop. The simulated DBM abundance, and hence pest problems, depend on planting regime, crop hygiene and biological control. A continuous supply of hosts, a low proportion of crop harvested and long residue times between harvest and replanting each exacerbate pest levels. Biological control provided by a larval parasitoid can reduce pest problems, but not eliminate them when climate is suitable for DBM and under certain planting practices. The classic Integrated Pest Management (IPM) method of insecticide application, when pest threshold is reached, proved effective and halved the number of insecticide sprays when compared with the typical practice of weekly insecticide application.

Thousand Island Lake (TIL) is a typical fragmented landscape and an ideal model to study ecological effects of fragmentation. Partial fragments of the mitochondrial cytochrome oxidase subunit I gene of 23 island populations of Dendrolimus punctatus in TIL were sequenced, 141 haplotypes being identified. The number of haplotypes increased significantly with the increase in island area and shape index, whereas no significant correlation was detected between three island attributes (area, shape and isolation) and haplotype diversity. However, the correlation with number of haplotypes was no longer significant when the ‘outlier’ island JSD (the largest island) was not included. Additionally, we found no significant relationship between geographic distance and genetic distance. Geographic isolation did not obstruct the gene flow among D. punctatus populations, which might be because of the high dispersal capacity of this pine moth. Fragmentation resulted in the conversion of large and continuous habitats into isolated, small and insular patches, which was the primary effect on the genetic diversity of D. punctatus in TIL. The conclusion to emphasize from our research is that habitat fragmentation reduced the biological genetic diversity to some extent, further demonstrating the importance of habitat continuity in biodiversity protection.

Superior and inferior rice grains have different weights and are located on the upper primary branch and lower secondary branches of the panicle, respectively. To study differences in germination vigour of these two types of grain, a number of factors were investigated from 0 to 48 h of germination. The present study demonstrated that in inferior grains the starch granule structure was looser at 0 h, with full water absorption at 48 h, while in superior grains the structure was tight and dense. Relative water content increased, and dry matter decreased, more rapidly in inferior grains than in superior ones. Abscisic acid and gibberellin levels, as well as α-amylase activity, also changed more rapidly in inferior grains, while soluble sugar content and amylase coding gene expression increased more rapidly in inferior than superior grains during early germination. The expression of OsGAMYB was higher in inferior grains at 24 h but higher in superior grains at 48 h. The phenotypic index of seedlings was higher in seedlings from superior grains at the two-leaf stage. However, the thousand-grain weight and yield per plant in superior and inferior plants showed no significant difference at harvest. The present study indicates that inferior grains germinate faster than superior ones in the early germination stage. Although inferior grains produced weaker seedlings, it is worthwhile using them in rice production due to their comparative yield potential over that of superior grains.

Cucumber powdery mildew is a destructive foliar disease caused by Podosphaera xanthii (formerly known as Sphaerotheca fuliginea) that substantially damages the yield and quality of crops. The control of this disease primarily involves the use of chemical pesticides that cause serious environmental problems. Currently, numerous studies have indicated that some plant extracts or products potentially have the ability to act as natural pesticides to control plant diseases. It has been reported that turmeric (Curcuma longa L.) and its extract can be used in agriculture due to their insecticidal and fungicidal properties. However, the most effective fungicidal component of this plant is still unknown. In the current study, the crude extract of C. longa L. was found to have a fungicidal effect against P. xanthii. Afterwards, eight fractions (Fr.1–Fr.8) were gradually separated from the crude extract by column chromatography. Fraction 1 had the highest fungicidal effect against this pathogen among the eight fractions. The active compound, (+)-(S)-ar-turmerone, was separated from Fr 1 by semi-preparative high-performance liquid chromatography and identified based on its 1H nuclear magnetic resonance (NMR) and 13C NMR spectrum data. The EC50 value of (+)-(S)-ar-turmerone was found to be 28.7 µg/ml. The compound also proved to have a curative effect. This is the first study to report that the compound (+)-(S)-ar-turmerone has an effect on controlling this disease. These results provide a basis for developing a new phytochemical fungicide from C. longa L. extract.

In vivo and in vitro trials were conducted to assess the effects of tributyrin (TB) supplementation on short-chain fatty acid (SFCA) concentrations, fibrolytic enzyme activity, nutrient digestibility and methanogenesis in adult sheep. Nine 12-month-old ruminally cannulated Small Tail ewes (initial body weight 55 ± 5.0 kg) without pregnancy were used for the in vitro trial. In vitro substrate made to offer TB at 0, 2, 4, 6 and 8 g/kg on a dry matter (DM) basis was incubated by ruminal microbes for 72 h at 39°C. Forty-five adult Small Tail ewes used for the in vivo trial were randomly assigned to five treatments with nine animals each for an 18-d period according to body weight (55 ± 5.0 kg). Total mixed ration fed to ewes was also used to offer TB at 0, 2, 4, 6 and 8 g/kg on a DM basis. The in vitro trial showed that TB supplementation linearly increased apparent digestibility of DM, crude protein, neutral detergent fibre and acid detergent fibre, and enhanced gas production and methane emissions. The in vivo trial showed that TB supplementation decreased DM intake, but enhanced ruminal fermentation efficiency. Both in vitro and in vivo trials showed that TB supplementation enhanced total SFCA concentrations and carboxymethyl cellulase activity. The results indicate that TB supplementation might exert advantage effects on rumen microbial metabolism, despite having an enhancing effect on methanogenesis.

Litter size has a great impact on the profit of swine producers. Uterine development is an important determinant of reproduction efficiency and could hence affect litter size. Chinese Erhualian pig is one of the most prolific breeds in the world, even though large phenotypic variation in litter size was observed within Erhualian sows. To dissect the genetic basis of the phenotypic variation, we herein conducted genome-wide association studies for total number born and number born alive (NBA) of Erhualian sows. In total, one significant single nucleotide polymorphism (SNP) (P<1.78e−06) and 11 suggestive SNPs (P<3.57e−05) were identified on 10 chromosomes, confirming seven previously reported quantitative trait loci (QTL) and uncovering six QTL for litter size or uterus length. One locus on Sus scrofa chromosome (SSC) 13 (79.28 to 90.43 Mb) harbored a cluster of suggestive SNPs associated with multiparous NBA. The SNP (rs81447100) within this region was confirmed to be significantly (P<0.05) associated with litter size in Erhualian (n=313), Sutai (n=173) and Yorkshire (n=488) populations. Retinol binding protein 2 and retinol binding protein 1 functionally related to the development of uterus were located in a region of 2 Mb around rs81447100. Moreover, four genes related to embryo implantation and development were also detected around other significant SNPs. Taken together, our findings provide a potential marker (rs81447100) for the genetic improvement of litter size not only in Chinese Erhualian pigs but also in European commercial pig breeds like Yorkshire, and would facilitate the final identification of causative variant(s) underlying the effect of SSC13 QTL on litter size.

Nitrogen (N) application and irrigation to winter wheat may decrease leaf temperature and enhance photosynthesis: as a result, more photosynthates will be allocated to the grains, resulting in higher grain yields. To investigate this hypothesis, a 2-year field study was conducted with three levels of N fertilizer application (no fertilizer, N0; 240 kg N/ha, N1; 360 kg N/ha, N2) and two different water regimes (rainfed with no irrigation, R; irrigation at the over-wintering, stem elongation and grain filling stages, W). The results show that both N application and supplemental irrigation significantly increased grain yield with increases in both grain number/m2 and the 1000-grain weight, viz., WN2>WN1>WN0>RN2>RN1>RN0. In addition, application of N under both water regimes significantly increased flag leaf area, above-ground biomass and single stem productivity and decreased leaf temperature, which led to an increase in net photosynthesis rates and ribulose bisphosphate (RuBP) carboxylase activity. Moreover, analysis of the chlorophyll α fluorescence transient showed that N fertilizer application and supplemental irrigation significantly increased electron donor and acceptor performance of the photosystem II reaction centre.

l-Arginine is a nutritionally essential amino acid for spermatogenesis and plays versatile roles in animal health and can be utilized as a potential agent to improve reproductive performance of boars under high ambient temperature. The present study aimed to determine whether dietary l-arginine could alleviate heat stress-induced infertility in boars. In all, 20 boars (PIC 1040; 248.59±3.84 kg BW and 407.65±6.40 days of age) were selected and randomly assigned to four groups (group 0.0%, basal diet; group 0.6%, 0.8% or 1.0%, basal diet added with 0.6%, 0.8% or 1.0% l-arginine (wt:wt), respectively.) The four diets were made isonitrogenous by addition of appropriate amounts of l-alanine. Boars were pre-fed the corresponding experimental diet for 42 days. Then, the semen characteristics and libido were accessed for 6 weeks during the hot summer period (25.5° to 33.0°C). Results show that dietary l-arginine remarkably improved sperm motility, normality, total sperm number and effective total sperm number. Also, dietary l-arginine improved semen antioxidant capacity, such as decrease of malondialdehyde and 8-Hydroxy-2'-deoxyguanosine content in sperm (P<0.05), increase of the ratio of glutathione and oxidized glutathione, total antioxidant capacity, glutathione peroxidase and catalase activities in seminal plasma (P<0.05). Most of mitochondria contained intact ultrastructure in l-arginine-supplemented group which also accompany with higher ATP content than the 0.0% group. The boars fed 0.8% l-arginine show increased levels of estradiol-17β and testosterone and exhibit improved libido performance than boars in the 0.0% group. Adding dietary l-arginine linearly increased (P=0.002) nitric oxide content (as l-arginine increased). The scrotal surface temperature in the 0.6%, 0.8% and 1.0% group were decreased by 0.9°C, 0.9°C and 0.4°C, respectively, compared with the 0.0% group. l-Arginine levels caused linear effect on semen quality and antioxidant capacity, also caused quadratic effect on libido performance. During the hot summer months, the predicted optimal l-arginine levels for best semen quality and antioxidant capacity was 0.8% to 1.0% and for best libido performance was 0.8%. It can be concluded that l-arginine can be used as an effective agent to alleviate heat stress-induced infertility of boar, and that 0.8% to 1.0% can be considered as the optimum dosage.

Genomic selection has become increasingly important in the breeding of animals and plants. The response variable is an important factor, influencing the accuracy of genomic selection. The de-regressed proof (DRP) based on traditional estimated breeding value (EBV) is commonly used as response variable. In the current study, simulated data from 16th QTL-MAS Workshop and real data from Chinese Holstein cattle were used to compare accuracy and bias of genomic prediction with two methods of calculating DRP. Our results with simulated data showed that the correlation between genomic EBV and true breeding value achieved using the Jairath method (DRP_J) was superior to that achieved using the Garrick method (DRP_G) for simulated trait 1 but the reverse was true for simulated trait 3, and these two methods performed comparably for simulated trait 2. For all three simulated traits, DRP_J yielded larger bias of genomic prediction. However, DRP_J outperformed DRP_G in both accuracy and unbiasedness for four milk production traits in Chinese Holstein. In the estimation of genomic breeding value using genomic BLUP model, two methods for weighting diagonal elements of incidence matrix associated with residual error were also compared. With increasing the proportion of genetic variance unexplained by markers, the accuracy of genomic prediction was decreased and the bias was increased. Weighting by the reliability of DRP produced accuracy comparable to the evaluation where the proportion of genetic variance unexplained by markers was considered, but with smaller bias in general.

The origins and phylogeny of different sheep breeds has been widely studied using polymorphisms within the mitochondrial hypervariable region. However, little is known about the mitochondrial DNA (mtDNA) content and phylogeny based on mtDNA protein-coding genes. In this study, we assessed the phylogeny and copy number of the mtDNA in eight indigenous (population size, n=184) and three introduced (n=66) sheep breeds in China based on five mitochondrial coding genes (COX1, COX2, ATP8, ATP6 and COX3). The mean haplotype and nucleotide diversities were 0.944 and 0.00322, respectively. We identified a correlation between the lineages distribution and the genetic distance, whereby Valley-type Tibetan sheep had a closer genetic relationship with introduced breeds (Dorper, Poll Dorset and Suffolk) than with other indigenous breeds. Similarly, the Median-joining profile of haplotypes revealed the distribution of clusters according to genetic differences. Moreover, copy number analysis based on the five mitochondrial coding genes was affected by the genetic distance combining with genetic phylogeny; we also identified obvious non-synonymous mutations in ATP6 between the different levels of copy number expressions. These results imply that differences in mitogenomic compositions resulting from geographical separation lead to differences in mitochondrial function.