Genetic diversity in oak populations under intensive management for fuelwood in the Sierra de Zongolica,Mexico

Firewood and charcoal are used on a daily basis both in rural areas and in cities. This type of energy is produced by one of the most ancient traditional methods, known as coppice, which harvest tree sprouts. There is controversy about its effects on forests: it preserves populations and tree cover of species used, but reduces density, inhibits sexual reproduction and generates genetic erosion. We inquired if it was possible to identify a loss of genetic diversity in oak populations traditionally used for charcoal by the Zongolica Nahuas in Veracruz state, Mexico. We studied populations of Quercus laurina, Quercus calophylla and Quercus rugosa in three different altitudes. Molecular analysis with eight nuclear codominant microsatellites was performed to determine the diversity, structure and gene flow of these species. Results for Q. laurina were Na = 8.458, I = 1.766, Ho = 0.679, polymorphism = 100%, Fis = 0.079, with intraindividual variation of 81.55%. For Q. calophylla: Na = 7.250, I = 1.563, Ho = 0.646, polymorphism = 91.67%, Fis = 0.083, with intraindividual variation of 83.80%. For Q. rugosa: Na = 6.958, I = 1.510, Ho = 0.574, polymorphism = 91.67%, Fis = 0.204, with intraindividual variation of 81.99%; this species shows signals of an early genetic isolation process. Our findings indicate that Quercus genetic diversity for the three species is high and comparable with oak species in Mexico and worldwide. We conclude that at the present, coppice is preserving a historical diversity in adult trees kept alive through sprouting. Nonetheless, problems with coppice systems elsewhere, unregulated harvesting and expansion of pine plantation in the region suggest that further studies, hand in hand with a landscape management approach that improve charcoal and firewood production, may be valuable for Sierra de Zongolica genetic biodiversity conservation.     

On computation of semiparametric maximum likelihood estimators with shape constraints

Large sample theory of semiparametric models based on maximum likelihood estimation (MLE) with shape constraint on the nonparametric component is well studied. Relatively less attention has been paid to the computational aspect of semiparametric MLE. The computation of semiparametric MLE based on existing approaches such as the expectation‐maximization (EM) algorithm can be computationally prohibitive when the missing rate is high. In this paper, we propose a computational framework for semiparametric MLE based on an inexact block coordinate ascent (BCA) algorithm. We show theoretically that the proposed algorithm converges. This computational framework can be applied to a wide range of data with different structures, such as panel count data, interval‐censored data, and degradation data, among others. Simulation studies demonstrate favorable performance compared with existing algorithms in terms of accuracy and speed. Two data sets are used to illustrate the proposed computational method. We further implement the proposed computational method in R package BCA1SG , available at CRAN.     

A penalized structural equation modeling method accounting for secondary phenotypes for variable selection on genetically regulated expression from PrediXcan for Alzheimer's disease

As the global burden of mental illness is estimated to become a severe issue in the near future, it demands the development of more effective treatments. Most psychiatric diseases are moderately to highly heritable and believed to involve many genes. Development of new treatment options demands more knowledge on the molecular basis of psychiatric diseases. Toward this end, we propose to develop new statistical methods with improved sensitivity and accuracy to identify disease‐related genes specialized for psychiatric diseases. The qualitative psychiatric diagnoses such as case control often suffer from high rates of misdiagnosis and oversimplify the disease phenotypes. Our proposed method utilizes endophenotypes, the quantitative traits hypothesized to underlie disease syndromes, to better characterize the heterogeneous phenotypes of psychiatric diseases. We employ the structural equation modeling using the liability‐index model to link multiple genetically regulated expressions from PrediXcan and the manifest variables including endophenotypes and case‐control status. The proposed method can be considered as a general method for multivariate regression, which is particularly helpful for psychiatric diseases. We derive penalized retrospective likelihood estimators to deal with the typical small sample size issue. Simulation results demonstrate the advantages of the proposed method and the real data analysis of Alzheimer's disease illustrates the practical utility of the techniques. Data used in preparation of this article were obtained from the Alzheimer's Disease Neuroimaging Initiative database.     

The interaction between host and host plant influences the oviposition and performance of a generalist ectoparasitoid

The successful development of parasitoids of herbivores depends on the quality of their host, which is often affected by the host plant. Therefore, a parasitoid’s oviposition decisions will directly depend on the host, but also on plant quality. Here, we investigated the direct effects of host species and the indirect effects of the host’s food plant on the oviposition decisions and performance of the gregarious ectoparasitoid Euplectrus platyhypenae Howard (Hymenoptera: Eulophidae). With a series of no‐choice experiments, we determined the oviposition and performance of the parasitoid on: (1) two caterpillar species, fall armyworm, Spodoptera frugiperda JE Smith (Lepidoptera: Noctuidae), and velvet armyworm, Spodoptera latifascia Walker, reared on maize (Zea mays L., Poaceae), (2) the same caterpillars reared on maize, bean (Phaseolus vulgaris L., Fabaceae), or squash (Cucurbita pepo L., Cucurbitaceae) leaves, and (3) S. latifascia caterpillars reared on leaves of wild and cultivated lima bean, Phaseolus lunatus L. All these insects and plants originate from Mesoamerica where they have coexisted for thousands of years in the traditional agricultural system known as Milpa in which maize, beans, and squash are planted together. We found that the preferred and best combination of host and host plant for parasitoid performance was S. frugiperda on maize. Parasitoids laid larger clutches, had higher survival, and more females and larger adults emerged from S. frugiperda reared on maize. However, when both caterpillar species were reared on squash, S. latifascia was the better host. Contrary to the literature, S. frugiperda was not able to develop on bean plants. Results from the lima bean experiment showed that parasitoid performance was best when S. latifascia was reared on leaves of cultivated compared to wild lima bean. These findings are discussed in the context of mixed cropping in which the ability of generalist parasitoids to switch among hosts and host plant species could be advantageous for pest management.     

