Exploring the Genetic Correlation Between Growth and Immunity Based on Summary Statistics of Genome-Wide Association Studies

Unraveling the genomic background of resilience based on variability in milk yield and milk production levels in North American Holstein cattle through genome-wide association study and Mendelian randomization analyses

Breeding more resilient animals will benefit the dairy cattle industry in the long term, especially as global climate changes become more severe. Previous studies have reported genetic parameters for various milk yield-based resilience indicators, but the underlying genomic background of these traits remain unknown. In this study, we conducted GWAS of 62,029 SNPs with 4 milk yield-based resilience indicators, including the weighted occurrence frequency (wfPert) and accumulated milk losses (dPert) of milk yield perturbations, and log-transformed variance (LnVar) and lag-1 autocorrelation (rauto) of daily yield residuals. These variables were previously derived from 5.6 million daily milk yield records from 21,350 lactations (parities 1-3) of 11,787 North American Holstein cows. The average daily milk yield (ADMY) throughout lactation was also included to compare the shared genetic background of resilience indicators with milk yield. The differential genetic background of these indicators was first revealed by the significant genomic regions identified and significantly enriched biological pathways of positional candidate genes, which confirmed the genetic difference among resilience indicators. Interestingly, the functional analyses of candidate genes suggested that the regulation of intestinal homeostasis is most likely affecting resilience derived based on variability in milk yield. Based on Mendelian randomization analyses of multiple instrumental SNPs, we further found an unfavorable causal association of ADMY with LnVar. In conclusion, the resilience indicators evaluated are genetically different traits, and there are causal associations of milk yield with some of the resilience indicators evaluated. In addition to providing biological insights into the molecular regulation mechanisms of resilience derived based on variability in milk yield, this study also indicates the need for developing selection indexes combining multiple indicator traits and taking into account their genetic relationship for breeding more resilient dairy cattle.

Pleiotropic architectures of porcine immune and growth trait pairs revealed by a self-product-based transcriptome method

Pleiotropy is an important biological phenomenon with complicated genetic architectures for multiple traits. To date, pleiotropy has been mainly identified by multi-trait genome-wide association studies, but this method has its disadvantages, and new developments for pleiotropy detection methods are needed. Here we define a novel metric, self-product, to measure individual-level co-variation of two traits, and develop a novel self-product-based transcriptome method to detect pleiotropic genes (PGs). Our method was tested using four immune-growth trait pairs and four immune-immune trait pairs in pigs. Comparative transcriptome analyses identified hundreds of candidate PGs related to eight trait pairs from two tails of self-product distribution. Gene Ontology enrichment analysis indicated that most of identified PGs were involved in immune- or growth-related biological processes. We established PG interaction networks to exhibit core genes shared by eight trait pairs, of which CCL5 and IL-10 genes were the hub genes. Genetic association analyses showed that SmaI-polymorphisms of CCL5 and IL-10 genes had significant associations with phenotypic co-variations of multiple trait pairs, indicating that the variants in pleiotropic genes were also pleiotropic variants. Taken together, the validity of our proposed method was preliminarily verified, and our findings provide new insights into the genetic basis of pleiotropic architectures of immune and growth trait pairs in pigs.

Direct and maternal genetic effects for preinflection point growth traits and humoral immunity in quail

Early growth traits in quails are considered as the growth performances before the inflection point which are genetically different from body weights (BW) at later stages. Moreover, in addition to growth performance, humoral immunity is moderately heritable and is considered in some breeding programs. However, estimating the direct genetic, particularly the maternal genetic correlations between growth and immunity in quail, are not studied sufficiently, which were the aims of the present study. The quails' BW were recorded at hatch (BW0) to 25 d of age with a 5-d interval and body weight gains (BWG) were measured as average growth performance of the birds in a 5-d period. Antibody titer against Newcastle disease virus (IgN) was measured through the hemagglutination inhibition (HI) test. For titration of anti-SRBC antibodies (IgY and IgM), a hemagglutination microtiter assay was used. In general, growth records in 4,181 birds and humoral immune responses in 1,023 birds were assigned to the study. The genetic parameters were estimated by single-trait analysis via Gibb's sampling. After finding the best model for each trait, multi-trait analysis was done to estimate the direct and maternal genetic correlations. Direct heritabilities (h²) were estimated to be moderate for BW (0.481-0.551) and BWG (0.524-0.557), while h² for immune responses were low (0.035-0.079). Maternal environmental effect (c²) was only significant for BW0, BW5, and BWG0-5. Maternal heritabilities (m²) for BW and BWG were all lower than corresponding h², ranging from 0.072 (BW25) to 0.098 (BW0). The m² for IgN (0.098) was more than 2.5 times greater than h² (0.040) for this trait. Direct (r_a) and maternal (r_m) genetic correlations between IgN-BW, IgY-BW, and IgY-BWG were negative, while ra and rm for IgM-BW, IgN-BWG, and IgM-BWG were positive. The r_a between humoral immune responses were low to moderate and rm was significant only for IgY-IgM (0.339). Given positive genetic correlations in BWG-IgN and BWG-IgM as well as positive genetic correlations between both IgN and IgM with IgY, it is suggested that including the BWG in the breeding programs would directly result in the improvement of the birds' growth performance. It would also contribute indirectly to the improvement of the birds' humoral immune responses.

