Single and longitudinal genome-wide association studies for dairy traits available in goats with three recorded lactations.

Milk yield and composition phenotypes are systematically recorded across several lactations in goats, but the majority of genome-wide association studies (GWAS) performed so far have rather ignored the longitudinal nature of such data. Here, we have used two different GWAS approaches to analyse data from three lactations recorded in Murciano-Granadina goats. In Analysis 1, independent GWAS have been carried out for each trait and lactation, while a single longitudinal GWAS, jointly considering all data, has been performed in Analysis 2. In both analyses, genome-wide significant QTL for lactose percentage on chromosome 2 (129.77-131.01 Mb) and for milk protein percentage on the chromosome 6 (74.8-94.6 Mb) casein gene cluster region were detected. In Analysis 1, several QTL were not replicated in all three lactations, possibly due to the existence of lactation-specific genetic determinants. In Analysis 2, we identified several genome-wide significant QTL related to milk yield and protein content that were not uncovered in Analysis 1. The increased number of QTL identified in Analysis 2 suggests that the longitudinal GWAS is particularly well suited for the genetic analysis of dairy traits. Moreover, our data confirm that variability within or close to the casein complex is the main genetic determinant of milk protein percentage in Murciano-Granadina goats.

Markers with low GenTrain scores can generate spurious signals in genome-wide scans for transmission ratio distortion.

Transmission ratio distortion (TRD) is the preferential transmission of one specific allele to offspring at the expense of the other. The existence of TRD is mostly explained by the segregation of genetic variants with deleterious effects on the developmental processes that go from the formation of gametes to fecundation and birth. A few years ago, a statistical methodology was implemented in order to detect TRD signals on a genome-wide scale as a first step toward uncovering the biological basis of TRD and reproductive success in domestic species. In the current work, we have analyzed the impact of SNP calling quality on the detection of TRD signals in a population of Murciano-Granadina goats. Seventeen bucks and their offspring (N = 288) were typed with the Goat SNP50 BeadChip, whereas the genotypes of the dams were lacking. Performance of a genome-wide scan revealed the existence of 36 SNPs showing significant evidence of TRD. When we calculated GenTrain scores for each of the SNPs, we observed that 25 SNPs showed scores below 0.8. The allele frequencies of these SNPs in the offspring were not correlated with the allele frequencies estimated in the dams with statistical methods, providing evidence that flawed SNP calling quality might lead to the detection of spurious TRD signals. We conclude that, when performing TRD scans, the GenTrain scores of markers should be taken into account to discriminate SNPs that are truly under TRD from those yielding spurious signals owing to technical problems.

A genome-wide analysis of copy number variation in Murciano-Granadina goats.

BACKGROUND: In this work, our aim was to generate a map of the copy number variations (CNV) segregating in a population of Murciano-Granadina goats, the most important dairy breed in Spain, and to ascertain the main biological functions of the genes that map to copy number variable regions. RESULTS: Using a dataset that comprised 1036 Murciano-Granadina goats genotyped with the Goat SNP50 BeadChip, we were able to detect 4617 and 7750 autosomal CNV with the PennCNV and QuantiSNP software, respectively. By applying the EnsembleCNV algorithm, these CNV were assembled into 1461 CNV regions (CNVR), of which 486 (33.3% of the total CNVR count) were consistently called by PennCNV and QuantiSNP and used in subsequent analyses. In this set of 486 CNVR, we identified 78 gain, 353 loss and 55 gain/loss events. The total length of all the CNVR (95.69 Mb) represented 3.9% of the goat autosomal genome (2466.19 Mb), whereas their size ranged from 2.0 kb to 11.1 Mb, with an average size of 196.89 kb. Functional annotation of the genes that overlapped with the CNVR revealed an enrichment of pathways related with olfactory transduction (fold-enrichment = 2.33, q-value = 1.61 × 10-10), ABC transporters (fold-enrichment = 5.27, q-value = 4.27 × 10-04) and bile secretion (fold-enrichment = 3.90, q-value = 5.70 × 10-03). CONCLUSIONS: A previous study reported that the average number of CNVR per goat breed was ~ 20 (978 CNVR/50 breeds), which is much smaller than the number we found here (486 CNVR). We attribute this difference to the fact that the previous study included multiple caprine breeds that were represented by small to moderate numbers of individuals. Given the low frequencies of CNV (in our study, the average frequency of CNV is 1.44%), such a design would probably underestimate the levels of the diversity of CNV at the within-breed level. We also observed that functions related with sensory perception, metabolism and embryo development are overrepresented in the set of genes that overlapped with CNV, and that these loci often belong to large multigene families with tens, hundreds or thousands of paralogous members, a feature that could favor the occurrence of duplications or deletions by non-allelic homologous recombination.