Characterization of binding kinetics of A2AR to Gαs protein by surface plasmon resonance

Because of their surface localization, G protein-coupled receptors (GPCRs) are often pharmaceutical targets as they respond to a variety of extracellular stimuli (e.g., light, hormones, small molecules) that may activate or inhibit a downstream signaling response. The adenosine A_2A receptor (A_2AR) is a well-characterized GPCR that is expressed widely throughout the human body, with over 10 crystal structures determined. Truncation of the A_2AR C-terminus is necessary for crystallization as this portion of the receptor is long and unstructured; however, previous work suggests shortening of the A_2AR C-terminus from 412 to 316 amino acids (A_2AΔ316R) ablates downstream signaling, as measured by cAMP production, to below that of constitutive full-length A_2AR levels. As cAMP production is downstream of the first activation event—coupling of G protein to its receptor—investigating that first step in activation is important in understanding how the truncation effects native GPCR function. Here, using purified receptor and Gα_s proteins, we characterize the association of A_2AR and A_2AΔ316R to Gα_s with and without GDP or GTPγs using surface plasmon resonance (SPR). Gα_s affinity for A_2AR was greatest for apo-Gα_s, moderately affected in the presence of GDP and nearly completely ablated by the addition of GTPγs. Truncation of the A_2AR C-terminus (A_2AΔ316R) decreased the affinity of the unliganded receptor for Gα_s by ～20%, suggesting small changes to binding can greatly impact downstream signaling.     

Linking Physiology To Ecological Function: Environmental Conditions Affect Performance And Size Of The Intertidal Kelp Hedophyllum Sessile (Laminariales,Phaeophyceae)

For autogenic ecosystem engineers, body size is an aspect of individual performance that has direct connections to community structure; yet the complex morphology of these species can make it difficult to draw clear connections between the environment and performance. We combined laboratory experiments and field surveys to test the hypothesis that individual body size was determined by disparate localized physiological responses to environmental conditions across the complex thallus of the intertidal kelp Hedophyllum sessile, a canopy‐forming physical ecosystem engineer. We documented substantial (> 40%) declines in whole‐thallus photosynthetic potential (as Maximum Quantum Yield, MQY) as a consequence of emersion, which were related to greater than 10‐fold increases in intra‐thallus MQY variability (as Coefficient of Variation). In laboratory experiments, desiccation and high light levels during emersion led to lasting impairment of photosynthetic potential and an immediate > 25% reduction in area due to tissue contraction, which was followed by complete loss of structural integrity after three days of submersion. Tissue exposed to desiccation and high light during emersion had higher nitrogen concentrations and lower phlorotannin concentrations than tissue in control treatments (on average 1.36 and 0.1x controls, respectively), suggesting that conditions during emersion have the potential to affect food quality for consumers. Our data indicate that the complex thallus morphology of H. sessile may be critical to this kelp’s ability to persist in the intertidal zone despite the physiological challenges of emersion and encourage a more nuanced view of the concept of “sub‐lethal stress” on the scale of the whole individual.     

Using RAD‐seq to develop sex‐linked markers in a haplodiplontic alga

For many taxa, including isomorphic haplodiplontic macroalgae, determining sex and ploidy is challenging, thereby limiting the scope of some population demographic and genetic studies. Here, we used double‐digest restriction site‐associated DNA sequencing (ddRAD‐seq) to identify sex‐linked molecular markers in the widespread red alga Agarophyton vermiculophyllum. In the ddRAD‐seq library, we included 10 female gametophytes, 10 male gametophytes, and 16 tetrasporophytes from one native and one non‐native site (N = 40 gametophytes and N = 32 tetrasporophytes total). We identified seven putatively female‐linked and 19 putatively male‐linked sequences. Four female‐ and eight male‐linked markers amplified in all three life cycle stages. Using one female‐ and one male‐linked marker that were sex‐specific, we developed a duplex PCR and tested the efficacy of this assay on a subset of thalli sampled at two sites in the non‐native range. We confirmed ploidy based on the visual observation of reproductive structures and previous microsatellite genotyping at 10 polymorphic loci. For 32 vegetative thalli, we were able to assign sex and confirm ploidy in these previously genotyped thalli. These markers will be integral to ongoing studies of A. vermiculophyllum invasion. We discuss the utility of RAD‐seq over other approaches previously used, such as RAPDs (random amplified polymorphic DNA), for future work designing sex‐linked markers in other haplodiplontic macroalgae for which genomes are lacking.