Phosphorus remobilization from rice flag leaves during grain filling: an RNA‐seq study

The physiology and molecular regulation of phosphorus (P) remobilization from vegetative tissues to grains during grain filling is poorly understood, despite the pivotal role it plays in the global P cycle. To test the hypothesis that a subset of genes involved in the P starvation response are involved in remobilization of P from flag leaves to developing grains, we conducted an RNA‐seq analysis of rice flag leaves during the preremobilization phase (6 DAA) and when the leaves were acting as a P source (15 DAA). Several genes that respond to phosphate starvation, including three purple acid phosphatases (OsPAP3, OsPAP9b and OsPAP10a), were significantly up‐regulated at 15 DAA, consistent with a role in remobilization of P from flag leaves during grain filling. A number of genes that have not been implicated in the phosphate starvation response, OsPAP26, SPX‐MFS1 (a putative P transporter) and SPX‐MFS2, also showed expression profiles consistent with involvement in P remobilization from senescing flag leaves. Metabolic pathway analysis using the KEGG system suggested plastid membrane lipid synthesis is a critical process during the P remobilization phase. In particular, the up‐regulation of OsPLDz2 and OsSQD2 at 15 DAA suggested phospholipids were being degraded and replaced by other lipids to enable continued cellular function while liberating P for export to developing grains. Three genes associated with RNA degradation that have not previously been implicated in the P starvation response also showed expression profiles consistent with a role in P mobilization from senescing flag leaves.     

Next generation semiconductor based sequencing of bitter taste receptor genes in different pig populations and association analysis using a selective DNA pool‐seq approach

Taste perception in animals affects feed intake and may influence production traits. In particular, bitter is sensed by receptors encoded by the family of TAS2R genes. In this research, using a DNA pool‐seq approach coupled with next generation semiconductor based target resequencing, we analysed nine porcine TAS2R genes (TAS2R1, TAS2R3, TAS2R4, TAS2R7, TAS2R9, TAS2R10, TAS2R16, TAS2R38 and TAS2R39) to identify variability and, at the same time, estimate single nucleotide polymorphism (SNP) allele frequencies in several populations and testing differences in an association analysis. Equimolar DNA pools were prepared for five pig breeds (Italian Duroc, Italian Landrace, Pietrain, Meishan and Casertana) and wild boars (5–10 individuals each) and for two groups of Italian Large White pigs with extreme and divergent back fat thickness (50 + 50 pigs). About 1.8 million reads were obtained by sequencing amplicons generated from these pools. A total of 125 SNPs were identified, of which 37 were missense mutations. Three of them (p.Ile53Phe and p.Trp85Leu in TAS2R4; p.Leu37Ser in TAS2R39) could have important effects on the function of these bitter taste receptors, based on in silico predictions. Variability in wild boars seems lower than that in domestic breeds potentially as a result of selective pressure in the wild towards defensive bitter taste perception. Three SNPs in TAS2R38 and TAS2R39 were significantly associated with back fat thickness. These results may be important to understand the complexity of taste perception and their associated effects that could be useful to develop nutrigenetic approaches in pig breeding and nutrition.     

Evolution and selection of Rhg1, a copy‐number variant nematode‐resistance locus

The soybean cyst nematode (SCN) resistance locus Rhg1 is a tandem repeat of a 31.2 kb unit of the soybean genome. Each 31.2‐kb unit contains four genes. One allele of Rhg1, Rhg1‐b, is responsible for protecting most US soybean production from SCN. Whole‐genome sequencing was performed, and PCR assays were developed to investigate allelic variation in sequence and copy number of the Rhg1 locus across a population of soybean germplasm accessions. Four distinct sequences of the 31.2‐kb repeat unit were identified, and some Rhg1 alleles carry up to three different types of repeat unit. The total number of copies of the repeat varies from 1 to 10 per haploid genome. Both copy number and sequence of the repeat correlate with the resistance phenotype, and the Rhg1 locus shows strong signatures of selection. Significant linkage disequilibrium in the genome outside the boundaries of the repeat allowed the Rhg1 genotype to be inferred using high‐density single nucleotide polymorphism genotyping of 15 996 accessions. Over 860 germplasm accessions were found likely to possess Rhg1 alleles. The regions surrounding the repeat show indications of non‐neutral evolution and high genetic variability in populations from different geographic locations, but without evidence of fixation of the resistant genotype. A compelling explanation of these results is that balancing selection is in operation at Rhg1.     