A genome-wide association analysis for body, udder, and leg conformation traits recorded in Murciano-Granadina goats.

Morphological traits are of great importance to dairy goat production given their effect on phenotypes of economic interest. However, their underlying genomic architecture has not yet been extensively characterized. Herein, we aimed to identify genomic regions associated with body, udder, and leg conformation traits recorded in 825 Murciano-Granadina goats. We genotyped this resource population using the GoatSNP50 BeadChip (Illumina Inc., San Diego, CA) and performed genome-wide association analyses using the GEMMA software. We found 2 genome-wide significant associations between markers rs268273468 [Capra hircus (CHI) 16:69617700] and rs268249346 (CHI 28:18321523) and medial suspensory ligament. In contrast, we did not detect any genome-wide significant associations for body and leg traits. Moreover, we found 12, 19, and 7 chromosome-wide significant associations for udder, body, and leg traits, respectively. Comparison of our data with previous studies revealed a low level of positional concordance between regions associated with morphological traits. In addition to technical factors, this lack of concordance could be due to a substantial level of genetic heterogeneity among breeds or to the strong polygenic background of morphological traits, which makes it difficult to detect genetic factors that have small phenotypic effects.

Patterns of homozygosity in insular and continental goat breeds.

BACKGROUND: Genetic isolation of breeds may result in a significant loss of diversity and have consequences on health and performance. In this study, we examined the effect of geographic isolation on caprine genetic diversity patterns by genotyping 480 individuals from 25 European and African breeds with the Goat SNP50 BeadChip and comparing patterns of homozygosity of insular and nearby continental breeds. RESULTS: Among the breeds analysed, number and total length of ROH varied considerably and depending on breeds, ROH could cover a substantial fraction of the genome (up to 1.6 Gb in Icelandic goats). When compared with their continental counterparts, goats from Iceland, Madagascar, La Palma and Ireland (Bilberry and Arran) displayed a significant increase in ROH coverage, ROH number and FROH values (P value < 0.05). Goats from Mediterranean islands represent a more complex case because certain populations displayed a significantly increased level of homozygosity (e.g. Girgentana) and others did not (e.g. Corse and Sarda). Correlations of number and total length of ROH for insular goat populations with the distance between islands and the nearest continental locations revealed an effect of extremely long distances on the patterns of homozygosity. CONCLUSIONS: These results indicate that the effects of insularization on the patterns of homozygosity are variable. Goats raised in Madagascar, Iceland, Ireland (Bilberry and Arran) and La Palma, show high levels of homozygosity, whereas those bred in Mediterranean islands display patterns of homozygosity that are similar to those found in continental populations. These results indicate that the diversity of insular goat populations is modulated by multiple factors such as geographic distribution, population size, demographic history, trading and breed management.

A genome-wide perspective about the diversity and demographic history of seven Spanish goat breeds.

BACKGROUND: The main goal of the current work was to infer the demographic history of seven Spanish goat breeds (Malagueña, Murciano-Granadina, Florida, Palmera, Mallorquina, Bermeya and Blanca de Rasquera) based on genome-wide diversity data generated with the Illumina Goat SNP50 BeadChip (population size, N = 176). Five additional populations from Europe (Saanen and Carpathian) and Africa (Tunisian, Djallonké and Sahel) were also included in this analysis (N = 80) for comparative purposes. RESULTS: Our results show that the genetic background of Spanish goats traces back mainly to European breeds although signs of North African admixture were detected in two Andalusian breeds (Malagueña and Murciano-Granadina). In general, observed and expected heterozygosities were quite similar across the seven Spanish goat breeds under analysis irrespective of their population size and conservation status. For the Mallorquina and Blanca de Rasquera breeds, which have suffered strong population declines during the past decades, we observed increased frequencies of large-sized (ROH), a finding that is consistent with recent inbreeding. In contrast, a substantial part of the genome of the Palmera goat breed comprised short ROH, which suggests a strong and ancient founder effect. CONCLUSIONS: Admixture with African goats, genetic drift and inbreeding have had different effects across the seven Spanish goat breeds analysed in the current work. This has generated distinct patterns of genome-wide diversity that provide new clues about the demographic history of these populations.

