Statistics to prioritize rare variants in family-based sequencing studies with disease subtypes.

Family-based sequencing studies are increasingly used to find rare genetic variants of high risk for disease traits with familial clustering. In some studies, families with multiple disease subtypes are collected and the exomes of affected relatives are sequenced for shared rare variants (RVs). Since different families can harbor different causal variants and each family harbors many RVs, tests to detect causal variants can have low power in this study design. Our goal is rather to prioritize shared variants for further investigation by, for example, pathway analyses or functional studies. The transmission-disequilibrium test prioritizes variants based on departures from Mendelian transmission in parent-child trios. Extending this idea to families, we propose methods to prioritize RVs shared in affected relatives with two disease subtypes, with one subtype more heritable than the other. Global approaches condition on a variant being observed in the study and assume a known probability of carrying a causal variant. In contrast, local approaches condition on a variant being observed in specific families to eliminate the carrier probability. Our simulation results indicate that global approaches are robust to misspecification of the carrier probability and prioritize more effectively than local approaches even when the carrier probability is misspecified.

Potential energy surfaces and dynamic properties via ab initio composite and density functional approaches.

Vibrational spectroscopy enables critical insight into the structural and dynamic properties of molecules. Presently, the majority of theoretical approaches to spectroscopy employ wavefunction-based ab initio or density functional methods that rely on the harmonic approximation. This approximation breaks down for large molecules with strongly anharmonic bonds or for molecules with large internuclear separations. An alternative to these methods involves generating molecular anharmonic potential energy surfaces (potentials) and using them to extrapolate the vibrational frequencies. This study examines the efficacy of density functional theory (DFT) and the correlation consistent Composite Approach (ccCA) in generating anharmonic frequencies from potentials of small main group molecules. Vibrational self-consistent field Theory (VSCF) and post-VSCF methods were used to calculate the fundamental frequencies of these molecules from their potentials. Functional choice, basis set selection, and mode-coupling are also examined as factors in influencing accuracy. The absolute deviations for the calculated frequencies using potentials at the ccCA level of theory were lower than the potentials at the DFT level. With DFT resulting in bending modes that are better described than those of ccCA, a multilevel DFT:ccCA approach where DFT potentials are used for single vibrational mode potentials and ccCA is used for vibrational mode-mode couplings can be utilized for larger polyatomic systems. The frequencies obtained with this multilevel approach using VCIPSI-PT2 were closer to experimental frequencies than the scaled harmonic frequencies, indicating the success of utilizing post-VSCF methods to generate more accurate representations of computed infrared spectra.

Interaction of Per- and Polyfluoroalkyl Substances with Estrogen Receptors in Rainbow Trout (Oncorhynchus mykiss): An In Silico Investigation.

Fresh water sources, including lakes, such as the Great Lakes, are some of the most important ecosystems in the world. Despite the importance of these lakes, there is increasing concern about the presence of per- and polyfluoroalkyl substances (PFAS)─among the most prevalent contaminants of our time─due to the ability of PFAS to bioaccumulate and persist in the environment, as well as to its linkages to detrimental human and animal health effects. In this study, PFAS exposure on rainbow trout (Oncorhynchus mykiss) is examined at the molecular level, focusing on the impact of PFAS binding on the alpha (α) and beta (ß) estrogen receptors (ERs) using molecular dynamics simulations, binding free energy calculations, and structural analysis. ERs are involved in fundamental physiological processes, including reproductive system development, muscle regeneration, and immunity. This study shows that PFAS binds to both the estrogen α and estrogen ß receptors, albeit via different binding modes, due to a modification of an amino acid in the binding site as a result of a reorientation of residues in the binding pocket. As ER overactivation can occur through environmental toxins and pollutants, this study provides insights into the influence of different types of PFAS on protein function.

Atypical speech production of multisyllabic words and phrases by children with developmental dyslexia.