Dynamic expression of N‐myc in mouse embryonic development using an enhanced green fluorescent protein reporter gene in the N‐myc locus

N‐myc belongs to the Myc oncogene family and plays an essential role in mammalian embryonic development. The expression of N‐myc is dynamically regulated during embryonic development; however, its expression pattern has not been well characterized due to the lack of a suitable animal model. In this paper, a genetically modified mouse model was generated in which the enhanced green fluorescent protein (EGFP) coding sequence was inserted into the N‐myc locus, so that endogenous N‐myc expression could be traced by the signal of EGFP. The EGFP signal in the transgenic mouse was confirmed to be consistent with the expression pattern of endogenous N‐myc by fluorescence microscopy and immunohistochemical staining. Furthermore, the spatial and temporal expression of EGFP was observed in the central and peripheral nervous system, heart, lung and kidney, given the known indispensable role of N‐myc in their formation. EGFP was also strongly detected in the liver, paranephros and the epithelium of the intestine. The EGFP signal can be used to trace N‐myc expression in this transgenic mouse model. N‐myc expression was observed in specific locations and cell lineages, and dynamically changed during embryonic development. The changing N‐myc expression pattern seen in mouse embryonic development and the animal model described in this paper provide important insights and a new tool to research N‐myc function.     

Development and evaluation of novel salt‐tolerant Eucalyptus trees by molecular breeding using an RNA‐Binding‐Protein gene derived from common ice plant (Mesembryanthemum crystallinum L.)

The breeding of plantation forestry trees for the possible afforestation of marginal land would be one approach to addressing global warming issues. Here, we developed novel transgenic Eucalyptus trees (Eucalyptus camaldulensis Dehnh.) harbouring an RNA‐Binding‐Protein (McRBP) gene derived from a halophyte plant, common ice plant (Mesembryanthemum crystallinum L.). We conducted screened‐house trials of the transgenic Eucalyptus using two different stringency salinity stress conditions to evaluate the plants’ acute and chronic salt stress tolerances. Treatment with 400 mM NaCl, as the high‐stringency salinity stress, resulted in soil electrical conductivity (EC) levels >20 mS/cm within 4 weeks. With the 400 mM NaCl treatment, >70% of the transgenic plants were intact, whereas >40% of the non‐transgenic plants were withered. Treatment with 70 mM NaCl, as the moderate‐stringency salinity stress, resulted in soil EC levels of approx. 9 mS/cm after 2 months, and these salinity levels were maintained for the next 4 months. All plants regardless of transgenic or non‐transgenic status survived the 70 mM NaCl treatment, but after 6‐month treatment the transgenic plants showed significantly higher growth and quantum yield of photosynthesis levels compared to the non‐transgenic plants. In addition, the salt accumulation in the leaves of the transgenic plants was 30% lower than that of non‐transgenic plants after 15‐week moderate salt stress treatment. There results suggest that McRBP expression in the transgenic Eucalyptus enhances their salt tolerance both acutely and chronically.     

Assessing forest vulnerability to climate warming using a process‐based model of tree growth: bad prospects for rear‐edges

Growth models can be used to assess forest vulnerability to climate warming. If global warming amplifies water deficit in drought‐prone areas, tree populations located at the driest and southernmost distribution limits (rear‐edges) should be particularly threatened. Here, we address these statements by analyzing and projecting growth responses to climate of three major tree species (silver fir, Abies alba; Scots pine, Pinus sylvestris; and mountain pine, Pinus uncinata) in mountainous areas of NE Spain. This region is subjected to Mediterranean continental conditions, it encompasses wide climatic, topographic and environmental gradients, and, more importantly, it includes rear‐edges of the continuous distributions of these tree species. We used tree‐ring width data from a network of 110 forests in combination with the process‐based Vaganov–Shashkin‐Lite growth model and climate–growth analyses to forecast changes in tree growth during the 21st century. Climatic projections were based on four ensembles CO₂ emission scenarios. Warm and dry conditions during the growing season constrain silver fir and Scots pine growth, particularly at the species rear‐edge. By contrast, growth of high‐elevation mountain pine forests is enhanced by climate warming. The emission scenario (RCP 8.5) corresponding to the most pronounced warming (+1.4 to 4.8 °C) forecasted mean growth reductions of ?10.7% and ?16.4% in silver fir and Scots pine, respectively, after 2050. This indicates that rising temperatures could amplify drought stress and thus constrain the growth of silver fir and Scots pine rear‐edge populations growing at xeric sites. Contrastingly, mountain pine growth is expected to increase by +12.5% due to a longer and warmer growing season. The projections of growth reduction in silver fir and Scots pine portend dieback and a contraction of their species distribution areas through potential local extinctions of the most vulnerable driest rear‐edge stands. Our modeling approach provides accessible tools to evaluate forest vulnerability to warmer conditions.     

