Reference genes for quantitative gene expression studies in multiple avian species

Postbiotic-based recombinant receptor activator of NF-κB ligand enhanced oral vaccine efficiency in chicken

Functional M cells are differentiated by receptor activator of NF-κB ligand (RANKL) and capture of luminal antigens to initiate immune responses. We aimed to use postbiotic-based recombinant chicken RANKL (cRANKL) to promote M cell differentiation and test the efficacy of oral vaccines. Chicks were divided into three groups that were administered phosphate-buffered saline (PBS), cell extracts of wild-type Lactococcus lactis subsp. lactis IL1403 (WT_CE), or cell extracts of recombinant L. lactis expressing cRANKL (cRANKL_CE). The expression of the M cell marker was measured, and the gut microbiome was profiled. The efficiency of the infectious bursal disease (IBD) vaccine was tested after 12 consecutive days of administering cRANKL_CE. The chickens that were administered cRANKL_CE (p = 0.038) had significantly higher Annexin A5 (ANXA5) mRNA expression levels than those in the PBS group (PBS vs. WT_CE, p = 0.657). In the gut microbiome analysis, no significant changes were observed. However, the relative abundance of Escherichia-Shigella was negatively correlated (r =  - 0.43, p = 0.019) with ANXA5 mRNA expression in Peyer's patches. cRANKL_CE/IBD (p = 0.018) had significantly higher IBD-specific faecal IgA levels than PBS/IBD (PBS/IBD vs. WT_CE/IBD, p = 0.217). Postbiotic-based recombinant cRANKL effectively improved the expression of M cell markers and the efficiency of oral vaccines. No significant changes were observed in the gut microbiome after administration of postbiotic-based recombinant cRANKL. This strategy can be used for the development of feed additives and adjuvants. KEY POINTS: • Postbiotic-based recombinant cRANKL enhanced the expression of ANXA5 in chicken. • The relative abundance of Escherichia-Shigella was negatively correlated with ANXA5 expression. • Postbiotic-based recombinant cRANKL effectively improved the efficiency of oral vaccine.

Selection of suitable reference genes for normalization of RT-qPCR in three tissues of Northern bobwhite (Colinus virginianus) infected with eyeworm (Oxyspirura petrowi)

BACKGROUND

The Northern bobwhite (Colinus virginianus) is an economically important, and popular game bird in North America. Northern bobwhites have experiencing declines of > 3.5% annually in recent decades due to several factors. The eyeworm Oxyspirura petrowi is a nematode parasite frequently found in the eyes of bobwhites. Although reported frequently in wild bobwhites, there is no research to understand the host-parasite mechanism. Hence, it is important to investigate mechanisms of eyeworm invasion and immune modulation in bobwhite. Cytokine gene expression using RT-PCR is widely used to identify the innate immune response of a host to an infection.

METHODOLOGY

In this study, we evaluated ten reference genes (HMBS, RPL19, RPL32, RPS7, RPS8, TATA, SDHA, YWHAZ, GAPDH, and ACTB) for their stability across three tissues (liver, spleen, and caecal tonsils) of control and O. petrowi infected Northern bobwhites. Primer efficiency and reference genes stability were assessed using GeNorm, NormFinder, and BestKeeper.

RESULTS

Expression of these reference genes with respect to O. petrowi infection in bobwhites showed RPL32 and HMBS were the most stable genes in the liver, HMBS and SDHA were the most stable genes in the spleen, and HMBS and YWHAZ were equally stable reference genes in the caecal tonsils.

CONCLUSION

Based on the geometric mean of all three analyses, our results indicate that the combination of RPL32 and HMBS for the liver, HMBS and SDHA for the spleen, and YWHAZ and HMBS for caecal tonsils might be used as reference genes for normalization in gene expression investigations on Northern bobwhites.

Developmental toxicity of pyriproxyfen induces changes in the ultrastructure of neural cells and in the process of skull ossification

Some studies relate the use of pyriproxyfen (PPF) in drinking water with damage to embryonic neurodevelopment, including a supposed association with cases of microcephaly. However, the effects on neural cells and skull ossification in embryos remain unclear. This study aims to investigate the effects of PPF on the structure and ultrastructure of brain cells and its influence on the skull ossification process during embryonic development. Chicken embryos, used as an experimental model, were exposed to concentrations of 0.01 and 10 mg/l PPF at E1. The findings demonstrated that PPF led to notable ultrastructural alterations such as reduced cilia and microvilli of ependymal cells and damage to mitochondria, endoplasmic reticulum, Golgi bodies, and cell membranes in neural cells. The frequency of changes and the degree of these cell damage between the forebrain and midbrain were similar. PPF induced a reduction in fox3 transcript levels, specific for differentiation of neurons, and a reduction in the NeuN protein content related to mature neurons and dendritic branches. PPF impacted the ossification process of the skull, as evidenced by the increase in the ossified area and the decrease in inter-bone spacing. In conclusion, this study highlights the ability of PPF to affect neurodevelopmental processes by inducing ultrastructural damage to neural cells, concomitant with a reduction in NeuN and fox3 expression. This detrimental impact coupled with deficiencies in skull ossification can prevent the proper growth and development of the brain.

Expression patterns of heat-shock genes during stopover and the trade-off between refueling and stress response in a passerine migrant

Migrating birds are often exposed to variable environments and face a multitude of stress exposures along their long-distance flights. During stopover refueling, migratory birds must balance the need to accumulate energy reserves to continue their migration with the need to respond to environmental and physiological stressors. We examined the gene expression patterns of different Heat Shock Proteins (HSPs) in migrating birds during stopover at different body condition states (lean vs. fat), to provide some first insights on the role of HSPs in bird migration and explore the concept of a trade-off between refueling and stress response. Our results showed upregulation of HSP expression at release that could be associated with muscle growth and increased cholesterol and lipid synthesis needed for birds to fuel their upcoming migration. On the other hand, during capture, upregulation of HSP5 could be attributed to physiological recovery from the non-stop endurance flight when crossing the Sahara Desert-Mediterranean Sea ecological barrier. All birds significantly increased their fuel loads up to 48% of lean body mass and we provide evidence for muscle rebuilding during stopover as flight muscle mass increased by 10%, highlighting the fact that stopover sites can play a major role in the physiological recovery of migrants.

Sarcolipin relates to fattening, but not sarco/endoplasmic reticulum Ca2+-ATPase uncoupling, in captive migratory gray catbirds

In order to complete their energetically demanding journeys, migratory birds undergo a suite of physiological changes to prepare for long-duration endurance flight, including hyperphagia, fat deposition, reliance on fat as a fuel source, and flight muscle hypertrophy. In mammalian muscle, SLN is a small regulatory protein which binds to sarco/endoplasmic reticulum Ca2+-ATPase (SERCA) and uncouples Ca2+ transport from ATP hydrolysis, increasing energy consumption, heat production, and cytosolic Ca2+ transients that signal for mitochondrial biogenesis, fatigue resistance and a shift to fatty acid oxidation. Using a photoperiod manipulation of captive gray catbirds (Dumetella carolinensis), we investigated whether SLN may play a role in coordinating the development of the migratory phenotype. In response to long-day photostimulation, catbirds demonstrated migratory restlessness and significant body fat stores, alongside higher SLN transcription while SERCA2 remained constant. SLN transcription was strongly correlated with h-FABP and PGC1α transcription, as well as fat mass. However, SLN was not significantly correlated with HOAD or CD36 transcripts or measurements of SERCA activity, SR membrane Ca2+ leak, Ca2+ uptake rates, pumping efficiency or mitochondrial biogenesis. Therefore, SLN may be involved in the process of storing fat and shifting to fat as a fuel, but the mechanism of its involvement remains unclear.

Varied impacts of social relationships on neuroendocrine state

