Genetic diversity of mtDNA D-loop sequences in four native Chinese chicken breeds

1. To explore the genetic diversity of Chinese indigenous chicken breeds, a 585 bp fragment of the mitochondrial DNA (mtDNA) region was sequenced in 102 birds from the Xichuan black-bone chicken, Yunyang black-bone chicken and Lushi chicken. In addition, 30 mtDNA D-loop sequences of Silkie fowls were downloaded from NCBI. The mtDNA D-loop sequence polymorphism and maternal origin of 4 chicken breeds were analysed in this study.

2. The results showed that a total of 33 mutation sites and 28 haplotypes were detected in the 4 chicken breeds. The haplotype diversity and nucleotide diversity of these 4 native breeds were 0.916 ± 0.014 and 0.012 ± 0.002, respectively. Three clusters were formed in 4 Chinese native chickens and 12 reference breeds. Both the Xichuan black-bone chicken and Yunyang black-bone chicken were grouped into one cluster. Four haplogroups (A, B, C and E) emerged in the median-joining network in these breeds.

3. It was concluded that these 4 Chinese chicken breeds had high genetic diversity. The phylogenetic tree and median network profiles showed that Chinese native chickens and its neighbouring countries had at least two maternal origins, one from Yunnan, China and another from Southeast Asia or its surrounding area.     

略阳乌鸡线粒体DNA控制区(D-loop)的遗传多样性分析

为了研究略阳乌鸡线粒体DNA控制区（mtDNA D-loop）的遗传多样性和起源,本研究对30只略阳乌鸡样品的mtDNA D-loop全序列进行了PCR扩增和测序,结合GenBank中公布的其他鸡的D-loop区序列,分析略阳乌鸡线粒体多态性及其起源。结果表明,略阳乌鸡mtDNA D-loop区全序列中,A、C、G、T平均含量分别为26.6%、26.6%、13.4%和33.4%,26个核苷酸多态位点均为转换位点,核苷酸多样度（Pi）为0.00705,单倍型变异度（Hd）为1.000,中性检验Tajima's D值为-0.47272。通过群体构建的系统进化树发现,略阳乌鸡样品在系统进化树上聚为4大分支。研究结果表明,略阳乌鸡群体内个体序列变异程度较大,遗传多样性丰富,揭示略阳乌鸡在遗传组成上具有4个母系来源。     

Genetic Diversity of Lueyang Black-bone Chicken Based on Mitochondrial DNA D-loop Region Sequence

This study was aimed to assess mitochondrial DNA D-loop sequence diversity and origin of Lueyang Black-bone chicken.The mtDNA D-loop sequence of 30 individuals from Lueyang Black-bone chicken were amplified by PCR and subsequently sequenced.The mtDNA D-loop sequence of other chicken were collected from GenBank and used as reference sequences to analyze the diversity and origin of Lueyang Black-bone chicken.The results revealed that the average values of base composition of A,C,G and T in mtDNA D-loop the sequence of Lueyang Black-bone chicken were 26.6%,26.6%,13.4% and 33.4%,respectively.26 nucleotide polymorphic sites were transition.The average nucleotide diversity (Pi) of the sites and haplotype diversity (Hd) were 0.00705 and 1.000,and the value of Tajima's D was -0.47272.Phylogenetic tree showed that samples were clusted in 4 clades. In the research, it could be concluded that the genetic diversity was relatively rich and wealthy and there were 4 maternal origins to Lueyang Black-bone chicken population.     

Genetic diversity of Bangladeshi native chickens based on complete sequence of mitochondrial DNA D-loop region

ABSTRACT

1. The aim of this study was to explore genetic diversity and possible origin of Bangladeshi (BD) native chickens. The complete mtDNA D-loop region was sequenced in 60 chickens representing five populations; naked neck, full feathered, Aseel, Hilly and autosomal dwarf. The 61 reference sequences representing different domestic chicken clades in China, India, Laos, Indonesia, Myanmar, and other Eurasian regions were included. The mtDNA D-loop sequence polymorphism and maternal origin of five BD populations were analysed.