The prevalent "core phonological deficit" model of dyslexia proposes that the reading and spelling difficulties characterizing affected children stem from prior developmental difficulties in processing speech sound structure, for example, perceiving and identifying syllable stress patterns, syllables, rhymes and phonemes. Yet spoken word production appears normal. This suggests an unexpected disconnect between speech input and speech output processes. Here we investigated the output side of this disconnect from a speech rhythm perspective by measuring the speech amplitude envelope (AE) of multisyllabic spoken phrases. The speech AE contains crucial information regarding stress patterns, speech rate, tonal contrasts and intonational information. We created a novel computerized speech copying task in which participants copied aloud familiar spoken targets like "Aladdin." Seventy-five children with and without dyslexia were tested, some of whom were also receiving an oral intervention designed to enhance multi-syllabic processing. Similarity of the child's productions to the target AE was computed using correlation and mutual information metrics. Similarity of pitch contour, another acoustic cue to speech rhythm, was used for control analyses. Children with dyslexia were significantly worse at producing the multi-syllabic targets as indexed by both similarity metrics for computing the AE. However, children with dyslexia were not different from control children in producing pitch contours. Accordingly, the spoken production of multisyllabic phrases by children with dyslexia is atypical regarding the AE. Children with dyslexia may not appear to listeners to exhibit speech production difficulties because their pitch contours are intact. RESEARCH HIGHLIGHTS: Speech production of syllable stress patterns is atypical in children with dyslexia. Children with dyslexia are significantly worse at producing the amplitude envelope of multi-syllabic targets compared to both age-matched and reading-level-matched control children. No group differences were found for pitch contour production between children with dyslexia and age-matched control children. It may be difficult to detect speech output problems in dyslexia as pitch contours are relatively accurate.

Unraveling the electronic structure of LuH, LuN, and LuNH: building blocks of new materials.

Advances in superconductor technology have been pursued for decades, moving towards room temperature models, such as a postulated nitrogen-doped lutetium hydride network. While experimental observations have been contradictory, insight into the building blocks of potential new superconductor materials can be gained theoretically, unravelling the fascinating electronic structure of these compounds at a molecular level. Here, the fundamental building blocks of lutetium materials (LuH, LuN, and LuNH) have been examined. The structures, spectroscopic constants for the ground and excited states, and the potential energy curves have been obtained for these species using complete active self-consistent field (CASSCF) and multireference configuration interaction with Davidson's correction (MRCI+Q) methods. For LuNH, the energetic properties of its isomers are determined. The bond dissociation energies of the three building blocks are calculated with the state-of-the-art f-block ab initio correlation consistent composite approach (f-ccCA) and the high accuracy extrapolated ab initio thermochemistry (HEAT) scheme. As well, an analysis of different formation pathways of LuNH has been provided.

Correlation consistent basis sets designed for density functional theory: Third-row atoms (Ga-Br).

The correlation consistent basis sets (cc-pVnZ with n = D, T, Q, 5) for the Ga-Br elements have been redesigned, tuning the sets for use for density functional approximations. Steps to redesign these basis sets for an improved correlation energy recovery and efficiency include truncation of higher angular momentum functions, recontraction of basis set coefficients, and reoptimization of basis set exponents. These redesigned basis sets are compared with conventional cc-pVnZ basis sets and other basis sets, which are, in principle, designed to achieve systematic improvement with respect to increasing basis set size. The convergence of atomic energies, bond lengths, bond dissociation energies, and enthalpies of formation to the Kohn-Sham limit is improved relative to other basis sets where convergence to the Kohn-Sham limit is typically not observed.

Cavity-modified molecular dipole switching dynamics.

Polaritonic states, which are formed by resonances between a molecular excitation and the photonic mode of a cavity, have a number of useful properties that offer new routes to control molecular photochemistry using electric fields. To provide a theoretical description of how polaritonic states affect the real-time electron dynamics in molecules, a new method is described where the effects of strong light-molecule coupling are implemented using real-time electronic structure theory. The coupling between the molecular electronic states and the cavity is described by the Pauli-Fierz Hamiltonian, and transitions between polaritonic states are induced via an external time-dependent electric field using time-dependent configuration interaction (TDCI) theory, producing quantum electrodynamics TDCI (QED-TDCI). This method is used to study laser-induced ultrafast charge transfer and dipole-switching dynamics of the LiCN molecule inside a cavity. The increase in cavity coupling strength is found to have a significant impact on the energies and transition dipole moments of the molecule-cavity system. The convergence of the polaritonic state energies as a function of the number of included electronic and photonic basis states is discussed.