Social relationships, affiliative social attachments, are important for many species. The best studied types of relationships are monogamous pair bonds. However, it remains unclear how generalizable models of pair bonding are across types of social attachments. Zebra finches are a fascinating system to explore the neurobiology of social relationships because they form various adult bonds with both same- and opposite-sex partners. To test whether different bonds are supported by a single brain network, we quantified individuals' neuroendocrine state after either 24 h or 2 weeks of co-housing with a novel same- or opposite-sex partner. We defined neuroendocrine state by the expression of 22 genes related to 4 major signaling pathways (dopamine, steroid, nonapeptide, and opioid) in six brain regions associated with affiliation or communication [nucleus accumbens (NAc), nucleus taeniae of the amygdala (TnA), medial preoptic area (POM), and periaqueductal gray (PAG), ventral tegmental area, and auditory cortex]. Overall, we found dissociable effects of social contexts (same- or opposite-sex partnerships) and duration of co-housing. Social bonding impacted the neuroendocrine state of four regions in males (NAc, TnA, POM, and PAG) and three regions in females (NAc, TnA, and POM). Monogamous pair bonding specifically appeared to impact male NAc. However, the patterns of gene expression in zebra finches were different than has previously been reported in mammals. Together, our results support the view that there are numerous mechanisms regulating social relationships and highlight the need to further our understanding of how social interactions shape social bonds.

Influence of Live Attenuated Salmonella Vaccines on Cecal Microbiome Composition and Microbiota Abundances in Young Broiler Chickens

Salmonellosis is a global food safety challenge caused by Salmonella, a gram-negative bacterium of zoonotic importance. Poultry is considered a major reservoir for the pathogen, and humans are exposed through consumption of raw or undercooked products derived from them. Prophylaxis of Salmonella in poultry farms generally mainly involves biosecurity measures, flock testing and culling, use of antibiotics, and vaccination programs. For decades, the use of antibiotics has been a common practice to limit poultry contamination with important pathogenic bacteria such as Salmonella at the farm level. However, due to an increasing prevalence of resistance, non-therapeutic use of antibiotics in animal production has been banned in many parts of the world. This has prompted the search for non-antimicrobial alternatives. Live vaccines are among the developed and currently used methods for Salmonella control. However, their mechanism of action, particularly the effect they might have on commensal gut microbiota, is not well understood. In this study, three different commercial live attenuated Salmonella vaccines (AviPro^® Salmonella Vac T, AviPro^® Salmonella DUO, and AviPro^® Salmonella Vac E) were used to orally vaccinate broiler chickens, and cecal contents were collected for microbiomes analysis by 16S rRNA next generation sequencing. Quantitative real-time PCR (qPCR) was used to study the cecal immune-related genes expression in the treatment groups, while Salmonella-specific antibodies were analyzed from sera and cecal extracts by enzyme-linked immunosorbent assay (ELISA). We show that vaccination with live attenuated Salmonella vaccines had a significant influence on the variability of the broiler cecal microbiota (p = 0.016). Furthermore, the vaccines AviPro^® Salmonella Vac T and AviPro^® Salmonella DUO, but not AviPro^® Salmonella Vac E, had a significant effect (p = 0.024) on microbiota composition. This suggests that the live vaccine type used can differently alter the microbiota profiles, driving the gut colonization resistance and immune responses to pathogenic bacteria, and might impact the overall chicken health and productivity. Further investigation is, however, required to confirm this.

Identification of stable reference genes for quantitative gene expression analysis in the duodenum of meat-type ducks

Quantitative polymerase chain reaction (qPCR) is an important method to detect gene expression at the molecular level. The selection of appropriate housekeeping genes is the key to accurately calculating the expression level of target genes and conducting gene function studies. In this study, the expression of eight candidate reference genes, glyceraldehyde-3-phosphate dehydrogenase (GAPDH), beta-actin (β-actin), 18S ribosomal RNA (18S rRNA), hydroxymethylbilane synthase (HMBS), hypoxanthine phosphoribosyltransferase 1 (HPRT1), TATA box binding protein (TBP), ribosomal protein L13 (RPL13), and tyrosine 3-monooxygenase/tryptophan 5-monooxygenase activation protein (YWHAZ), in the duodenal epithelial tissue of 42-day-old meat-type ducks were detected using qPCR. Furthermore, their expression stability was analyzed using the geNorm, NormFinder, and BestKeeper programs. The results indicated that HMBS and YWHAZ were the most stably expressed genes. All three programs indicated that the expression of 18S rRNA was the least stable, making it unsuitable for the study of gene expression in meat-type duck tissues. This study provides stable reference genes for gene expression analysis and contributes to further studies on the gene function of meat-type ducks.

Genomic, biochemical and expressional properties reveal strong conservation of the CLCA2 gene in birds and mammals

Recent studies have revealed the dynamic and complex evolution of CLCA1 gene homologues in and between mammals and birds with a particularly high diversity in mammals. In contrast, CLCA2 has only been found as a single copy gene in mammals, to date. Furthermore, CLCA2 has only been investigated in few mammalian species but not in birds. Here, we established core genomic, protein biochemical and expressional properties of CLCA2 in several bird species and compared them with mammalian CLCA2. Chicken, turkey, quail and ostrich CLCA2 were compared to their mammalian orthologues using in silico, biochemical and expressional analyses. CLCA2 was found highly conserved not only at the level of genomic and exon architecture but also in terms of the canonical CLCA2 protein domain organization. The putatively prototypical galline CLCA2 (gCLCA2) was cloned and immunoblotting as well as immunofluorescence analyses of heterologously expressed gCLCA2 revealed protein cleavage, glycosylation patterns and anchoring in the plasma membrane similar to those of most mammalian CLCA2 orthologues. Immunohistochemistry found highly conserved CLCA2 expression in epidermal keratinocytes in all birds and mammals investigated. Our results suggest a highly conserved and likely evolutionarily indispensable role of CLCA2 in keratinocyte function. Its high degree of conservation on the genomic, biochemical and expressional levels stands in contrast to the dynamic structural complexities and proposed functional diversifications between mammalian and avian CLCA1 homologues, insinuating a significant degree of negative selection of CLCA2 orthologues among birds and mammals. Finally, and again in contrast to CLCA1, the high conservation of CLCA2 makes it a strong candidate for studying basic properties of the functionally still widely unresolved CLCA gene family.

Embryonic transcriptome unravels mechanisms and pathways underlying embryonic development with respect to muscle growth, egg production, and plumage formation in native and broiler chickens

Background: Muscle development, egg production, and plumage colors are different between native and broiler chickens. The study was designed to investigate why improved Aseel (PD4) is colorful, stronger, and grew slowly compared with the control broiler (CB).

Methods: A microarray was conducted using the 7th-day embryo (7EB) and 18th-day thigh muscle (18TM) of improved Aseel and broiler, respectively. Also, we have selected 24 Gallus gallus candidate reference genes from NCBI, and total RNA was isolated from the broiler, improved Aseel embryo tissues, and their expression profiles were studied by real-time quantitative PCR (qPCR). Furthermore, microarray data were validated with qPCR using improved Aseel and broiler embryo tissues.

Results: In the differential transcripts screening, all the transcripts obtained by microarray of slow and fast growth groups were screened by fold change ≥ 1 and false discovery rate (FDR) ≤ 0.05. In total, 8,069 transcripts were differentially expressed between the 7EB and 18TM of PD4 compared to the CB. A further analysis showed that a high number of transcripts are differentially regulated in the 7EB of PD4 (6,896) and fewer transcripts are differentially regulated (1,173) in the 18TM of PD4 compared to the CB. On the 7th^- and 18th-day PD4 embryos, 3,890, 3,006, 745, and 428 transcripts were up- and downregulated, respectively. The commonly up- and downregulated transcripts are 91 and 44 between the 7th- and 18th-day of embryos. In addition, the best housekeeping gene was identified. Furthermore, we validated the differentially expressed genes (DEGs) related to muscle growth, myostatin signaling and development, and fatty acid metabolism genes in PD4 and CB embryo tissues by qPCR, and the results correlated with microarray expression data.

Conclusion: Our study identified DEGs that regulate the myostatin signaling and differentiation pathway; glycolysis and gluconeogenesis; fatty acid metabolism; Jak-STAT, mTOR, and TGF-β signaling pathways; tryptophan metabolism; and PI3K-Akt signaling pathways in PD4. The results revealed that the gene expression architecture is present in the improved Aseel exhibiting embryo growth that will help improve muscle development, differentiation, egg production, protein synthesis, and plumage formation in PD4 native chickens. Our findings may be used as a model for improving the growth in Aseel as well as optimizing the growth in the broiler.

Molecular Surveillance for Lymphoproliferative Disease Virus and Reticuloendotheliosis Virus in Rio Grande Wild Turkeys (Meleagris gallopavo intermedia) in Texas, USA

