Restriction site-associated DNA sequencing, genotyping error estimation and de novo assembly optimization for population genetic inference.

Restriction site-associated DNA sequencing (RADseq) provides researchers with the ability to record genetic polymorphism across thousands of loci for nonmodel organisms, potentially revolutionizing the field of molecular ecology. However, as with other genotyping methods, RADseq is prone to a number of sources of error that may have consequential effects for population genetic inferences, and these have received only limited attention in terms of the estimation and reporting of genotyping error rates. Here we use individual sample replicates, under the expectation of identical genotypes, to quantify genotyping error in the absence of a reference genome. We then use sample replicates to (i) optimize de novo assembly parameters within the program Stacks, by minimizing error and maximizing the retrieval of informative loci; and (ii) quantify error rates for loci, alleles and single-nucleotide polymorphisms. As an empirical example, we use a double-digest RAD data set of a nonmodel plant species, Berberis alpina, collected from high-altitude mountains in Mexico.

Pseudohypoaldosteronisms, report on a 10-patient series.

BACKGROUND: Type 1 pseudohypoaldosteronism (PHA1) is a salt-wasting syndrome caused by mineralocorticoid resistance. Autosomal recessive and dominant hereditary forms are caused by Epithelial Na Channel and Mineralocorticoid Receptor mutation respectively, while secondary PHA1 is usually associated with urological problems. METHODS: Ten patients were studied in four French pediatric units in order to characterize PHA1 spectrum in infants. Patients were selected by chart review. Genetic, clinical and biochemistry data were collected and analyzed. RESULTS: Autosomal recessive PHA1 (n = 3) was diagnosed at 6 and 7 days of life in three patients presenting with severe hyperkalaemia and weight loss. After 8 months, 3 and 5 years on follow-up, neurological development and longitudinal growth was normal with high sodium supplementation. Autosomal dominant PHA1 (n = 4) was revealed at 15, 19, 22 and 30 days of life because of failure to thrive. At 8 months, 3 and 21 years of age, longitudinal growth was normal in three patients who were given salt supplementation; no significant catch-up growth was obtained in the last patient at 20 months of age. Secondary PHA1 (n = 3) was diagnosed at 11, 26 days and 5 months of life concomitantly with acute pyelonephritis in three children with either renal hypoplasia, urinary duplication or bilateral megaureter. The outcome was favourable and salt supplementation was discontinued after 3, 11 and 13 months. CONCLUSIONS: PHA1 should be suspected in case of severe hyperkalemia and weight loss in infants and need careful management. Pathogenesis of secondary PHA1 is still challenging and further studies are mandatory to highlight the link between infection, developing urinary tract and pseudohypoaldosteronism.

Subgroup analyses in randomised controlled trials: cohort study on trial protocols and journal publications

OBJECTIVE: To investigate the planning of subgroup analyses in protocols of randomised controlled trials and the agreement with corresponding full journal publications. DESIGN: Cohort of protocols of randomised controlled trial and subsequent full journal publications. SETTING: Six research ethics committees in Switzerland, Germany, and Canada. DATA SOURCES: 894 protocols of randomised controlled trial involving patients approved by participating research ethics committees between 2000 and 2003 and 515 subsequent full journal publications. RESULTS: Of 894 protocols of randomised controlled trials, 252 (28.2%) included one or more planned subgroup analyses. Of those, 17 (6.7%) provided a clear hypothesis for at least one subgroup analysis, 10 (4.0%) anticipated the direction of a subgroup effect, and 87 (34.5%) planned a statistical test for interaction. Industry sponsored trials more often planned subgroup analyses compared with investigator sponsored trials (195/551 (35.4%) v 57/343 (16.6%), P<0.001). Of 515 identified journal publications, 246 (47.8%) reported at least one subgroup analysis. In 81 (32.9%) of the 246 publications reporting subgroup analyses, authors stated that subgroup analyses were prespecified, but this was not supported by 28 (34.6%) corresponding protocols. In 86 publications, authors claimed a subgroup effect, but only 36 (41.9%) corresponding protocols reported a planned subgroup analysis. CONCLUSIONS: Subgroup analyses are insufficiently described in the protocols of randomised controlled trials submitted to research ethics committees, and investigators rarely specify the anticipated direction of subgroup effects. More than one third of statements in publications of randomised controlled trials about subgroup prespecification had no documentation in the corresponding protocols. Definitive judgments regarding credibility of claimed subgroup effects are not possible without access to protocols and analysis plans of randomised controlled trials.

