Efficient Bayesian Sample Size Calculation for Designing a Clinical Trial with Multi‐Cluster Outcome Data

Health care utilization and outcome studies call for hierarchical approaches. The objectives were to predict major complications following percutaneous coronary interventions by health providers, and to compare Bayesian and non‐Bayesian sample size calculation methods. The hierarchical data structure consisted of: (1) Strata: PGY4, PGY7, and physician assistant as providers with varied experiences; (2) Clusters: k_s providers per stratum; (3) Individuals: n_s patients reviewed by each provider. The main outcome event illustrated was mortality modeled by a Bayesian beta‐binomial model. Pilot information and assumptions were utilized to elicit beta prior distributions. Sample size calculations were based on the approximated average length, fixed at 1%, of 95% posterior intervals of the mean event rate parameter. Necessary sample sizes by both non‐Bayesian and Bayesian methods were compared. We demonstrated that the developed Bayesian methods can be efficient and may require fewer subjects to satisfy the same length criterion.     

Inference for meta‐analysis with a suspected temporal trend

There is sometimes a clear evidence of a strong secular trend in the treatment effect of studies included in a meta‐analysis. In such cases, estimating the present‐day treatment effect by meta‐regression is both reasonable and straightforward. We however consider the more common situation where a secular trend is suspected, but is not strongly statistically significant. Typically, this lack of significance is due to the small number of studies included in the analysis, so that a meta‐regression could give wild point estimates. We introduce an empirical Bayes meta‐analysis methodology, which shrinks the secular trend toward zero. This has the effect that treatment effects are adjusted for trend, but where the evidence from data is weak, wild results are not obtained. We explore several frequentist approaches and a fully Bayesian method is also implemented. A measure of trend analogous to I² is described, and exact significance tests for trend are given. Our preferred method is one based on penalized or h‐likelihood, which is computationally simple, and allows invariance of predictions to the (arbitrary) choice of time origin. We suggest that a trendless standard random effects meta‐analysis should routinely be supplemented with an h‐likelihood analysis as a sensitivity analysis.     

Structures of reduced and ligand‐bound nitric oxide reductase provide insights into functional differences in respiratory enzymes

Nitric oxide reductase (NOR) catalyzes the generation of nitrous oxide (N₂O) via the reductive coupling of two nitric oxide (NO) molecules at a heme/non‐heme Fe center. We report herein on the structures of the reduced and ligand‐bound forms of cytochrome c‐dependent NOR (cNOR) from Pseudomonas aeruginosa at a resolution of 2.3–2.7 Å, to elucidate structure‐function relationships in NOR, and compare them to those of cytochrome c oxidase (CCO) that is evolutionarily related to NOR. Comprehensive crystallographic refinement of the CO‐bound form of cNOR suggested that a total of four atoms can be accommodated at the binuclear center. Consistent with this, binding of bulky acetaldoxime (CH₃‐CH=N‐OH) to the binuclear center of cNOR was confirmed by the structural analysis. Active site reduction and ligand binding in cNOR induced only ～0.5 Å increase in the heme/non‐heme Fe distance, but no significant structural change in the protein. The highly localized structural change is consistent with the lack of proton‐pumping activity in cNOR, because redox‐coupled conformational changes are thought to be crucial for proton pumping in CCO. It also permits the rapid decomposition of cytotoxic NO in denitrification. In addition, the shorter heme/non‐heme Fe distance even in the bulky ligand‐bound form of cNOR (～4.5 Å) than the heme/Cu distance in CCO (～5 Å) suggests the ability of NOR to maintain two NO molecules within a short distance in the confined space of the active site, thereby facilitating N‐N coupling to produce a hyponitrite intermediate for the generation of N₂O. Proteins 2014; 82:1258–1271. © 2013 Wiley Periodicals, Inc.     

High‐performance liquid chromatographic enantioseparation of N‐aryl‐thioureidoalkylphosphonates and thiourylenedi(alkylphosphonates) on polysaccharide‐based chiral stationary phases

