Sex modifies genetic effects on residual variance in urinary calcium excretion in rat (Rattus norvegicus)

Genetic architecture of trait variance in craniofacial morphology

The genetic architecture of trait variance has long been of interest in genetics and evolution. One of the earliest attempts to understand this architecture was presented in Lerner's Genetic Homeostasis (1954). Lerner proposed that heterozygotes should be better able to tolerate environmental perturbations because of functional differences between the alleles at a given locus, with each allele optimal for slightly different environments. This greater robustness to environmental variance, he argued, would result in smaller trait variance for heterozygotes. The evidence for Lerner's hypothesis has been inconclusive. To address this question using modern genomic methods, we mapped loci associated with differences in trait variance (vQTL) on 1,101 individuals from the F34 of an advanced intercross between LG/J and SM/J mice. We also mapped epistatic interactions for these vQTL in order to understand the influence of epistasis for the architecture of trait variance. We did not find evidence supporting Lerner's hypothesis, that heterozygotes tend to have smaller trait variances than homozygotes. We further show that the effects of most mapped loci on trait variance are produced by epistasis affecting trait means and that those epistatic effects account for about a half of the differences in genotypic-specific trait variances. Finally, we propose a model where the different interactions between the additive and dominance effects of the vQTL and their epistatic partners can explain Lerner's original observations but can also be extended to include other conditions where heterozygotes are not the least variable genotype.

Network pharmacology and experimental validation to elucidate the pharmacological mechanisms of Bushen Huashi decoction against kidney stones

INTRODUCTION

Kidney stone disease (KS) is a complicated disease with an increasing global incidence. It was shown that Bushen Huashi decoction (BSHS) is a classic Chinese medicine formula that has therapeutic benefits for patients with KS. However, its pharmacological profile and mechanism of action are yet to be elucidated.

METHODS

The present study used a network pharmacology approach to characterize the mechanism by which BSHS affects KS. Compounds were retrieved from corresponding databases, and active compounds were selected based on their oral bioavailability (≥30) and drug-likeness index (≥0.18). BSHS potential proteins were obtained from the Traditional Chinese Medicine Systems Pharmacology (TCMSP) database, whereas KS potential genes were obtained from GeneCards and OMIM, TTD, and DisGeNET. Gene ontology and pathway enrichment analysis were used to determine potential pathways associated with genes. The ingredients of BSHS extract were identified by the ultra-high-performance liquid chromatography coupled with quadrupole orbitrap mass spectrometry (UHPLC-Q/Orbitrap MS). The network pharmacology analyses predicted the potential underlying action mechanisms of BSHS on KS, which were further validated experimentally in the rat model of calcium oxalate kidney stones.

RESULTS

Our study found that BSHS reduced renal crystal deposition and improved renal function in ethylene glycol(EG)+ammonium chloride(AC)-induced rats, and also reversed oxidative stress levels and inhibited renal tubular epithelial cell apoptosis in rats. BSHS upregulated protein and mRNA expression of E2, ESR1, ESR2, BCL2, NRF2, and HO-1 in EG+AC-induced rat kidney while downregulating BAX protein and mRNA expression, consistent with the network pharmacology results.

DISCUSSION

This study provides evidence that BSHS plays a critical role in anti-KS via regulation of E2/ESR1/2, NRF2/HO-1, and BCL2/BAX signaling pathways, indicating that BSHS is a candidate herbal drug for further investigation in treating KS.

Evaluation of the phenotypic and genomic background of variability based on litter size of Large White pigs

BACKGROUND

The genetic background of trait variability has captured the interest of ecologists and animal breeders because the genes that control it could be involved in buffering various environmental effects. Phenotypic variability of a given trait can be assessed by studying the heterogeneity of the residual variance, and the quantitative trait loci (QTL) that are involved in the control of this variability are described as variance QTL (vQTL). This study focuses on litter size (total number born, TNB) and its variability in a Large White pig population. The variability of TNB was evaluated either using a simple method, i.e. analysis of the log-transformed variance of residuals (LnVar), or the more complex double hierarchical generalized linear model (DHGLM). We also performed a single-SNP (single nucleotide polymorphism) genome-wide association study (GWAS). To our knowledge, this is only the second study that reports vQTL for litter size in pigs and the first one that shows GWAS results when using two methods to evaluate variability of TNB: LnVar and DHGLM.