Reticuloendotheliosis virus (REV) and lymphoproliferative disease virus (LPDV) are avian retroviruses that can cause neoplastic disease and present with similar pathologies. Lymphoproliferative disease virus has been reported in the Eastern US and states bordering Texas, USA, but has not been previously detected within the state. In a prior study, we detected REV in native Rio Grande Wild Turkeys (Meleagris gallopavo intermedia) and an Eastern Wild Turkey (Meleagris gallopavo silvestris) originating from West Virginia. Given LPDV detection in states bordering Texas and our finding of an REV-positive Eastern Wild Turkey imported from a LPDV endemic region, we sought to determine LPDV prevalence in Texas and continue surveillance for REV. During 2018-20, dried blood spots from 373 individual Rio Grande Wild Turkeys from 20 different counties were tested for the presence of proviral REV or LPDV DNA. In affected counties, approximately 4% of individuals were infected with REV (7/197) or LPDV (10/273) and one bird was coinfected with both viruses. Phylogenetic analysis indicated a close relationship of the LPDV isolates to variants from other Southern and Central states. This study provides molecular evidence of LPDV in Texas, and continued surveillance is necessary to determine the potential effects of the virus on reproductive success, coinfections, and overall health of Wild Turkey populations.

Investigation of chicken housekeeping genes using next-generation sequencing data

Accurate normalization of the gene expression assays, using housekeeping genes (HKGs), is critically necessary. To do so, selection of a proper set of HKGs for a specific experiment is of great importance. Despite many studies, there is no consensus about the suitable set of HKGs for implementing in the quantitative real-time PCR analyses of chicken tissues. A limited number of HKGs have been widely used. However, wide utilization of a little number of HKGs for all tissues is challenging. The emergence of high-throughput gene expression RNA-seq data has enabled the simultaneous comparison of the stability of multiple HKGs. Therefore, employing the average coefficient of variations of at least three datasets per tissue, we sorted all reliably expressed genes (REGs; with FPKM ≥ 1 in at least one sample) and introduced the top 10 most suitable and stable reference genes for each of the 16 chicken tissues. We evaluated the consistency of the results of five tissues using the same methodology on other datasets. Furthermore, we assessed 96 previously widely used HKGs (WU-HKGs) in order to challenge the accuracy of the previous studies. The New Tuxedo software suite was used for the main analyses. The results revealed novel, different sets of reference genes for each of the tissues with 17 common genes among the top 10 genes lists of 16 tissues. The results did disprove the suitability of WU-HKGs such as Actb, Ldha, Scd, B2m, and Hprt1 for any of the tissues examined. On the contrary, a total of 6, 13, 14, 23, and 32 validated housekeeping genes (V-HKGs) were discovered as the most stable and suitable reference genes for muscle, spleen, liver, heart, and kidney tissues, respectively. Although we identified a few new HKGs usable for multiple tissues, the selection of suitable HKGs is required to be tissue specific. The newly introduced reference genes from the present study, despite lacking experimental validation, will be able to contribute to the more accurate normalization for future expression analysis of chicken genes.

Identification of suitable reference genes for normalization of quantitative real-time PCR-based gene expression in chicken (Gallus gallus)

The recent availability of genome information greatly facilitates the fundamental research on chicken. In different organs, gene expression patterns can provide clues to understanding the biological functions. For rapid and accurate quantification of gene expression, quantitative real-time PCR (qPCR) has become one of the most widely used methods. However, the success of qPCR data normalization depends on the use of a suitable reference gene and a single reference gene is not universally suitable for all the experiments. Therefore, reference gene validation is a crucial step for different organ tissues of chicken where suitable reference genes for qPCR analysis in varieties of tissues have not been investigated exhaustively so far. In this study, we have selected 30 Gallus gallus candidate reference genes from NCBI, amplified and studied their expression profiles by qPCR in different organ tissues (breast muscle, thigh muscle, heart, liver, spleen, gizzard, and bursa) of chicken. The result showed that, for breast muscle HSP10 and RPL23, thigh muscle RPL14 and RPL13, liver ALB and HSP70, spleen ALB and GAPDH, heart CYCS and TUBA8B, gizzard RPL5 and 18S rRNA, and bursa EEF1A1 and PGK2 are most stable genes respectively. The results also showed that for different organ tissues, individual or a combination of reference genes should be selected for data normalization. In this study, we have identified and validated 30 reference genes in seven different organ tissues to provide accurate transcript normalization and quantification, which can be useful for gene expression studies in other avian species.

Clarifying expression patterns by renal lesion using transcriptome analysis and vanin-1 as a potential novel biomarker for renal injury in chickens

Bird death is often caused by renal lesions induced by chemicals. The avian kidney has a renal portal system with significant blood flow that is sensitive to many chemicals. However, early avian biomarkers for kidney injury are yet to be identified. This study aimed to identify novel renal biomarkers. Acute kidney injury (AKI) can be divided into acute interstitial nephritis (AIN) and acute tubular necrosis (ATN). A chicken model of kidney damage was created by an injection of diclofenac or cisplatin, which caused either AIN or ATN, respectively. Microarray analysis was performed to profile the gene expression patterns in the chickens with nephropathy. A gene enrichment analysis suggested that the genes related to responses to external stimuli showed expression changes in both AIN and ATN. However, hierarchical clustering analyses suggested that gene expression patterns differed between AIN and ATN, and the number of biomarkers relating to renal damage was low. To identify early biomarkers for nephropathy, we focused on genes that were induced at various levels of renal damage. The gene, vanin-1 (VNN1) was highly induced in the early stages of renal damage. A quantitative real-time PCR analysis supported this finding. These results suggest VNN1 could be a useful early biomarker of kidney injury in avian species.

The expression of the gastrin/cholecystokinin (GAST/CCK) family and their receptors (CCKAR/CCKBR) in the chicken changes in response to quantitative restriction and reveals a functional role of CCK in the crop

Gastrin and cholecystokinin peptides bind a common G-protein coupled receptor, cholecystokinin receptor B (CCKBR) whilst cholecystokinin receptor A (CCKAR) is preferentially bound by CCK. Gastrin and cholecystokinin mediate signalling from the gastrointestinal tract to regulate appetite and digestive function. In this study, expression of the cholecystokinin/gastrin family and distribution of their receptors expression was measured to understand the target organs for the peptides and how expression responds to changes in food intake. We confirmed the restricted expression of gastrin in the antrum and the abundant expression of cholecystokinin in the hypothalamus. The expression of gastrin in the antrum was significantly elevated in broiler breeders when released from feed restriction. CCKBR was most abundant in the hypothalamus and proventriculus. CCKAR was most abundant in the pancreas and crop, more than tenfold greater than the gastrointestinal tract. Cholecystokinin expression in the pancreas increased after removal of food restriction. CCKAR in the gastrointestinal tract peaks around the distal ileum, distal to the peak of cholecystokinin expression. There was virtually no cholecystokinin expression in the caecum but CCKAR expression was high. The CCKAR expression in the crop was unexpected, supporting a role of cholecystokinin in mediating crop emptying which was supported by the observation of in-vitro contraction after cholecystokinin administration. The response to changes in food intake and the expression pattern of the cholecystokinin/gastrin family and their receptors will stimulate and inform new hypotheses on their role in growth in poultry.

Expression of vimentin, TPI and MAT2A in human dermal microvascular endothelial cells during angiogenesis in vitro

Introduction

In vitro assays of angiogenesis face immense problems considering their reproducibility based on the inhomogeneous characters of endothelial cells (ECs). It is necessary to detect influencing factors, which affect the angiogenic potency of ECs.

Objective

This study aimed to analyse expression profiles of vimentin (VIM), triosephosphate isomerase (TPI) and adenosylmethionine synthetase isoform type–2 (MAT2A) during the whole angiogenic cascade in vitro. Furthermore, the impact of knocking down vimentin (VIM) on angiogenesis in vitro was evaluated, while monitoring TPI and MAT2A expression.

Methods

A long–term cultivation and angiogenic stimulation of human dermal microvascular ECs was performed. Cells were characterized via VEGFR–1 and VEGFR–2 expression and a shRNA–mediated knockdown of VIM was performed. The process of angiogenesis in vitro was quantified via morphological staging and mRNA–and protein–levels of all proteins were analysed.

Results

While native cells ran through the angiogenic cascade chronologically, knockdown cells only entered beginning stages of angiogenesis and died eventually. Cell cultures showing a higher VEGFR–1 expression survived exclusively and displayed an upregulation of MAT2A and TPI expression. Native cells highly expressed VIM in early stages, MAT2A mainly in the beginning and TPI during the course of angiogenesis in vitro.