The first successful enantioseparation of representative O,O‐diphenyl‐N‐arylthioureidoalkylphosphonates, (±)‐Ptc‐Val^P(OPh)₂ & (±)‐Ptc‐Leu^P(OPh)₂ and thiourylenedi(isobutyl phosphonate), Tcm[Val^P(OPh)₂]₂ on analytical and semipreparative scale was achieved by high‐performance liquid chromatography using polysaccharide‐based chiral stationary phases (CPs). Atc‐AA^P(OPh)₂ was obtained using modified tricomponent condensations of the corresponding aldehydes, N‐arylthiourea and triphenyl phosphite whereas Tcm[Val^P(OPh)₂]₂ by the condensations of aldehydes, thiourea, and triphenyl phosphite. The prepared, racemic (±)‐Atc‐AA^P(OPh)₂ [(±)‐Ptc‐Val^P(OPh)₂, (±)‐Ptc‐Leu^P(OPh)₂, (±)‐Ptc‐Pgly^P(OPh)₂ and (±)‐Ntc‐Pgly^P(OPh)₂] and racemic (±)‐Tcm[AA^P(OPh)₂]₂ [(±)‐Tcm[Nva^P(OPh)₂]₂ & (±)‐Tcm[Val^P(OPh)₂]₂] were adequately characterized and used for chromatographic separations on high‐performance liquid chromatography–chiral stationary phases. The best results were obtained for (±)‐Ptc‐Val^P(OPh)₂, (±)‐Ptc‐Leu^P(OPh)₂ and (±)‐Tcm[Val^P(OPh)₂]₂.     

WFABC: a Wright–Fisher ABC‐based approach for inferring effective population sizes and selection coefficients from time‐sampled data

With novel developments in sequencing technologies, time‐sampled data are becoming more available and accessible. Naturally, there have been efforts in parallel to infer population genetic parameters from these data sets. Here, we compare and analyse four recent approaches based on the Wright–Fisher model for inferring selection coefficients (s) given effective population size (N_e), with simulated temporal data sets. Furthermore, we demonstrate the advantage of a recently proposed approximate Bayesian computation (ABC)‐based method that is able to correctly infer genomewide average N_e from time‐serial data, which is then set as a prior for inferring per‐site selection coefficients accurately and precisely. We implement this ABC method in a new software and apply it to a classical time‐serial data set of the medionigra genotype in the moth Panaxia dominula. We show that a recessive lethal model is the best explanation for the observed variation in allele frequency by implementing an estimator of the dominance ratio (h).     

Cumulative Damage Survival Models for the Randomized Complete Block Design

A continuous time discrete state cumulative damage process {X(t), t ≥ 0} is considered, based on a non‐homogeneous Poisson hit‐count process and discrete distribution of damage per hit, which can be negative binomial, Neyman type A, Polya‐Aeppli or Lagrangian Poisson. Intensity functions considered for the Poisson process comprise a flexible three‐parameter family. The survival function is S(t) = P(X(t) ≤ L) where L is fixed. Individual variation is accounted for within the construction for the initial damage distribution {P(X(0) = x) | x = 0, 1, …,}. This distribution has an essential cut‐off before x = L and the distribution of L – X(0) may be considered a tolerance distribution. A multivariate extension appropriate for the randomized complete block design is developed by constructing dependence in the initial damage distributions. Our multivariate model is applied (via maximum likelihood) to litter‐matched tumorigenesis data for rats. The litter effect accounts for 5.9 percent of the variance of the individual effect. Cumulative damage hazard functions are compared to nonparametric hazard functions and to hazard functions obtained from the PVF‐Weibull frailty model. The cumulative damage model has greater dimensionality for interpretation compared to other models, owing principally to the intensity function part of the model.     

Investigation of a sugar N‐formyltransferase from the plant pathogen Pantoea ananatis

Pantoea ananatis is a Gram‐negative bacterium first recognized in 1928 as the causative agent of pineapple rot in the Philippines. Since then various strains of the organism have been implicated in the devastation of agriculturally important crops. Some strains, however, have been shown to function as non‐pathogenic plant growth promoting organisms. To date, the factors that determine pathogenicity or lack thereof between the various strains are not well understood. All P. ananatis strains contain lipopolysaccharides, which differ with respect to the identities of their associated sugars. Given our research interest on the presence of the unusual sugar, 4‐formamido‐4,6‐dideoxy‐d ‐glucose, found on the lipopolysaccharides of Campylobacter jejuni and Francisella tularensis, we were curious as to whether other bacteria have the appropriate biosynthetic machinery to produce these unique carbohydrates. Four enzymes are typically required for their biosynthesis: a thymidylyltransferase, a 4,6‐dehydratase, an aminotransferase, and an N‐formyltransferase. Here, we report that the gene SAMN03097714_1080 from the P. ananatis strain NFR11 does, indeed, encode for an N‐formyltransferase, hereafter referred to as PA1080c. Our kinetic analysis demonstrates that PA1080c displays classical Michaelis–Menten kinetics with dTDP‐4‐amino‐4,6‐dideoxy‐d ‐glucose as the substrate and N¹⁰‐formyltetrahydrofolate as the carbon source. In addition, the X‐ray structure of PA1080c, determined to 1.7 Å resolution, shows that the enzyme adopts the molecular architecture observed for other sugar N‐formyltransferases. Analysis of the P. ananatis NFR11 genome suggests that the three other enzymes necessary for N‐formylated sugar biosynthesis are also present. Intriguingly, those strains of P. ananatis that are non‐pathogenic apparently do not contain these genes.     

