Microarray and proteomic analysis of the human alcoholic brain

The superior frontal cortex (SFC) is selectively damaged in chronic alcohol abuse, with localized neuronal loss and tissue atrophy. Regions such as motor cortex show little neuronal loss except in severe co-morbidity (liver cirrhosis or WKS). Altered gene expression was found in microarray comparisons of alcoholic and control SFC samples [1]. We used Western blots and proteomic analysis to identify the proteins that also show differential expression. Tissue was obtained at autopsy under informed, written consent from uncomplicated alcoholics and age- and sex-matched controls. Alcoholics had consumed 80 g ethanol/day chronically (often, 200 g/day for 20 y). Controls either abstained or were social drinkers ( 20 g/day). All subjects had pathological confirmation of liver and brain diagnosis; none had been polydrug abusers. Samples were homogenized in water and clarified by brief centrifugation (1000g, 3 min) before storage at –80°C. For proteomics the thawed suspensions were centrifuged (15000g, 50 min) to prepare soluble fractions. Aliquots were pooled from SFC samples from the 5 chronic alcoholics and 5 matched controls used in the previous microarray study [1]. 2-Dimensional electrophoresis was performed in triplicate using 18 cm format pH 4–7 and pH 6–11 immobilized pH gradients for firstdimension isoelectric focusing. Following second-dimension SDS-PAGE the proteins were fluorescently stained and the images collected by densitometry. 182 proteins differed by 2-fold between cases and controls. 141 showed lower expression in alcoholics, 33 higher, and 8 were new or had disappeared. To date 63 proteins have been identified using MALDI-MS and MS-MS. Western blots were performed on uncentrifuged individual samples from 76 subjects (controls, uncomplicated alcoholics and cirrhotic alcoholics). A common standard was run on every gel. After transfer, immunolabeling, and densitometry, the intensities of the unknown bands were compared to those of the standards. We focused on proteins from transcripts that showed clear differences in a series of microarray studies, classified into common sets including Regulators of G-protein Signaling and Myelin-associated proteins. The preponderantly lower level of differentially expressed proteins in alcoholics parallels the microarray mRNA analysis in the same samples. We found that mRNA and protein expression do not frequently correspond; this may help identify pathogenic processes acting at the level of transcription, translation, or post-translationally.

Seawater carbonate chemistry and protein sports of barnicle Balanus amphitrite during experiments, 2011

The majority of benthic marine invertebrates have a complex life cycle, during which the pelagic larvae select a suitable substrate, attach to it, and then metamorphose into benthic adults. Anthropogenic ocean acidification (OA) is postulated to affect larval metamorphic success through an altered protein expression pattern (proteome structure) and post-translational modifications. To test this hypothesis, larvae of an economically and ecologically important barnacle species Balanus amphitrite, were cultured from nauplius to the cyprid stage in the present (control) and in the projected elevated concentrations of CO2 for the year 2100 (the OA treatment). Cyprid response to OA was analyzed at the total proteome level as well as two protein post-translational modification (phosphorylation and glycosylation) levels using a 2-DE based proteomic approach. The cyprid proteome showed OA-driven changes. Proteins that were differentially up or down regulated by OA come from three major groups, namely those related to energy-metabolism, respiration, and molecular chaperones, illustrating a potential strategy that the barnacle larvae may employ to tolerate OA stress. The differentially expressed proteins were tentatively identified as OA-responsive, effectively creating unique protein expression signatures for OA scenario of 2100. This study showed the promise of using a sentinel and non-model species to examine the impact of OA at the proteome level.

Proteomic profile of acute myeloid leukaemia: A review update

Proteome analysis is a complex and dynamic process that encompasses several analytical platforms that include protein sequencing, structural or expression proteomics, protein modification, sub-cellular protein localization, protein-protein interaction and biological functional proteomics. In fact, expression proteomics is extensively applied in a majority of biomarker detection studies because it provides a detailed overview of differentially expressed proteins in cellular pathways and disease processes. Proteomics are also effective and dynamic in protein-protein interactions and cross-talks between interacting molecules of the cell. Proteomics has evolved into a crucial tool used to investigate the biochemical changes that possibly lead to development of cancer biomarkers. This review draws attention to the progress and advancements in cancer proteomics technology with the aim of simplifying the understanding of the mechanisms underlying the disease and to contribute to detection of biomarkers in addition to the development of novel treatments. Given that proteome is a dynamic entity of cellular functions in health and disease, it is capable of reflecting the immediate environmental state of cells and tissues as shown in this review. The review shows the possibility of elucidating the pathophysiology of acute myeloid leukaemia (AML) through proteome expressions, thus confirming the viability of proteome analysis in profiling AML.