2. A total of 35 polymorphic sites, and 21 haplotypes were detected in 60 mtDNA D-loop sequences. The haplotype and nucleotide diversity of the five populations were 0.921 ± 0.018 and 0.0061 ± 0.0019, respectively. Both mtDNA network and phylogenetic analysis indicated four clades (four haplogroups) in BD populations (21 haplotypes) along with 61 reference haplotypes. Clade E contained the most individuals (20) and haplotypes (11) of BD chickens, followed by clade D (17, 6), clade C (12, 2) and clade F (11, 2), respectively.

3. The higher number of unique haplotypes found in Yunnan, China, suggested that the origin of BD chickens was in this region. The haplotypes from different haplogroups were introduced in Bangladeshi chickens from India, China and Myanmar. The phylogenetic tree showed a close relationship of BD chickens with the clusters from India, China, Myanmar and Laos, and indicated the dispersion of BD chickens from these sources. The phylogenetic information revealed high genetic diversity of BD chickens because of their origin from different lineages with high genetic variation and distance, which was determined from four cluster and neighbour-joining trees.

4. In conclusion, BD populations had high genetic diversity. The mtDNA network profiles and phylogenetic trees showed multiple maternal origins of BD chickens from India, China, Myanmar and Laos.     

Genetic Diversity and Evolution of Danzhou Chicken Assessed with Mitochondrial Complete DNA D-loop Region Sequencing

The objective of this study was to determine the genetic diversity and evolution of Danzhou chicken.The complete mitochondrial DNA (mtDNA) D-loop regions of 36 Danzhou chickens were amplified,sequenced and analyzed.The sequencing reads were compared with the complete mtDNA D-loop sequence of several relative strains of chicken annotated in GenBank,and analyzed by bioinformatics methods.The genetic diversity and its evolutionary relationship in Danzhou chicken were analyzed.The results showed that the lengths of PCR products at the D-loop region were 1 210 bp,with 59.9% being A+T and 40.1% as C+G.The variable regions were 167-1 215 bp,and the high variable regions were mainly 167-367 bp.A total of 20 variable sites that defined 6 haplotypes were identified.The average haplotype diversity (Hd) and average number of nucleotide difference (k) were 0.571 and 6.449,respectively,the nucleotide diversity (Pi) was 0.00537,and the Tajima's D value of neutrality test was 1.61643.6 haplotypes could be grouped to 3 haplogroups (A,B and C) as determined by phylogenic analysis,with B clade,as the most abundant population.It concluded that the genetic diversity and haplotype diversity of Danzhou chicken were relatively low.Phylogenetic tree showed that the genetic composition of Danzhou chicken came from 3 maternal ancestors,Gallus gallus spadiceus,Gallus gallus bankiva and Gallus gallus jabouillei were potential ancestors.There was few influence of exotic lineage detected,which indicated that Danzhou chicken was a relatively conserved breed.     

应用微卫星标记对12个中国地方乌骨鸡品种遗传多样性的研究

利用27个微卫星DNA标记对12个地方乌骨鸡品种进行遗传多样性分析，通过计算等位基因频率、多态信息含量、杂合度、有效等位基因数和遗传距离，评估各品种内遗传变异和各品种间遗传相关，并根据遗传距离对12个鸡品种进行了聚类分析。结果表明，各乌骨鸡群体的遗传多样性较为丰富，并具有较高的选择潜力。平均杂合度最高的是竹乡鸡，为0．670；最低的是江山乌骨鸡，为0．581。27个微卫星标记中4个为中度多态座位，1个低度多态座位，22个为高度多态座位；LEI0234与LEI0192分别检测到了10．3与12．2个等位基因。以标准遗传距离为准，遗传距离最近的为沐川乌骨鸡与略阳鸡，为0．1002，而乌蒙乌骨鸡与略阳鸡最远，为0．2546。通过UPG—MA法聚类后，略阳鸡、沐川乌骨鸡、兴文乌骨鸡、盐津乌骨鸡、竹乡鸡首先聚为1类，江山乌骨鸡、郧阳白羽乌鸡、余干乌骨鸡聚为1类；金湖乌凤鸡、丝绒乌骨鸡和丝毛乌骨鸡聚为第3类；乌蒙乌骨鸡独为一类。     

海南白唇竹叶青线粒体基因组结构与特征分析