Conclusion

VIM knockdown led to a deceleration of angiogenesis in vitro and knockdown cells displayed expressional changes in TPI and MAT2A. Cell populations with a higher number of stalk cells emerged as being more stable against manipulations and native expression profiles provided an indication of VIM and MAT2A being relevant predominantly in beginning stages and TPI during the whole angiogenic cascade in vitro.

Evolutionarily conserved properties of CLCA proteins 1, 3 and 4, as revealed by phylogenetic and biochemical studies in avian homologues

Species-specific diversities are particular features of mammalian chloride channel regulator, calcium activated (CLCA) genes. In contrast to four complex gene clusters in mammals, only two CLCA genes appear to exist in chickens. CLCA2 is conserved in both, while only the galline CLCA1 (gCLCA1) displays close genetic distance to mammalian clusters 1, 3 and 4. In this study, sequence analyses and biochemical characterizations revealed that gCLCA1 as a putative avian prototype shares common protein domains and processing features with all mammalian CLCA homologues. It has a transmembrane (TM) domain in the carboxy terminal region and its mRNA and protein were detected in the alimentary canal, where the protein was localized in the apical membrane of enterocytes, similar to CLCA4. Both mammals and birds seem to have at least one TM domain containing CLCA protein with complex glycosylation in the apical membrane of enterocytes. However, some characteristic features of mammalian CLCA1 and 3 including entire protein secretion and expression in cell types other than enterocytes seem to be dispensable for chicken. Phylogenetic analyses including twelve bird species revealed that avian CLCA1 and mammalian CLCA3 form clades separate from a major branch containing mammalian CLCA1 and 4. Overall, our data suggest that gCLCA1 and mammalian CLCA clusters 1, 3 and 4 stem from a common ancestor which underwent complex gene diversification in mammals but not in birds.

Baseline Gene Expression Levels in Falkland-Malvinas Island Penguins: Towards a New Monitoring Paradigm

Health diagnostics of wildlife have historically relied on the evaluation of select serum biomarkers and the identification of a contaminant or pathogen burden within specific tissues as an indicator of a level of insult. However, these approaches fail to measure the physiological reaction of the individual to stressors, thus limiting the scope of interpretation. Gene-based health diagnostics provide an opportunity for an alternate, whole-system, or holistic assessment of health, not only in individuals or populations but potentially in ecosystems. Seabirds are among the most threatened marine taxonomic groups in the world, with ~25% of this species currently listed as threatened or considered of special concern; among seabirds, the penguins (Family Spheniscidae) are the most threatened seabird Family. We used gene expression to develop baseline physiological indices for wild penguins in the Falkland-Malvinas Islands, and captive zoo penguins. We identified the almost complete statistical separation of penguin groups (gentoo Detroit Zoo, gentoo Falkland-Malvinas Islands, rockhopper Detroit Zoo, and rockhopper Falkland-Malvinas Islands) based on gene expression profiles. Implementation of long-term longitudinal studies would allow for the assessment of temporal increases or decreases of select transcripts and would facilitate interpretation of the drivers of change.

Meta-transcriptomic analysis of the virome and microbiome of the invasive Indian myna (Acridotheres tristis) in Australia

Invasive species exert a serious impact on native fauna and flora and have become the target of eradication and management efforts worldwide. Invasive avian species can also be important pathogen reservoirs, although their viromes and microbiomes have rarely been studied. As one of the top 100 invasive pest species globally, the expansion of Indian mynas (Acridotheres tristis) into peri-urban and rural environments, in conjunction with increasing free-ranging avian agricultural practices, may increase the risk of microbial pathogens jumping species boundaries. Herein, we used a meta-transcriptomic approach to explore the microbes present in brain, liver and large intestine of 16 invasive Indian myna birds in Sydney, Australia. From this, we discovered seven novel viruses from the families Adenoviridae, Caliciviridae, Flaviviridae, Parvoviridae and Picornaviridae. Interestingly, each of the novel viruses identified shared less than 80% genomic similarity with their closest relatives from other avian species, indicative of a lack of detectable virus transmission between invasive mynas to native or domestic species. Of note, we also identified two coccidian protozoa, Isospora superbusi and Isospora greineri, from the liver and gut tissues of mynas. Overall, these data demonstrate that invasive mynas can harbor a diversity of viruses and other microorganisms such that ongoing pathogen surveillance in this species is warranted.

Identification of reliable reference genes for expression studies in the magnum of laying hens housed in cage and cage-free systems

Stress factors during poultry production can evoke changes in gene transcription and protein synthesis in the hen oviduct and could affect the internal and external egg quality. Studies of relative gene expression require the identification of the most stable reference genes for the quantitative polymerase chain reaction (qPCR) to investigate the reproductive tissues' response in laying hens kept in different production systems. The objective of this study was to determine the most stable reference genes of the magnum tissues of laying hens housed in two different production systems. Hy-Line Brown hens were reared under the same sanitary conditions until 15 weeks of age. Later on, they were transferred into two different production systems, conventional cage (CC) and cage free (CF), up to 82 weeks of age. At 50 and 60 weeks, a total of six hens from each production system were euthanized, and magnum samples were collected. The qPCR was used to determine the RNA transcription level of five reference genes, ACTB, 18S, GAPDH, MSX2 and HMBS. These genes were evaluated for transcript expression in magnum tissues by NormFinder, BestKeeper, geNorm and RefFinder software. The results indicated that the most stable gene in the CF housing system was HMBS in three of the algorithms and in the CC housing system was the 18S, and the best combination of reference genes was HMBS/GAPDH in CF and 18S/HMBS in CC. In conclusion, HMBS, 18S and GAPDH genes could be used together as reference genes for the normalization of the magnum tissues transcript expression of hens in CF and CC housing systems.

Molecular characterization of Gyps africanus (African white-backed vulture) organic anion transporter 1 and 2 expressed in the kidney

Gyps species have been previously shown to be highly sensitive to the toxic effects of diclofenac, when present in their food sources as drug residues following use as a veterinary medicine. Vultures exposed to diclofenac soon become depressed and die with signs of severe visceral gout and renal damage on necropsy. The molecular mechanism behind toxicity and renal excretion of uric acid is still poorly understood. With the clinical pictures suggesting renal uric acid excretion as the target site for toxicity, as a first step the following study was undertaken to determine the uric acid excretory pathways present in the African white-backed vulture (Gyps africanus) (AWB), one of the species susceptible to toxicity. Using transcriptome analysis, immunohistochemistry and functional predictions, we demonstrated that AWB makes use of the organic anion transporter 2 (OAT2) for their uric acid excretion. RT-qPCR analysis subsequently demonstrated relatively similar expression of the OAT2 transporter in the vulture and chicken. Lastly docking analysis, predicted that the non-steroidal drugs induce their toxicity through an allosteric binding.

Epigenetic potential affects immune gene expression in house sparrows

Epigenetic mechanisms may play a central role in mediating phenotypic plasticity, especially during range expansions, when populations face a suite of novel environmental conditions. Individuals may differ in their epigenetic potential (EP; their capacity for epigenetic modifications of gene expression), which may affect their ability to colonize new areas. One form of EP, the number of CpG sites, is higher in introduced house sparrows (Passer domesticus) than in native birds in the promoter region of a microbial surveillance gene, Toll-like Receptor 4 (TLR4), which may allow invading birds to fine-tune their immune responses to unfamiliar parasites. Here, we compared TLR4 gene expression from whole blood, liver and spleen in house sparrows with different EP, first challenging some birds with lipopolysaccharide (LPS), to increase gene expression by simulating a natural infection. We expected that high EP would predict high inducibility and reversibility of TLR4 expression in the blood of birds treated with LPS, but we did not make directional predictions regarding organs, as we could not repeatedly sample these tissues. We found that EP was predictive of TLR4 expression in all tissues. Birds with high EP expressed more TLR4 in the blood than individuals with low EP, regardless of treatment with LPS. Only females with high EP exhibited reversibility in gene expression. Further, the effect of EP varied between sexes and among tissues. Together, these data support EP as one regulator of TLR4 expression.

Transversal gene expression panel to evaluate intestinal health in broiler chickens in different challenging conditions