Quantitative Proteomic Analysis Reveals Metabolic Alterations, Calcium Dysregulation, And Increased Expression Of Extracellular Matrix Proteins In Laminin α2 Chain-deficient Muscle.

Congenital muscular dystrophy with laminin α2 chain deficiency (MDC1A) is one of the most severe forms of muscular disease and is characterized by severe muscle weakness and delayed motor milestones. The genetic basis of MDC1A is well known, yet the secondary mechanisms ultimately leading to muscle degeneration and subsequent connective tissue infiltration are not fully understood. In order to obtain new insights into the molecular mechanisms underlying MDC1A, we performed a comparative proteomic analysis of affected muscles (diaphragm and gastrocnemius) from laminin α2 chain-deficient dy(3K)/dy(3K) mice, using multidimensional protein identification technology combined with tandem mass tags. Out of the approximately 700 identified proteins, 113 and 101 proteins, respectively, were differentially expressed in the diseased gastrocnemius and diaphragm muscles compared with normal muscles. A large portion of these proteins are involved in different metabolic processes, bind calcium, or are expressed in the extracellular matrix. Our findings suggest that metabolic alterations and calcium dysregulation could be novel mechanisms that underlie MDC1A and might be targets that should be explored for therapy. Also, detailed knowledge of the composition of fibrotic tissue, rich in extracellular matrix proteins, in laminin α2 chain-deficient muscle might help in the design of future anti-fibrotic treatments. All MS data have been deposited in the ProteomeXchange with identifier PXD000978 (http://proteomecentral.proteomexchange.org/dataset/PXD000978).

Genes that code for T cell signaling proteins establish transcriptional regulatory networks during thymus ontogeny

Gene expression profiling by cDNA microarrays during murine thymus ontogeny has contributed to dissecting the large-scale molecular genetics of T cell maturation. Gene profiling, although useful for characterizing the thymus developmental phases and identifying the differentially expressed genes, does not permit the determination of possible interactions between genes. In order to reconstruct genetic interactions, on RNA level, within thymocyte differentiation, a pair of microarrays containing a total of 1,576 cDNA sequences derived from the IMAGE MTB library was applied on samples of developing thymuses (14-17 days of gestation). The data were analyzed using the GeneNetwork program. Genes that were previously identified as differentially expressed during thymus ontogeny showed their relationships with several other genes. The present method provided the detection of gene nodes coding for proteins implicated in the calcium signaling pathway, such as Prrg2 and Stxbp3, and in protein transport toward the cell membrane, such as Gosr2. The results demonstrate the feasibility of reconstructing networks based on cDNA microarray gene expression determinations, contributing to a clearer understanding of the complex interactions between genes involved in thymus/thymocyte development.

p53-induced protein with a death domain (PIDD) isoforms differentially activate nuclear factor-kappaB and caspase-2 in response to genotoxic stress

Cells respond to DNA damage in a complex way and the fate of damaged cells depends on the balance between pro- and antiapoptotic signals. This is of crucial importance in cancer as genotoxic stress is implied both in oncogenesis and in classical tumor therapies. p53-induced protein with a death domain (PIDD), initially described as a p53-inducible gene, is one of the molecular switches able to activate a survival or apoptotic program. Two isoforms of PIDD, PIDD (isoform 1) and LRDD (isoform 2), have already been reported and we describe here a third isoform. These three isoforms are differentially expressed in tissues and cell lines. Genotoxic stress only affects PIDD isoform 3 mRNA levels, whereas isoforms 1 and 2 mRNA levels remain unchanged. All isoforms are capable of activating nuclear factor-kappaB in response to genotoxic stress, but only isoform 1 interacts with RIP-associated ICH-1/CED-3 homologous protein with a death domain and activates caspase-2. Isoform 2 counteracts the pro-apoptotic function of isoform 1, whereas isoform 3 enhances it. Thus, the differential splicing of PIDD mRNA leads to the formation of at least three proteins with antagonizing/agonizing functions, thereby regulating cell fate in response to DNA damage

Formin proteins in normal tissues and cancer

The actin cytoskeleton is a dynamic structure that determines cell shape. Actin turnover is mandatory for migration in normal and malignant cells. In epithelial cancers invasion is frequently accompanied by epithelial to mesenchymal transition (EMT). In EMT, cancer cells acquire a migratory phenotype through transcriptional reprogramming. EMT requires substantial re-organization of actin. During the past decade, new actin regulating proteins have been discovered. Among these are members of the formin family. To study formin expression in tissues and cells, antibodies for detection of formin proteins FMNL1 (Formin-like protein 1), FMNL2 (Formin-like protein 2) and FHOD1 (Formin homology 2 domain containing protein 1) were used. The expression of formins was characterized in normal tissues and selected cancers using immunohistochemistry. The functional roles of formins were studied in cancer cell lines. We found that FMNL2 is widely expressed. It is a filopodial component in cultured melanoma cells. In clinical melanoma, FMNL2 expression has prognostic significance. FHOD1 is a formin expressed in mesenchymal cell types. FHOD1 expression is increased in oral squamous cell carcinoma (SCC) EMT. Importantly, FHOD1 participates in invasion of cultured oral SCC cells. FMNL1 expression is low in normal epithelia, but high in leukocytes and smooth muscle cells. Expression of FMNL1 can be found in carcinoma; we detected FMNL1 expressing cells in basal type of breast cancer. Our results indicate that formins are differentially expressed in normal tissues and that their expression may shift in cancer. Functionally FMNL2 and FHOD1 participate in processes related to cancer progression. Studying formins is increasingly important since they are potential drug targets.

Specific alterations in plasma proteins during depressed, manic, and euthymic states of bipolar disorder

Bipolar disorder (BD) is a common psychiatric mood disorder affecting more than 1-2% of the general population of different European countries. Unfortunately, there is no objective laboratory-based test to aid BD diagnosis or monitor its progression, and little is known about the molecular basis of BD. Here, we performed a comparative proteomic study to identify differentially expressed plasma proteins in various BD mood states (depressed BD, manic BD, and euthymic BD) relative to healthy controls. A total of 10 euthymic BD, 20 depressed BD, 15 manic BD, and 20 demographically matched healthy control subjects were recruited. Seven high-abundance proteins were immunodepleted in plasma samples from the 4 experimental groups, which were then subjected to proteome-wide expression profiling by two-dimensional electrophoresis and matrix-assisted laser desorption/ionization-time-of-flight/time-of-flight tandem mass spectrometry. Proteomic results were validated by immunoblotting and bioinformatically analyzed using MetaCore. From a total of 32 proteins identified with 1.5-fold changes in expression compared with healthy controls, 16 proteins were perturbed in BD independent of mood state, while 16 proteins were specifically associated with particular BD mood states. Two mood-independent differential proteins, apolipoprotein (Apo) A1 and Apo L1, suggest that BD pathophysiology may be associated with early perturbations in lipid metabolism. Moreover, down-regulation of one mood-dependent protein, carbonic anhydrase 1 (CA-1), suggests it may be involved in the pathophysiology of depressive episodes in BD. Thus, BD pathophysiology may be associated with early perturbations in lipid metabolism that are independent of mood state, while CA-1 may be involved in the pathophysiology of depressive episodes.

Germin and germin-like proteins: evolution, structure, and function

Germin and germin-like proteins (GLPs) are encoded by a family of genes found in all plants. They are part of the cupin superfamily of biochemically diverse proteins, a superfamily that has a conserved tertiary structure, though with limited similarity in primary sequence. The subgroups of GLPs have different enzyme functions that include the two hydrogen peroxide-generating enzymes, oxalate oxidase (OxO) and superoxide dismutase. This review summarizes the sequence and structural details of GLPs and also discusses their evolutionary progression, particularly their amplification in gene number during the evolution of the land plants. In terms of function, the GLPs are known to be differentially expressed during specific periods of plant growth and development, a pattern of evolutionary subfunctionalization. They are also implicated in the response of plants to biotic (viruses, bacteria, mycorrhizae, fungi, insects, nematodes, and parasitic plants) and abiotic (salt, heat/cold, drought, nutrient, and metal) stress. Most detailed data come from studies of fungal pathogenesis in cereals. This involvement with the protection of plants from environmental stress of various types has led to numerous plant breeding studies that have found links between GLPs and QTLs for disease and stress resistance. In addition the OxO enzyme has considerable commercial significance, based principally on its use in the medical diagnosis of oxalate concentration in plasma and urine. Finally, this review provides information on the nutritional importance of these proteins in the human diet, as several members are known to be allergenic, a feature related to their thermal stability and evolutionary connection to the seed storage proteins, also members of the cupin superfamily.

Structural synaptic elements are differentially regulated in superior temporal cortex of schizophrenia patients

Inaccurate wiring and synaptic pathology appear to be major hallmarks of schizophrenia. A variety of gene products involved in synaptic neurotransmission and receptor signaling are differentially expressed in brains of schizophrenia patients. However, synaptic pathology may also develop by improper expression of intra- and extra-cellular structural elements weakening synaptic stability. Therefore, we have investigated transcription of these elements in the left superior temporal gyrus of 10 schizophrenia patients and 10 healthy controls by genome-wide microarrays (Illumina). Fourteen up-regulated and 22 downregulated genes encoding structural elements were chosen from the lists of differentially regulated genes for further qRT-PCR analysis. Almost all genes confirmed by this method were downregulated. Their gene products belonged to vesicle-associated proteins, that is, synaptotagmin 6 and syntaxin 12, to cytoskeletal proteins, like myosin 6, pleckstrin, or to proteins of the extracellular matrix, such as collagens, or laminin C3. Our results underline the pivotal roles of structural genes that control formation and stabilization of pre- and post-synaptic elements or influence axon guidance in schizophrenia. The glial origin of collagen or laminin highlights the close interrelationship between neurons and glial cells in establishment and maintenance of synaptic strength and plasticity. It is hypothesized that abnormal expression of these and related genes has a major impact on the pathophysiology of schizophrenia.

SAGE analysis highlights the putative role of underexpression of ribosomal proteins in GH-secreting pituitary adenomas

Background: Although the molecular pathogenesis of pituitary adenomas has been assessed by several different techniques, it still remains partially unclear. Ribosomal proteins (RPs) have been recently related to human tumorigenesis, but they have not yet been evaluated in pituitary tumorigenesis. Objective: The aim of this study was to introduce serial analysis of gene expression (SAGE), a high-throughput method, in pituitary research in order to compare differential gene expression. Methods: Two SAGE cDNA libraries were constructed, one using a pool of mRNA obtained from five GH-secreting pituitary tumors and another from three normal pituitaries. Genes differentially expressed between the libraries were further validated by real-time PCR in 22 GH-secreting pituitary tumors and in 15 normal pituitaries. Results: Computer-generated genomic analysis tools identified 13 722 and 14 993 exclusive genes in normal and adenoma libraries respectively. Both shared 6497 genes, 2188 were underexpressed and 4309 overexpressed in tumoral library. In adenoma library, 33 genes encoding RPs were underexpressed. Among these, RPSA, RPS3, RPS14, and RPS29 were validated by real-time PCR. Conclusion: We report the first SAGE library from normal pituitary tissue and GH-secreting pituitary tumor, which provide quantitative assessment of cellular transcriptome. We also validated some downregulated genes encoding RPs. Altogether, the present data suggest that the underexpression of the studied RP genes possibly collaborates directly or indirectly with other genes to modify cell cycle arrest, DNA repair, and apoptosis, leading to an environment that might have a putative role in the tumorigenesis, introducing new perspectives for further studies on molecular genesis of somatotrophinomas.

Novel functional profiling approach combining reverse phase protein microarrays and human 3-D ex vivo tissue cultures: expression of apoptosis-related proteins in human colon cancer

Cancer is caused by a complex pattern of molecular perturbations. To understand the biology of cancer, it is thus important to look at the activation state of key proteins and signaling networks. The limited amount of available sample material from patients and the complexity of protein expression patterns make the use of traditional protein analysis methods particularly difficult. In addition, the only approach that is currently available for performing functional studies is the use of serial biopsies, which is limited by ethical constraints and patient acceptance. The goal of this work was to establish a 3-D ex vivo culture technique in combination with reverse-phase protein microarrays (RPPM) as a novel experimental tool for use in cancer research. The RPPM platform allows the parallel profiling of large numbers of protein analytes to determine their relative abundance and activation level. Cancer tissue and the respective corresponding normal tissue controls from patients with colorectal cancer were cultured ex vivo. At various time points, the cultured samples were processed into lysates and analyzed on RPPM to assess the expression of carcinoembryonic antigen (CEA) and 24 proteins involved in the regulation of apoptosis. The methodology displayed good robustness and low system noise. As a proof of concept, CEA expression was significantly higher in tumor compared with normal tissue (p<0.0001). The caspase 9 expression signal was lower in tumor tissue than in normal tissue (p<0.001). Cleaved Caspase 8 (p=0.014), Bad (p=0.007), Bim (p=0.007), p73 (p=0.005), PARP (p<0.001), and cleaved PARP (p=0.007) were differentially expressed in normal liver and normal colon tissue. We demonstrate here the feasibility of using RPPM technology with 3-D ex vivo cultured samples. This approach is useful for investigating complex patterns of protein expression and modification over time. It should allow functional proteomics in patient samples with various applications such as pharmacodynamic analyses in drug development.

Diversity and phylogeny of gephyrin: Tissue-specific splice variants, gene structure, and sequence similarities to molybdenum cofactor-synthesizing and cytoskeleton-associated proteins

Gephyrin is essential for both the postsynaptic localization of inhibitory neurotransmitter receptors in the central nervous system and the biosynthesis of the molybdenum cofactor (Moco) in different peripheral organs. Several alternatively spliced gephyrin transcripts have been identified in rat brain that differ in their 5′ coding regions. Here, we describe gephyrin splice variants that are differentially expressed in non-neuronal tissues and different regions of the adult mouse brain. Analysis of the murine gephyrin gene indicates a highly mosaic organization, with eight of its 29 exons corresponding to the alternatively spliced regions identified by cDNA sequencing. The N- and C-terminal domains of gephyrin encoded by exons 3–7 and 16–29, respectively, display sequence similarities to bacterial, invertebrate, and plant proteins involved in Moco biosynthesis, whereas the central exons 8, 13, and 14 encode motifs that may mediate oligomerization and tubulin binding. Our data are consistent with gephyrin having evolved from a Moco biosynthetic protein by insertion of protein interaction sequences.

Comparative expressed-sequence-tag analysis of differential gene expression profiles in PC-12 cells before and after nerve growth factor treatment.

Nerve growth factor-induced differentiation of adrenal chromaffin PC-12 cells to a neuronal phenotype involves alterations in gene expression and represents a model system to study neuronal differentiation. We have used the expressed-sequence-tag approach to identify approximately 600 differentially expressed mRNAs in untreated and nerve growth factor-treated PC-12 cells that encode proteins with diverse structural and biochemical functions. Many of these mRNAs encode proteins belonging to cellular pathways not previously known to be regulated by nerve growth factor. Comparative expressed-sequence-tag analysis provides a basis for surveying global changes in gene-expression patterns in response to biological signals at an unprecedented scale, is a powerful tool for identifying potential interactions between different cellular pathways, and allows the gene-expression profiles of individual genes belonging to a particular pathway to be followed.

Acute phase proteins and biomarkers for health in chickens

Acute phase proteins (APPs) are proteins synthesised predominantly in the liver, whose plasma concentrations increase (positive APP) or decrease (negative APP) as a result of infection, inflammation, trauma and tissue injury. They also change as a result of the introduction of immunogens such as bacterial lipopolysaccharide (LPS), turpentine and vaccination. While publications on APPs in chickens are numerous, the limited availability of anti-sera and commercial ELISAs has resulted in a lot of information on only a few APPs. Disease is a threat to the poultry industry, as pathogens have the potential to evolve, spread and cause rapid onset of disease that is detrimental to the welfare of birds. Low level, sub-acute disease with non-specific, often undiagnosed causes can greatly affect bird health and growth and impact greatly on productivity and profitability. Developing and validating methods to measure and characterise APPs in chickens will allow these proteins to be used diagnostically for monitoring flock health. Using immune parameters such as APPs that correlate with disease resistance or improvements in production and welfare will allow the use of APPs as selection parameters for breeding to be evaluated. For APPs to be useful parameters on which to evaluate chicken health, information on normal APP concentrations is required. Ceruloplasmin (Cp) and PIT54 concentrations were found to be much lower in healthy birds form commercial production farms than the reported normal values obtained from the literature. These APPs were found to be significantly higher in culled birds from a commercial farm and Cp, PIT54 and ovotransferrin (Ovt) were significantly higher in birds classified as having obvious gait defects. Using quantitative shotgun proteomics to identify the differentially abundant proteins between three pools: highly acute phase (HAP), acute phase (AP) and non-acute phase (NAP), generated data from which a selection of proteins, based on the fold difference between the three pools was made. These proteins were targeted on a individual samples alongside proteins known to be APPs in chickens or other species: serum amyloid A (SAA), C-reactive protein (CRP), Ovt, apolipoprotein A-I (apo-AI), transthyretin (Ttn), haemopexin (Hpx) and PIT54. Together with immunoassay data for SAA, Ovt, alpha-1-acid glycoprotein (AGP) and Cp the results of this research reveal that SAA is the only major APP in chickens. Ovotransferrin and AGP behave as moderate APPs while PIT54 and Cp are minor APPs. Haemopexin was not significantly different between the three acute phase groups. Apolipoprotein AI and Ttn were significantly lower in the HAP and AP groups and as such can be classed as negative APPs. In an effort to identify CRP, multiple anti-sera cross reacting with CRP from other species were used and a phosphorylcholine column known to affinity purify CRP were used. Enriched fractions containing low molecular weight proteins, elutions from the affinity column together with HAP, AP and NAP pooled samples were applied to a Q-Exactive Hybrid Quadrupole–Orbitrap mass spectrometer (Thermo Scientific) for Shotgun analysis and CRP was not identified. It would appear that CRP is not present as a plasma protein constitutively or during an APR in chickens and as such is not an APP in this species. Of the proteins targeted as possible novel biomarkers of the APR in chickens mannan binding lectin associated serine protease-2, α-2-HS-glycoprotein (fetuin) and major facilitator superfamily domain-containing protein 10 were reduced in abundance in the HAP group, behaving as negative biomarkers. Myeloid protein and putative ISG(12)2 were positively associated with the acute phase being significantly higher in the HAP and AP groups. The protein cathepsin D was significantly higher in both HAP and AP compared to the NAP indicating that of all the proteins targeted, this appears to have the most potential as a biomarker of the acute phase, as it was significantly increased in the AP as well as the HAP group. To evaluate APPs and investigate biomarkers of intestinal health, a study using re-used poultry litter was undertaken. The introduction of litter at 12 days of age did not significantly increase any APPs measured using immunoassays and quantitative proteomics at 3, 6 and 10 days post introduction. While no APP was found to be significantly different between the challenged and control groups at anytime point, the APPs AGP, SAA and Hpx did increase over time in all birds. The protein apolipoprotein AIV (apo-AIV) was targeted as a possible APP and because of its reported role in controlling satiety. An ELISA was developed, successfully validated and used to measure apo-AIV in this study. While no significant differences in apo-AIV plasma concentrations between challenged and control groups were identified apo-AIV plasma concentrations did change significantly between certain time points in challenged and control groups. Apoliporotein AIV does not appear to behave as an APP in chickens, as it was not significantly different between acute phase groups. The actin associated proteins villin and gelsolin were investigated as possible biomarkers of intestinal health. Villin was found not to be present in the plasma of chickens and as such not a biomarker target. Gelsolin was found not to be differentially expressed during the acute phase or as a result of intestinal challenge. Finally a proteomic approach was undertaken to investigate gastrocnemius tendon (GT) rupture in broiler chickens with a view of elucidating to and identify proteins associated with risk of rupture. A number of proteins were found to be differentially expressed between tendon pools and further work would enable further detailing of these findings. In conclusion this work has made a number of novel findings and addressed a number of data poor areas. The area of chicken APPs research has stagnated over the last 15 years with publications becoming repetitive and reliant on a small number of immunoassays. This work has sought to characterise the classic APPs in chickens, and use a quantitative proteomic approach to measure and categorise them. This method was also used to take a fresh approach to biomarker identification for both the APR and intestinal health. The development and validation of assays for Ovt and apo-AIV and the shotgun data mean that these proteins can be further characterised in chickens with a view of applying their measurement to diagnostics and selective breeding programs.