miR‐148a inhibits early relapsed colorectal cancers and the secretion of VEGF by indirectly targeting HIF‐1α under non‐hypoxia/hypoxia conditions

Vascular endothelial growth factor (VEGF) is correlated with angiogenesis and early relapse of colorectal cancer (CRC). This study investigated the role of miR‐148a in the regulation of VEGF/angiogenesis and early relapse of CRC. We established a stable clone with miR‐148a expression in HCT116 and HT29 cell lines and created a hypoxic condition by using CoCl₂ to determine the underlying mechanism of miR‐148a. The effects of miR‐148a on the phosphoryl‐ERK (pERK)/hypoxia‐inducible factor‐1α (HIF‐1α)/VEGF pathway were evaluated through Western blotting and the inhibitory effect of miR‐148a on angiogenesis was demonstrated through a tube formation assay. Sixty‐three CRC tissues (28 early relapse and 35 non‐early relapse) were analysed to assess the relationship between miR‐148a and HIF‐1α/VEGF. The protein expression of pERK/HIF‐1α/VEGF in HCT116 and HT29 cells was significantly decreased by miR‐148a (all P < 0.05). The protein expression of VEGF/HIF‐1α was strongly inversely associated with the expression of miR‐148a in the 63 CRC tissue samples (all P < 0.05). Tube formation assay demonstrated that miR‐148a significantly obliterated angiogenesis. miR‐148a suppresses VEGF through down‐regulation of the pERK/HIF‐1α/VEGF pathway and might lead to the inhibition of angiogenesis; miR‐148a down‐regulation increased the early relapse rate of CRC. This demonstrates that miR‐148a is a potential diagnostic and therapeutic target.     

ABC inference of multi‐population divergence with admixture from unphased population genomic data

Rapidly developing sequencing technologies and declining costs have made it possible to collect genome‐scale data from population‐level samples in nonmodel systems. Inferential tools for historical demography given these data sets are, at present, underdeveloped. In particular, approximate Bayesian computation (ABC) has yet to be widely embraced by researchers generating these data. Here, we demonstrate the promise of ABC for analysis of the large data sets that are now attainable from nonmodel taxa through current genomic sequencing technologies. We develop and test an ABC framework for model selection and parameter estimation, given histories of three‐population divergence with admixture. We then explore different sampling regimes to illustrate how sampling more loci, longer loci or more individuals affects the quality of model selection and parameter estimation in this ABC framework. Our results show that inferences improved substantially with increases in the number and/or length of sequenced loci, while less benefit was gained by sampling large numbers of individuals. Optimal sampling strategies given our inferential models included at least 2000 loci, each approximately 2 kb in length, sampled from five diploid individuals per population, although specific strategies are model and question dependent. We tested our ABC approach through simulation‐based cross‐validations and illustrate its application using previously analysed data from the oak gall wasp, Biorhiza pallida.     

Phylogeography and diversification of an Amazonian understorey hummingbird: paraphyly and evidence for widespread cryptic speciation in the Plio‐Pleistocene