There is a high interest on gut health in poultry with special focus on consequences of the intestinal diseases, such as coccidiosis and C. perfringens-induced necrotic enteritis (NE). We developed a custom gene expression panel, which could provide a snapshot of gene expression variation under challenging conditions. Ileum gene expression studies were performed through high throughput reverse transcription quantitative real-time polymerase chain reaction. A deep review on the bibliography was done and genes related to intestinal health were selected for barrier function, immune response, oxidation, digestive hormones, nutrient transport, and metabolism. The panel was firstly tested by using a nutritional/Clostridium perfringens model of intestinal barrier failure (induced using commercial reused litter and wheat-based diets without exogenous supplementation of enzymes) and the consistency of results was evaluated by another experiment under a coccidiosis challenge (orally gavaged with a commercial coccidiosis vaccine, 90× vaccine dose). Growth traits and intestinal morphological analysis were performed to check the gut barrier failure occurrence. Results of ileum gene expression showed a higher expression in genes involved in barrier function and nutrient transport in chickens raised in healthy conditions, while genes involved in immune response presented higher expression in C.perfringens-challenged birds. On the other hand, the Eimeria challenge also altered the expression of genes related to barrier function and metabolism, and increased the expression of genes related to immune response and oxidative stress. The panel developed in the current study gives us an overview of genes and pathways involved in broiler response to pathogen challenge. It also allows us to deep into the study of differences in gene expression pattern and magnitude of responses under either a coccidial vaccine or a NE.

Acute restraint stress does not alter corticosteroid receptors or 11β-hydroxysteroid dehydrogenase gene expression at hypothalamic-pituitary-adrenal axis regulatory sites in captive male white-crowned sparrows (Zonotrichia leucophrys gambelii)

Capture-restraint is often used to investigate the acute hypothalamic-pituitary-adrenal axis (HPA) response to stress in wild and captive animals through the production of glucocorticoids. Although this approach is useful for understanding changes in glucocorticoids, it overlooks potential changes in the complex regulatory systems associated with the glucocorticoid response, including genomic receptors, steroid metabolizing enzymes, carrier proteins, and downstream target proteins (e.g. gonadotropin-inhibitory hormone; GnIH). The present study in captive male white-crowned sparrows (Zonotrichia leucophrys) tests the hypothesis that corticosteroid receptors (mineralocorticoid - MR and glucocorticoid - GR), 11β-hydroxysteroid dehydrogenase 1 (11βHSD1) and 2 (11βHSD2), corticosteroid binding globulin (CBG), and GnIH undergo rapid changes in expression to mediate the glucocorticoid response to acute stress. To determine dynamic changes in gene mRNA expression in the hippocampus, hypothalamus, pituitary gland, and liver, birds were sampled within 3 min of entering the room and after 10, 30, and 60 min of capture restraint stress in a cloth bag. Restraint stress handling increased CBG and decreased GnIH mRNA expression in the liver and hypothalamus, respectively. MR, GR, 11βHSD1, and 11βHSD2 mRNA expression in the brain, pituitary gland, and liver did not change. No correlations were found between gene expression and baseline or stress-induced plasma corticosterone levels. No rapid changes of MR, GR, 11βHSD1, and 11βHSD2 mRNA expression during a standardized acute restraint protocol suggests that tissue level sensitivity may remain constant during acute stressors. However, the observed rise in CBG mRNA expression could act to facilitate transport to target tissues or buffer the rise in circulating glucocorticoids. Further studies on tissue specific sensitivity are warranted.

Selection of reference genes for quantitative real-time PCR normalization in European quail tissues

Coturniculture has been standing out as an industrial poultry activity in several countries around the world because of the several adaptive advantages of quails. Research that considers the analysis of gene expression can enhance this activity. This study aimed to analyze the stability of reference genes (RGs) in different tissues of quails (both males and females) for the recommendation of use in gene expression studies by the quantitative reverse transcription-polymerase chain reaction (RT-qPCR). The expression stability of ten RGs (ACTA1, ACTB, B2M, GAPDH, HMBS, SDHA, HPRT1, MRPS27, MRPS30, and RPL5) was analyzed in four tissues (breast muscle, abdominal fat, liver, and intestine), and assessed using the statistical tools geNorm, NormFinder, comparative ΔC_q method, and BestKeeper. The HPRT1 gene was the most stable in all quail tissues tested, followed by MRPS27 and MRPS30 in breast muscle, B2M and RPL5 in abdominal fat, HMBS and B2M in the liver, and RPL5 and HMBS in the intestine. These results may help studies using RT-qPCR assays to assess quail tissues from both sexes because they provide data on the most stable genes, which should be tested as candidate RGs for other experimental conditions.

Determination of novel reference genes for improving gene expression data normalization in selected canine reproductive tissues - a multistudy analysis

Background

Real time RT-PCR (qPCR) is a useful and powerful tool for quantitative measurement of gene expression. The proper choice of internal standards such as reference genes is crucial for correct data evaluation. In female dogs, as in other species, the reproductive tract is continuously undergoing hormonal and cycle stage-dependent morphological changes, which are associated with altered gene expression. However, there have been few attempts published so far targeted to the dog aimed at determining optimal reference genes for the reproductive organs. Most of these approaches relied on genes previously described in other species. Large-scale transcriptome-based experiments are promising tools for defining potential candidate reference genes, but were never considered in this context in canine research.

Results

Here, using available microarray and RNA-seq datasets derived from reproductive organs (corpus luteum, placenta, healthy and diseased uteri) of dogs, we have performed multistudy analysis to identify the most stably expressed genes for expression studies, in each tissue separately and collectively for different tissues. The stability of newly identified reference genes (EIF4H, KDELR2, KDM4A and PTK2) has been determined and ranked relative to previously used reference genes, i.e., GAPDH, β-actin and cyclophillin A/PPIA, using RefFinder and NormFinder algorithms. Finally, expression of selected target genes (luteal IL-1b and MHCII, placental COX2 and VEGFA, and uterine IGF2 and LHR) was re-evaluated and normalized. All proposed candidate reference genes were more stable, ranked higher and introduced less variation than previously used genes.

Conclusions

Based on our analyses, we recommend applying KDM4A and PTK2 for normalization of gene expression in the canine CL and placenta. The inclusion of a third reference gene, EIF4H, is suggested for healthy uteri. With this, the interpretation of qPCR data will be more reliable, allowing better understanding of canine reproductive physiology.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12917-020-02635-6.

Varying Expression of Mu and Kappa Opioid Receptors in Cockatiels (Nymphicus hollandicus) and Domestic Pigeons (Columba livia domestica)

Avian species have varying analgesic responses to opioid drugs. Some of this variability could be due to extrinsic factors such as administration route or dose. However, intrinsic factors such as gene expression or polymorphic differences in opioid receptors may be important components.

Sex-and Region-Dependent Expression of the Autism-Linked ADNP Correlates with Social- and Speech-Related Genes in the Canary Brain

The activity-dependent neuroprotective protein (ADNP) syndrome is an autistic-like disorder, instigated by mutations in ADNP. This syndrome is characterized by developmental delays, impairments in speech, motor function, abnormal hearing, and intellectual disabilities. In the Adnp-haploinsufficient mouse model, many of these impediments are evident, appearing in a sex-dependent manner. In zebra finch songbird (ZF; Taeniopygia guttata), an animal model used for song/language studies, ADNP mRNA most robust expression is observed in the cerebrum of young males, potentially corroborating with male ZF exclusive singing behavior and developed cerebral song system. Herein, we report a similar sex-dependent ADNP expression profile, with the highest expression in the cerebrum (qRT-PCR) in the brain of another songbird, the domesticated canary (Serinus canaria domestica). Additional analyses for the mRNA transcripts of the ADNP regulator, vasoactive intestinal peptide (VIP), sister gene ADNP2, and speech-related Forkhead box protein P2 (FoxP2) revealed multiple sex and brain region-dependent positive correlations between the genes (including ADNP). Parallel transcript expression patterns for FoxP2 and VIP were observed alongside specific FoxP2 increase in males compared with females as well as VIP/ADNP2 correlations. In spatial view, a sexually independent extensive form of expression was found for ADNP in the canary cerebrum (RNA in situ hybridization). The songbird cerebral mesopallium area stood out as a potentially high-expressing ADNP tissue, further strengthening the association of ADNP with sense integration and auditory memory formation, previously implicated in mouse and human.

Identification of Reference Genes in Chicken Intraepithelial Lymphocyte Natural Killer Cells Infected with Very-virulent Infectious Bursal Disease Virus

