How many biological replicates are needed in an RNA-seq experiment and which differential expression tool should you use?

RNA-seq is now the technology of choice for genome-wide differential gene expression experiments, but it is not clear how many biological replicates are needed to ensure valid biological interpretation of the results or which statistical tools are best for analyzing the data. An RNA-seq experiment with 48 biological replicates in each of two conditions was performed to answer these questions and provide guidelines for experimental design. With three biological replicates, nine of the 11 tools evaluated found only 20%–40% of the significantly differentially expressed (SDE) genes identified with the full set of 42 clean replicates. This rises to >85% for the subset of SDE genes changing in expression by more than fourfold. To achieve >85% for all SDE genes regardless of fold change requires more than 20 biological replicates. The same nine tools successfully control their false discovery rate at ≲5% for all numbers of replicates, while the remaining two tools fail to control their FDR adequately, particularly for low numbers of replicates. For future RNA-seq experiments, these results suggest that at least six biological replicates should be used, rising to at least 12 when it is important to identify SDE genes for all fold changes. If fewer than 12 replicates are used, a superior combination of true positive and false positive performances makes edgeR and DESeq2 the leading tools. For higher replicate numbers, minimizing false positives is more important and DESeq marginally outperforms the other tools.     

Sonic hedgehog signaling from the urethral epithelium controls external genital development

External genital development begins with formation of paired genital swellings, which develop into the genital tubercle. Proximodistal outgrowth and axial patterning of the genital tubercle are coordinated to give rise to the penis or clitoris. The genital tubercle consists of lateral plate mesoderm, surface ectoderm, and endodermal urethral epithelium derived from the urogenital sinus. We have investigated the molecular control of external genital development in the mouse embryo. Previous work has shown that the genital tubercle has polarizing activity, but the precise location of this activity within the tubercle is unknown. We reasoned that if the tubercle itself is patterned by a specialized signaling region, then polarizing activity may be restricted to a subset of cells. Transplantation of urethral epithelium, but not genital mesenchyme, to chick limbs results in mirror-image duplication of the digits. Moreover, when grafted to chick limbs, the urethral plate orchestrates morphogenetic movements normally associated with external genital development. Signaling activity is therefore restricted to urethral plate cells. Before and during normal genital tubercle outgrowth, urethral plate epithelium expresses Sonic hedgehog (Shh). In mice with a targeted deletion of Shh, external genitalia are absent. Genital swellings are initiated, but outgrowth is not maintained. In the absence of Shh signaling, Fgf8, Bmp2, Bmp4, Fgf10, and Wnt5a are downregulated, and apoptosis is enhanced in the genitalia. These results identify the urethral epithelium as a signaling center of the genital tubercle, and demonstrate that Shh from the urethral epithelium is required for outgrowth, patterning, and cell survival in the developing external genitalia.     

Environmental conditions at arrival to the wintering grounds and during spring migration affect population dynamics of barn swallows Hirundo rustica breeding in Northern Italy

Several populations of long-distance migratory birds are currently suffering steep demographic declines. The identification of the causes of such declines is difficult because population changes may be driven by events occurring in distant geographical areas during different phases of the annual life-cycle of migrants. Furthermore, wintering areas and migration routes of populations of small-sized species are still largely unknown, with few exceptions. In this paper we identified the critical phases of the annual life-cycle that most influence the population dynamics of a small passerine, the Barn Swallow Hirundo rustica. We used information on temporal dynamics of a population breeding in Northern Italy, whose wintering range and timing of migration have been recently described by miniaturised tracking dataloggers. Our results indicated that primary productivity in the wintering grounds in the month when most individuals arrive from autumn migration and primary productivity in an area that is probably a stopover site during spring migration, influenced population dynamics more than habitat conditions at the breeding grounds. By using annual variation in primary productivity at the wintering grounds and stopover sites as predictors, we replicated the observed interannual population changes with great accuracy. However, the steep decline recently suffered by the population could be replicated only by including a constant annual decline in the model, suggesting that changes in primary productivity only predicted the interannual variation around the long-term trend. Our study therefore suggests the existence of critical periods during wintering and migration that may have large impact on population fluctuations of migrant birds.     

Improving Maternal Care through a State-Wide Health Insurance Program: A Cost and Cost-Effectiveness Study in Rural Nigeria

Background

While the Nigerian government has made progress towards the Millennium Development Goals, further investments are needed to achieve the targets of post-2015 Sustainable Development Goals, including Universal Health Coverage. Economic evaluations of innovative interventions can help inform investment decisions in resource-constrained settings. We aim to assess the cost and cost-effectiveness of maternal care provided within the new Kwara State Health Insurance program (KSHI) in rural Nigeria.

Methods and Findings