RESULTS

Based on LnVar, three candidate vQTL regions were detected, on Sus scrofa chromosomes (SSC) 1, 7, and 18, which comprised 18 SNPs. Based on the DHGLM, three candidate vQTL regions were detected, i.e. two on SSC7 and one on SSC11, which comprised 32 SNPs. Only one candidate vQTL region overlapped between the two methods, on SSC7, which also contained the most significant SNP. Within this vQTL region, two candidate genes were identified, ADGRF1, which is involved in neurodevelopment of the brain, and ADGRF5, which is involved in the function of the respiratory system and in vascularization. The correlation between estimated breeding values based on the two methods was 0.86. Three-fold cross-validation indicated that DHGLM yielded EBV that were much more accurate and had better prediction of missing observations than LnVar.

CONCLUSIONS

The results indicated that the LnVar and DHGLM methods resulted in genetically different traits. Based on their validation, we recommend the use of DHGLM over the simpler method of log-transformed variance of residuals. These conclusions can be useful for future studies on the evaluation of the variability of any trait in any species.

Variance heterogeneity genome-wide mapping for cadmium in bread wheat reveals novel genomic loci and epistatic interactions

Genome-wide association mapping identifies quantitative trait loci (QTL) that influence the mean differences between the marker genotypes for a given trait. While most loci influence the mean value of a trait, certain loci, known as variance heterogeneity QTL (vQTL) determine the variability of the trait instead of the mean trait value (mQTL). In the present study, we performed a variance heterogeneity genome-wide association study (vGWAS) for grain cadmium (Cd) concentration in bread wheat. We used double generalized linear model and hierarchical generalized linear model to identify vQTL associated with grain Cd. We identified novel vQTL regions on chromosomes 2A and 2B that contribute to the Cd variation and loci that affect both mean and variance heterogeneity (mvQTL) on chromosome 5A. In addition, our results demonstrated the presence of epistatic interactions between vQTL and mvQTL, which could explain variance heterogeneity. Overall, we provide novel insights into the genetic architecture of grain Cd concentration and report the first application of vGWAS in wheat. Moreover, our findings indicated that epistasis is an important mechanism underlying natural variation for grain Cd concentration.

The Evolution of Variance Control

Genetically identical individuals can be phenotypically variable, even in constant environmental conditions. The ubiquity of this phenomenon, known as 'intra-genotypic variability', is increasingly evident and the relevant mechanistic underpinnings are beginning to be understood. In parallel, theory has delineated a number of formal expectations for contexts in which such a feature would be adaptive. Here, we review empirical evidence across biological systems and theoretical expectations, including nonlinear averaging and bet hedging. We synthesize existing results to illustrate the dependence of selection outcomes both on trait characteristics, features of environmental variability, and species' demographic context. We conclude by discussing ways to bridge the gap between empirical evidence of intra-genotypic variability, studies demonstrating its genetic component, and evidence that it is adaptive.

'Fat's chances': Loci for phenotypic dispersion in plasma leptin in mouse models of diabetes mellitus

Background

Leptin, a critical mediator of feeding, metabolism and diabetes, is expressed on an incidental basis according to satiety. The genetic regulation of leptin should similarly be episodic.

Methodology

Data from three mouse cohorts hosted by the Jackson Laboratory– 402 (174F, 228M) F₂ Dilute Brown non-Agouti (DBA/2)×DU6i intercrosses, 142 Non Obese Diabetic (NOD/ShiLtJ×(NOD/ShiLtJ×129S1/SvImJ.H2^g7) N₂ backcross females, and 204 male Nonobese Nondiabetic (NON)×New Zealand Obese (NZO/HlLtJ) reciprocal backcrosses–were used to test for loci associated with absolute residuals in plasma leptin and arcsin-transformed percent fat (‘phenotypic dispersion’; PD_pLep and PD_AFP). Individual data from 1,780 mice from 43 inbred strains was also used to estimate genetic variances and covariances for dispersion in each trait.

Principal findings

Several loci for PD_pLep were detected, including possibly syntenic Chr 17 loci, but there was only a single position on Chr 6 for PD_AFP. Coding SNP in genes linked to the consensus Chr 17 PD_pLep locus occurred in immunological and cancer genes, genes linked to diabetes and energy regulation, post-transcriptional processors and vomeronasal variants. There was evidence of intersexual differences in the genetic architecture of PD_pLep. PD_pLep had moderate heritability (hs2=0.29) and PD_AFP low heritability (hs2=0.12); dispersion in these traits was highly genetically correlated r = 0.8).