Genomic patterns of homozygosity and inbreeding depression in Murciano-Granadina goats.

BACKGROUND: Inbreeding depression can adversely affect traits related to fitness, reproduction and productive performance. Although current research suggests that inbreeding levels are generally low in most goat breeds, the impact of inbreeding depression on phenotypes of economic interest has only been investigated in a few studies based on genealogical data. RESULTS: We genotyped 1040 goats with the Goat SNP50 BeadChip. This information was used to estimate different molecular inbreeding coefficients and characterise runs of homozygosity and homozygosity patterns. We detected 38 genomic regions with increased homozygosity as well as 8 ROH hotspots mapping to chromosomes 1, 2, 4, 6, 14, 16 and 17. Eight hundred seventeen goats with available records for dairy traits were analysed to evaluate the potential consequences of inbreeding depression on milk phenotypes. Four regions on chromosomes 8 and 25 were significantly associated with inbreeding depression for the natural logarithm of the somatic cell count. Notably, these regions contain several genes related with immunity, such as SYK, IL27, CCL19 and CCL21. Moreover, one region on chromosome 2 was significantly associated with inbreeding depression for milk yield. CONCLUSIONS: Although genomic inbreeding levels are low in Murciano-Granadina goats, significant evidence of inbreeding depression for the logarithm of the somatic cell count, a phenotype closely associated with udder health and milk yield, have been detected in this population. Minimising inbreeding would be expected to augment economic gain by increasing milk yield and reducing the incidence of mastitis, which is one of the main causes of dairy goat culling.

Effects of α s1-casein (CSN1S1) and κ-casein (CSN3) genotypes on milk coagulation properties in Murciano-Granadina goats.

The effects of the caprine α s1-casein (CSN1S1) polymorphisms on milk quality and cheese yield have been widely studied in French and Italian goat breeds. Much less is known about the consequences of κ-casein (CSN3) genotype on the technological and coagulation properties of goat milk. In the current study, we have performed an association analysis between polymorphisms at the goat CSN1S1 and CSN3 genes and milk coagulation (rennet coagulation time, curdling rate and curd firmness) and technological (time to cutting of curd and cheese yield) properties. In this analysis, we have included 193 records from 74 Murciano-Granadina goats (with genotypes constituted by different combinations of alleles B, E and F of the gene CSN1S1 and alleles A and B of the gene CSN3) distributed in three herds, which were collected bimonthly during a whole lactation. Data analysis, using a linear mixed model for repeated observations, revealed significant associations between CSN1S1 genotypes and the rate of the curdling process. In this way, milk from EE goats had a significantly higher curdling rate than milk from BB individuals (P<0·05). Contrary to previous experiments performed in French breeds, cheese yield was not significantly different in BB, EE and EF goats. Moreover, we have shown that CSN3 genotype has a significant effect on the rennet coagulation time (BB>AB, P<0·05) but not on cheese yield. No interaction between the CSN1S1 and CSN3 genotypes was observed.

Genetic variation at the goat hormone-sensitive lipase ( LIPE) gene and its association with milk yield and composition.

Hormone-sensitive lipase (LIPE) plays a fundamental role in the regulation of energy balance by releasing free fatty acids from adipose triacylglycerol stores. These fatty acids can be subsequently transferred to other body compartments to be oxidized or employed in other biochemical reactions. This enzymic function is particularly important in lactating animals because the synthesis of milk components involves the mobilization of lipid depots to satisfy the large energy demands of the mammary gland. In the current study, we partially sequenced the goat LIPE gene in several individuals. In doing so, we identified two synonymous polymorphisms at exons 2 (c.327C>A>T, triallelic polymorphism) and 3 (c.558C>T). Moreover, we found a mis-sense polymorphism at exon 6 (c.1162G>T) that involves an alanine to serine substitution at position 388. Analysis with Polyphen and Panther softwares revealed that this amino acid replacement is expected to be neutral. Performance of an association analysis with a variety of milk traits revealed that goat LIPE genotype has highly suggestive effects on milk yield (P=0.0032) as well as on C18:3 n-6g (P=0.0051), trans-10 cis-12 CLA (P=0.007) and C12:0 (P=0.0084) milk contents. These associations are concordant with the preference of LIPE to selectively mobilize medium-chain and unsaturated fatty acids.

