Principal component approach in variance component estimation for international sire evaluation

Background

The dairy cattle breeding industry is a highly globalized business, which needs internationally comparable and reliable breeding values of sires. The international Bull Evaluation Service, Interbull, was established in 1983 to respond to this need. Currently, Interbull performs multiple-trait across country evaluations (MACE) for several traits and breeds in dairy cattle and provides international breeding values to its member countries. Estimating parameters for MACE is challenging since the structure of datasets and conventional use of multiple-trait models easily result in over-parameterized genetic covariance matrices. The number of parameters to be estimated can be reduced by taking into account only the leading principal components of the traits considered. For MACE, this is readily implemented in a random regression model.

Methods

This article compares two principal component approaches to estimate variance components for MACE using real datasets. The methods tested were a REML approach that directly estimates the genetic principal components (direct PC) and the so-called bottom-up REML approach (bottom-up PC), in which traits are sequentially added to the analysis and the statistically significant genetic principal components are retained. Furthermore, this article evaluates the utility of the bottom-up PC approach to determine the appropriate rank of the (co)variance matrix.

Results

Our study demonstrates the usefulness of both approaches and shows that they can be applied to large multi-country models considering all concerned countries simultaneously. These strategies can thus replace the current practice of estimating the covariance components required through a series of analyses involving selected subsets of traits. Our results support the importance of using the appropriate rank in the genetic (co)variance matrix. Using too low a rank resulted in biased parameter estimates, whereas too high a rank did not result in bias, but increased standard errors of the estimates and notably the computing time.

Conclusions

In terms of estimation''s accuracy, both principal component approaches performed equally well and permitted the use of more parsimonious models through random regression MACE. The advantage of the bottom-up PC approach is that it does not need any previous knowledge on the rank. However, with a predetermined rank, the direct PC approach needs less computing time than the bottom-up PC.     

Vitamin D Binding Protein Genotype Is Associated with Serum 25-Hydroxyvitamin D and PTH Concentrations,as Well as Bone Health in Children and Adolescents in Finland

Vitamin D binding protein (DBP)/group-specific component (Gc), correlates positively with serum vitamin D metabolites, and phenotype influences serum 25-hydroxyvitamin D (S-25(OH)D) concentration. The protein isoform has been associated with decreased bone mineral density (BMD) and increased fracture risk. We examined the role of GC genotypes in S-25(OH)D status and BMD in 231 Finnish children and adolescents aged 7−19 yr. BMD was measured with DXA from lumbar spine (LS), total hip, and whole body, and for 175 subjects, radial volumetric BMD was measured with pQCT. Background characteristic and total dietary intakes of vitamin D and calcium were collected. The concentrations of 25(OH)D, parathyroid hormone (PTH), calcium and other markers of calcium homeostasis were determined from blood and urine. Genotyping was based on single-nucleotide polymorphism (rs4588) in the GC gene. The genotype distribution was: GC 1/1 68%, GC 1/2 26% and GC 2/2 6%. A significant difference emerged in 25(OH)D and PTH concentrations between the genotypes, (p = 0.001 and 0.028 respectively, ANCOVA). There was also a linear trend in: Gc 2/2 had the lowest 25(OH)D and PTH concentrations (p = 0.025 and 0.012, respectively). Total hip bone mineral content was associated with GC genotype (BMC) (p = 0.05, ANCOVA) in boys. In regression analysis, after adjusting for relevant covariates, GC genotype was associated with LS BMC and strength and strain index (SSI) Z-score in both genders, and LS BMD in boys. In conclusion, the present study demonstrates the association between GC genotypes and S-25(OH)D and PTH concentrations. The results show the influence of DBP genetic variation on bone mass accrual in adolescence.     

Predictors for pathologically confirmed aortitis after resection of the ascending aorta: A 12-year Danish nationwide population-based cross-sectional study

Introduction  

Assessing the prevalence of, and predictors for, pathologically-confirmed inflammation of the aorta in Denmark, using a nationwide population-based study design.     

Alarm pheromones and chemical communication in nymphs of the tropical bed bug Cimex hemipterus (Hemiptera: Cimicidae)

The recent resurge of bed bug infestations (Cimex spp.; Cimicidae) and their resistance to commonly used pesticides calls for alternative methods of control. Pheromones play an important role in environmentally sustainable methods for the management of many pest insects and may therefore be applicable for the control of bed bugs. The tropical bed bug, Cimex hemipterus, is a temporary ectoparasite on humans and causes severe discomfort. Compared to the common bed bug, Cimex lectularius, little is known about the chemical signalling and pheromone-based behaviour of the tropical species. Here, we show that the antennal morphology and volatile emission of C. hemipterus closely resembles those of C. lectularius and we test their behavioural responses to conspecific odour emissions. Two major volatiles are emitted by male, female and nymph C. hemipterus under stress, (E)-2-hexenal and (E)-2-octenal. Notably, nymph emissions show contrasting ratios of these compounds to adults and are further characterized by the addition of 4-oxo-(E)-2-hexenal and 4-oxo-(E)-2-octenal. The discovery of this nymph pheromone in C. hemipterus is potentially the cause of a repellent effect observed in the bio-tests, where nymph odours induce a significantly stronger repellent reaction in conspecifics than adult odours. Our results suggest that pheromone-based pest control methods developed for C. lectularius could be applicable to C. hemipterus, with the unique nymph blend showing promising practical properties.     

