Isolation and characterization of sexual dimorphism genes expressed in chicken embryonic gonads

Characterization of heterozygosity-rich regions in Italian and worldwide goat breeds

Heterozygosity-rich regions (HRR) are genomic regions of high heterozygosity, which may harbor loci related to key functional traits such as immune response, survival rate, fertility, and other fitness traits. This study considered 30 Italian and 19 worldwide goat breeds genotyped with the Illumina GoatSNP50k BeadChip. The aim of the work was to study inter-breed relationships and HRR patterns using Sliding Window (SW) and Consecutive Runs (CR) detection methods. Genetic relationships highlighted a clear separation between non-European and European breeds, as well as the north-south geographic cline within the latter. The Pearson correlation coefficients between the descriptive HRR parameters obtained with the SW and CR methods were higher than 0.9. A total of 166 HRR islands were detected. CHI1, CHI11, CHI12 and CHI18 were the chromosomes harboring the highest number of HRR islands. The genes annotated in the islands were linked to various factors such as productive, reproductive, immune, and environmental adaptation mechanisms. Notably, the Montecristo feral goat showed the highest number of HRR islands despite the high level of inbreeding, underlining potential balancing selection events characterizing its evolutionary history. Identifying a species-specific HRR pattern could provide a clearer view of the mechanisms regulating the genome modelling following anthropogenic selection combined with environmental interaction.

Selection of reference genes for RT-qPCR analysis in developing chicken embryonic ovary

BACKGROUND

Normalization of the expression profiling of target genes, in a tissue-specific manner and under different experimental conditions, requires stably expressed gene(s) to be used as internal reference(s). However, to study the molecular regulation of oocyte meiosis initiation during ovary development in chicken embryos, stable reference gene(s) still need to be compared and confirmed.

METHODS AND RESULTS

Six candidate genes previously used as internal references for the chicken embryo (Actb, Cvh, Dazl, Eef1a, Gapdh and Rpl15) were chosen, and their expression profiles in left ovaries dissected at five chicken embryonic days (E12.5, E15.5, E17.5, E18.5 and E20.5) were evaluated, respectively. Separately, GeNorm, NormFinder, BestKeeper and Comparative ΔCt methods were used to assess the stability of candidate reference genes, and all results were combined to give the final rank by RefFinder. All methods identified that Eef1a and Rpl15 were the two most stable internal reference genes, whereas Cvh is the most unstable one. Moreover, expression levels of three marker genes for chicken oocyte meiosis entry (Stra8, Scp3 and Dmc1) were normalized, based on Eef1a, Rpl15, or their combinations, respectively.

CONCLUSION

Our findings provide the most suitable internal reference genes (Eef1a and Rpl15), to investigate further molecular regulation of ovary development and oocyte meiosis initiation in chicken embryos.

Keep Garfagnina alive. An integrated study on patterns of homozygosity, genomic inbreeding, admixture and breed traceability of the Italian Garfagnina goat breed