Preferential assembly of epithelial sodium channel (ENaC) subunits in Xenopus oocytes: role of furin-mediated endogenous proteolysis.

The epithelial sodium channel (ENaC) is preferentially assembled into heteromeric alphabetagamma complexes. The alpha and gamma (not beta) subunits undergo proteolytic cleavage by endogenous furin-like activity correlating with increased ENaC function. We identified full-length subunits and their fragments at the cell surface, as well as in the intracellular pool, for all homo- and heteromeric combinations (alpha, beta, gamma, alphabeta, alphagamma, betagamma, and alphabetagamma). We assayed corresponding channel function as amiloride-sensitive sodium transport (I(Na)). We varied furin-mediated proteolysis by mutating the P1 site in alpha and/or gamma subunit furin consensus cleavage sites (alpha(mut) and gamma(mut)). Our findings were as follows. (i) The beta subunit alone is not transported to the cell surface nor cleaved upon assembly with the alpha and/or gamma subunits. (ii) The alpha subunit alone (or in combination with beta and/or gamma) is efficiently transported to the cell surface; a surface-expressed 65-kDa alpha ENaC fragment is undetected in alpha(mut)betagamma, and I(Na) is decreased by 60%. (iii) The gamma subunit alone does not appear at the cell surface; gamma co-expressed with alpha reaches the surface but is not detectably cleaved; and gamma in alphabetagamma complexes appears mainly as a 76-kDa species in the surface pool. Although basal I(Na) of alphabetagamma(mut) was similar to alphabetagamma, gamma(mut) was not detectably cleaved at the cell surface. Thus, furin-mediated cleavage is not essential for participation of alpha and gamma in alphabetagamma heteromers. Basal I(Na) is reduced by preventing furin-mediated cleavage of the alpha, but not gamma, subunits. Residual current in the absence of furin-mediated proteolysis may be due to non-furin endogenous proteases.

Airway surface liquid volume regulation determines different airway phenotypes in liddle compared with betaENaC-overexpressing mice.

Studies in cystic fibrosis patients and mice overexpressing the epithelial Na(+) channel beta-subunit (betaENaC-Tg) suggest that raised airway Na(+) transport and airway surface liquid (ASL) depletion are central to the pathogenesis of cystic fibrosis lung disease. However, patients or mice with Liddle gain-of-function betaENaC mutations exhibit hypertension but no lung disease. To investigate this apparent paradox, we compared the airway phenotype (nasal versus tracheal) of Liddle with CFTR-null, betaENaC-Tg, and double mutant mice. In mouse nasal epithelium, the region that functionally mimics human airways, high levels of CFTR expression inhibited Liddle epithelial Nat channel (ENaC) hyperfunction. Conversely, in mouse trachea, low levels of CFTR failed to suppress Liddle ENaC hyperfunction. Indeed, Na(+) transport measured in Ussing chambers ("flooded" conditions) was raised in both Liddle and betaENaC-Tg mice. Because enhanced Na(+) transport did not correlate with lung disease in these mutant mice, measurements in tracheal cultures under physiologic "thin film" conditions and in vivo were performed. Regulation of ASL volume and ENaC-mediated Na(+) absorption were intact in Liddle but defective in betaENaC-Tg mice. We conclude that the capacity to regulate Na(+) transport and ASL volume, not absolute Na(+) transport rates in Ussing chambers, is the key physiologic function protecting airways from dehydration-induced lung disease.

Genetic dissection of sodium and potassium transport along the aldosterone-sensitive distal nephron: Importance in the control of blood pressure and hypertension.

In this review, we discuss genetic evidence supporting Guyton's hypothesis stating that blood pressure control is critically depending on fluid handling by the kidney. The review is focused on the genetic dissection of sodium and potassium transport in the distal nephron and the collecting duct that are the most important sites for the control of sodium and potassium balance by aldosterone and angiotensin II. Thanks to the study of Mendelian forms of hypertension and their corresponding transgenic mouse models, three main classes of diuretic receptors (furosemide, thiazide, amiloride) and the main components of the aldosterone- and angiotensin-dependent signaling pathways were molecularly identified over the past 20years. This will allow to design rational strategies for the treatment of hypertension and for the development of the next generation of diuretics.