Lineal Discrimination of Horses and Mules. A Sympatric Case from Arauca, Colombia.

This is the first morphological comparative study between local horses and mules from Arauca, Colombia. It was realized to compare morphological traits between both species by analysing 15 adult mules (7 males and 8 females) and 150 adult horses (137 males and 13 females), with an age interval from 2 to 22 years. Data consisted of 24 different body quantitative traits which can explain the body conformation: thoracic circumference, body length (BL), thoracic depth and width, withers height (WH), sternum height, shoulders width, chest width, forelimb cannon perimeter and length, head length and width, skull length and width, face length and width, ear length and width, loin height, croup height (CrH), width and length, dock height (DoH), and hock height. Heart girth circumference, body length, withers height, croup height, and dock height were the most discriminative traits, showing statistical differences between species. The formula is X = (BL × 0.402) + (WH × 0.323) + (CrH × 0.352) + (DoH × 0.384). A value of X > 184.5 assigns with total certainty that a skeleton belongs to a horse, and if X < 174.0, it is a mule. The proposed formula has a 100% specificity but a 71.4% sensibility for mules and an 84.4% for horses into the rank of 174.0-184.5. Therefore, results demonstrate that some postcranial anatomical elements of Equus could give enough information for a bone differentiation between horses and mules, at least in animals from the Araucan region, but the main interest is that it reflects the possibility to differentiate morphometrically both species from bone remains when horses and mules were sympatric.

Analyzing the genomic and transcriptomic architecture of milk traits in Murciano-Granadina goats.

BACKGROUND: In this study, we aimed to investigate the molecular basis of lactation as well as to identify the genetic factors that influence milk yield and composition in goats. To achieve these two goals, we have analyzed how the mRNA profile of the mammary gland changes in seven Murciano-Granadina goats at each of three different time points, i.e. 78 d (T1, early lactation), 216 d (T2, late lactation) and 285 d (T3, dry period) after parturition. Moreover, we have performed a genome-wide association study (GWAS) for seven dairy traits recorded in the 1st lactation of 822 Murciano-Granadina goats. RESULTS: The expression profiles of the mammary gland in the early (T1) and late (T2) lactation were quite similar (42 differentially expressed genes), while strong transcriptomic differences (more than one thousand differentially expressed genes) were observed between the lactating (T1/T2) and non-lactating (T3) mammary glands. A large number of differentially expressed genes were involved in pathways related with the biosynthesis of amino acids, cholesterol, triglycerides and steroids as well as with glycerophospholipid metabolism, adipocytokine signaling, lipid binding, regulation of ion transmembrane transport, calcium ion binding, metalloendopeptidase activity and complement and coagulation cascades. With regard to the second goal of the study, the performance of the GWAS allowed us to detect 24 quantitative trait loci (QTLs), including three genome-wide significant associations: QTL1 (chromosome 2, 130.72-131.01 Mb) for lactose percentage, QTL6 (chromosome 6, 78.90-93.48 Mb) for protein percentage and QTL17 (chromosome 17, 11.20 Mb) for both protein and dry matter percentages. Interestingly, QTL6 shows positional coincidence with the casein genes, which encode 80% of milk proteins. CONCLUSIONS: The abrogation of lactation involves dramatic changes in the expression of genes participating in a broad array of physiological processes such as protein, lipid and carbohydrate metabolism, calcium homeostasis, cell death and tissue remodeling, as well as immunity. We also conclude that genetic variation at the casein genes has a major impact on the milk protein content of Murciano-Granadina goats.

Donkey genomes provide new insights into domestication and selection for coat color.