The objective of this study was to investigate the genetic diversity of the Garfagnina (GRF) goat, a breed that currently risks extinction. For this purpose, 48 goats were genotyped with the Illumina CaprineSNP50 BeadChip and analyzed together with 214 goats belonging to 9 other Italian breeds (~25 goats/breed), whose genotypes were available from the AdaptMap project [Argentata (ARG), Bionda dell'Adamello (BIO), Ciociara Grigia (CCG), Di Teramo (DIT), Garganica (GAR), Girgentana (GGT), Orobica (ORO), Valdostana (VAL) and Valpassiria (VSS)]. Comparative analyses were conducted on i) runs of homozygosity (ROH), ii) admixture ancestries and iii) the accuracy of breed traceability via discriminant analysis on principal components (DAPC) based on cross-validation. ROH analyses was used to assess the genetic diversity of GRF, while admixture and DAPC to evaluate its relationship to the other breeds. For GRF, common ROH (more than 45% in GRF samples) was detected on CHR 12 at, roughly 50.25-50.94Mbp (ARS1 assembly), which spans the CENPJ (centromere protein) and IL17D (interleukin 17D) genes. The same area of common ROH was also present in DIT, while a broader region (~49.25-51.94Mbp) was shared among the ARG, CCG, and GGT. Admixture analysis revealed a small region of common ancestry from GRF shared by BIO, VSS, ARG and CCG breeds. The DAPC model yielded 100% assignment success for GRF. Overall, our results support the identification of GRF as a distinct native Italian goat breed. This work can contribute to planning conservation programmes to save GRF from extinction and will improve the understanding of the socio-agro-economic factors related with the farming of GRF.

Transcriptomic and Epigenetic Preservation of Genetic Sex Identity in Estrogen-feminized Male Chicken Embryonic Gonads

Whereas in ovo exposure of genetically male (ZZ) chicken embryos to exogenous estrogens temporarily feminizes gonads at the time of hatching, the morphologically ovarian ZZ-gonads (FemZZs for feminized ZZ gonads) are masculinized back to testes within 1 year. To identify the feminization-resistant "memory" of genetic male sex, FemZZs showing varying degrees of feminization were subjected to transcriptomic, DNA methylome, and immunofluorescence analyses. Protein-coding genes were classified based on their relative mRNA expression across normal ZZ-testes, genetically female (ZW) ovaries, and FemZZs. We identified a group of 25 genes that were strongly expressed in both ZZ-testes and FemZZs but dramatically suppressed in ZW-ovaries. Interestingly, 84% (21/25) of these feminization-resistant testicular marker genes, including the DMRT1 master masculinizing gene, were located in chromosome Z. Expression of representative marker genes of germline cells (eg, DAZL or DDX4/VASA) was stronger in FemZZs than normal ZZ-testes or ZW-ovaries. We also identified 231 repetitive sequences (RSs) that were strongly expressed in both ZZ-testes and FemZZs, but these RSs were not enriched in chromosome Z. Although 94% (165/176) of RSs exclusively expressed in ZW-ovaries were located in chromosome W, no feminization-inducible RS was detected in FemZZs. DNA methylome analysis distinguished FemZZs from normal ZZ- and ZW-gonads. Immunofluorescence analysis of FemZZ gonads revealed expression of DMRT1 protein in medullary SOX9+ somatic cells and apparent germline cell populations in both medulla and cortex. Taken together, our study provides evidence that both somatic and germline cell populations in morphologically feminized FemZZs maintain significant transcriptomic and epigenetic memories of genetic sex.

Cardiorespiratory and thermal responses to hypercapnia in chickens exposed to CO2 during embryonic development

The concentration of CO₂ in the environment surrounding the embryo impacts development and may also influence the cardiorespiratory responses after hatching. Therefore, we aimed to evaluate the cardiorespiratory and thermal responses to hypercapnia in chicks that were exposed to CO₂ during embryonic development, i.e., incubation. Embryos were incubated without and with a gradual increase in CO₂ concentration up to 1 % during the first ten days of incubation. Ten-day-old chicks (males and females) were again acutely exposed to hypercapnia (7 % CO₂), or to room air (normocapnia) and pulmonary ventilation, arterial pH and blood gases, arterial blood pressure and heart rate, body temperature (Tb) and oxygen consumption (V⋅O₂) were measured. Compared to control animals, male chicks incubated with 1 % CO₂ presented an attenuated ventilatory response to hypercapnia (P < 0.05), whereas no difference was found in the hypercapnic ventilatory response in both female chick groups (0 % vs 1 % CO₂ incubation). Hypercapnia induced bradycardia in all groups (P < 0.001). The CO₂ exposure during incubation did not alter the cardiovascular responses to hypercapnia in post-hatch animals. There were no significant effects of incubation treatment (0 % vs 1 % CO₂) or sex in the mean arterial pressure, Tb, and V⋅O₂ of animals in normocapnia and hypercapnia. As for the V⋅_E/V⋅O₂, hypercapnia caused an increase in both groups (P < 0.05), regardless of incubation treatment. In conclusion, among cardiorespiratory and metabolic variables, the ventilatory response to hypercapnia can be attenuated by pre-exposure to 1 % CO₂ during embryonic development, especially in male chicks up to 10 days.

Genome-wide patterns of homozygosity provide clues about the population history and adaptation of goats

Background

Patterns of homozygosity can be influenced by several factors, such as demography, recombination, and selection. Using the goat SNP50 BeadChip, we genotyped 3171 goats belonging to 117 populations with a worldwide distribution. Our objectives were to characterize the number and length of runs of homozygosity (ROH) and to detect ROH hotspots in order to gain new insights into the consequences of neutral and selection processes on the genome-wide homozygosity patterns of goats.

Results

The proportion of the goat genome covered by ROH is, in general, less than 15% with an inverse relationship between ROH length and frequency i.e. short ROH (< 3 Mb) are the most frequent ones. Our data also indicate that ~ 60% of the breeds display low F_ROH coefficients (< 0.10), while ~ 30 and ~ 10% of the goat populations show moderate (0.10 < F_ROH < 0.20) or high (> 0.20) F_ROH values. For populations from Asia, the average number of ROH is smaller and their coverage is lower in goats from the Near East than in goats from Central Asia, which is consistent with the role of the Fertile Crescent as the primary centre of goat domestication. We also observed that local breeds with small population sizes tend to have a larger fraction of the genome covered by ROH compared to breeds with tens or hundreds of thousands of individuals. Five regions on three goat chromosomes i.e. 11, 12 and 18, contain ROH hotspots that overlap with signatures of selection.

Conclusions

Patterns of homozygosity (average number of ROH of 77 and genome coverage of 248 Mb; F_ROH < 0.15) are similar in goats from different geographic areas. The increased homozygosity in local breeds is the consequence of their small population size and geographic isolation as well as of founder effects and recent inbreeding. The existence of three ROH hotspots that co-localize with signatures of selection demonstrates that selection has also played an important role in increasing the homozygosity of specific regions in the goat genome. Finally, most of the goat breeds analysed in this work display low levels of homozygosity, which is favourable for their genetic management and viability.

Electronic supplementary material

The online version of this article (10.1186/s12711-018-0424-8) contains supplementary material, which is available to authorized users.

Temperature alterations during embryogenesis have a sex-dependent influence on growth properties and muscle metabolism of day-old chicks and 35-day-old broilers

Broiler eggs were either incubated at 37.8°C during the whole incubation period (control), or at higher (38.8°C, group H) and lower temperatures (36.8°C, group L) from embryonic day (ED) 7 up to ED 10 (ED 7 to 10) or from ED 10 up to ED 13 (ED 10 to 13). Before and after this temperature treatment the eggs were incubated at 37.8°C. The day-old chicks were weighted, sexed and fed up to day 35. On days 1 and 35 samples were taken from the breast and leg muscles for analyzing of the mitochondrial respiratory activity (MRA) and from the breast muscles for analysis of the cross-sectional areas (CSA) and the glycogen phosphorylase (GP), phosphofructokinase (PFK), lactate dehydrogenase (LDH), citrate synthase (CS) and cytochrome oxidase (COX) activities. Statistical analysis showed that treatment (control, group H, group L), sex and their interaction, but not the treatment period (ED 7 to 10; ED 10 to 13), significantly influenced the results. Group H chicks had lower (P⩽0.05) body and heart weights but higher (P⩽0.05) liver weights, CSA values, leg MRA as well as PFK, LDH, CS, GP and COX activities compared with the group L chicks. The results of the control chicks differ (P⩽0.05) from those of the group H (body, heart weight, COX), the group L chicks (liver weight, PFK, LDH, CS, GP) or the birds of both other groups (CSA). The group H broiler had higher (P⩽0.05) body and leg weights as well as LDH, CS, COX and GP activities than the group L broilers. The BWs and the LDH and GP results of the control broiler differ (P⩽0.05) from those of both other groups or from the results of the group H (CS) and group L broiler (COX). Female broilers had lower (P⩽0.05) body, breast and leg weights, but higher (P⩽0.05) CSA, LDH, CS and GP activities than the male animals. Analysis of treatment×sex interaction showed that group H hens had higher (P⩽0.05) body and breast weights, LDH and GP activities compared with the group L hens, whereas in the male broiler no effect of the interaction could be found, except for the lower (P⩽0.05) CSA values in the group H than group L cocks. The treatment effects are probably due to altered embryonic activity and related molecular mechanisms. The sex-related differences in the broiler indicate that these alterations already occur in the embryos and chicks, but become significant with the sexual dimorphism after hatch.

Sexual Dimorphism in the Early Embryogenesis in Zebra Finches

Sex-specific gene expression before the onset of gonadogensis has been documented in embryos of mammals and chickens. In several mammalian species, differences in gene expression are accompanied by faster growth of pre-implantation male embryos. Here we asked whether avian embryos before gonadal differentiation are also sex-dimorphic in size and what genes regulate their growth. We used captive zebra finches (Taeniopygia guttata) whose freshly laid eggs were artificially incubated for 36–40 hours. Analyses controlling for the exact time of incubation of 81 embryos revealed that males were larger than females in terms of Hamburger and Hamilton stage and number of somites. Expression of 15 genes involved in cell cycle regulation, growth, metabolic activity, steroidogenic pathway and stress modulation were measured using RT-PCR in 5 male and 5 female embryos incubated for exactly 36 h. We found that in the presence of equal levels of the growth hormone itself, the faster growth of male embryos is most likely achieved by the overexpression of the growth hormone receptor gene and three other genes responsible for cell cycle regulation and metabolism, all of them located on the Z chromosome. Autosomal genes did not show sex-specific expression, except for the steroidogenic factor 1 which was expressed only in female embryos. To our knowledge this is the first report of sexual size dimorphism before gonadogenesis in birds. The finding suggests that faster growth of early male embryos is conserved through the mammalian and bird phyla, irrespective of their differential sex chromosome systems.

Gene expression profile during testicular development in patients with SRY-negative 46,XX testicular disorder of sex development

OBJECTIVE

To elucidate alternative pathways in testicular development, we attempted to clarify the genetic characteristics of SRY-negative XX testes.

MATERIALS AND METHODS

We previously reported 5 cases of SRY-negative 46,XX testicular disorders of sex development and demonstrated that coordinated expression of genes such as SOX9, SOX3, and DAX1 was associated with testicular development. We performed a case-control study between the aforementioned boy with 46,XX testicular disorders of sex development and an age-matched patient with hydrocele testis (46,XY). During their consecutive surgeries, testicular biopsy specimens were obtained. Genes with differential expression compared with XY testis were identified using polymerase chain reaction (PCR)-based subtractive hybridization and sequencing. For validation of differential gene expression, real-time RT-PCR was performed using gene-specific primers. The distribution of candidate proteins in the testicular tissue was clarified by immunohistochemistry in human and rodent specimens. Moreover, in vitro inhibitory assays were performed.

RESULTS

We identified 13 upregulated and 7 downregulated genes in XX testis. Among the candidate genes, we focused on ROCK1 (Rho-associated, coiled-coil protein kinase 1) in the upregulated gene group, because high expression in XX testis was validated by real-time RT-PCR. ROCK1 protein was detected in germ cells, Leydig cells, and Sertoli cells by immunohistochemistry. Moreover, the addition of specific ROCK1 inhibitor to Sertoli cells decreased SOX9 gene expression.

CONCLUSION

On the basis of in vitro inhibitory assay, it is suggested that ROCK1 phosphorylates and activates SOX9 in Sertoli cells. Testes formation might be initiated by an alternative signaling pathway attributed to ROCK1, not SRY, activation in XX testes.