We used a decision analytic model to simulate a cohort of pregnant women. The primary outcome is the incremental cost effectiveness ratio (ICER) of the KSHI scenario compared to the current standard of care. Intervention cost from a healthcare provider perspective included service delivery costs and above-service level costs; these were evaluated in a participating hospital and using financial records from the managing organisations, respectively. Standard of care costs from a provider perspective were derived from the literature using an ingredient approach. We generated 95% credibility intervals around the primary outcome through probabilistic sensitivity analysis (PSA) based on a Monte Carlo simulation. We conducted one-way sensitivity analyses across key model parameters and assessed the sensitivity of our results to the performance of the base case separately through a scenario analysis. Finally, we assessed the sustainability and feasibility of this program’s scale up within the State’s healthcare financing structure through a budget impact analysis. The KSHI scenario results in a health benefit to patients at a higher cost compared to the base case. The mean ICER (US$46.4/disability-adjusted life year averted) is considered very cost-effective compared to a willingness-to-pay threshold of one gross domestic product per capita (Nigeria, US$ 2012, 2,730). Our conclusion was robust to uncertainty in parameters estimates (PSA: median US$49.1, 95% credible interval 21.9–152.3), during one-way sensitivity analyses, and when cost, quality, cost and utilization parameters of the base case scenario were changed. The sustainability of this program’s scale up by the State is dependent on further investments in healthcare.

Conclusions

This study provides evidence that the investment made by the KSHI program in rural Nigeria is likely to have been cost-effective; however, further healthcare investments are needed for this program to be successfully expanded within Kwara State. Policy makers should consider supporting financial initiatives to reduce maternal mortality tackling both supply and demand issues in the access to care.     

Lung hypoplasia and neonatal death in Fgf9-null mice identify this gene as an essential regulator of lung mesenchyme

Mammalian lung develops as an evagination of ventral gut endoderm into the underlying mesenchyme. Iterative epithelial branching, regulated by the surrounding mesenchyme, generates an elaborate network of airways from the initial lung bud. Fibroblast growth factors (FGFs) often mediate epithelial-mesenchymal interactions and mesenchymal Fgf10 is essential for epithelial branching in the developing lung. However, no FGF has been shown to regulate lung mesenchyme. In embryonic lung, Fgf9 is detected in airway epithelium and visceral pleura at E10.5, but is restricted to the pleura by E12.5. We report that mice homozygous for a targeted disruption of Fgf9 exhibit lung hypoplasia and early postnatal death. Fgf9(-/-) lungs exhibit reduced mesenchyme and decreased branching of airways, but show significant distal airspace formation and pneumocyte differentiation. Our results suggest that Fgf9 affects lung size by stimulating mesenchymal proliferation. The reduction in the amount of mesenchyme in Fgf9(-/-) lungs limits expression of mesenchymal Fgf10. We suggest a model whereby FGF9 signaling from the epithelium and reciprocal FGF10 signaling from the mesenchyme coordinately regulate epithelial airway branching and organ size during lung embryogenesis.     

The classical complement cascade mediates CNS synapse elimination

During development, the formation of mature neural circuits requires the selective elimination of inappropriate synaptic connections. Here we show that C1q, the initiating protein in the classical complement cascade, is expressed by postnatal neurons in response to immature astrocytes and is localized to synapses throughout the postnatal CNS and retina. Mice deficient in complement protein C1q or the downstream complement protein C3 exhibit large sustained defects in CNS synapse elimination, as shown by the failure of anatomical refinement of retinogeniculate connections and the retention of excess retinal innervation by lateral geniculate neurons. Neuronal C1q is normally downregulated in the adult CNS; however, in a mouse model of glaucoma, C1q becomes upregulated and synaptically relocalized in the adult retina early in the disease. These findings support a model in which unwanted synapses are tagged by complement for elimination and suggest that complement-mediated synapse elimination may become aberrantly reactivated in neurodegenerative disease.     

Impact of definitive therapy with beta-lactam monotherapy or combination with an aminoglycoside or a quinolone for Pseudomonas aeruginosa bacteremia

Background

Bacteremia by Pseudomonas aeruginosa represents one severe infection. It is not clear whether beta-lactam monotherapy leads to similar rates of treatment success compared to combinations of beta-lactams with aminoglycosides or quinolones.

Methods

Retrospective cohort study from 3 tertiary hospitals (2 in Greece and 1 in Italy). Pseudomonas aeruginosa isolates were susceptible to a beta-lactam and an aminoglycoside or a quinolone. Patients received appropriate therapy for at least 48 hours. Primary outcome of interest was treatment success in patients with definitive beta-lactam combination therapy compared to monotherapy. Secondary outcomes were treatment success keeping the same empirical and definitive regimen, mortality, and toxicity.

Results

Out of 92 bacteremias there were 54 evaluable episodes for the primary outcome (20 received monotherapy). Treatment success was higher with combination therapy (85%) compared to beta-lactam monotherapy (65%), however not statistically significantly [Odds ratio (OR) 3.1; 95% Confidence Interval (CI) 0.69–14.7, p = 0.1]. Very long (>2 months) hospitalisation before bacteremia was the only factor independently associated with treatment success (OR 0.73; 95% CI 0.01–0.95, p = 0.046), however this result entailed few episodes. All-cause mortality did not differ significantly between combination therapy [6/31 (19%)] and monotherapy [8/19 (42%)], p = 0.11. Only Charlson comorbidity index was associated with excess mortality (p = 0.03).

Conclusion

Our study, in accordance with previous ones, indicates that the choice between monotherapy and combination therapy may not affect treatment success significantly. However, our study does not have statistical power to identify small or moderate differences. A large randomized controlled trial evaluating this issue is justified.