Genetic correlation and causal relationships between cardio-metabolic traits and lung function impairment

BACKGROUND

Associations of low lung function with features of poor cardio-metabolic health have been reported. It is, however, unclear whether these co-morbidities reflect causal associations, shared genetic heritability or are confounded by environmental factors.

METHODS

We performed three analyses: (1) cardio-metabolic health to lung function association tests in Northern Finland Birth cohort 1966, (2) cross-trait linkage disequilibrium score regression (LDSC) to compare genetic backgrounds and (3) Mendelian randomisation (MR) analysis to assess the causal effect of cardio-metabolic traits and disease on lung function, and vice versa (bidirectional MR). Genetic associations were obtained from the UK Biobank data or published large-scale genome-wide association studies (N > 82,000).

RESULTS

We observed a negative genetic correlation between lung function and cardio-metabolic traits and diseases. In Mendelian Randomisation analysis (MR), we found associations between type 2 diabetes (T2D) instruments and forced vital capacity (FVC) as well as FEV1/FVC. Body mass index (BMI) instruments were associated to all lung function traits and C-reactive protein (CRP) instruments to FVC. These genetic associations provide evidence for a causal effect of cardio-metabolic traits on lung function. Multivariable MR suggested independence of these causal effects from other tested cardio-metabolic traits and diseases. Analysis of lung function specific SNPs revealed a potential causal effect of FEV1/FVC on blood pressure.

CONCLUSIONS

The present study overcomes many limitations of observational studies by using Mendelian Randomisation. We provide evidence for an independent causal effect of T2D, CRP and BMI on lung function with some of the T2D effect on lung function being attributed to inflammatory mechanisms. Furthermore, this analysis suggests a potential causal effect of FEV1/FVC on blood pressure. Our detailed analysis of the interplay between cardio-metabolic traits and impaired lung function provides the opportunity to improve the quality of existing intervention strategies.

Cutting-edge genetics in obsessive-compulsive disorder

This article reviews recent advances in the genetics of obsessive-compulsive disorder (OCD). We cover work on the following: genome-wide association studies, whole-exome sequencing studies, copy number variation studies, gene expression, polygenic risk scores, gene–environment interaction, experimental animal systems, human cell models, imaging genetics, pharmacogenetics, and studies of endophenotypes. Findings from this work underscore the notion that the genetic architecture of OCD is highly complex and shared with other neuropsychiatric disorders. Also, the latest evidence points to the participation of gene networks involved in synaptic transmission, neurodevelopment, and the immune and inflammatory systems in this disorder. We conclude by highlighting that further study of the genetic architecture of OCD, a great part of which remains to be elucidated, could benefit the development of diagnostic and therapeutic approaches based on the biological basis of the disorder. Studies to date revealed that OCD is not a simple homogeneous entity, but rather that the underlying biological pathways are variable and heterogenous. We can expect that translation from bench to bedside, through continuous effort and collaborative work, will ultimately transform our understanding of what causes OCD and thus how best to treat it.

The schizophrenia risk locus in SLC39A8 alters brain metal transport and plasma glycosylation

A common missense variant in SLC39A8 is convincingly associated with schizophrenia and several additional phenotypes. Homozygous loss-of-function mutations in SLC39A8 result in undetectable serum manganese (Mn) and a Congenital Disorder of Glycosylation (CDG) due to the exquisite sensitivity of glycosyltransferases to Mn concentration. Here, we identified several Mn-related changes in human carriers of the common SLC39A8 missense allele. Analysis of structural brain MRI scans showed a dose-dependent change in the ratio of T2w to T1w signal in several regions. Comprehensive trace element analysis confirmed a specific reduction of only serum Mn, and plasma protein N-glycome profiling revealed reduced complexity and branching. N-glycome profiling from two individuals with SLC39A8-CDG showed similar but more severe alterations in branching that improved with Mn supplementation, suggesting that the common variant exists on a spectrum of hypofunction with potential for reversibility. Characterizing the functional impact of this variant will enhance our understanding of schizophrenia pathogenesis and identify novel therapeutic targets and biomarkers.

Searching for parent-of-origin effects on cardiometabolic traits in imprinted genomic regions