Conclusions

Greater genetic variance for dispersion in plasma leptin, a physiological trait, may reflect its more ephemeral nature compared to body fat, an accrued progressive character. Genetic effects on incidental phenotypes such as leptin might be effectively characterized with randomization-detection methodologies in addition to classical approaches, helping identify incipient or borderline cases or providing new therapeutic targets.

Genetic Effects on Dispersion in Urinary Albumin and Creatinine in Three House Mouse (Mus musculus) Cohorts

Conventionally, quantitative genetics concerns the heredity of trait means, but there is growing evidence for the existence of architectures in which certain alleles cause random variance in phenotype, termed ‘phenotypic dispersion’ (PD) or ‘variance QTL’ (vQTL), including in physiological traits like disease signs. However, the structure of this phenomenon is still poorly known. PD for urinary albumin (PD_UAlb) and creatinine (PD_UCrea) was mapped using curated data from two nearly genetically identical F₂ mouse (Mus musculus) cohorts (383 male F₂ C57BL/6J×A/J (97 SNP) and 207 male F₂ C57BL/6J×A/J ApoE knockout mice (144 SNP)) and a related mapping cohort (340 male F₂ DBA/2J×C57BL/6J (83 SNP, 8 microsatellites)). PD_UAlb was associated with markers in regions of Chr 1 (5-64 megabases (MB); 141-158 MB), 3 (∼113 MB), 8 (37-68 MB), 14 (92-117 MB) and 17 (14-24 MB) with several positions and quantitative architectures in common between the two C57BL/6J×A/J cohorts, most of which had a negative dominant construction. One locus for PD_UCrea was detected on Chr 19 (57 MB) in the C57BL/6J×A/J ApoE^−/− cohort. The large number of negative dominant loci for albuminuria dispersion relative to conventional quantitative trait loci suggests that the development of albuminuria may be largely genetically dynamic and that randomization in this development is detrimental.

A sibling method for identifying vQTLs

The propensity of a trait to vary within a population may have evolutionary, ecological, or clinical significance. In the present study we deploy sibling models to offer a novel and unbiased way to ascertain loci associated with the extent to which phenotypes vary (variance-controlling quantitative trait loci, or vQTLs). Previous methods for vQTL-mapping either exclude genetically related individuals or treat genetic relatedness among individuals as a complicating factor addressed by adjusting estimates for non-independence in phenotypes. The present method uses genetic relatedness as a tool to obtain unbiased estimates of variance effects rather than as a nuisance. The family-based approach, which utilizes random variation between siblings in minor allele counts at a locus, also allows controls for parental genotype, mean effects, and non-linear (dominance) effects that may spuriously appear to generate variation. Simulations show that the approach performs equally well as two existing methods (squared Z-score and DGLM) in controlling type I error rates when there is no unobserved confounding, and performs significantly better than these methods in the presence of small degrees of confounding. Using height and BMI as empirical applications, we investigate SNPs that alter within-family variation in height and BMI, as well as pathways that appear to be enriched. One significant SNP for BMI variability, in the MAST4 gene, replicated. Pathway analysis revealed one gene set, encoding members of several signaling pathways related to gap junction function, which appears significantly enriched for associations with within-family height variation in both datasets (while not enriched in analysis of mean levels). We recommend approximating laboratory random assignment of genotype using family data and more careful attention to the possible conflation of mean and variance effects.

Genome-wide association study reveals novel loci for litter size and its variability in a Large White pig population

Background