Current knowledge about the evolutionary history of donkeys is still incomplete due to the lack of archeological and whole-genome diversity data. To fill this gap, we have de novo assembled a chromosome-level reference genome of one male Dezhou donkey and analyzed the genomes of 126 domestic donkeys and seven wild asses. Population genomics analyses indicate that donkeys were domesticated in Africa and conclusively show reduced levels of Y chromosome variability and discordant paternal and maternal histories, possibly reflecting the consequences of reproductive management. We also investigate the genetic basis of coat color. While wild asses show diluted gray pigmentation (Dun phenotype), domestic donkeys display non-diluted black or chestnut coat colors (non-Dun) that were probably established during domestication. Here, we show that the non-Dun phenotype is caused by a 1 bp deletion downstream of the TBX3 gene, which decreases the expression of this gene and its inhibitory effect on pigment deposition.

Analysing the Expression of Eight Clock Genes in Five Tissues From Fasting and Fed Sows.

In a previous study, we observed that circadian clock genes are differentially expressed in the skeletal muscle of fasting and fed sows. The goal of the current work was to investigate if these genes are also differentially expressed in tissues containing the central (hypothalamus) and peripheral (duodenum, dorsal fat, muscle, and liver) clocks. As animal material, we used 12 sows that fasted 12 h before slaughtering (T0) and 12 sows that were fed ad libitum 7 h prior slaughtering (T2). Tissue samples were collected immediately after slaughter and total RNA was subsequently extracted. The expression of the ARNTL, BHLHE40, CRY2, NPAS2, NR1D1, PER1, PER2, and SIK1 genes was measured by quantitative reverse transcription PCR. The numbers of clock genes showing differential expression before and after feeding varied depending on the tissue i.e., four in dorsal fat and duodenum, six in skeletal muscle, and seven in the liver. In contrast, none of the eight analysed genes displayed a significant differential expression in hypothalamus, the tissue where the central clock resides. This result supports that the differential expression of clock genes in the four tissues mentioned above is probably induced by nutrition and not by the central clock entrained by light. Moreover, we have observed that the NPAS2 and ARNTL genes display positive log2(FC) values in the five tissues under analysis, whilst the CRY2, PER1 (except dorsal fat) and PER2 (except hypothalamus) genes generally show negative log2(FC) values. Such result might be explained by the existence of a negative feedback loop between the ARNTL/NPAS2 and CRY/PER genes. Collectively, these results support that nutrition plays an important role in modulating the timing of porcine peripheral circadian clocks. Such regulation could be essential for coordinating the subsequent metabolic response to nutrient supply.

Measuring and modeling for the assessment of the genetic background behind cognitive processes in donkeys.

New productive niches can offer new commercial perspectives linked to donkeys' products and human therapeutic or leisure applications. However, no assessment for selection criteria has been carried out yet. First, we assessed the animal inherent features and environmental factors that may potentially influence several cognitive processes in donkeys. Then, we aimed at describing a practical methodology to quantify such cognitive processes, seeking their inclusion in breeding and conservation programmes, through a multifactorial linear model. Sixteen cognitive process-related traits were scored on a problem-solving test in a sample of 300 Andalusian donkeys for three consecutive years from 2013 to 2015. The linear model assessed the influence and interactions of four environmental factors, sex as an animal-inherent factor, age as a covariable, and the interactions between these factors. Analyses of variance were performed with GLM procedure of SPSS Statistics for Windows, Version 24.0 software to assess the relative importance of each factor. All traits were significantly (P<0.05) affected by all factors in the model except for sex that was not significant for some of the cognitive processes, and stimulus which was not significant (P<0.05) for all of them except for the coping style related ones. The interaction between all factors within the model was non-significant (P<0.05) for almost all cognitive processes. The development of complex multifactorial models to study cognitive processes may counteract the inherent variability in behavior genetics and the estimation and prediction of related breeding parameters, key for the implementation of successful conservation programmes in apparently functionally misplaced endangered breeds.

Joint QTL mapping and gene expression analysis identify positional candidate genes influencing pork quality traits.

Meat quality traits have an increasing importance in the pig industry because of their strong impact on consumer acceptance. Herewith, we have combined phenotypic and microarray expression data to map loci with potential effects on five meat quality traits recorded in the longissimus dorsi (LD) and gluteus medius (GM) muscles of 350 Duroc pigs, i.e. pH at 24 hours post-mortem (pH24), electric conductivity (CE) and muscle redness (a*), lightness (L*) and yellowness (b*). We have found significant genome-wide associations for CE of LD on SSC4 (~104 Mb), SSC5 (~15 Mb) and SSC13 (~137 Mb), while several additional regions were significantly associated with meat quality traits at the chromosome-wide level. There was a low positional concordance between the associations found for LD and GM traits, a feature that reflects the existence of differences in the genetic determinism of meat quality phenotypes in these two muscles. The performance of an eQTL search for SNPs mapping to the regions associated with meat quality traits demonstrated that the GM a* SSC3 and pH24 SSC17 QTL display positional concordance with cis-eQTL regulating the expression of several genes with a potential role on muscle metabolism.