Salmonella transiently reside in luminal neutrophils in the inflamed gut

Background

Enteric pathogens need to grow efficiently in the gut lumen in order to cause disease and ensure transmission. The interior of the gut forms a complex environment comprising the mucosal surface area and the inner gut lumen with epithelial cell debris and food particles. Recruitment of neutrophils to the intestinal lumen is a hallmark of non-typhoidal Salmonella enterica infections in humans. Here, we analyzed the interaction of gut luminal neutrophils with S. enterica serovar Typhimurium (S. Tm) in a mouse colitis model.

Results

Upon S. Tm^wt infection, neutrophils transmigrate across the mucosa into the intestinal lumen. We detected a majority of pathogens associated with luminal neutrophils 20 hours after infection. Neutrophils are viable and actively engulf S. Tm, as demonstrated by live microscopy. Using S. Tm mutant strains defective in tissue invasion we show that pathogens are mostly taken up in the gut lumen at the epithelial barrier by luminal neutrophils. In these luminal neutrophils, S. Tm induces expression of genes typically required for its intracellular lifestyle such as siderophore production iroBCDE and the Salmonella pathogenicity island 2 encoded type three secretion system (TTSS-2). This shows that S. Tm at least transiently survives and responds to engulfment by gut luminal neutrophils. Gentamicin protection experiments suggest that the life-span of luminal neutrophils is limited and that S. Tm is subsequently released into the gut lumen. This “fast cycling” through the intracellular compartment of gut luminal neutrophils would explain the high fraction of TTSS-2 and iroBCDE expressing intra- and extracellular bacteria in the lumen of the infected gut.

Conclusion

In conclusion, live neutrophils recruited during acute S. Tm colitis engulf pathogens in the gut lumen and may thus actively engage in shaping the environment of pathogens and commensals in the inflamed gut.     

Genome-Wide Analysis of DNA Methylation Differences in Muscle and Fat from Monozygotic Twins Discordant for Type 2 Diabetes

Background

Monozygotic twins discordant for type 2 diabetes constitute an ideal model to study environmental contributions to type 2 diabetic traits. We aimed to examine whether global DNA methylation differences exist in major glucose metabolic tissues from these twins.

Methodology/Principal Findings

Skeletal muscle (n = 11 pairs) and subcutaneous adipose tissue (n = 5 pairs) biopsies were collected from 53–80 year-old monozygotic twin pairs discordant for type 2 diabetes. DNA methylation was measured by microarrays at 26,850 cytosine-guanine dinucleotide (CpG) sites in the promoters of 14,279 genes. Bisulfite sequencing was applied to validate array data and to quantify methylation of intergenic repetitive DNA sequences. The overall intra-pair variation in DNA methylation was large in repetitive (LINE1, D4Z4 and NBL2) regions compared to gene promoters (standard deviation of intra-pair differences: 10% points vs. 4% points, P<0.001). Increased variation of LINE1 sequence methylation was associated with more phenotypic dissimilarity measured as body mass index (r = 0.77, P = 0.007) and 2-hour plasma glucose (r = 0.66, P = 0.03) whereas the variation in promoter methylation did not associate with phenotypic differences. Validated methylation changes were identified in the promoters of known type 2 diabetes-related genes, including PPARGC1A in muscle (13.9±6.2% vs. 9.0±4.5%, P = 0.03) and HNF4A in adipose tissue (75.2±3.8% vs. 70.5±3.7%, P<0.001) which had increased methylation in type 2 diabetic individuals. A hypothesis-free genome-wide exploration of differential methylation without correction for multiple testing identified 789 and 1,458 CpG sites in skeletal muscle and adipose tissue, respectively. These methylation changes only reached some percentage points, and few sites passed correction for multiple testing.

Conclusions/Significance

Our study suggests that likely acquired DNA methylation changes in skeletal muscle or adipose tissue gene promoters are quantitatively small between type 2 diabetic and non-diabetic twins. The importance of methylation changes in candidate genes such as PPARGC1A and HNF4A should be examined further by replication in larger samples.     

Kainate induces various domain closures in AMPA and kainate receptors