白唇竹叶青蛇(Trimeresurus albolabris)是亚洲竹叶青毒蛇之一,主要分布在东亚和东南亚。为研究海南岛产白唇竹叶青蛇线粒体基因,测定了采自海南文昌的个体线粒体基因组全序列,并与近缘物种进行比较。结果显示,其环状基因组全长17222 bp,包含13个蛋白质编码基因、22个转移RNA基因、2个核糖体RNA基因、2个控制区和1个L-链复制起源。总碱基组成为33.4%A,27.0%T,27.3%C,12.3%G,除8个tRNA基因和ND 6亚基基因在轻链上编码的外,其余基因均分布在重链上。此序列与其他产地的白唇竹叶青蛇相似度为99.61%,与其他竹叶青蛇相似度大于88%,与原矛头蝮蛇相似度大于85%。     

A Genetic Diversity Analysis of Mitochondrial DNA D-loop Region in Four High Quality Chicken Breeds

This study was conducted to elucidate the genetic diversity of mitochondrial DNA (mtDNA) D-loop region in Qingyuan partridge chicken group 1,Qingyuan partridge chicken group 2,Yangshan chicken and Qingyuan Yellow feather black-bone chicken.The specific primers were designed according to mtDNA D-loop region of Gullus gullus spadiceus (accession No.:NC_007235.1) in GenBank.The sequence was analyzed after PCR amplification and sequencing,and the haplotype number,polymorphism number,haplotype diversity,nucleotide diversity and nucleotide mean difference were counted.The evolution divergence among breeds was calculated by Mega 5.10 software,and the phylogenetic tree was constructed.The results showed that the length of mtDNA D-loop region in four high quality chicken breeds was 591 bp,and 549 bp were used for subsequent analysis.The content of A,T,C and G were 27.2% to 27.3%,30.1% to 30.4%,29.5% to 29.8% and 12.8% to 12.9%,respectively,and the average content of G+C was 42.5%.There were 92 polymorphic sites which contained 14 singleton variable sites and 78 parsimony informative sites,and the percentage of transitions and transversions were 89.13% (82/92) and 10.87% (10/92),respectively.The haplotype diversity ranged from 0.682 to 0.835,and the nucleotide diversity ranged from 0.00849 to 0.01167.There were 32 haplotypes in all sequences,which could be divided into clades A,B,C and E,however,most of the individuals belonged to clades B (51.2%) and E (37.6%).The phylogenetic tree results showed that four high quality chicken breeds could be classified as 4 branches which were consistent with the haplotypes classification results.The results indicated that the four high quality chicken populations from Qingyuan had relatively high haplotype and nucleotide diversity and likely shared two or more common maternal lineages.     

基于线粒体控制区的鸡不同杂交组合遗传多样性研究

旨在探讨鸡不同杂交组合线粒体控制区（mtDNA D-loop区）的遗传多样性和单倍型特性。选取固始鸡和隐性白羽鸡及其正、反交F1代、藏鸡以及F2代等6个群体共387个个体的mtDNA D-loop区进行测序,分析其遗传规律和单倍型特性,并与不同红色原鸡亚种进行聚类,分析其母系起源。结果显示,6个群体D-loop区全序列大小为1 231 bp,共检测到28个多态位点和1个C碱基缺失,共构成19种单倍型,分为A、B、C和E 4个单倍型群,其中,固始鸡和反交F1代主要为A、C单倍型,固始鸡A、C单倍型比例分别为53.42%和46.58%,反交F1代A、C单倍型比例分别为50.75%和49.25%;隐性白羽鸡、正交F1代和F2代优势单倍型均为E单倍型,占比分别为48.89%、48.84%和50.00%。6个鸡群体单倍型多样度（Hd）在0.496~0.729之间,核苷酸多样度（Pi）在0.003 40~0.005 41之间,Hd值和Pi值最大的均为正交F1代,其次为隐性白羽鸡和F2代,固始鸡和反交F1代群体遗传多样性接近。聚类分析显示,A、B单倍型群与滇南亚种交叉聚为一枝;E单倍型群与印度亚种交叉聚为一枝;C单倍型群与印度亚种、指名亚种、印尼亚种以及滇南亚种聚为一枝。结果提示,mtDNA D-loop区遵循严格的母系遗传,后代的遗传多样性和单倍型比例与其母本基本一致;我国家鸡群体具有多个红色原鸡母系起源,且主要起源于原鸡滇南亚种。     

Identification of key microRNAs affecting melanogenesis of breast muscle in Muchuan black-boned chickens by RNA sequencing