Due to the limitations in the range of antibodies recognising avian viruses, quantitative real-time PCR (RT-qPCR) is still the most widely used method to evaluate the expression of immunologically related genes in avian viruses. The objective of this study was to identify suitable reference genes for mRNA expression analysis in chicken intraepithelial lymphocyte natural killer (IEL-NK) cells after infection with very-virulent infectious bursal disease virus (vvIBDV). Fifteen potential reference genes were selected based on the references available. The coefficient of variation percentage (CV%) and average count of these 15 genes were determined by NanoString technology for control and infected samples. The M and V values for shortlisted reference genes (ACTB, GAPDH, HMBS, HPRT1, SDHA, TUBB1 and YWHAZ) were calculated using geNorm and NormFinder. GAPDH, YWHAZ and HMBS were the most stably expressed genes. The expression levels of three innate immune response related target genes, CASP8, IL22 and TLR3, agreed in the NanoString and RNA sequencing (RNA-Seq) results using one or two reference genes for normalisation (not HMBS). In conclusion, GAPDH and YWHAZ could be used as reference genes for the normalisation of chicken IEL-NK cell gene responses to infection with vvIBDV.

Standard guidelines for the publication of telomere qPCR results in evolutionary ecology

Telomere length has been used as a proxy of fitness, aging and lifespan in vertebrates. In the last decade, dozens of articles reporting on telomere dynamics in the fields of ecology and evolution have been published for a wide range of taxa. With this growing interest, it is necessary to ensure the accuracy and reproducibility of telomere length measurement techniques. Real-time quantitative PCR (qPCR) is routinely applied to measure relative telomere length. However, this technique is highly sensitive to several methodological variables and the optimization of qPCR telomere assays remains highly variable between studies. Therefore, standardized guidelines are required to enable the optimization of robust protocols, and to help in judging the validity of the presented results. This review provides an overview of preanalytical and analytical factors that can lead to qPCR inconsistencies and biases, including: (a) sample type, collection and storage; (b) DNA extraction, storage and quality; (c) qPCR primers, laboratory reagents, and assay conditions; and (d) data analysis. We propose a minimum level of information for publication of qPCR telomere assays in evolutionary ecology considering the methodological pitfalls and sources of error. This review highlights the complexity of the optimization and validation of qPCR for telomere measurement per se, demonstrating the importance of transparency and clarity of reporting methodological details required for reliable, reproducible and comparable qPCR telomere assays. We encourage efforts to implement standardized protocols that ensure the rigour and quality of telomere dynamics studies.

Assessing the effect of starch digestion characteristics on ileal brake activation in broiler chickens

The objective of this research was to evaluate activation of the ileal brake in broiler chickens using diets containing semi-purified wheat (WS; rapidly and highly digested) and pea (PS; slowly and poorly digested) starch. Diets were formulated to contain six WS:PS ratios (100:0, 80:20, 60:40, 40:60, 20:80, 0:100) and each starch ratio was fed to 236 Ross 308 male broilers housed in 4 litter floor pens. At 28 d of age, the effect of PS concentration was assessed on starch digestion, digestive tract morphology, and digesta pH and short-chain fatty acid (SCFA) concentration. Glucagon-like peptide-1 (GLP-1) and peptide tyrosine-tyrosine (PYY) status were assessed in serum (ELISA) and via gene expression in jejunal and ileal tissue (proglucagon for GLP-1). Data were analyzed using regression analyses, and significance was accepted at P ≤ 0.05. Increasing dietary PS resulted in reduced starch digestibility in the small intestine, but had no effect in the colon. Crop content pH responded quadratically to PS level with an estimated minimum at 55% PS. Total SCFA increased linearly in the crop with PS level, but changed in a quadratic fashion in the ileum (estimated maximum at 62% PS). Ceacal SCFA concentrations were highest for the 80 and 100% PS levels. The relative empty weight (crop, small intestine, colon), length (small intestine) and content (crop jejunum, Ileum) of digestive tract sections increased linearly with increasing PS concentration. Dietary treatment did not affect serum GLP-1 or PYY or small intestine transcript abundance. In conclusion, feeding PS increased the presence of L-cell activators (starch, SCFA) and increased trophic development and content of the digestive tract, suggestive of L-cell activation. However, no direct evidence of ileal brake activation was found by measuring venous blood levels of GLP-1 or PYY or corresponding gene expression in small intestine tissue.

Influence of heat stress on reference genes stability in heart and liver of two chickens genotypes

Introduction

Real-time polymerase chain reaction (RT-qPCR) is an important tool for analyzing gene expression. However, before analyzing the expression of target genes, it is crucial to normalize the reference genes, in order to find the most stable gene to be used as an endogenous control. A gene that remains stable in all samples under different treatments is considered a suitable normalizer. In this sense, we aimed to identify stable reference genes for normalization of target genes in the heart and liver tissues from two genetically divergent groups of chickens (Cobb 500^® commercial line and Peloco backyard chickens) under comfort and acute heat stress environmental conditions. Eight reference genes (ACTB, HPRT1, RPL5, EEF1, MRPS27, MRPS30, TFRC and LDHA) were analyzed for expression stability. The samples were obtained from 24 chickens, 12 from the backyard Peloco and 12 from the Cobb 500® line, exposed to two environmental conditions (comfort and heat stress). Comfort temperature was 23°C and heat stress temperature was 39.5°C for one hour. Subsequently, the animals were euthanized, and heart and liver tissue fragments were collected for RNA extraction and amplification. To determine the stability rate of gene expression, three different statistical algorithms were applied: BestKeeper, geNorm and NormFinder, and to obtain an aggregated stability list, the RankAgregg package of R software was used.

Results

The most stable genes using BestKeeper tool, including the two factors (genetic group and environmental condition), were LDHA, RPL5 and MRPS27 for heart tissue, and TFRC, RPL5 and EEF1 for liver tissue. Applying geNorm algorithm, the best reference genes were RPL5, EEF1 and MRPS30 for heart tissue and LDHA, EEF1 and RPL5 for liver. Using the NormFinder algorithm, the best normalizer genes were EEF1, RPL5 and LDHA in heart, and EEF1, RPL5 and ACTB in liver tissue. In the overall ranking obtained by RankAggreg package, considering the three algorithms, the RPL5, EEF1 and LDHA genes were the most stable for heart tissue, whereas RPL5, EEF1 and ACTB were the most stable for liver tissue.

Conclusion

According to the RankAggreg tool classification based on the three different algorithms (BestKeeper, geNorm and NormFinder), the most stable genes were RPL5, EEF1 and LDHA for heart tissue and RPL5, EEF1 and ACTB for liver tissue of chickens subjected to comfort and acute heat stress environmental conditions. However, the best reference genes may vary depending on the experimental conditions of each study, such as different breeds, environmental stressors, and tissues analyzed. Therefore, the need to perform priori studies to assay the best reference genes at the outset of each study is emphasized.

Avian interspecific differences in VKOR activity and inhibition: Insights from amino acid sequence and mRNA expression ratio of VKORC1 and VKORC1L1

Worldwide use of anticoagulant rodenticides (ARs) for rodents control has frequently led to secondary poisoning of non-target animals, especially raptors. In order to suggest some factors that may help considering the mechanism of the incidents, this study focused on the avian vitamin K 2, 3-epoxide reductase (VKOR) that is the target protein of ARs. We addressed the interspecific differences in VKOR activity and inhibition related to amino acid sequence and mRNA expression of VKORC1 and VKORC1-like1 (VKORC1L1). Poultry have been considered to be more tolerant to ARs than mammals. However, VKOR activity of owls, hawks, falcon and surprisingly, canaries, was lower and inhibited by warfarin more easily than that of chickens and turkeys. The amino acid sequence of VKORC1 and VKORC1L1 implied that the value of Ki for VKOR activity to ARs could depend on the amino acid at position 140 in the TYX warfarin-binding motif in VKORC1, and other amino acid mutations in VKORC1L1. The mRNA expression ratio of VKORC1:VKORC1L1 differed between turkey (8:1) and chicken (2:3) liver. VKORC1L1 has been reported to be resistant to warfarin compared to VKORC1. Hence, both the Ki of specific VKORC1 and VKORC1L1, and the mRNA expression ratio would cause avian interspecific difference of the VKOR inhibition. Our study also suggested the high inhibition of VKOR activities in raptors and surprisingly that in canaries as well. These factors are the most likely to contribute to the high sensitivity to ARs found in raptors.