In many traits, not only individual trait levels are under genetic control, but also the variation around that level. In other words, genotypes do not only differ in mean, but also in (residual) variation around the genotypic mean. New statistical methods facilitate gaining knowledge on the genetic architecture of complex traits such as phenotypic variability. Here we study litter size (total number born) and its variation in a Large White pig population using a Double Hierarchical Generalized Linear model, and perform a genome-wide association study using a Bayesian method.

Results

In total, 10 significant single nucleotide polymorphisms (SNPs) were detected for total number born (TNB) and 9 SNPs for variability of TNB (varTNB). Those SNPs explained 0.83 % of genetic variance in TNB and 1.44 % in varTNB. The most significant SNP for TNB was detected on Sus scrofa chromosome (SSC) 11. A possible candidate gene for TNB is ENOX1, which is involved in cell growth and survival. On SSC7, two possible candidate genes for varTNB are located. The first gene is coding a swine heat shock protein 90 (HSPCB = Hsp90), which is a well-studied gene stabilizing morphological traits in Drosophila and Arabidopsis. The second gene is VEGFA, which is activated in angiogenesis and vasculogenesis in the fetus. Furthermore, the genetic correlation between additive genetic effects on TNB and on its variation was 0.49. This indicates that the current selection to increase TNB will also increase the varTNB.

Conclusions

To the best of our knowledge, this is the first study reporting SNPs associated with variation of a trait in pigs. Detected genomic regions associated with varTNB can be used in genomic selection to decrease varTNB, which is highly desirable to avoid very small or very large litters in pigs. However, the percentage of variance explained by those regions was small. The SNPs detected in this study can be used as indication for regions in the Sus scrofa genome involved in maintaining low variability of litter size, but further studies are needed to identify the causative loci.

Unifying genetic canalization, genetic constraint, and genotype-by-environment interaction: QTL by genomic background by environment interaction of flowering time in Boechera stricta

Natural populations exhibit substantial variation in quantitative traits. A quantitative trait is typically defined by its mean and variance, and to date most genetic mapping studies focus on loci altering trait means but not (co)variances. For single traits, the control of trait variance across genetic backgrounds is referred to as genetic canalization. With multiple traits, the genetic covariance among different traits in the same environment indicates the magnitude of potential genetic constraint, while genotype-by-environment interaction (GxE) concerns the same trait across different environments. While some have suggested that these three attributes of quantitative traits are different views of similar concepts, it is not yet clear, however, whether they have the same underlying genetic mechanism. Here, we detect quantitative trait loci (QTL) influencing the (co)variance of phenological traits in six distinct environments in Boechera stricta, a close relative of Arabidopsis. We identified nFT as the QTL altering the magnitude of phenological trait canalization, genetic constraint, and GxE. Both the magnitude and direction of nFT's canalization effects depend on the environment, and to our knowledge, this reversibility of canalization across environments has not been reported previously. nFT's effects on trait covariance structure (genetic constraint and GxE) likely result from the variable and reversible canalization effects across different traits and environments, which can be explained by the interaction among nFT, genomic backgrounds, and environmental stimuli. This view is supported by experiments demonstrating significant nFT by genomic background epistatic interactions affecting phenological traits and expression of the candidate gene, FT. In contrast to the well-known canalization gene Hsp90, the case of nFT may exemplify an alternative mechanism: Our results suggest that (at least in traits with major signal integrators such as flowering time) genetic canalization, genetic constraint, and GxE may have related genetic mechanisms resulting from interactions among major QTL, genomic backgrounds, and environments.

1,25(OH)₂D₃ induces a mineralization defect and loss of bone mineral density in genetic hypercalciuric stone-forming rats