Straight‐billed Hermit Phaethornis bourcieri inhabits the understorey of upland terra firme forest throughout most of the Amazon basin. Currently, two allopatric taxa regarded as subspecies are recognized: P. b. bourcieri and P. b. major. However, the validity, interspecific limits and evolutionary history of these taxa are not yet fully elucidated. We use molecular characters to propose a phylogenetic hypothesis for populations and taxa grouped under Phaethornis bourcieri. Our results showed that P. bourcieri is part of the ‘Ametrornis’ clade, along with Phaethornis philippii and Phaethornis koepckeae, and that the subspecies major is more closely related to the latter two species than to populations grouped under nominate bourcieri. Our phylogenetic hypotheses recovered three main reciprocally monophyletic clades under nominate bourcieri separated by the lower Negro River and the Branco River or the Branco–Negro interfluve (clades B and C) and the upper Amazon (Solimões) or lower Marañon/Ucayali Rivers (clades C and D). Based on multi‐locus phylogeographic and population genetics approaches, we show that P. b. major is best treated as a separate species, and that P. b. bourcieri probably includes more than one evolutionary species, whose limits remain uncertain. The diversification of the ‘Ametrornis’ clade (P. bourcieri, P. philippii and P. koepckeae) is centred in the Amazon and appears to be closely linked to the formation of the modern Amazon drainage during the Plio‐Pleistocene.     

Process‐based models not always better than empirical models for simulating budburst of Norway spruce and birch in Europe

Budburst models have mainly been developed to capture the processes of individual trees, and vary in their complexity and plant physiological realism. We evaluated how well eleven models capture the variation in budburst of birch and Norway spruce in Germany, Austria, the United Kingdom and Finland. The comparison was based on the models performance in relation to their underlying physiological assumptions with four different calibration schemes. The models were not able to accurately simulate the timing of budburst. In general the models overestimated the temperature effect, thereby the timing of budburst was simulated too early in the United Kingdom and too late in Finland. Among the better performing models were three models based on the growing degree day concept, with or without day length or chilling, and an empirical model based on spring temperatures. These models were also the models least influenced by the calibration data. For birch the best calibration scheme was based on multiple sites in either Germany or Europe, and for Norway spruce the best scheme included multiple sites in Germany or cold years of all sites. Most model and calibration combinations indicated greater bias with higher spring temperatures, mostly simulating earlier than observed budburst.     

A single‐nucleotide polymorphism of miR‐146a and psoriasis: an association and functional study

Epidermal growth factor receptor (EGFR), which is overexpressed in psoriatic lesions, has been proven to contribute to the hyperproliferation of keratinocytes in psoriasis. Single nucleotide polymorphisms (SNPs) involved in miRNAs that can regulate the expression of EGFR could potentially influence the development of psoriasis. The present study investigated the association between a functional SNP of rs2910164 in miR‐146a and the risk of psoriasis in the Chinese Han population. A total of 521 Han Chinese patients with psoriasis and 582 healthy controls were recruited in this study. The miR‐146a rs2910164 SNP was genotyped by polymerase chain reaction‐restriction fragment length polymorphism. Overall, a significantly increased risk of psoriasis was associated with the rs2910164 miR‐146a CG and GG genotypes (adjusted OR, 1.38; 95% CI, 1.06–1.80). Furthermore, the rs2910164G allele in miR‐146a attenuated its inhibitory regulation on the expression of EGFR as well as the proliferation of human keratinocytes, and lowered the level of miR‐146a in the psoriatic lesions. These findings indicate that the rs2910164G allele in miR‐146a weakens its suppression on the proliferation of keratinocytes probably through the decreased inhibition of the target gene, EGFR, which may account for the increased risk of psoriasis in this study population.     

Empirical calibrated radiocarbon sampler: a tool for incorporating radiocarbon‐date and calibration error into Bayesian phylogenetic analyses of ancient DNA

Studies of DNA from ancient samples provide a valuable opportunity to gain insight into past evolutionary and demographic processes. Bayesian phylogenetic methods can estimate evolutionary rates and timescales from ancient DNA sequences, with the ages of the samples acting as calibrations for the molecular clock. Sample ages are often estimated using radiocarbon dating, but the associated measurement error is rarely taken into account. In addition, the total uncertainty quantified by converting radiocarbon dates to calendar dates is typically ignored. Here, we present a tool for incorporating both of these sources of uncertainty into Bayesian phylogenetic analyses of ancient DNA. This empirical calibrated radiocarbon sampler (ECRS) integrates the age uncertainty for each ancient sequence over the calibrated probability density function estimated for its radiocarbon date and associated error. We use the ECRS to analyse three ancient DNA data sets. Accounting for radiocarbon‐dating and calibration error appeared to have little impact on estimates of evolutionary rates and related parameters for these data sets. However, analyses of other data sets, particularly those with few or only very old radiocarbon dates, might be more sensitive to using artificially precise sample ages and should benefit from use of the ECRS.