APOE ε4 carrier status modifies plasma p-tau181 concentrations in cognitively healthy super-seniors.

INTRODUCTION: This study investigates the effect of apolipoprotein E (APOE) genotype on neurology plasma biomarkers in cognitively healthy Super-Seniors. METHODS: Three hundred seventy plasma specimens from Super-Senior participants ≥ 85 years old, who have never been diagnosed with dementia, cancer, diabetes, cardiovascular, or major pulmonary disease, were analyzed on the Quanterix Simoa HD-X analyzer using commercial Neurology 4-plex E and phosphorylated tau (p-tau)181 assays. RESULTS: Eighty (22%) participants were APOE ε4 carriers and 290 (73%) were non-carriers. No significant differences were found between APOE ε4 carriers and non-carriers regarding age, sex, or Mini-Mental State Examination scores. In APOE ε4 carriers, plasma amyloid beta 42/40 was lower and p-tau181 and glial fibrillary acidic protein were higher compared to non-APOE ε4 carriers. After adjusting for demographic variables, p-tau181 was the only biomarker to remain significantly associated with APOE ε4 carrier status. DISCUSSION: APOE ε4 genotype modifies plasma p-tau181 concentration in seniors resilient to age-related clinical disease, suggesting that some Super-Seniors may have Alzheimer's disease pathology without progressing to cognitive decline. HIGHLIGHTS: Healthy seniors enable identification of associations that may be masked by disease. Plasma phosphorylated tau (p-tau)181 concentrations associate with apolipoprotein E (APOE) ε4 carriership in healthy seniors. APOE should be accounted for when interpreting p-tau181, regardless of disease.

Solid-Phase-Supported Chemoenzymatic Synthesis and Analysis of Chondroitin Sulfate Proteoglycan Glycopeptides.

Proteoglycans (PGs), consisting of glycosaminoglycans (GAGs) linked with the core protein through a tetrasaccharide linkage region, play roles in many important biological events. The chemical synthesis of PG glycopeptides is extremely challenging. In this work, the enzymes required for synthesis of chondroitin sulfate (CS) PG (CSPG) have been expressed and the suitable sequence of enzymatic reactions has been established. To expedite CSPG synthesis, the peptide acceptor was immobilized on solid phase and the glycan units were directly installed enzymatically onto the peptide. Subsequent enzymatic chain elongation and sulfation led to the successful synthesis of CSPG glycopeptides. The CS dodecasaccharide glycopeptide was the longest homogeneous CS glycopeptide synthesized to date. The enzymatic synthesis was much more efficient than the chemical synthesis of the corresponding CS glycopeptides, which could reduce the total number of synthetic steps by 80 %. The structures of the CS glycopeptides were confirmed by mass spectrometry analysis and NMR studies. In addition, the interactions between the CS glycopeptides and cathepsin G were studied. The sulfation of glycan chain was found to be important for binding with cathepsin G. This efficient chemoenzymatic strategy opens new avenues to investigate the structures and functions of PGs.

A genome-wide gene-environment interaction study of breast cancer risk for women of European ancestry.

BACKGROUND: Genome-wide studies of gene-environment interactions (G×E) may identify variants associated with disease risk in conjunction with lifestyle/environmental exposures. We conducted a genome-wide G×E analysis of ~ 7.6 million common variants and seven lifestyle/environmental risk factors for breast cancer risk overall and for estrogen receptor positive (ER +) breast cancer. METHODS: Analyses were conducted using 72,285 breast cancer cases and 80,354 controls of European ancestry from the Breast Cancer Association Consortium. Gene-environment interactions were evaluated using standard unconditional logistic regression models and likelihood ratio tests for breast cancer risk overall and for ER + breast cancer. Bayesian False Discovery Probability was employed to assess the noteworthiness of each SNP-risk factor pairs. RESULTS: Assuming a 1 × 10-5 prior probability of a true association for each SNP-risk factor pairs and a Bayesian False Discovery Probability < 15%, we identified two independent SNP-risk factor pairs: rs80018847(9p13)-LINGO2 and adult height in association with overall breast cancer risk (ORint = 0.94, 95% CI 0.92-0.96), and rs4770552(13q12)-SPATA13 and age at menarche for ER + breast cancer risk (ORint = 0.91, 95% CI 0.88-0.94). CONCLUSIONS: Overall, the contribution of G×E interactions to the heritability of breast cancer is very small. At the population level, multiplicative G×E interactions do not make an important contribution to risk prediction in breast cancer.