Genetic hypercalciuric stone-forming (GHS) rats, bred to maximize urine (u) calcium (Ca) excretion, demonstrate increased intestinal Ca absorption, increased bone Ca resorption, and reduced renal Ca reabsorption, all leading to elevated uCa compared to the parental Sprague-Dawley (SD) rats. GHS rats have increased numbers of vitamin D receptors (VDRs) at each site, with normal levels of 1,25(OH)₂D₃ (1,25D), suggesting their VDR is undersaturated with 1,25D. We have shown that 1,25D induces a greater increase in uCa in GHS than SD rats. To examine the effect of the increased VDR on the osseous response to 1,25D, we fed GHS and SD rats an ample Ca diet and injected either 1,25D [low dose (LD) 12.5 or high dose (HD) 25 ng/100 g body weight/day] or vehicle (veh) daily for 16 days. Femoral areal bone mineral density (aBMD, by DEXA) was decreased in GHS+LD and GHS+HD relative to GHS+veh, while there was no effect on SD. Vertebral aBMD was lower in GHS compared to SD and further decreased in GHS+HD. Both femoral and L6 vertebral volumetric BMD (by μCT) were lower in GHS and further reduced by HD. Histomorphometry indicated a decreased osteoclast number in GHS+HD compared to GHS+veh or SD+HD. In tibiae, GHS+HD trabecular thickness and number increased, with a 12-fold increase in osteoid volume but only a threefold increase in bone volume. Bone formation rate was decreased in GHS+HD relative to GHS+veh, confirming the mineralization defect. The loss of BMD and the mineralization defect in GHS rats contribute to increased hypercalciuria; if these effects persist, they would result in decreased bone strength, making these bones more fracture-prone. The enhanced effect of 1,25D in GHS rats indicates that the increased VDRs are biologically active.

Persistence of 1,25D-induced hypercalciuria in alendronate-treated genetic hypercalciuric stone-forming rats fed a low-calcium diet

Genetic hypercalciuric stone-forming (GHS) rats demonstrate increased intestinal Ca absorption, increased bone resorption, and reduced renal tubular Ca reabsorption leading to hypercalciuria and all form kidney stones. GHS have increased vitamin D receptors (VDR) at these sites of Ca transport. Injection of 1,25(OH)2D3 (1,25D) leads to a greater increase in urine (u)Ca in GHS than in control Sprague-Dawley (SD), possibly due to the additional VDR. In GHS the increased uCa persists on a low-Ca diet (LCD) suggesting enhanced bone resorption. We tested the hypothesis that LCD, coupled to inhibition of bone resorption by alendronate (alen), would eliminate the enhanced 1,25D-induced hypercalciuria in GHS. SD and GHS were fed LCD and half were injected daily with 1,25D. After 8 days all were also given alen until euthanasia at day 16. At 8 days, 1,25D increased uCa in SD and to a greater extent in GHS. At 16 days, alen eliminated the 1,25D-induced increase in uCa in SD. However, in GHS alen decreased, but did not eliminate, the 1,25D-induced hypercalciuria, suggesting maximal alen cannot completely prevent the 1,25D-induced bone resorption in GHS, perhaps due to increased VDR. There was no consistent effect on mRNA expression of renal transcellular or paracellular Ca transporters. Urine CaP and CaOx supersaturation (SS) increased with 1,25D alone in both SD and GHS. Alen eliminated the increase in CaP SS in SD but not in GHS. If these results are confirmed in humans with IH, the use of bisphosphonates, such as alen, may not prevent the decreased bone density observed in these patients.

Additive, epistatic, and environmental effects through the lens of expression variability QTL in a twin cohort

The expression of a gene can vary across individuals in the general population, as well as between monozygotic twins. This variable expression is assumed to be due to the influence of both genetic and nongenetic factors. Yet little evidence supporting this assumption has been obtained from empirical data. In this study, we used expression data from a large twin cohort to investigate the influences of genetic and nongenetic factors on variable gene expression. We focused on a set of expression variability QTL (evQTL)--i.e., genetic loci associated with the variance, as opposed to the mean, of gene expression. We identified evQTL for 99, 56, and 79 genes in lymphoblastoid cell lines, skin, and fat, respectively. The differences in gene expression, measured by the relative mean difference (RMD), tended to be larger between pairs of dizygotic (DZ) twins than between pairs of monozygotic (MZ) twins, showing that genetic background influenced the expression variability. Furthermore, a more profound RMD was observed between pairs of MZ twins whose genotypes were associated with greater expression variability than the RMD found between pairs of MZ twins whose genotypes were associated with smaller expression variability. This suggests that nongenetic (e.g., environmental) factors contribute to the variable expression. Lastly, we demonstrated that the formation of evQTL is likely due to partial linkages between eQTL SNPs that are additively associated with the mean of gene expression; in most cases, no epistatic effect is involved. Our findings have implications for understanding divergent sources of gene expression variability.