Nitrous oxide emissions during establishment of eight alternative cellulosic bioenergy cropping systems in the North Central United States

Greenhouse gas (GHG) emissions from soils are a key sustainability metric of cropping systems. During crop establishment, disruptive land‐use change is known to be a critical, but under reported period, for determining GHG emissions. We measured soil N₂O emissions and potential environmental drivers of these fluxes from a three‐year establishment‐phase bioenergy cropping systems experiment replicated in southcentral Wisconsin (ARL) and southwestern Michigan (KBS). Cropping systems treatments were annual monocultures (continuous corn, corn–soybean–canola rotation), perennial monocultures (switchgrass, miscanthus, and poplar), and perennial polycultures (native grass mixture, early successional community, and restored prairie) all grown using best management practices specific to the system. Cumulative three‐year N₂O emissions from annuals were 142% higher than from perennials, with fertilized perennials 190% higher than unfertilized perennials. Emissions ranged from 3.1 to 19.1 kg N₂O‐N ha^?1 yr^?1 for the annuals with continuous corn > corn–soybean–canola rotation and 1.1 to 6.3 kg N₂O‐N ha^?1 yr^?1 for perennials. Nitrous oxide peak fluxes typically were associated with precipitation events that closely followed fertilization. Bayesian modeling of N₂O fluxes based on measured environmental factors explained 33% of variability across all systems. Models trained on single systems performed well in most monocultures (e.g., R² = 0.52 for poplar) but notably worse in polycultures (e.g., R² = 0.17 for early successional, R² = 0.06 for restored prairie), indicating that simulation models that include N₂O emissions should be parameterized specific to particular plant communities. Our results indicate that perennial bioenergy crops in their establishment phase emit less N₂O than annual crops, especially when not fertilized. These findings should be considered further alongside yield and other metrics contributing to important ecosystem services.     

Efficient super‐resolution volumetric imaging by radial fluctuation Bayesian analysis light‐sheet microscopy

Various computational super‐resolution methods are available based on the analysis of fluorescence fluctuation behind acquired frames. However, dilemmas often exist in the balance of fluorophore characteristics, computation cost, and achievable resolution. Here we present an approach that uses a super‐resolution radial fluctuations (SRRF) image to guide the Bayesian analysis of fluorophore blinking and bleaching (3B) events, allowing greatly accelerated localization of overlapping fluorophores with high accuracy. This radial fluctuation Bayesian analysis (RFBA) approach is also extended to three dimensions for the first time and combined with light‐sheet fluorescence microscopy, to achieve super‐resolution volumetric imaging of thick samples densely labeled with common fluorophores. For example, a 700‐nm thin Bessel plane illumination is developed to optically section the Drosophila brain, providing a high‐contrast 3D image of rhythmic neurons. RFBA analyzes 30 serial volumes to reconstruct a super‐resolved 3D image at 4‐times higher resolutions (~70 and 170 nm), and precisely resolve the axon terminals. The computation is over 2‐orders faster than conventional 3B analysis microscopy. The capability of RFBA is also verified through dual‐color imaging of cell nucleus in live Drosophila brain. The spatial co‐localization patterns of the nuclear envelope and DNA in a neuron deep inside the brain can be precisely extracted by our approach.     

Melatonin Levels and Low‐Frequency Magnetic Fields in Humans and Rats: New Insights From a Bayesian Logistic Regression