CCNE1 and survival of patients with tubo-ovarian high-grade serous carcinoma: An Ovarian Tumor Tissue Analysis consortium study.

BACKGROUND: Cyclin E1 (CCNE1) is a potential predictive marker and therapeutic target in tubo-ovarian high-grade serous carcinoma (HGSC). Smaller studies have revealed unfavorable associations for CCNE1 amplification and CCNE1 overexpression with survival, but to date no large-scale, histotype-specific validation has been performed. The hypothesis was that high-level amplification of CCNE1 and CCNE1 overexpression, as well as a combination of the two, are linked to shorter overall survival in HGSC. METHODS: Within the Ovarian Tumor Tissue Analysis consortium, amplification status and protein level in 3029 HGSC cases and mRNA expression in 2419 samples were investigated. RESULTS: High-level amplification (>8 copies by chromogenic in situ hybridization) was found in 8.6% of HGSC and overexpression (>60% with at least 5% demonstrating strong intensity by immunohistochemistry) was found in 22.4%. CCNE1 high-level amplification and overexpression both were linked to shorter overall survival in multivariate survival analysis adjusted for age and stage, with hazard stratification by study (hazard ratio [HR], 1.26; 95% CI, 1.08-1.47, p = .034, and HR, 1.18; 95% CI, 1.05-1.32, p = .015, respectively). This was also true for cases with combined high-level amplification/overexpression (HR, 1.26; 95% CI, 1.09-1.47, p = .033). CCNE1 mRNA expression was not associated with overall survival (HR, 1.00 per 1-SD increase; 95% CI, 0.94-1.06; p = .58). CCNE1 high-level amplification is mutually exclusive with the presence of germline BRCA1/2 pathogenic variants and shows an inverse association to RB1 loss. CONCLUSION: This study provides large-scale validation that CCNE1 high-level amplification is associated with shorter survival, supporting its utility as a prognostic biomarker in HGSC.

Increased FOXJ1 protein expression is associated with improved overall survival in high-grade serous ovarian carcinoma: an Ovarian Tumor Tissue Analysis Consortium Study.

BACKGROUND: Recently, we showed a >60% difference in 5-year survival for patients with tubo-ovarian high-grade serous carcinoma (HGSC) when stratified by a 101-gene mRNA expression prognostic signature. Given the varied patient outcomes, this study aimed to translate prognostic mRNA markers into protein expression assays by immunohistochemistry and validate their survival association in HGSC. METHODS: Two prognostic genes, FOXJ1 and GMNN, were selected based on high-quality antibodies, correlation with protein expression and variation in immunohistochemical scores in a preliminary cohort (n = 134 and n = 80, respectively). Six thousand four hundred and thirty-four (FOXJ1) and 5470 (GMNN) formalin-fixed, paraffin-embedded ovarian neoplasms (4634 and 4185 HGSC, respectively) represented on tissue microarrays from the Ovarian Tumor Tissue Analysis consortium underwent immunohistochemical staining and scoring, then univariate and multivariate survival analysis. RESULTS: Consistent with mRNA, FOXJ1 protein expression exhibited a linear, increasing association with improved overall survival in HGSC patients. Women with >50% expression had the most favourable outcomes (HR = 0.78, 95% CI 0.67-0.91, p < 0.0001). GMNN protein expression was not significantly associated with overall HSGC patient survival. However, HGSCs with >35% GMNN expression showed a trend for better outcomes, though this was not significant. CONCLUSION: We provide foundational evidence for the prognostic value of FOXJ1 in HGSC, validating the prior mRNA-based prognostic association by immunohistochemistry.