Tissue-specific expression profiles and positive selection analysis in the tree swallow (Tachycineta bicolor) using a de novo transcriptome assembly

Tree swallows (Tachycineta bicolor) are one of the most commonly studied wild birds in North America. They have advanced numerous research areas, including life history, physiology, and organismal responses to global change; however, transcriptomic resources are scarce. To further advance the utility of this system for biologists across disciplines, we generated a transcriptome for the tree swallow using six tissues (brain, blood, ovary, spleen, liver, and muscle) collected from breeding females. We de novo assembled 207,739 transcripts, which we aligned to 14,717 high confidence protein-coding genes. We then characterized each tissue with regard to its unique genes and processes and applied this transcriptome to two fundamental questions in evolutionary biology and endocrinology. First, we analyzed 3,015 single-copy orthologs and identified 46 genes under positive selection in the tree swallow lineage, including those with putative links to adaptations in this species. Second, we analyzed tissue-specific expression patterns of genes involved in sex steroidogenesis and processing. Enzymes capable of synthesizing these behaviorally relevant hormones were largely limited to the ovary, whereas steroid binding genes were found in nearly all other tissues, highlighting the potential for local regulation of sex steroid-mediated traits. These analyses provide new insights into potential sources of phenotypic variation in a free-living female bird and advance our understanding of fundamental questions in evolutionary and organismal biology.

Differential Song Deficits after Lentivirus-Mediated Knockdown of FoxP1, FoxP2, or FoxP4 in Area X of Juvenile Zebra Finches

Mutations in the transcription factors FOXP1 and FOXP2 are associated with speech impairments. FOXP1 is additionally linked to cognitive deficits, as is FOXP4. These FoxP proteins are highly conserved in vertebrates and expressed in comparable brain regions, including the striatum. In male zebra finches, experimental manipulation of FoxP2 in Area X, a striatal song nucleus essential for vocal production learning, affects song development, adult song production, dendritic spine density, and dopamine-regulated synaptic transmission of striatal neurons. We previously showed that, in the majority of Area X neurons FoxP1, FoxP2, and FoxP4 are coexpressed, can dimerize and multimerize with each other and differentially regulate the expression of target genes. These findings raise the possibility that FoxP1, FoxP2, and FoxP4 (FoxP1/2/4) affect neural function differently and in turn vocal learning. To address this directly, we downregulated FoxP1 or FoxP4 in Area X of juvenile zebra finches and compared the resulting song phenotypes with the previously described inaccurate and incomplete song learning after FoxP2 knockdown. We found that experimental downregulation of FoxP1 and FoxP4 led to impaired song learning with partly similar features as those reported for FoxP2 knockdowns. However, there were also specific differences between the groups, leading us to suggest that specific features of the song are differentially impacted by developmental manipulations of FoxP1/2/4 expression in Area X.SIGNIFICANCE STATEMENT We compared the effects of experimentally reduced expression of the transcription factors FoxP1, FoxP2, and FoxP4 in a striatal song nucleus, Area X, on vocal production learning in juvenile male zebra finches. We show, for the first time, that these temporally and spatially precise manipulations of the three FoxPs affect spectral and temporal song features differentially. This is important because it raises the possibility that the different FoxPs control different aspects of vocal learning through combinatorial gene expression or by acting in different microcircuits within Area X. These results are consistent with the deleterious effects of human FOXP1 and FOXP2 mutations on speech and language and add FOXP4 as a possible candidate gene for vocal disorders.

PPAR expression, muscle size and metabolic rates across the gray catbird's annual cycle are greatest in preparation for fall migration

Phenotypic flexibility across the annual cycle allows birds to adjust to fluctuating ecological demands. Varying energetic demands associated with time of year have been demonstrated to drive metabolic and muscle plasticity in birds, but it remains unclear what molecular mechanisms control this flexibility. We sampled gray catbirds at five stages across their annual cycle: tropical overwintering (January), northward spring (late) migration (early May), breeding (mid June), the fall pre-migratory period (early August) and southward fall (early) migration (end September). Across the catbird's annual cycle, cold-induced metabolic rate (V̇O2summit) was highest during migration and lowest during tropical wintering. Flight muscles exhibited significant hypertrophy and/or hyperplasia during fall migratory periods compared with breeding and the fall pre-migratory period. Changes in heart mass were driven by the tropical wintering stage, when heart mass was lowest. Mitochondrial content of the heart and pectoralis remained constant across the annual cycle as quantified by aerobic enzyme activities (CS, CCO), as did lipid catabolic capacity (HOAD). In the pectoralis, transcription factors PPARα, PPARδ and ERRβ, coactivators PGC-1α and β, and genes encoding proteins associated with fat uptake (FABPpm, Plin3) were unexpectedly upregulated in the tropical wintering stage, whereas those involved in fatty acid oxidation (ATGL, LPL, MCAD) were downregulated, suggesting a preference for fat storage over utilization. Transcription factors and coactivators were synchronously upregulated during pre-migration and fall migration periods in the pectoralis but not the heart, suggesting that these pathways are important in preparation for and during early migration to initiate changes to phenotypes that facilitate long-distance migration.

Acute social isolation alters neurogenomic state in songbird forebrain

Prolonged social isolation has negative effects on brain and behavior in humans and other social organisms, but neural mechanisms leading to these effects are not understood. Here we tested the hypothesis that even brief periods of social isolation can alter gene expression and DNA methylation in higher cognitive centers of the brain, focusing on the auditory/associative forebrain of the highly social zebra finch. Using RNA sequencing, we first identified genes that individually increase or decrease expression after isolation and observed general repression of gene sets annotated for neurotrophin pathways and axonal guidance functions. We then pursued 4 genes of large effect size: EGR1 and BDNF (decreased by isolation) and FKBP5 and UTS2B (increased). By in situ hybridization, each gene responded in different cell subsets, arguing against a single cellular mechanism. To test whether effects were specific to the social component of the isolation experience, we compared gene expression in birds isolated either alone or with a single familiar partner. Partner inclusion ameliorated the effect of solo isolation on EGR1 and BDNF, but not on FKBP5 and UTS2B nor on circulating corticosterone. By bisulfite sequencing analysis of auditory forebrain DNA, isolation caused changes in methylation of a subset of differentially expressed genes, including BDNF. Thus, social isolation has rapid consequences on gene activity in a higher integrative center of the brain, triggering epigenetic mechanisms that may influence processing of ongoing experience.

Armoured exogenous internal control for real-time PCR diagnosis of avian influenza

An exogenous "armoured" PCR internal control (IC) short RNA was analyzed in conjunction with real-time RT-PCR method for diagnosis of avian influenza. The resistance to nucleases and increased physical stability of the IC was ensured using branched polyethyleneimine (PEI) which was in complex with IC-RNA. The option to add the IC directly to pathological material suspensions allows measurement of the nucleic acids extraction efficiency. Stability of armoured RNA-IC during storage and tissue suspension preparation was shown. The advantage of exogenous "armoured" IC was demonstrated in the experiment with AIV genome detection by qPCR in samples from different species of wild birds. The exogenous IC gave reproducible homogeneous Ct values in all tests.

Selection of best reference genes for qRT-PCR analysis of human neural stem cells preconditioned with hypoxia or baicalein-enriched fraction extracted from Oroxylum indicum medicinal plant

Whilst the potential of neural stem cell (NSC)-based treatment is recognized worldwide and seems to offer a promising therapeutic option for stroke treatments, there is currently no full understanding regarding the effects of hypoxic and baicalein-enriched fraction (BEF) preconditioning approaches on the therapeutic potential of these cells for stroke. The potential of preconditioned NSC can be determined based on the expression of several key neuroprotective genes using qRT-PCR technique. However, prior to that, it is imperative and extremely important to carefully select reference gene(s) for accurate qRT-PCR data normalization to avoid error in data interpretation. This study aimed to evaluate the stability of ten candidate reference genes via comprehensive analysis using three algorithms software: geNorm, NormFinder and BestKeeper. Our results revealed that HPRT1 and RPL13A were the most reliable reference genes for BEF-preconditioned NSCs, but ironically, HPRT1 was ranked as the least stable reference gene for hypoxic-preconditioned NSCs. On the other hand, RPLP1 and RPL13A were selected as the most stably expressed pair of reference genes for hypoxic-preconditioned NSCs. In conclusion, this study has pointed out the importance of identifying valid reference genes and has presented the first significant validation on best reference genes recommended for qRT-PCR study involves NSC preconditioned with hypoxia or with BEF extracted from Oroxylum indicum medicinal plant.