The present analysis revisits the impact of extremely low‐frequency magnetic fields (ELF‐MF) on melatonin (MLT) levels in human and rat subjects using both a parametric and non‐parametric approach. In this analysis, we use 62 studies from review articles. The parametric approach consists of a Bayesian logistic regression (LR) analysis and the non‐parametric approach consists of a Support Vector analysis, both of which are robust against spurious/false results. Both approaches reveal a unique well‐ordered pattern, and show that human and rat studies are consistent with each other once the MF strength is restricted to cover the same range (with B ? 50 μT). In addition, the data reveal that chronic exposure (longer than ～22 days) to ELF‐MF appears to decrease MLT levels only when the MF strength is below a threshold of ~30 μT (), i.e., when the man‐made ELF‐MF intensity is below that of the static geomagnetic field. Studies reporting an association between ELF‐MF and changes to MLT levels and the opposite (no association with ELF‐MF) can be reconciled under a single framework. Bioelectromagnetics. 2019;40:539–552. © 2019 Bioelectromagnetics Society.     

A comparison of single‐sample effective size estimators using empirical toad (Bufo calamita) population data: genetic compensation and population size‐genetic diversity correlations

The accuracy and precision of four single‐sample estimators of effective population size, N_e (heterozygote excess, linkage disequilibrium, Bayesian partial likelihood and sibship analysis) were compared using empirical data (microsatellite genotypes) from multiple natterjack toad (Bufo calamita) populations in Britain (n = 16) and elsewhere in Europe (n = 10). Census size data were available for the British populations. Because toads have overlapping generations, all of these methods estimated the number of effective breeders N_b rather than N_e. The heterozygote excess method only provided results, without confidence limits, for nine of the British populations. Linkage disequilibrium gave estimates for 10 British populations, but only six had finite confidence limits. The Bayesian and sibship methods both produced estimates with finite confidence limits for all the populations. Although the Bayesian method was the most precise, on most criteria (insensitivity to locus number, correlation with other effective and census size estimates and correlation with genetic diversity) the sibship method performed best. The results also provided evidence of genetic compensation in natterjack toads, and highlighted how the relationship between effective size and genetic diversity can vary as a function of geographical scale.     

Electronic Effect of Substituents on the DNA Intercalation of Ruthenium(II)?Polypyridyl Complexes

Five ruthenium(II) complexes, i.e., [Ru(bpy)₂(TIP)]²⁺ (bpy=2,2′‐bipyridine; TIP=2‐thiophenimidazo[4,5‐f] [1,10]phenanthroline; 1 ), [Ru(bpy)₂(5‐NTIP)]²⁺ (5‐NTIP=2‐(5‐nitrothiophen)imidazo[4,5‐f] [1,10]phenanthroline; 2 ), [Ru(bpy)₂(5‐MOTIP)]²⁺ (5‐MOTIP=2‐(5‐methoxythiophen)imidazo[4,5‐f] [1,10]phenanthroline; 3 ), [Ru(bpy)₂(5‐BTIP)]²⁺ (5‐BTIP=2‐(5‐bromothiophen)imidazo[4,5‐f] [1,10]phenanthroline; 4 ), and [Ru(bpy)₂(4‐BTIP)]²⁺ (4‐BTIP=2‐(4‐bromothiophen)imidazo[4,5‐f] [1,10]phenanthroline; 5 ), were synthesized and characterized by elemental analysis and UV/VIS, IR, and ¹H‐NMR spectroscopic methods. The photophysical and DNA‐binding properties were investigated by means of UV and fluorescence spectroscopic methods and viscosity measurements, respectively. The results suggest that all five complexes can bind to CT‐DNA with various binding strength. Complexes 2 and 3 showed the strongest and the weakest binding affinity, respectively, among these five complexes. Due to the substituent position of the Br‐atom in the ligand, complex 5 interacted stronger with CT‐DNA than complex 4 . The binding affinities of the complexes decreased in the order 2, 5, 4, 1 , and 3 .     

Using stable hydrogen isotopes (δ2H) and ring recoveries to trace natal origins in a Eurasian passerine with a migratory divide