Thermochemistry of per- and polyfluoroalkyl substances.

The determination of gas phase thermochemical properties of per- and polyfluoroalkyl substances (PFAS) is central to understanding the long-range transport behavior of PFAS in the atmosphere. Prior gas-phase studies have reported the properties of perfluorinated sulfonic acid (PFOS) and perfluorinated octanoic acid (PFOA). Here, this study reports the gas phase enthalpies of formation of short- and long-chain PFAS and their precursor molecules determined using density functional theory (DFT) and ab initio approaches. Two density functionals, two ab initio methods and an empirical method were used to compute enthalpies of formation with the total atomization approach and an isogyric reaction. The performance of the computational methods employed in this work were validated against the experimental enthalpies of linear alkanoic acids and perfluoroalkanes. The gas-phase determinations will be useful for future studies of PFAS in the atmosphere, and the methodological choices will be helpful in the study of other PFAS.

Microglia specific deletion of miR-155 in Alzheimer's disease mouse models reduces amyloid-ß pathology but causes hyperexcitability and seizures.

Alzheimer's Disease (AD) is characterized by the accumulation of extracellular amyloid-ß (Aß) as well as CNS and systemic inflammation. Microglia, the myeloid cells resident in the CNS, use microRNAs to rapidly respond to inflammatory signals. MicroRNAs (miRNAs) modulate inflammatory responses in microglia, and miRNA profiles are altered in Alzheimer's disease (AD) patients. Expression of the pro-inflammatory miRNA, miR-155, is increased in the AD brain. However, the role of miR-155 in AD pathogenesis is not well-understood. We hypothesized that miR-155 participates in AD pathophysiology by regulating microglia internalization and degradation of Aß. We used CX3CR1CreER/+ to drive-inducible, microglia-specific deletion of floxed miR-155 alleles in two AD mouse models. Microglia-specific inducible deletion of miR-155 in microglia increased anti-inflammatory gene expression while reducing insoluble Aß1-42 and plaque area. Yet, microglia-specific miR-155 deletion led to early-onset hyperexcitability, recurring spontaneous seizures, and seizure-related mortality. The mechanism behind hyperexcitability involved microglia-mediated synaptic pruning as miR-155 deletion altered microglia internalization of synaptic material. These data identify miR-155 as a novel modulator of microglia Aß internalization and synaptic pruning, influencing synaptic homeostasis in the setting of AD pathology.

Profiling the immune landscape in mucinous ovarian carcinoma.

OBJECTIVE: Mucinous ovarian carcinoma (MOC) is a rare histotype of ovarian cancer, with low response rates to standard chemotherapy, and very poor survival for patients diagnosed at advanced stage. There is a limited understanding of the MOC immune landscape, and consequently whether immune checkpoint inhibitors could be considered for a subset of patients. METHODS: We performed multicolor immunohistochemistry (IHC) and immunofluorescence (IF) on tissue microarrays in a cohort of 126 MOC patients. Cell densities were calculated in the epithelial and stromal components for tumor-associated macrophages (CD68+/PD-L1+, CD68+/PD-L1-), T cells (CD3+/CD8-, CD3+/CD8+), putative T-regulatory cells (Tregs, FOXP3+), B cells (CD20+/CD79A+), plasma cells (CD20-/CD79a+), and PD-L1+ and PD-1+ cells, and compared these values with clinical factors. Univariate and multivariable Cox Proportional Hazards assessed overall survival. Unsupervised k-means clustering identified patient subsets with common patterns of immune cell infiltration. RESULTS: Mean densities of PD1+ cells, PD-L1- macrophages, CD4+ and CD8+ T cells, and FOXP3+ Tregs were higher in the stroma compared to the epithelium. Tumors from advanced (Stage III/IV) MOC had greater epithelial infiltration of PD-L1- macrophages, and fewer PD-L1+ macrophages compared with Stage I/II cancers (p = 0.004 and p = 0.014 respectively). Patients with high epithelial density of FOXP3+ cells, CD8+/FOXP3+ cells, or PD-L1- macrophages, had poorer survival, and high epithelial CD79a + plasma cells conferred better survival, all upon univariate analysis only. Clustering showed that most MOC (86%) had an immune depleted (cold) phenotype, with only a small proportion (11/76,14%) considered immune inflamed (hot) based on T cell and PD-L1 infiltrates. CONCLUSION: In summary, MOCs are mostly immunogenically 'cold', suggesting they may have limited response to current immunotherapies.