ABSTRACT

1. Melanin content is considered an important indicator of meat quality in black-boned chickens, which have a high market value. To understand the complex physiological processes underlying muscle melanogenesis in this chicken, differentially expressed miRNAs (DEMs) were detected between black muscle (BM) and white muscle (WM) of chickens using high-throughput sequencing technology. Six small RNA libraries were constructed, and more than 16.75 million clean reads were obtained for each library.

2. A total of 582 known miRNAs and 65 novel miRNAs were identified from the six chicken sequence libraries. A total of 19 DEMs were identified between the two groups, of which nine were upregulated and 10 were downregulated. Furthermore, the DEMs were predicted to target 572 genes.

3. Certain DEMs (such as miR-204, miR-133b, and miR-12 229-3p) and their target genes may play an important role in muscle melanogenesis of chickens. These findings provide a foundation for clarifying the miRNA regulatory mechanisms involved in muscle pigmentation in avian species.     

Genome‐wide re‐sequencing and transcriptome analysis reveal candidate genes associated with the pendulous comb phenotype in domestic chickens

To determine the causative variations associated with two chicken comb phenotypes, pendulous comb (PC) or upright comb (UC), two pooled genomic DNA samples from PC and UC chickens were re‐sequenced by Next‐Generation Sequencer, and genome‐wide Single nucleotide polymorphisms (SNPs) were detected. Using three selective sweep approaches, F_ST, θπ, and Tajima's D, with top 5% window values serving as the threshold, a total of 84 positively selective genes (PSGs) were identified. There were no SNPs in exons of the PSGs with significant differences in allele frequencies between the two comb phenotype groups. Then, 515 differentially expressed genes (DEGs) between the PC and UC were identified by RNA‐seq. Three genes including CD36 (CD36 molecule), ADAMTSL3 (ADAMTS‐like 3), and AOX1 (aldehyde oxidases 1) are overlapped between PSGs and DEGs. After genotyping seven candidate SNPs in the regulatory regions of the three overlapping genes in 120 chickens from two other breeds, two variants (rs14607046 and rs731818051) in the regulatory regions of AOX1 and ADAMTSL3 were found to have significant differences in allele frequency between the PC and UC, suggesting that the two variants may be causative mutations for PC. Overall, our study shed light on the genetic basis underlying the PC phenotype in chickens.     

4种鼢鼠线粒体DNA D-loop区全序列分析

采用PCR测序技术对采自甘肃境内的4种鼢鼠22个个体的线粒体DAND-loop区全序列进行了测定。结果表明：鼢鼠线粒体DNAD-loop序列长度为893、894、895、896或899bp,核苷酸位点突变类型有5种,即转换、颠换、插入、缺失及转换与颠换共存。碱基转换以T〈=〉C形式为主。其中T、C、A、G4种核苷酸的平均比例分别为34．1％、24．3％、30．1％、11．5％,A＋T含量（64．2％）明显高于G＋C含量（35．8％）。这4种鼢鼠22条线粒体DNAD-loop区发现20种单倍型,单倍型比例为81．82％,说明中国鼢鼠mtDNA遗传多态性很丰富。从系统发育树分析,4种鼢鼠明显聚为两类,推测鼢鼠可能有两个起源。     

A serological and virological survey for evidence of infection with Newcastle disease virus in Australian chicken farms

OBJECTIVE: To determine the prevalence and distribution of antibodies to Newcastle disease virus on Australian chicken farms and to determine the pathotype and relationships of the Newcastle disease viruses present on those farms. DESIGN: A cross-sectional survey of 753 commercial chicken farms. PROCEDURE: The survey comprised a detailed questionnaire and collection of venous blood samples. The titre of antibodies to Newcastle disease virus was determined by haemagglutination inhibition. Virus isolation was conducted from cloacal and tracheal swabs taken from chickens in serologically positive flocks. Virus isolates were pathotyped on the basis of the deduced Fusion protein cleavage site determined by nucleotide sequencing of a 265 bp region of the genome in the region of the cleavage site. RESULTS: Antibody evidence of Newcastle disease virus infection was found on 300 of the 753 surveyed farms throughout all 11 geographic regions of the survey. The highest prevalence occurred in the Sydney basin, New South Wales and Victoria east regions. Antibody titres were also highest in the regions where serologically positive flocks were most prevalent. The 259 virus isolates revealed nine different RNA sequences. Of the nine virus groups isolated, the most common group W was identical in sequence to the V4 vaccine strain. Five of the other groups had novel RNA sequences in the region of the F protein cleavage site. CONCLUSIONS: Antibodies to Newcastle disease virus are highly prevalent in the Australian chicken flock but all identified strains were avirulent in nature.     