Despite recent advances in technology, it remains difficult to connect breeding and non‐breeding areas of populations of migratory organisms due to the challenges of year‐round tracking. Here, we used the Eurasian reed warbler Acrocephalus scirpaceus, a passerine with a pronounced migratory divide to demonstrate the promise of integrating several sources of information within the Bayesian modelling framework for the study of migratory connectivity. To this end, we combined data from stable hydrogen isotope ratios (δ²H) of feathers, ring recoveries, and the geographic delineation of sub‐populations on either side of the migratory divide. Feather δ²H measurements from local juvenile birds sampled across the breeding range tightly correlated with amount‐weighted mean annual precipitation δ²H values predicted for the natal sites. Predicted natal origins of birds intercepted en route in the Mediterranean region largely differed among the five stopover sites. Thanks to the different migratory pathways used by different breeding populations and the existence of a migratory divide, we were able to effectively narrow the assigned regions of origin. Our results show that spatial resolution of likelihood‐based assignments of geographic origins based on δ²H measurements may improve significantly when prior probabilities derived from population‐specific migratory directions are included. Integrating information from stable isotopes, ring recoveries, geolocators and other sources within the Bayesian modelling framework will provide an extremely useful toolbox for the study of animal movements in the future.     

Phoronid phylogenetics (Brachiopoda; Phoronata): evidence from morphological cladistics,small and large subunit rDNA sequences,and mitochondrial cox1

A matrix of 24 morphodevelopmental characters and an alignment of small subunit (SSU) and large subunit (LSU) rDNA nuclear and cox1 mitochondrial gene sequences (～4500 sites) were compiled from up to 12 phoronids including most named taxa, but probably constituting only a portion of worldwide diversity. Morphological data were analysed by weighted parsimony; sequence data by maximum and Bayesian likelihood, both with Phoronis ovalis as the local outgroup. Morphological and sequence‐based phylogenies were similar, but not fully congruent. Phoronid rDNAs were almost free from mutational saturation, but cox1 showed strong saturation unless distant outgroups and P. ovalis were omitted, suggesting that many phoronid divergences are old (≥100 Myr). rDNA divergence between named phoronid taxa is generally substantial, but Phoronopsis harmeri (from Vladivostock) and Phoronopsis viridis (from California) are genetically close enough to be conspecific. In another alignment, of 24 taxa, phoronid rDNAs were combined with data from brachiopods and distant (molluscan) outgroups. The relative ages of divergence between phoronids and their brachiopod sister‐groups, of the split between the P. ovalis and non‐ovalis lineages, and of other phoronid splits, were estimated from this alignment with a Bayesian lognormal uncorrelated molecular clock model. Although confidence limits (95% highest probability density) are wide, the results are compatible with an Early Cambrian split between phoronids and brachiopods and with the Upper Devonian latest age suggested for the P. ovalis/non‐ovalis split by the putative phoronid ichnofossil, Talpina. Most other ingroup splits appear to be ～50–200 Myr old. Inclusion of phoronids with brachiopods in the crown clade pan‐Brachiopoda suggests that a distinctive metamorphosis and absence of mineralization are ancestral phoronid apomorphies. Worldwide diversity and possible associations between character‐states and life‐history attributes deserve comprehensive further study.     

Epigenetic patterns within the haplotype phased fig (Ficus carica L.) genome

Due to DNA heterozygosity and repeat content, assembly of non‐model plant genomes is challenging. Herein, we report a high‐quality genome reference of one of the oldest known domesticated species, fig (Ficus carica L.), using Pacific Biosciences single‐molecule, real‐time sequencing. The fig genome is ~333 Mbp in size, of which 80% has been anchored to 13 chromosomes. Genome‐wide analysis of N⁶‐methyladenine and N⁴‐methylcytosine revealed high methylation levels in both genes and transposable elements, and a prevalence of methylated over non‐methylated genes. Furthermore, the characterization of N⁶‐methyladenine sites led to the identification of ANHGA, a species‐specific motif, which is prevalent for both genes and transposable elements. Finally, exploiting the contiguity of the 13 pseudomolecules, we identified 13 putative centromeric regions. The high‐quality reference genome and the characterization of methylation profiles, provides an important resource for both fig breeding and for fundamental research into the relationship between epigenetic changes and phenotype, using fig as a model species.     

New insight into the identification and molecular phylogeny of dagger nematodes of the genus Xiphinema (Nematoda: Longidoridae) with description of two new species