Pregnancy-induced chromatin remodeling in the breast of postmenopausal women.

Early pregnancy and multiparity are known to reduce the risk of women to develop breast cancer at menopause. Based on the knowledge that the differentiation of the breast induced by the hormones of pregnancy plays a major role in this protection, this work was performed with the purpose of identifying what differentiation-associated molecular changes persist in the breast until menopause. Core needle biopsies (CNB) obtained from the breast of 42 nulliparous (NP) and 71 parous (P) postmenopausal women were analyzed in morphology, immunocytochemistry and gene expression. Whereas in the NP breast, nuclei of epithelial cells were large and euchromatic, in the P breast they were small and hyperchromatic, showing strong methylation of histone 3 at lysine 9 and 27. Transcriptomic analysis performed using Affymetrix HG_U133 oligonucleotide arrays revealed that in CNB of the P breast, there were 267 upregulated probesets that comprised genes controlling chromatin organization, transcription regulation, splicing machinery, mRNA processing and noncoding elements including XIST. We concluded that the differentiation process induced by pregnancy is centered in chromatin remodeling and in the mRNA processing reactome, both of which emerge as important regulatory pathways. These are indicative of a safeguard step that maintains the fidelity of the transcription process, becoming the ultimate mechanism mediating the protection of the breast conferred by full-term pregnancy.

Recurrent deletions and reciprocal duplications of 10q11.21q11.23 including CHAT and SLC18A3 are likely mediated by complex low-copy repeats.

We report 24 unrelated individuals with deletions and 17 additional cases with duplications at 10q11.21q21.1 identified by chromosomal microarray analysis. The rearrangements range in size from 0.3 to 12 Mb. Nineteen of the deletions and eight duplications are flanked by large, directly oriented segmental duplications of >98% sequence identity, suggesting that nonallelic homologous recombination (NAHR) caused these genomic rearrangements. Nine individuals with deletions and five with duplications have additional copy number changes. Detailed clinical evaluation of 20 patients with deletions revealed variable clinical features, with developmental delay (DD) and/or intellectual disability (ID) as the only features common to a majority of individuals. We suggest that some of the other features present in more than one patient with deletion, including hypotonia, sleep apnea, chronic constipation, gastroesophageal and vesicoureteral refluxes, epilepsy, ataxia, dysphagia, nystagmus, and ptosis may result from deletion of the CHAT gene, encoding choline acetyltransferase, and the SLC18A3 gene, mapping in the first intron of CHAT and encoding vesicular acetylcholine transporter. The phenotypic diversity and presence of the deletion in apparently normal carrier parents suggest that subjects carrying 10q11.21q11.23 deletions may exhibit variable phenotypic expressivity and incomplete penetrance influenced by additional genetic and nongenetic modifiers.

ENaC-mediated alveolar fluid clearance and lung fluid balance depend on the channel-activating protease 1.

Sodium transport via epithelial sodium channels (ENaC) expressed in alveolar epithelial cells (AEC) provides the driving force for removal of fluid from the alveolar space. The membrane-bound channel-activating protease 1 (CAP1/Prss8) activates ENaC in vitro in various expression systems. To study the role of CAP1/Prss8 in alveolar sodium transport and lung fluid balance in vivo, we generated mice lacking CAP1/Prss8 in the alveolar epithelium using conditional Cre-loxP-mediated recombination. Deficiency of CAP1/Prss8 in AEC induced in vitro a 40% decrease in ENaC-mediated sodium currents. Sodium-driven alveolar fluid clearance (AFC) was reduced in CAP1/Prss8-deficient mice, due to a 48% decrease in amiloride-sensitive clearance, and was less sensitive to beta(2)-agonist treatment. Intra-alveolar treatment with neutrophil elastase, a soluble serine protease activating ENaC at the cell surface, fully restored basal AFC and the stimulation by beta(2)-agonists. Finally, acute volume-overload increased alveolar lining fluid volume in CAP1/Prss8-deficient mice. This study reveals that CAP1 plays a crucial role in the regulation of ENaC-mediated alveolar sodium and water transport and in mouse lung fluid balance.