不同类型肉鸡线粒体单倍型与遗传起源研究

旨在探讨鸡线粒体D-loop区单倍型特性以及与生长速度的相关性,并分析其遗传起源。选取不同类型肉鸡品种（配套系）6个,测定其生产性能,并对6个群体共计314个个体线粒体D-loop区全长进行测序,分析其单倍型特性;同时与不同红色原鸡亚种进行聚类,分析其母系起源。结果显示,6个群体D-loop区全序列共检测到37个突变位点,构成40种单倍型,分为A、B、C和E 4个单倍型群,其中AA肉鸡、罗斯308、禽雁麻鸡和裕禾1号肉鸡主要为E单倍型,占比分别为85.92%、50.00%、100.00%和70.21%;园丰麻鸡2号和港丰瑶黑麻鸡E单倍型占比相对较低,分别为15.00%和22.39%,园丰麻鸡2号B单倍型含量最高,占比66.67%,港丰瑶黑麻鸡C单倍型含量最高,占比40.30%。相关性分析显示,初生重与E单倍型比例之间呈极显著正相关（P<0.01）,与A、B和C单倍型比例之间均呈负相关;E单倍型比例与公母平均体重约1.8 kg时日龄之间呈显著负相关（P<0.05）,而与公母平均体重约1.8 kg时饲料转化比之间呈极显著负相关（P<0.01）。聚类分析显示,A、B单倍型所有个体均与原鸡滇南亚种聚为一枝;E单倍型所有个体均与原鸡印度亚种聚为一枝;C单倍型所有个体与原鸡印度亚种、滇南亚种、指名亚种以及印尼亚种交叉聚为一类。结果提示,线粒体单倍型与肉鸡生长速度之间相关显著,E单倍型与肉鸡生长速度具有较强的正相关;我国家鸡群体母系起源丰富,但主要起源于原鸡滇南亚种。     

Two fiber genes of nearly equal lengths are a common and distinctive feature of Fowl adenovirus C members

The complete genome or the genome region containing the two fiber genes of two reference strains and one field isolate representing both serotypes of Fowl adenovirus C were sequenced. Two fiber genes were revealed in the genomes of all three isolates. Fiber-1 and fiber-2 genes of several Fowl adenovirus C isolates were sequenced as well. Both serotypes 4 and 10 have two fiber genes. The genome region containing the fiber gene was also sequenced for the reference strain of Fowl adenovirus B. Just one fiber gene was revealed in this strain. Predicted amino acid sequences were compared to already published fiber sequences of different adenovirus isolates and one amino acid substitution within fiber-2 was detected in all of the Fowl adenovirus C isolates that were isolated from chickens with hepatitis-hydropericardium syndrome in comparison to apathogenic isolates. Phylogenetic analyses provided insights about the evolution of fiber genes in avian adenoviruses and their genetic relationships.     

基于线粒体DNA控制区(mtDNA D-loop)序列分析贵妃鸡的遗传多样性

采用PCR产物直接测序技术对30只贵妃鸡样品的线粒体DNA控制区(mtDNA D-loop)第Ⅰ高变区序列进行了分析。结果表明,在所分析的D-loop区部分序列中(520 bp),A、G、C、T平均含量分别为26.8%、13.0%、31.1%和29.1%;共发现13个核苷酸多态位点,均为转换位点,未检测到插入/缺失和颠换,核苷酸多样度(Pi)为0.0059,单倍型变异度(Hd)为0.538,中性检验Tajima’s D值为-0.67065。通过群体构建的NJ聚类图分子系统树发现,贵妃鸡起源于红原鸡。     

Phylogenetic analysis of Myanmar indigenous chickens using mitochondrial D-loop sequence reveals their characteristics as a genetic resource