The genus Xiphinema constitutes a large group of about 260 species of plant‐ectoparasitic nematodes. The group is polyphagous and distributed almost worldwide. Some of the species of this genus damage agricultural crops by direct feeding on root cells as well as by transmitting nepoviruses. Species discrimination in Xiphinema is complicated by phenotypic plasticity leading to potential misidentification. We conducted nematode surveys in cultivated and natural environments in Spain from 2009 to 2012, from which we identified 20 populations of Xiphinema species morphologically close to the virus‐vector nematode species Xiphinema diversicaudatum, three apomictic populations tentatively identified as species from the complex Xiphinema aceri‐pyrenaicum group, and one population morphologically different from all others that is characterized by a female tail elongate to conical and absence of uterine differentiation. We developed comparative multivariate analyses for these related species by using morphological and morphometrical features together with molecular data from nuclear ribosomal DNA genes [D2‐D3 expansion segments of large ribosomal subunit 28S, internal transcribed spacer 1 (ITS1), and partial small ribosomal subunit (18S)]. The results of multivariate, molecular, and phylogenetic analysis confirmed the morphological hypotheses and allowed the delimitation and discrimination of two new species in the genus described herein as Xiphinema baetica sp. nov. and Xiphinema turdetanensis sp. nov. , and ten known species: Xiphinema adenohystherum, Xiphinema belmontense, Xiphinema cohni, Xiphinema coxi europaeum, Xiphinema gersoni, Xiphinema hispidum, Xiphinema italiae, Xiphinema lupini, Xiphinema nuragicum, and Xiphinema turcicum. Multivariate analyses based on quantitative and qualitative characters and phylogenetic relationships of Xiphinema spp. based on the three molecular ribosomal markers resulted in a partial consensus of these species grouping as nematode populations were maintained for the majority of morphospecies groups (e.g. morphospecies groups 5 and 6), but not in some others (e.g. position of Xiphinema granatum), demonstrating the usefulness of these analyses for helping in the diagnosis and identification of Xiphinema spp. The clade topology of phylogenetic trees of D2‐D3 and partial 18S regions in this study were congruent in supporting the polyphyletic status of some characters, such as the female tail shape and the degree of development of the genital system in species with both genital branches equally developed. This is the most complete phylogenetic study for Xiphinema non‐americanum‐group species. Agreement between phylogenetic trees and some morphological characters (uterine spines, pseudo‐Z organ, and tail shape) was tested by reconstruction of their histories on rDNA‐based trees using parsimony and Bayesian approaches. Thus, integrative taxonomy, based on the combination of multivariate, molecular analyses with morphology, constitutes a new insight into the identification of Xiphinema species. © 2013 The Linnean Society of London     

Micro‐ and macroheterogeneity of N‐glycosylation yields size and charge isoforms of human sex hormone binding globulin circulating in serum

Human sex hormone binding globulin (hSHBG) is a serum glycoprotein central to the transport and targeted delivery of sex hormones to steroid‐sensitive tissues. Several molecular mechanisms of action of hSHBG, including the function of its attached glycans remain unknown. Here, we perform a detailed site‐specific characterization of the N‐ and O‐linked glycosylation of serum‐derived hSHBG. MS‐driven glycoproteomics and glycomics combined with exoglycosidase treatment were used in a bottom‐up and top‐down manner to determine glycosylation sites, site‐specific occupancies and monosaccharide compositions, detailed glycan structures, and the higher level arrangement of glycans on intact hSHBG. It was found that serum‐derived hSHBG is N‐glycosylated at Asn³⁵¹ and Asn³⁶⁷ with average molar occupancies of 85.1 and 95.3%, respectively. Both sites are occupied by the same six sialylated and partly core fucosylated bi‐ and triantennary N‐Glycoforms with lactosamine‐type antennas of the form (±NeuAcα6)Galβ4GlcNAc. N‐Glycoforms of Asn³⁶⁷ were slightly more branched and core fucosylated than Asn³⁵¹ N‐glycoforms due probably to a more surface‐exposed glycosylation site. The N‐terminal Thr⁷ was fully occupied by the two O‐linked glycans NeuAcα3Galβ3(NeuAcα6)GalNAc (where NeuAc is N‐acetylneuraminic acid and GalNAc is N‐acetylgalactosamine) and NeuAcα3Galβ3GalNAc in a 1:6 molar ratio. Electrophoretic analysis of intact hSHBG revealed size and charge heterogeneity of the isoforms circulating in blood serum. Interestingly, the size and charge heterogeneity were shown to originate predominantly from differential Asn³⁵¹ glycan occupancies and N‐glycan sialylation that may modulate the hSHBG activity. To date, this work represents the most detailed structural map of the heterogeneous hSHBG glycosylation, which is a prerequisite for investigating the functional aspects of the hSHBG glycans.     

Skew‐normal/independent linear mixed models for censored responses with applications to HIV viral loads