Binding of Per- and Polyfluoroalkyl Substances (PFAS) to the PPARγ/RXRα-DNA Complex.

Nuclear receptors are the fundamental building blocks of gene expression regulation and the focus of many drug targets. While binding to DNA, nuclear receptors act as transcription factors, governing a multitude of functions in the human body. Peroxisome proliferator-activator receptor γ (PPARγ) and the retinoid X receptor α (RXRα) form heterodimers with unique properties and have a primordial role in insulin sensitization. This PPARγ/RXRα heterodimer has been shown to be impacted by per- and polyfluoroalkyl substances (PFAS) and linked to a variety of significant health conditions in humans. Herein, a selection of the most common PFAS (legacy and emerging) was studied utilizing molecular dynamics simulations for PPARγ/RXRα. The local and global structural effects of PFAS binding on the known ligand binding pockets of PPARγ and RXRα as well as the DNA binding domain (DBD) of RXRα were inspected. The binding free energies were predicted computationally and were compared between the different binding pockets. In addition, two electronic structure approaches were utilized to model the interaction of PFAS within the DNA binding domain, density functional theory (DFT) and domain-based pair natural orbital coupled cluster with perturbative triples (DLPNO-CCSD(T)) approaches, with implicit solvation. Residue decomposition and hydrogen-bonding analysis were also performed, detailing the role of prominent residues in molecular recognition. The role of l-carnitine is explored as a potential in vivo remediation strategy for PFAS interaction with the PPARγ/RXRα heterodimer. In this work, it was found that PFAS can bind and act as agonists for all of the investigated pockets. For the first time in the literature, PFAS are postulated to bind to the DNA binding domain in a nonspecific manner. In addition, for the PPARγ ligand binding domain, l-carnitine shows promise in replacing smaller PFAS from the pocket.

Multireference Wavefunction-Based Investigation of the Ground and Excited States of LrF and LrO.

Complete active space self-consistent field (CASSCF) and multireference configuration interaction with Davidson correction (MRCI+Q) calculations have been carried out for lawrencium fluoride (LrF) and lawrencium oxide (LrO) molecules, detailing 19 and 20 electronic states for LrF and LrO, respectively. For LrF, two dissociation channels were considered, Lr(2P)+F(2P) and Lr(2D)+F(2P). However, due to the more complex electronic manifold of LrO, three dissociation channels were computed: Lr(2P)+O(3P), Lr(2D)+O(3P), and Lr(2P)+O(1D). In addition, equilibrium bond lengths, harmonic vibrational frequencies ωe, anharmonicity constants ωeχe, ΔG1/2 values, and excitation energies Te for the ground and several excited electronic states were calculated for both molecules, for the first time. Bond dissociation energies (BDEs) were calculated for LrF and LrO using several different levels of theory: unrestricted coupled-cluster with single, double, and perturbative triple excitations (UCCSD(T)), density functional theory (B3LYP, TPSS, M06-L, and PBE), and the correlation-consistent composite approach developed for f-elements (f-ccCA).

Tailoring light-induced charge transfer and intersystem crossing in FeCO using time-dependent spin-orbit configuration interaction.

Real-time (RT) electronic structure methods provide a natural framework for describing light-matter interactions in arbitrary time-dependent electromagnetic fields (EMF). Optically induced excited state transitions are of particular interest, which require tuned EMF to drive population transfer to and from the specific state(s) of interest. Intersystem crossing, or spin-flip, may be driven through shaped EMF or laser pulses. These transitions can result in long-lived "spin-trapped" excited states, which are especially useful for materials requiring charge separation or protracted excited state lifetimes. Time-dependent configuration interaction (TDCI) is unique among RT methods in that it may be implemented in a basis of eigenstates, allowing for rapid propagation of the time-dependent Schrödinger equation. The recent spin-orbit TDCI (TD-SOCI) enables a real-time description of spin-flip dynamics in an arbitrary EMF and, therefore, provides an ideal framework for rational pulse design. The present study explores the mechanism of multiple spin-flip pathways for a model transition metal complex, FeCO, using shaped pulses designed to drive controlled intersystem crossing and charge transfer. These results show that extremely tunable excited state dynamics can be achieved by considering the dipole transition matrix elements between the states of interest.