Myanmar indigenous chickens play important roles in food, entertainment, and farm business for the people of Myanmar. In this study, complete mitochondrial D-loop sequences (1232 bp) were analyzed using 176 chickens, including three indigenous breeds, two fighting cock populations, and three indigenous populations to elucidate genetic diversity and accomplish a phylogenetic analysis of Myanmar indigenous chickens. The average haplotype and nucleotide diversities were 0.948 ± 0.009 and 0.00814 ± 0.00024, respectively, exhibiting high genetic diversity of Myanmar indigenous chickens. Sixty-four haplotypes were classified as seven haplogroups, with the majority being haplogroup F. The breeds and populations except Inbinwa had multiple maternal haplogroups, suggesting that they experienced no recent purifying selection and bottleneck events. All breeds and populations examined shared haplogroup F. When 232 sequences belonging to haplogroup F (79 from Myanmar and 153 deposited sequences from other Asian countries/region) were analyzed together, the highest genetic diversity was observed in Myanmar indigenous chickens. Furthermore, Myanmar indigenous chickens and red junglefowls were observed in the center of the star-like median-joining network of 37 F-haplotypes, suggesting that Myanmar is one of the origins of haplogroup F. These findings revealed the unique genetic characteristic of Myanmar indigenous chickens as important genetic resources.     

Use of high‐density SNP data to identify patterns of diversity and signatures of selection in broiler chickens

The development of broiler chickens over the last 70 years has been accompanied by large phenotypic changes, so that the resulting genomic signatures of selection should be detectable by current statistical techniques with sufficiently dense genetic markers. Using two approaches, this study analysed high‐density SNP data from a broiler chicken line to detect low‐diversity genomic regions characteristic of past selection. Seven regions with zero diversity were identified across the genome. Most of these were very small and did not contain many genes. In addition, fifteen regions were identified with diversity increasing asymptotically from a low level. These regions were larger and thus generally included more genes. Several candidate genes for broiler traits were found within these ‘regression regions’, including IGF1, GPD2 and MTNR1AI. The results suggest that the identification of zero‐diversity regions is too restrictive for characterizing regions under selection, but that regions showing patterns of diversity along the chromosome that are consistent with selective sweeps contain a number of genes that are functional candidates for involvement in broiler development. Many regions identified in this study overlap or are close to regions identified in layer chicken populations, possibly due to their shared precommercialization history or to shared selection pressures between broilers and layers.     

Complete genome sequence analysis of chicken astrovirus isolate from India

Objective

Chicken astroviruses have been known to cause severe disease in chickens leading to increased mortality and “white chicks” condition. Here we aim to characterize the causative agent of visceral gout suspected for astrovirus infection in broiler breeder chickens.

Methods

Total RNA isolated from allantoic fluid of SPF embryo passaged with infected chicken sample was sequenced by whole genome shotgun sequencing using ion-torrent PGM platform. The sequence was analysed for the presence of coding and non-coding features, its similarity with reported isolates and epitope analysis of capsid structural protein.

Results

The consensus length of 7513 bp genome sequence of Indian isolate of chicken astrovirus was obtained after assembly of 14,121 high quality reads. The genome was comprised of 13 bp 5′-UTR, three open reading frames (ORFs) including ORF1a encoding serine protease, ORF1b encoding RNA dependent RNA polymerase (RdRp) and ORF2 encoding capsid protein, and 298 bp of 3′-UTR which harboured two corona virus stem loop II like “s2m” motifs and a poly A stretch of 19 nucleotides. The genetic analysis of CAstV/INDIA/ANAND/2016 suggested highest sequence similarity of 86.94% with the chicken astrovirus isolate CAstV/GA2011 followed by 84.76% with CAstV/4175 and 74.48%% with CAstV/Poland/G059/2014 isolates. The capsid structural protein of CAstV/INDIA/ANAND/2016 showed 84.67% similarity with chicken astrovirus isolate CAstV/GA2011, 81.06% with CAstV/4175 and 41.18% with CAstV/Poland/G059/2014 isolates. However, the capsid protein sequence showed high degree of sequence identity at nucleotide level (98.64-99.32%) and at amino acids level (97.74–98.69%) with reported sequences of Indian isolates suggesting their common origin and limited sequence divergence. The epitope analysis by SVMTriP identified two unique epitopes in our isolate, seven shared epitopes among Indian isolates and two shared epitopes among all isolates except Poland isolate which carried all distinct epitopes.