Differential distribution of Y-chromosome haplotypes in Swiss and Southern European goat breeds.

The analysis of Y-chromosome variation has provided valuable clues about the paternal history of domestic animal populations. The main goal of the current work was to characterize Y-chromosome diversity in 31 goat populations from Central Eastern (Switzerland and Romania) and Southern Europe (Spain and Italy) as well as in reference populations from Africa and the Near East. Towards this end, we have genotyped seven single nucleotide polymorphisms (SNPs), mapping to the SRY, ZFY, AMELY and DDX3Y Y-linked loci, in 275 bucks from 31 populations. We have observed a low level of variability in the goat Y-chromosome, with just five haplotypes segregating in the whole set of populations. We have also found that Swiss bucks carry exclusively Y1 haplotypes (Y1A: 24%, Y1B1: 15%, Y1B2: 43% and Y1C: 18%), while in Italian and Spanish bucks Y2A is the most abundant haplotype (77%). Interestingly, in Carpathian goats from Romania the Y2A haplotype is also frequent (42%). The high Y-chromosome differentiation between Swiss and Italian/Spanish breeds might be due to the post-domestication spread of two different Near Eastern genetic stocks through the Danubian and Mediterranean corridors. Historical gene flow between Southern European and Northern African goats might have also contributed to generate such pattern of genetic differentiation.

Detecting the existence of gene flow between Spanish and North African goats through a coalescent approach.

Human-driven migrations are one of the main processes shaping the genetic diversity and population structure of domestic species. However, their magnitude and direction have been rarely analysed in a statistical framework. We aimed to estimate the impact of migration on the population structure of Spanish and African goats. To achieve this goal, we analysed a dataset of 1,472 individuals typed with 23 microsatellites. Population structure of African and Spanish goats was moderate (mean FST = 0.07), with the exception of the Canarian and South African breeds that displayed a significant differentiation when compared to goats from North Africa and Nigeria. Measurement of gene flow with Migrate-n and IMa coalescent genealogy samplers supported the existence of a bidirectional gene flow between African and Spanish goats. Moreover, IMa estimates of the effective number of migrants were remarkably lower than those calculated with Migrate-n and classical approaches. Such discrepancies suggest that recent divergence, rather than extensive gene flow, is the main cause of the weak population structure observed in caprine breeds.

Effect of donkey seminal plasma on sperm movement and sperm-polymorphonuclear neutrophils attachment in vitro.

To evaluate the effect of seminal plasma in endometrial inflammation in donkeys, samples from fresh pure, fresh diluted and frozen-thawed semen of three different jackasses were co-incubated in water bath at 37°C with uterine Jennie's secretions collected 6h after artificial insemination with frozen-thawed donkey semen. Individual sperm movement parameters using the computerised sperm analysis system (CASA) and sperm-polymorphonuclear neutrophils (sperm-PMN) attachment observed in Diff-Quick stained smears were evaluated at 0, 1, 2, 3 and 4h of co-incubation. Controls consisted of incubating diluted or frozen-thawed sperm in the absence of uterine secretions. For data analyses, a repeated measures ANOVA was performed with incubation time as intra-subject factor and with treatment and donkey as inter-subject factor, followed by a post-hoc Bonferroni's test. Greater values (P<0.05) of sperm-PMN percentages and a loss of progressive motility were observed in frozen-thawed semen compared with pure and diluted fresh semen samples throughout the incubation time. In addition, the presence of seminal plasma in fresh and diluted semen samples reduced the inflammatory response of polymorphonuclear neutrophils produced after insemination by suppressing the sperm-PMN attachment in vitro. Motility sperm parameters analysed by CASA were also less affected than those in frozen-thawed semen samples. In conclusion, seminal plasma in jennies appears to have a modulation on the endometrial response after artificial insemination with frozen-thawed donkey semen. As a result, spermatozoa with the greater motility characteristics are selected.