Penalized Logistic Regression Analysis for Genetic Association Studies of Binary Phenotypes.

INTRODUCTION: Increasingly, logistic regression methods for genetic association studies of binary phenotypes must be able to accommodate data sparsity, which arises from unbalanced case-control ratios and/or rare genetic variants. Sparseness leads to maximum likelihood estimators (MLEs) of log-OR parameters that are biased away from their null value of zero and tests with inflated type 1 errors. Different penalized-likelihood methods have been developed to mitigate sparse-data bias. We study penalized logistic regression using a class of log-F priors indexed by a shrinkage parameter m to shrink the biased MLE towards zero. For a given m, log-F-penalized logistic regression may be easily implemented using data augmentation and standard software. METHOD: We propose a two-step approach to the analysis of a genetic association study: first, a set of variants that show evidence of association with the trait is used to estimate m; and second, the estimated m is used for log-F-penalized logistic regression analyses of all variants using data augmentation with standard software. Our estimate of m is the maximizer of a marginal likelihood obtained by integrating the latent log-ORs out of the joint distribution of the parameters and observed data. We consider two approximate approaches to maximizing the marginal likelihood: (i) a Monte Carlo EM algorithm (MCEM) and (ii) a Laplace approximation (LA) to each integral, followed by derivative-free optimization of the approximation. RESULTS: We evaluate the statistical properties of our proposed two-step method and compared its performance to other shrinkage methods by a simulation study. Our simulation studies suggest that the proposed log-F-penalized approach has lower bias and mean squared error than other methods considered. We also illustrate the approach on data from a study of genetic associations with "super senior" cases and middle aged controls. DISCUSSION/CONCLUSION: We have proposed a method for single rare variant analysis with binary phenotypes by logistic regression penalized by log-F priors. Our method has the advantage of being easily extended to correct for confounding due to population structure and genetic relatedness through a data augmentation approach.

Neural sampling of the speech signal at different timescales by children with dyslexia.

Phonological difficulties characterize individuals with dyslexia across languages. Currently debated is whether these difficulties arise from atypical neural sampling of (or entrainment to) auditory information in speech at slow rates (<10 Hz, related to speech rhythm), faster rates, or neither. MEG studies with adults suggest that atypical sampling in dyslexia affects faster modulations in the neurophysiological gamma band, related to phoneme-level representation. However, dyslexic adults have had years of reduced experience in converting graphemes to phonemes, which could itself cause atypical gamma-band activity. The present study was designed to identify specific linguistic timescales at which English children with dyslexia may show atypical entrainment. Adopting a developmental focus, we hypothesized that children with dyslexia would show atypical entrainment to the prosodic and syllable-level information that is exaggerated in infant-directed speech and carried primarily by amplitude modulations <10 Hz. MEG was recorded in a naturalistic story-listening paradigm. The modulation bands related to different types of linguistic information were derived directly from the speech materials, and lagged coherence at multiple temporal rates spanning 0.9-40 Hz was computed. Group differences in lagged speech-brain coherence between children with dyslexia and control children were most marked in neurophysiological bands corresponding to stress and syllable-level information (<5 Hz in our materials), and phoneme-level information (12-40 Hz). Functional connectivity analyses showed network differences between groups in both hemispheres, with dyslexic children showing significantly reduced global network efficiency. Global network efficiency correlated with dyslexic children's oral language development and with control children's reading development. These developmental data suggest that dyslexia is characterized by atypical neural sampling of auditory information at slower rates. They also throw new light on the nature of the gamma band temporal sampling differences reported in MEG dyslexia studies with adults.