Ionotropic glutamate receptors are key players in fast excitatory synaptic transmission within the central nervous system. These receptors have been divided into three subfamilies: the N-methyl-d-aspartic acid (NMDA), 2-amino-3-(3-hydroxy-5-methyl-4-isoxazolyl)propionic acid (AMPA) and kainate receptors. Kainate has previously been crystallized with the ligand binding domain (LBD) of AMPA receptors (GluA2 and GluA4) and kainate receptors (GluK1 and GluK2). Here, we report the structures of the kainate receptor GluK3 LBD in complex with kainate and GluK1 LBD in complex with kainate in the absence of glycerol. Kainate introduces a conformational change in GluK3 LBD comparable to that of GluK2, but different from the conformational changes induced in GluA2 and GluK1. Compared to their domain closures in a glutamate bound state, GluA2 and GluK1 become more open and kainate induces a domain closure of 60% and 62%, respectively, relative to glutamate (100%). In GluK2 and GluK3 with kainate, the domain closure is 88% and 83%, respectively. In previously determined structures of GluK1 LBD in complex with kainate, glycerol is present in the binding site where it bridges interlobe residues and thus, might contribute to the large domain opening. However, the structure of GluK1 LBD with kainate in the absence of glycerol confirms that the observed domain closure is not an artifact of crystallization conditions. Comparison of the LBD structures with glutamate and kainate reveals that contacts are lost upon binding of kainate in the three kainate receptors, which is in contrast to the AMPA receptors where similar contacts are seen. It was revealed by patch clamp electrophysiology studies that kainate is a partial agonist at GluK1 with 36% efficacy compared to glutamate, which is in between the published efficacies of kainate at GluK2 and AMPA receptors. The ranking of efficacies seems to correlate with LBD domain closures.     

Ccdc94 Protects Cells from Ionizing Radiation by Inhibiting the Expression of p53

DNA double-strand breaks (DSBs) represent one of the most deleterious forms of DNA damage to a cell. In cancer therapy, induction of cell death by DNA DSBs by ionizing radiation (IR) and certain chemotherapies is thought to mediate the successful elimination of cancer cells. However, cancer cells often evolve to evade the cytotoxicity induced by DNA DSBs, thereby forming the basis for treatment resistance. As such, a better understanding of the DSB DNA damage response (DSB-DDR) pathway will facilitate the design of more effective strategies to overcome chemo- and radioresistance. To identify novel mechanisms that protect cells from the cytotoxic effects of DNA DSBs, we performed a forward genetic screen in zebrafish for recessive mutations that enhance the IR-induced apoptotic response. Here, we describe radiosensitizing mutation 7 (rs7), which causes a severe sensitivity of zebrafish embryonic neurons to IR-induced apoptosis and is required for the proper development of the central nervous system. The rs7 mutation disrupts the coding sequence of ccdc94, a highly conserved gene that has no previous links to the DSB-DDR pathway. We demonstrate that Ccdc94 is a functional member of the Prp19 complex and that genetic knockdown of core members of this complex causes increased sensitivity to IR-induced apoptosis. We further show that Ccdc94 and the Prp19 complex protect cells from IR-induced apoptosis by repressing the expression of p53 mRNA. In summary, we have identified a new gene regulating a dosage-sensitive response to DNA DSBs during embryonic development. Future studies in human cancer cells will determine whether pharmacological inactivation of CCDC94 reduces the threshold of the cancer cell apoptotic response.     

Structural and functional characterization of disulfide isoforms of the human IgG2 subclass

In the accompanying report ( Wypych, J., Li, M., Guo, A., Zhang, Z., Martinez, T., Allen, M. J., Fodor, S., Kelner, D. N., Flynn, G. C., Liu, Y. D., Bondarenko, P. V., Ricci, M. S., Dillon, T. M., and Balland, A. (2008) J. Biol. Chem. 283, 16194-16205 ), we have identified that the human IgG2 subclass exists as an ensemble of distinct isoforms, designated IgG2-A, -B, and -A/B, which differ by the disulfide connectivity at the hinge region. In this report, we studied the structural and functional properties of the IgG2 disulfide isoforms and compared them to IgG1. Human monoclonal IgG1 and IgG2 antibodies were designed with identical antigen binding regions, specific to interleukin-1 cell surface receptor type 1. In vitro biological activity measurements showed an increased activity of the IgG1 relative to the IgG2 in blocking interleukin-1beta ligand from binding to the receptor, suggesting that some of the IgG2 isoforms had lower activity. Under reduction-oxidation conditions, the IgG2 disulfide isoforms converted to IgG2-A when 1 m guanidine was used, whereas IgG2-B was enriched in the absence of guanidine. The relative potency of the antibodies in cell-based assays was: IgG1 > IgG2-A > IgG2 > IgG2-B. This difference correlated with an increased hydrodynamic radius of IgG2-A relative to IgG2-B, as shown by biophysical characterization. The enrichment of disulfide isoforms and activity studies were extended to additional IgG2 monoclonal antibodies with various antigen targets. All IgG2 antibodies displayed the same disulfide conversion, but only a subset showed activity differences between IgG2-A and IgG2-B. Additionally, the distribution of isoforms was influenced by the light chain type, with IgG2lambda composed mostly of IgG2-A. Based on crystal structure analysis, we propose that IgG2 disulfide exchange is caused by the close proximity of several cysteine residues at the hinge and the reactivity of tandem cysteines within the hinge. Furthermore, the IgG2 isoforms were shown to interconvert in whole blood or a "blood-like" environment, thereby suggesting that the in vivo activity of human IgG2 may be dependent on the distribution of isoforms.