alphaENaC-mediated lithium absorption promotes nephrogenic diabetes insipidus.

Lithium-induced nephrogenic diabetes insipidus (NDI) is accompanied by polyuria, downregulation of aquaporin 2 (AQP2), and cellular remodeling of the collecting duct (CD). The amiloride-sensitive epithelial sodium channel (ENaC) is a likely candidate for lithium entry. Here, we subjected transgenic mice lacking αENaC specifically in the CD (knockout [KO] mice) and littermate controls to chronic lithium treatment. In contrast to control mice, KO mice did not markedly increase their water intake. Furthermore, KO mice did not demonstrate the polyuria and reduction in urine osmolality induced by lithium treatment in the control mice. Lithium treatment reduced AQP2 protein levels in the cortex/outer medulla and inner medulla (IM) of control mice but only partially reduced AQP2 levels in the IM of KO mice. Furthermore, lithium induced expression of H(+)-ATPase in the IM of control mice but not KO mice. In conclusion, the absence of functional ENaC in the CD protects mice from lithium-induced NDI. These data support the hypothesis that ENaC-mediated lithium entry into the CD principal cells contributes to the pathogenesis of lithium-induced NDI.

Sodium and potassium balance depends on αENaC expression in connecting tubule.

Mutations in α, β, or γ subunits of the epithelial sodium channel (ENaC) can downregulate ENaC activity and cause a severe salt-losing syndrome with hyperkalemia and metabolic acidosis, designated pseudohypoaldosteronism type 1 in humans. In contrast, mice with selective inactivation of αENaC in the collecting duct (CD) maintain sodium and potassium balance, suggesting that the late distal convoluted tubule (DCT2) and/or the connecting tubule (CNT) participates in sodium homeostasis. To investigate the relative importance of ENaC-mediated sodium absorption in the CNT, we used Cre-lox technology to generate mice lacking αENaC in the aquaporin 2-expressing CNT and CD. Western blot analysis of microdissected cortical CD (CCD) and CNT revealed absence of αENaC in the CCD and weak αENaC expression in the CNT. These mice exhibited a significantly higher urinary sodium excretion, a lower urine osmolality, and an increased urine volume compared with control mice. Furthermore, serum sodium was lower and potassium levels were higher in the genetically modified mice. With dietary sodium restriction, these mice experienced significant weight loss, increased urinary sodium excretion, and hyperkalemia. Plasma aldosterone levels were significantly elevated under both standard and sodium-restricted diets. In summary, αENaC expression within the CNT/CD is crucial for sodium and potassium homeostasis and causes signs and symptoms of pseudohypoaldosteronism type 1 if missing.

Mapping Bias Overestimates Reference Allele Frequencies at the HLA Genes in the 1000 Genomes Project Phase I Data.

Next-generation sequencing (NGS) technologies have become the standard for data generation in studies of population genomics, as the 1000 Genomes Project (1000G). However, these techniques are known to be problematic when applied to highly polymorphic genomic regions, such as the human leukocyte antigen (HLA) genes. Because accurate genotype calls and allele frequency estimations are crucial to population genomics analyses, it is important to assess the reliability of NGS data. Here, we evaluate the reliability of genotype calls and allele frequency estimates of the single-nucleotide polymorphisms (SNPs) reported by 1000G (phase I) at five HLA genes (HLA-A, -B, -C, -DRB1, and -DQB1). We take advantage of the availability of HLA Sanger sequencing of 930 of the 1092 1000G samples and use this as a gold standard to benchmark the 1000G data. We document that 18.6% of SNP genotype calls in HLA genes are incorrect and that allele frequencies are estimated with an error greater than ±0.1 at approximately 25% of the SNPs in HLA genes. We found a bias toward overestimation of reference allele frequency for the 1000G data, indicating mapping bias is an important cause of error in frequency estimation in this dataset. We provide a list of sites that have poor allele frequency estimates and discuss the outcomes of including those sites in different kinds of analyses. Because the HLA region is the most polymorphic in the human genome, our results provide insights into the challenges of using of NGS data at other genomic regions of high diversity.