Cardiometabolic traits pose a major global public health burden. Large-scale genome-wide association studies (GWAS) have identified multiple loci accounting for up to 30% of the genetic variance in complex traits such as cardiometabolic traits. However, the contribution of parent-of-origin effects (POEs) to complex traits has been largely ignored in GWAS. Family-based studies enable the assessment of POEs in genetic association analyses. We investigated POEs on a range of complex traits in 3 family-based studies. The discovery phase was carried out in large pedigrees from the Kibbutzim Family Study (n = 901 individuals) and in 872 parent-offspring trios from the Jerusalem Perinatal Study. Focusing on imprinted genomic regions, we examined parent-specific associations with 12 complex traits (i.e., body-size, blood pressure, lipids), mostly cardiometabolic risk traits. Forty five of the 11,967 SNPs initially found to have POE were evaluated for replication (p value < 1 × 10-4) in Framingham Heart Study families (max n = 8000 individuals). Three common variants yielded evidence of POE in the meta-analysis. Two variants, located on chr6 in the HLA region, showed a paternal effect on height (rs1042136: βpaternal = -0.023, p value = 1.5 × 10-8 and rs1431403: βpaternal = -0.011, p value = 5.4 × 10-6). The corresponding maternally-derived effects were statistically nonsignificant. The variant rs9332053, located on chr13 in RCBTB2 gene, demonstrated a maternal effect on hip circumference (βmaternal = -4.24, p value = 9.6 × 10-6; βpaternal = 1.29, p value = 0.23). These findings provide evidence for the utility of incorporating POEs into association studies of cardiometabolic traits, especially anthropometric traits. The study highlights the benefits of using family-based data for deciphering the genetic architecture of complex traits.

Changes in the Composition of the Gut Microbiota and the Blood Transcriptome in Preterm Infants at Less than 29 Weeks Gestation Diagnosed with Bronchopulmonary Dysplasia

Bronchopulmonary dysplasia (BPD) is a common chronic lung condition in preterm infants that results in abnormal lung development and leads to considerable morbidity and mortality, making BPD one of the most common complications of preterm birth. We employed RNA sequencing and 16S rRNA gene sequencing to profile gene expression in blood and the composition of the fecal microbiota in infants born at <29 weeks gestational age and diagnosed with BPD in comparison to those of preterm infants that were not diagnosed with BPD. 16S rRNA gene sequencing, performed longitudinally on 255 fecal samples collected from 50 infants in the first months of life, identified significant differences in the relative levels of abundance of Klebsiella, Salmonella, Escherichia/Shigella, and Bifidobacterium in the BPD infants in a manner that was birth mode dependent. Transcriptome sequencing (RNA-Seq) analysis revealed that more than 400 genes were upregulated in infants with BPD. Genes upregulated in BPD infants were significantly enriched for functions related to red blood cell development and oxygen transport, while several immune-related pathways were downregulated. We also identified a gene expression signature consistent with an enrichment of immunosuppressive CD71+ early erythroid cells in infants with BPD. Intriguingly, genes that were correlated in their expression with the relative abundances of specific taxa in the microbiota were significantly enriched for roles in the immune system, suggesting that changes in the microbiota might influence immune gene expression systemically.IMPORTANCE Bronchopulmonary dysplasia (BPD) is a serious inflammatory condition of the lung and is the most common complication associated with preterm birth. A large body of evidence now suggests that the gut microbiota can influence immunity and inflammation systemically; however, the role of the gut microbiota in BPD has not been evaluated to date. Here, we report that there are significant differences in the gut microbiota of infants born at <29 weeks gestation and subsequently diagnosed with BPD, which are particularly pronounced when infants are stratified by birth mode. We also show that erythroid and immune gene expression levels are significantly altered in BPD infants. Interestingly, we identified an association between the composition of the microbiota and immune gene expression in blood in early life. Together, these findings suggest that the composition of the microbiota may influence the risk of developing BPD and, more generally, may shape systemic immune gene expression.

Cell Metabolism Control Through O-GlcNAcylation of STAT5: A Full or Empty Fuel Tank Makes a Big Difference for Cancer Cell Growth and Survival

O-GlcNAcylation is a post-translational modification that influences tyrosine phosphorylation in healthy and malignant cells. O-GlcNAc is a product of the hexosamine biosynthetic pathway, a side pathway of glucose metabolism. It is essential for cell survival and proper gene regulation, mirroring the metabolic status of a cell. STAT3 and STAT5 proteins are essential transcription factors that can act in a mutational context-dependent manner as oncogenes or tumor suppressors. They regulate gene expression for vital processes such as cell differentiation, survival, or growth, and are also critically involved in metabolic control. The role of STAT3/5 proteins in metabolic processes is partly independent of their transcriptional regulatory role, but is still poorly understood. Interestingly, STAT3 and STAT5 are modified by O-GlcNAc in response to the metabolic status of the cell. Here, we discuss and summarize evidence of O-GlcNAcylation-regulating STAT function, focusing in particular on hyperactive STAT5A transplant studies in the hematopoietic system. We emphasize that a single O-GlcNAc modification is essential to promote development of neoplastic cell growth through enhancing STAT5A tyrosine phosphorylation. Inhibition of O-GlcNAcylation of STAT5A on threonine 92 lowers tyrosine phosphorylation of oncogenic STAT5A and ablates malignant transformation. We conclude on strategies for new therapeutic options to block O-GlcNAcylation in combination with tyrosine kinase inhibitors to target neoplastic cancer cell growth and survival.