Often in biomedical studies, the routine use of linear mixed‐effects models (based on Gaussian assumptions) can be questionable when the longitudinal responses are skewed in nature. Skew‐normal/elliptical models are widely used in those situations. Often, those skewed responses might also be subjected to some upper and lower quantification limits (QLs; viz., longitudinal viral‐load measures in HIV studies), beyond which they are not measurable. In this paper, we develop a Bayesian analysis of censored linear mixed models replacing the Gaussian assumptions with skew‐normal/independent (SNI) distributions. The SNI is an attractive class of asymmetric heavy‐tailed distributions that includes the skew‐normal, skew‐t, skew‐slash, and skew‐contaminated normal distributions as special cases. The proposed model provides flexibility in capturing the effects of skewness and heavy tail for responses that are either left‐ or right‐censored. For our analysis, we adopt a Bayesian framework and develop a Markov chain Monte Carlo algorithm to carry out the posterior analyses. The marginal likelihood is tractable, and utilized to compute not only some Bayesian model selection measures but also case‐deletion influence diagnostics based on the Kullback–Leibler divergence. The newly developed procedures are illustrated with a simulation study as well as an HIV case study involving analysis of longitudinal viral loads.     

Clade support measures and their adequacy

In addition to hypothesis optimality, the evaluation of clade (group, edge, split, node) support is an important aspect of phylogenetic analysis. Here we clarify the logical relationship between support and optimality and formulate adequacy conditions for support measures. Support, S, and optimality, O, are both empirical knowledge claims about the strength of hypotheses, h₁, h₂, …h_n, in relation to evidence, e, given background knowledge, b. Whereas optimality refers to the absolute strength of hypotheses, support refers to the relative strength of hypotheses. Consequently, support and optimality are logically related such that they vary in direct proportion to each other, S(h | e,b) ∝ O(h | e,b). Furthermore, in order for a support measure to be objective it must quantify support as a function of explanatory power. For example, Goodman–Bremer support and ratio of explanatory power (REP) support satisfy the adequacy requirement S(h | e,b) ∝ O(h | e,b) and calculate support as a function of explanatory power. As such, these are adequate measures of objective support. The equivalent measures for statistical optimality criteria are the likelihood ratio (or log‐likelihood difference) and likelihood difference support measures for maximum likelihood and the posterior probability ratio and posterior probability difference support measures for Bayesian inference. These statistical support measures satisfy the adequacy requirement S(h | e,b) ∝ O(h | e,b) and to that extent are internally consistent; however, they do not quantify support as a function of explanatory power and therefore are not measures of objective support. Neither the relative fit difference (RFD; relative GB support) nor any of the parsimony (bootstrap and jackknife character resampling) or statistical [bootstrap character resampling, Markov chain Monte Carlo (MCMC) clade frequencies] support measures based on clade frequencies satisfy the adequacy condition S(h | e,b) ∝ O(h | e,b) or calculate support as a function of explanatory power. As such, they are not adequate support measures. © The Willi Hennig Society 2008.     

Low‐Cost N,N′‐Bicarbazole‐Based Dopant‐Free Hole‐Transporting Materials for Large‐Area Perovskite Solar Cells

Despite the recent unprecedented development of efficient dopant‐free hole transporting materials (HTMs) for high‐performance perovskite solar cells (PSCs) on small‐area devices (≤0.1 cm²), low‐cost dopant‐free HTMs for large‐area PSCs (≥1 cm²) with high power conversion efficiencies (PCEs) have rarely been reported. Herein, two novel HTMs, 3,3′,6,6′ (or 2,2′,7,7′)‐tetrakis(N,N′‐di‐p‐methoxyphenylamine)‐N,N′‐bicarbazole (3,6 BCz‐OMeTAD or 2,7 BCz‐OMeTAD), are synthesized via an extremely simple route from very cheap raw materials. Owing to their excellent film‐forming abilities and matching energy levels, 3,6 BCz‐OMeTAD and 2,7 BCz‐OMeTAD can be successfully employed as a perfect ultrathin (≈30 nm) hole transporting layer in large‐area PSCs up to 1 cm². The 3,6 BCz‐OMeTAD and 2,7 BCz‐OMeTAD based large‐area PSCs show highest PCEs up to 17.0% and 17.6%, respectively. More importantly, high performance large‐area PSCs based on 2,7 BCz‐OMeTAD retain 90% of the initial efficiency after 2000 h storage in an ambient environment without encapsulation.