Integrative analysis of APOL3 gene CNV for adult cattle stature

Copy number variations (CNVs) have been identified as another important structural variation of genome. In recent years, a large amount of CNVRs have been identified in humans and animals. However, association and dosage effects studies of CNVs are very limited. Apolipoprotein L3 (APOL3) gene plays a central role in modulating gene transcription and is located within a CNVR that encompasses quantitative trait locis (QTLs) for economic traits like meat quality. Herein, we analyzed the CNV polymorphism of APOL3 in 421 individuals from five distinct cattle breeds, and then correlated their genotypes with growth traits. Association analysis revealed that the APOL3 CNV was significantly associated with hip height and cannon circumference of Xianan (XN) cattle (P < .01), and visibly associated with body slanting length and hucklebone width of Pinan (PN) cattle (P < .05). Overall, the data provide evidence for the functional role of APOL3 CNV and a basis for future applications in cattle breeding.

Robust reference gene design and validation for expression studies in the large milkweed bug, Oncopeltus fasciatus, upon cardiac glycoside stress

Despite its history as a developmental and evolutionary model organism, gene expression analysis in the large milkweed bug, Oncopeltus fasciatus, has rarely been explored using quantitative real-time PCR. The strength of this method depends greatly on the endogenous controls used for normalization, which are lacking for the milkweed bug system. Here, to fill in this gap in our knowledge, we validated the stability of a set of ten candidate reference genes identified from the O. fasciatus transcriptome, and did so upon exposure to a dietary toxin, a cardiac glycoside, and across four different exposure periods. To increase robustness against gDNA contaminants, genome resources were used to design intron-bridging primers. A comprehensive stability validation by the Bestkeeper, Normfinder, geNorm and comparative ΔCt methods identified ef1a and tubulin as the most stable genes across treatments and time points, whereas 18S rRNA was the most unstable. However, accounting for the temporal scale indicated that time point confined normalizers might enable higher quantification accuracy for treatment comparison. Overall this study demonstrates: (i) a robust RT-qPCR primer design approach is possible for non-model organisms where genome annotation is often incomplete, and (ii) the importance of detailed reference gene stability exploration in multifactorial experimental designs.

Identification of reliable reference genes for quantitative real-time PCR in ovary and uterus of laying hens under heat stress

The main stage in real-time quantitative PCR is a quantification of gene transcriptomes, in which suitable use of reliable reference genes is critical to normalize accurately. To determine the most stable reference genes in laying hens under heat stress, from a panel of nine typical candidate reference genes, the mRNA transcript of ACTB, HMBS, HPRT1, RPL13, RPL32, 18SrRNA, TBP, TFRC, and YWHAZ was evaluated in the ovary and uterus of both control and heat-stress groups of laying hens. Forty 23-week-old White Leghorn laying hens were housed in two rooms. The control (n = 20) and heat-stress (n = 20) groups were maintained at 21-23 °C and 36-38 °C for 8 weeks respectively. Analysis of this set of genes was done with BestKeeper, geNorm, and NormFinder software programs to find the most stable ones. Candidate reference genes ranked in the uterus of heat-stress and control groups of hens included YWHAZ, HPRT1, HMBS, RPL13, TFRC, ACTB, TBP, RPL32, and 18SrRNA; those in the ovary were YWHAZ, HPRT1, TFRC, HMBS, RPL13, TBP, RPL32, ACTB, and 18SrRNA. The overall results indicated that the most stable genes are YWHAZ, HPRT1, HMBS, RPL13, TFRC, TBP, ACTB, RPL32, and 18SrRNA respectively. In addition, the combination of YWHAZ, HPRT1, and HMBS is suggested as the most stable reference group of genes for more accurate quantitative data normalization in the ovarian and uterine tissues of laying hens under control and heat stress conditions. Lay summary Heat stress influences the expression of many genes in the reproductive tissues of birds. Accurate evaluation of these changes via real-time quantitative PCR depends on the determination of reliable reference genes. In this study, nine candidate housekeeping genes were evaluated, and the most stable were YWHAZ, HPRT1, HMBS, RPL13, TFRC, TBP, ACTB, RPL32, and 18SrRNA.

Identification of Reference Genes for Quantitative Gene Expression Studies in Three Tissues of Japanese Quail

RT-qPCR is the gold standard for candidate gene expression analysis. However, the interpretation of RT-qPCR results depends on the proper use of internal controls, i.e., reference genes. Japanese quail is an agronomic species also used as a laboratory model, but little is known about RT-qPCR reference genes for this species. Thus, we investigated 10 putative reference genes (ACTB, GAPDH, PGK1, RPS7, RPS8, RPL19, RPL32, SDHA, TBP and YWHAZ) in three different female and male quail tissues (liver, brain and pectoral muscle). Gene expression stability was evaluated with three different algorithms: geNorm, NormFinder and BestKeeper. For each tissue, a suitable set of reference genes was defined and validated by a differential analysis of gene expression between females and males (CCNH in brain and RPL19 in pectoral muscle). Collectively, our study led to the identification of suitable reference genes in liver, brain and pectoral muscle for Japanese quail, along with recommendations for the identification of reference gene sets for this species.

Survey of Reticuloendotheliosis Virus in Wild Turkeys (Meleagris gallopavo) in Texas, USA

Reticuloendotheliosis virus (REV) is an immunosuppressive and sometimes oncogenic avian retrovirus that establishes lifelong infection in a wide range of avian species. REV-infected wild birds roaming near at-risk captive flocks, such as is the case for the highly endangered Attwater's Prairie Chicken (APC; Tympanuchus cupido attwateri), could act as a reservoir for viral transmission. In wild birds, prevalence rates of REV are low and appearance of associated disease is uncommon. During 2016-17, nearly half of all captive adult APC mortality at Fossil Rim Wildlife Center captive breeding facility in Glen Rose, Texas, US was attributed to REV infection. The unusually high REV prevalence rate prompted us to survey for this virus in wild galliforms throughout the region. From 2016-17, 393 blood samples collected from two subspecies of Wild Turkeys (Meleagris gallopavo) were tested for REV proviral DNA through amplification of the viral 3' long terminal repeat and segments of the viral pol gene. In REV-affected counties, 5% (5/98) of native Rio Grande Wild Turkeys (Meleagris gallopavo intermedia) were identified as REV-positive. In addition, we detected REV in one of 62 Eastern Wild Turkeys (Meleagris gallopavo silvestris) that had been imported during conservation efforts. To better determine protective measures, continued surveillance, including collection and genetic analysis of REV-infected samples, is necessary to identify sources of REV outbreaks in captive APC flocks.

In vitro study on role of σB protein in avian reovirus pathogenesis

Avian reoviruses, members of Orthoreovirus genus was known to cause diseases like tenosynovitis, runting-stunting syndrome in chickens. Among eight structural proteins, the proteins of S-class are mainly associated with viral arthritis but the significance of σB protein in arthritis is not established till date. In this infection pathological condition together with infection of joints often leads to arthritis because joints consists of cartilage which forms lubricating surface between two bones, and has limited metabolic, replicative and repair capacity. To establish the role of σB protein in arthritis, an in-vitro microarray study was conducted consisting four groups viz. virus infected and control; pDsRed-Express-N1-σB and empty pDs-Red transfected, CEF cells. With cut-off value as FC ≥2, p value <0.05, 6709 and 4026 numbers of DEGs in virus and σB, respectively were identified. The Ingenuity Pathway Analysis gave an idea about the involvement of σB protein in "osteoarthritis pathway", which was activated with z-score with 3.151. The pathway "Role of IL-17A in arthritis pathway" was also enriched with -log (p-value) 1.64. Among total 122 genes involved in osteoarthritis pathway, 28 upregulated and 11 downregulated DEGs were common to both virus and σB treated cells. Moreover, 14 upregulated and 7 downregulated were unique in σB transfected cells. Using qRT-PCR for IL-1B, BMP2, SMAD1, SPP1 genes, the microarray data was validated. We concluded that during ARV infection σB protein, if not fully partially leads to molecular alteration of various genes of host orchestrating the different molecular pattern in joints, leading to tenosynovitis syndrome.