1,25(OH)₂D₃-enhanced hypercalciuria in genetic hypercalciuric stone-forming rats fed a low-calcium diet

The inbred genetic hypercalciuric stone-forming (GHS) rats exhibit many features of human idiopathic hypercalciuria and have elevated levels of vitamin D receptors (VDR) in calcium (Ca)-transporting organs. On a normal-Ca diet, 1,25(OH)2D3 (1,25D) increases urine (U) Ca to a greater extent in GHS than in controls [Sprague-Dawley (SD)]. The additional UCa may result from an increase in intestinal Ca absorption and/or bone resorption. To determine the source, we asked whether 1,25D would increase UCa in GHS fed a low-Ca (0.02%) diet (LCD). With 1,25D, UCa in SD increased from 1.2 ± 0.1 to 9.3 ± 0.9 mg/day and increased more in GHS from 4.7 ± 0.3 to 21.5 ± 0.9 mg/day (P < 0.001). In GHS rats on LCD with or without 1,25D, UCa far exceeded daily Ca intake (2.6 mg/day). While the greater excess in UCa in GHS rats must be derived from bone mineral, there may also be a 1,25D-mediated decrease in renal tubular Ca reabsorption. RNA expression of the components of renal Ca transport indicated that 1,25D administration results in a suppression of klotho, an activator of the renal Ca reabsorption channel TRPV5, in both SD and GHS rats. This fall in klotho would decrease tubular reabsorption of the 1,25D-induced bone Ca release. Thus, the greater increase in UCa with 1,25D in GHS fed LCD strongly suggests that the additional UCa results from an increase in bone resorption, likely due to the increased number of VDR in the GHS rat bone cells, with a possible component of decreased renal tubular calcium reabsorption.

Effects of Sex on Intra-Individual Variance in Urinary Solutes in Stone-Formers Collected from a Single Clinical Laboratory

Background/Aims

Our work in a rodent model of urinary calcium suggests genetic and gender effects on increased residual variability in urine chemistries. Based on these findings, we hypothesized that sex would similarly be associated with residual variation in human urine solutes. Sex-related effects on residuals might affect the establishment of physiological baselines and error in medical assays.

Methods

We tested the effects of sex on residual variation in urine chemistry by estimating coefficients of variation (CV) for urinary solutes in paired sequential 24-h urines (≤72 hour interval) in 6,758 females and 9,024 males aged 16–80 submitted to a clinical laboratory.

Results

Females had higher CVs than males for urinary phosphorus overall at the False Discovery Rate (P<0.01). There was no effect of sex on CV for calcium (P>0.3). Males had higher CVs for citrate (P<0.01) from ages 16–45 and females higher CVs for citrate (P<0.01) from ages 56–80, suggesting effects of an extant oestral cycle on residual variance.

Conclusions

Our findings indicate the effects of sex on residual variance of the excretion of urinary solutes including phosphorus and citrate; differences in CV by sex might reflect dietary lability, differences in the fidelity of reporting or genetic differentiation in renal solute consistency. Such an effect could complicate medical analysis by the addition of random error to phenotypic assays. Renal analysis might require explicit incorporation of heterogeneity among factorial effects, and for sex in particular.

Weight, age and coefficients of variation in renal solute excretion

BACKGROUND

Homoscedasticity (constant variance over axes or among statistical factors) is an integral assumption of most statistical analyses. However, a number of empirical studies in model organisms and humans demonstrate significant differences in residual variance (that component of phenotype unexplained by known factors) or intra-individual variation among genotypes. Our work suggests that renal traits may be particularly susceptible to randomization by genetic and non-genetic factors, including endogenous variables like age and weight.

METHODS

We tested associations between age, weight and intra-individual variation in urinary calcium, citrate, chloride, creatinine, potassium, magnesium, sodium, ammonium, oxalate, phosphorus, sulfate, uric acid and urea nitrogen in 9,024 male and 6,758 female kidney stone patients. Coefficients of variation (CVs) were calculated for each individual for each solute from paired 24-hour urines. Analysis of CVs was corrected for inter-measurement collection variance in creatinine and urine volume. CVs for sodium and urea nitrogen were included to correct for dietary salt and protein.

RESULTS

Age was positively associated with individual CVs for calcium and negatively associated with CVs for potassium, ammonium and phosphorus (p(FDR) < 0.01). Weight was associated with CVs for creatinine, magnesium and uric acid, and negatively associated with CVs for calcium, potassium and oxalate (p(FDR) < 0.05).

CONCLUSION

Intra-individual variation changes over age and weight axes for numerous urinary solutes. Changing residual variance over age and weight could cause bias in the detection or estimation of genetic or environmental effects. New methodologies may need to account for such residual unpredictability, especially in diverse collections.

Increased biological response to 1,25(OH)(2)D(3) in genetic hypercalciuric stone-forming rats

Genetic hypercalciuric stone-forming (GHS) rats, bred to maximize urine (U) calcium (Ca) excretion, have increased intestinal Ca absorption and bone Ca resorption and reduced renal Ca reabsorption, leading to increased UCa compared with the Sprague-Dawley (SD) rats. GHS rats have increased vitamin D receptors (VDR) at each of these sites, with normal levels of 1,25(OH)(2)D(3) (1,25D), indicating that their VDR is undersaturated with 1,25D. We tested the hypothesis that 1,25D would induce a greater increase in UCa in GHS rats by feeding both strains ample Ca and injecting 1,25D (25 ng · 100 g body wt(-1) · day(-1)) or vehicle for 16 days. With 1,25D, UCa in SD increased from 1.7 ± 0.3 mg/day to 24.4 ± 1.2 (Δ = 22.4 ± 1.5) and increased more in GHS from 10.5 ± 0.7 to 41.9 ± 0.7 (Δ = 29.8 ± 1.8; P = 0.003). To determine the mechanism of the greater increase in UCa in GHS rats, we measured kidney RNA expression of components of renal Ca transport. Expression of transient receptor potential vanilloid (TRPV)5 and calbindin D(28K) were increased similarly in SD + 1,25D and GHS + 1,25D. The Na(+)/Ca(2+) exchanger (NCX1) was increased in GHS + 1,25D. Klotho was decreased in SD + 1,25D and GHS + 1,25D. TRPV6 was increased in SD + 1,25D and increased further in GHS + 1,25D. Claudin 14, 16, and 19, Na/K/2Cl transporter (NKCC2), and secretory K channel (ROMK) did not differ between SD + 1,25D and GHS + 1,25D. Increased UCa with 1,25D in GHS exceeded that of SD, indicating that the increased VDR in GHS induces a greater biological response. This increase in UCa, which must come from the intestine and/or bone, must exceed any effect of 1,25D on TRPV6 or NCX1-mediated renal Ca reabsorption.

Recent developments in statistical methods for detecting genetic loci affecting phenotypic variability

A number of recent works have introduced statistical methods for detecting genetic loci that affect phenotypic variability, which we refer to as variability-controlling quantitative trait loci (vQTL). These are genetic variants whose allelic state predicts how much phenotype values will vary about their expected means. Such loci are of great potential interest in both human and non-human genetic studies, one reason being that a detected vQTL could represent a previously undetected interaction with other genes or environmental factors. The simultaneous publication of these new methods in different journals has in many cases precluded opportunity for comparison. We survey some of these methods, the respective trade-offs they imply, and the connections between them. The methods fall into three main groups: classical non-parametric, fully parametric, and semi-parametric two-stage approximations. Choosing between alternatives involves balancing the need for robustness, flexibility, and speed. For each method, we identify important assumptions and limitations, including those of practical importance, such as their scope for including covariates and random effects. We show in simulations that both parametric methods and their semi-parametric approximations can give elevated false positive rates when they ignore mean-variance relationships intrinsic to the data generation process. We conclude that choice of method depends on the trait distribution, the need to include non-genetic covariates, and the population size and structure, coupled with a critical evaluation of how these fit with the assumptions of the statistical model.