Gene-based association tests in family samples using GWAS summary statistics.

Genome-wide association studies (GWAS) have led to rapid growth in detecting genetic variants associated with various phenotypes. Owing to a great number of publicly accessible GWAS summary statistics, and the difficulty in obtaining individual-level genotype data, many existing gene-based association tests have been adapted to require only GWAS summary statistics rather than individual-level data. However, these association tests are restricted to unrelated individuals and thus do not apply to family samples directly. Moreover, due to its flexibility and effectiveness, the linear mixed model has been increasingly utilized in GWAS to handle correlated data, such as family samples. However, it remains unknown how to perform gene-based association tests in family samples using the GWAS summary statistics estimated from the linear mixed model. In this study, we show that, when family size is negligible compared to the total sample size, the diagonal block structure of the kinship matrix makes it possible to approximate the correlation matrix of marginal Z scores by linkage disequilibrium matrix. Based on this result, current methods utilizing summary statistics for unrelated individuals can be directly applied to family data without any modifications. Our simulation results demonstrate that this proposed strategy controls the type 1 error rate well in various situations. Finally, we exemplify the usefulness of the proposed approach with a dental caries GWAS data set.

The lipid droplet assembly complex consists of seipin and four accessory factors in budding yeast.

Seipin, a crucial protein for cellular lipid droplet (LD) assembly, oligomerizes at the interface between the endoplasmic reticulum and LDs to facilitate neutral lipid packaging. Using proximity labeling, we identified four proteins-Ldo45, Ldo16, Tgl4, and Pln1-that are recruited to the vicinity of yeast seipin, the Sei1-Ldb16 complex, exclusively when seipin function is intact, hence termed seipin accessory factors. Localization studies identified Tgl4 at the endoplasmic reticulum-LD contact site, in contrast to Ldo45, Ldo16, and Pln1 at the LD surface. Cells with compromised seipin function resulted in uneven distribution of these proteins with aberrant LDs, supporting a central role of seipin in orchestrating their association with the LD. Overexpression of any seipin accessory factor causes LD aggregation and affects a subset of LD protein distribution, highlighting the importance of their stoichiometry. Although single factor mutations show minor LD morphology changes, the combined mutations have additive effects. Lastly, we present evidence that seipin accessory factors assemble and interact with seipin in the absence of neutral lipids and undergo dynamical rearrangements during LD formation induction, with Ldo45 acting as a central hub recruiting other factors to interact with the seipin complex.

Maize lipid droplet-associated protein 2 is recruited by a virus to enhance viral multiplication and infection through regulating cellular fatty acid metabolism.

Pathogen infection induces massive reprogramming of host primary metabolism. Lipid and fatty acid (FA) metabolism is generally disrupted by pathogens and co-opted for their proliferation. Lipid droplets (LDs) that play important roles in regulating cellular lipid metabolism are utilized by a variety of pathogens in mammalian cells. However, the function of LDs during pathogenic infection in plants remains unknown. We show here that infection by rice black streaked dwarf virus (RBSDV) affects the lipid metabolism of maize, which causes elevated accumulation of C18 polyunsaturated fatty acids (PUFAs) leading to viral proliferation and symptom development. The overexpression of one of the two novel LD-associated proteins (LDAPs) of maize (ZmLDAP1 and ZmLDAP2) induces LD clustering. The core capsid protein P8 of RBSDV interacts with ZmLDAP2 and prevents its degradation through the ubiquitin-proteasome system mediated by a UBX domain-containing protein, PUX10. In addition, silencing of ZmLDAP2 downregulates the expression of FA desaturase genes in maize, leading to a decrease in C18 PUFAs levels and suppression of RBSDV accumulation. Our findings reveal that plant virus may recruit LDAP to regulate cellular FA metabolism to promote viral multiplication and infection. These results expand the knowledge of LD functions and viral infection mechanisms in plants.

Leveraging LD eigenvalue regression to improve the estimation of SNP heritability and confounding inflation.

Heritability is a fundamental concept in genetic studies, measuring the genetic contribution to complex traits and bringing insights about disease mechanisms. The advance of high-throughput technologies has provided many resources for heritability estimation. Linkage disequilibrium (LD) score regression (LDSC) estimates both heritability and confounding biases, such as cryptic relatedness and population stratification, among single-nucleotide polymorphisms (SNPs) by using only summary statistics released from genome-wide association studies. However, only partial information in the LD matrix is utilized in LDSC, leading to loss in precision. In this study, we propose LD eigenvalue regression (LDER), an extension of LDSC, by making full use of the LD information. Compared to state-of-the-art heritability estimating methods, LDER provides more accurate estimates of SNP heritability and better distinguishes the inflation caused by polygenicity and confounding effects. We demonstrate the advantages of LDER both theoretically and with extensive simulations. We applied LDER to 814 complex traits from UK Biobank, and LDER identified 363 significantly heritable phenotypes, among which 97 were not identified by LDSC.

The ER-Localized Transmembrane Protein EPG-3/VMP1 Regulates SERCA Activity to Control ER-Isolation Membrane Contacts for Autophagosome Formation.

During autophagosome formation in mammalian cells, isolation membranes (IMs; autophagosome precursors) dynamically contact the ER. Here, we demonstrated that the ER-localized metazoan-specific autophagy protein EPG-3/VMP1 controls ER-IM contacts. Loss of VMP1 causes stable association of IMs with the ER, thus blocking autophagosome formation. Interaction of WIPI2 with the ULK1/FIP200 complex and PI(3)P contributes to the formation of ER-IM contacts, and these interactions are enhanced by VMP1 depletion. VMP1 controls contact formation by promoting SERCA (sarco[endo]plasmic reticulum calcium ATPase) activity. VMP1 interacts with SERCA and prevents formation of the SERCA/PLN/SLN inhibitory complex. VMP1 also modulates ER contacts with lipid droplets, mitochondria, and endosomes. These ER contacts are greatly elevated by the SERCA inhibitor thapsigargin. Calmodulin acts as a sensor/effector to modulate the ER contacts mediated by VMP1/SERCA. Our study provides mechanistic insights into the establishment and disassociation of ER-IM contacts and reveals that VMP1 modulates SERCA activity to control ER contacts.

The E-Syt3 cleavage and traffic uncovers the primordial cisterna, a new organelle that mothers the lipid droplets in the adipocyte.

Extended synaptotagmins are endoplasmic reticulum proteins consisting of an SMP domain and multiple C2 domains that bind phospholipids and Ca2+ . E-Syts create contact junctions between the ER and plasma membrane (PM) to facilitate the exchange of glycerophospholipids between the apposed membranes. We find in the differentiating adipocyte that the E-Syt3 carboxyl domain is cleaved by a multi-step mechanism that includes removing the C2C domain. Confocal and live-cell time-lapse studies show that truncated E-Syt3ΔC2C, as well as endogenous E-Syt3 and the coat protein PLIN1, target the LDs from an annular, single giant ER cisterna. Inhibition of the proteasome blocks the proteolytic cleavage of Esyt3 and E-Syt3ΔC2C and causes the E-Syt3ΔC2C retention in the giant cisterna. The Esyt3 and PLIN1 distributions and LDs biogenesis show that the primordial cisterna, as we call it, is the birth and nurturing site of LDs in the adipocyte. Isoproterenol-induced lipolysis results in loss of cytoplasmic LDs and reappearance of the primordial cisterna. Electron microscopy and 3D-electron tomography studies show that the primordial cisterna consists of a tightly packed network of varicose tubules with extensively blistered membranes. Rounds of homotypic fusions from nascent to mature LDs play a central role in LD growth. The knockdown of E-Syt3 inhibits LD biogenesis. The identification of the primordial cisterna, an organelle that substitutes the randomly scattered ER foci that mother the LDs in non-adipose cells, sets the stage for a better understanding of LD biogenesis in the adipocyte.

Identification of eQTLs using different sets of single nucleotide polymorphisms associated with carcass and body composition traits in pigs.

BACKGROUND: Mapping expression quantitative trait loci (eQTLs) in skeletal muscle tissue in pigs is crucial for understanding the relationship between genetic variation and phenotypic expression of carcass traits in meat animals. Therefore, the primary objective of this study was to evaluate the impact of different sets of single nucleotide polymorphisms (SNP), including scenarios removing SNPs pruned for linkage disequilibrium (LD) and SNPs derived from SNP chip arrays and RNA-seq data from liver, brain, and skeletal muscle tissues, on the identification of eQTLs in the Longissimus lumborum tissue, associated with carcass and body composition traits in Large White pigs. The SNPs identified from muscle mRNA were combined with SNPs identified in the brain and liver tissue transcriptomes, as well as SNPs from the GGP Porcine 50 K SNP chip array. Cis- and trans-eQTLs were identified based on the skeletal muscle gene expression level, followed by functional genomic analyses and statistical associations with carcass and body composition traits in Large White pigs. RESULTS: The number of cis- and trans-eQTLs identified across different sets of SNPs (scenarios) ranged from 261 to 2,539 and from 29 to 13,721, respectively. Furthermore, 6,180 genes were modulated by eQTLs in at least one of the scenarios evaluated. The eQTLs identified were not significantly associated with carcass and body composition traits but were significantly enriched for many traits in the "Meat and Carcass" type QTL. The scenarios with the highest number of cis- (n = 304) and trans- (n = 5,993) modulated genes were the unpruned and LD-pruned SNP set scenarios identified from the muscle transcriptome. These genes include 84 transcription factor coding genes. CONCLUSIONS: After LD pruning, the set of SNPs identified based on the transcriptome of the skeletal muscle tissue of pigs resulted in the highest number of genes modulated by eQTLs. Most eQTLs are of the trans type and are associated with genes influencing complex traits in pigs, such as transcription factors and enhancers. Furthermore, the incorporation of SNPs from other genomic regions to the set of SNPs identified in the porcine skeletal muscle transcriptome contributed to the identification of eQTLs that had not been identified based on the porcine skeletal muscle transcriptome alone.

Molecular genetic diversity and linkage disequilibrium structure of the Egyptian faba bean using Single Primer Enrichment Technology (SPET).

Faba bean is an important legume crop. The genetic diversity among faba bean genotypes is very important for the genetic improvement of target traits. A set of 128 fab bean genotypes that are originally from Egypt were used in this study to investigate the genetic diversity and population structure. The 128 genotypes were genotyped using the Single Primer Enrichment Technology (SPET) by which a set of 6759 SNP markers were generated after filtration. The SNP markers were distributed on all chromosomes with a range extending from 822 (Chr. 6) to 1872 (Chr.1). The SNP markers had wide ranges of polymorphic information content (PIC), gene diversity (GD), and minor allele frequency. The analysis of population structure divided the Egyptian faba bean population into five subpopulations. Considerable genetic distance was found among all genotypes, ranging from 0.1 to 0.4. The highly divergent genotype was highlighted in this study and the genetic distance among genotypes ranged from 0.1 and 0.6. Moreover, the structure of linkage disequilibrium was studied, and the analysis revealed a low level of LD in the Egyptian faba bean population. A slow LD decay at the genomic and chromosomal levels was observed. Interestingly, the distribution of haplotype blocks was presented in each chromosome and the number of haplotype block ranged from 65 (Chr. 4) to 156 (Chr. 1). Migration and genetic drift are the main reasons for the low LD in the Egyptian faba bean population. The results of this study shed light on the possibility of the genetic improvement of faba bean crop in Egypt and conducting genetic association analyses to identify candidate genes associated with target traits (e.g. protein content, grain yield, etc.) in this panel.

Effective population size in field pea.

BACKGROUND: Effective population size (Ne) is a pivotal parameter in population genetics as it can provide information on the rate of inbreeding and the contemporary status of genetic diversity in breeding populations. The population with smaller Ne can lead to faster inbreeding, with little potential for genetic gain making selections ineffective. The importance of Ne has become increasingly recognized in plant breeding, which can help breeders monitor and enhance the genetic variability or redesign their selection protocols. Here, we present the first Ne estimates based on linkage disequilibrium (LD) in the pea genome. RESULTS: We calculated and compared Ne using SNP markers from North Dakota State University (NDSU) modern breeding lines and United States Department of Agriculture (USDA) diversity panel. The extent of LD was highly variable not only between populations but also among different regions and chromosomes of the genome. Overall, NDSU had a higher and longer-range LD than the USDA that could extend up to 500 Kb, with a genome-wide average r2 of 0.57 (vs 0.34), likely due to its lower recombination rates and the selection background. The estimated Ne for the USDA was nearly three-fold higher (Ne = 174) than NDSU (Ne = 64), which can be confounded by a high degree of population structure due to the selfing nature of pea. CONCLUSIONS: Our results provided insights into the genetic diversity of the germplasm studied, which can guide plant breeders to actively monitor Ne in successive cycles of breeding to sustain viability of the breeding efforts in the long term.

Comparability of comparative toxicity: insect sensitivity to imidacloprid reveals huge variations across species but also within species.

Pesticides have been identified as major drivers of insect biodiversity loss. Thus, the study of their effects on non-pest insect species has attracted a lot of attention in recent decades. In general toxicology, the 'gold standard' to assess the toxicity of a substance is to measure mass-specific LD50 (i.e. median lethal dose per unit body mass). In entomology, reviews attempting to compare these data across all available studies are lacking. To fill this gap in knowledge, we performed a systematic review of the lethality of imidacloprid for adult insects. Imidacloprid is possibly the most extensively studied insecticide in recent times, yet we found that little is comparable across studies, owing to both methodological divergence and missing estimates of body mass. By accounting for body mass whenever possible, we show how imidacloprid sensitivity spans across an apparent range of approximately six orders of magnitude across insect species. Very high variability within species can also be observed owing to differences in exposure methods and observation time. We suggest that a more comparable and comprehensive approach has both biological and economic relevance. Ultimately, this would help to identify differences that could direct research towards preventing non-target species from being negatively affected.

Bora-Diaza-Indacene Based Fluorescent Probes for Simultaneous Visualisation of Lipid Droplets and Endoplasmic Reticulum.

Organelle selective fluorescent probes, especially those capable of concurrent detection of specific organelles, are of benefit to the research community in delineating the interplay between various organelles and the impact of such interaction in maintaining cellular homeostasis and its disruption in the diseased state. Although very useful, such probes are synthetically challenging to design due to the stringent lipophilicity requirement posed by different organelles, and hence, the lack of such probes being reported so far. This work details the synthesis, photophysical properties, and cellular imaging studies of two bora-diaza-indacene based fluorescent probes that can specifically and simultaneously visualise lipid droplets and endoplasmic reticulum; two organelles suggested having close interactions and implicated in stress-induced cellular dysfunction and disease progression.

RNA-seq analysis revealed the pathogenicity of Vibrio vulnificus to American eel (Anguilla rostrata) and the strategy of host anti-V. vulnificus infection.

Vibrio vulnificus is a commonly pathogenic bacterium in cultivated eels, but its pathogenicity to American eel (Anguilla rostrata) and the molecular mechanism of host anti-V. vulnificus infection remains uncertain. In this study, American eels were infected with different dose of V. vulnificus to determine the LD50. Then, bacterial load in the liver and kidney histopathology were assessed post the LD50 of V. vulnificus infection. Additionally, gene expressions of 18 immune related genes in the liver, spleen and kidney were detected. Furthermore, transcriptome sequencing and enrichment of differentially expressed genes (DEGs) were analyzed in the eel spleens between pre-infection (Con_0), post-36 h (Vv_36), and post-60 h (Vv_60) infection. The results showed that LD50 of V. vulnificus to American eels was determined to be 5.0 × 105 cfu/g body weight, and the bacterial load peaked at 24 and 12 h post the infection (hpi) in the kidney and liver, respectively. The histopathology was highlighted by necrotic hepatocytes and splenic cells, congestion blood vessels in liver and spleen, atrophied glomeruli and vacuolization of renal tubular epithelial cells. The results of RT-PCR revealed that 18 host immune-related genes showed significantly up or downregulated expression post-infection compare to that of pre-infection. Finally, results of the RNA-seq revealed 16 DEGs play essential role to the immunosuppression in American eels, and the protein-protein interactions shed light on the widespread upregulation GEGs related to metabolism and immune response maintained the host cell homeostasis post the V. vulnificus infection, shedding new light on our understanding of the V. vulnificus pathogenesis towards understudied American eel and the host anti-V. vulnificus infection strategies in gene transcript.

Transcriptome analysis of host anti-Aeromonas hydrophila infection revealed the pathogenicity of A. hydrophila to American eels (Anguilla rostrata).

Aeromonas hydrophila is a commonly pathogenic bacterium in cultivated eels, but its pathogenicity to American eel (Anguilla rostrata) and the molecular mechanism of host anti-A. hydrophila infection remains uncertain. In this study, LD50 of A. hydrophila to American eels was determined and bacterial load in the liver and kidney of eels was assessed post 2.56 doses of LD50 of A. hydrophila infection. The results showed that the LD50 of A. hydrophila to American eels was determined to be 3.9 × 105 cfu/g body weight (7.8 × 106 cfu/fish), and the bacterial load peaked at 36 h post the infection (hpi) in the liver. Then, the histopathology was highlighted by congestion in splenic blood vessels, atrophied glomeruli, and necrotic hepatocytes. Additionally, the results of qRT-PCR revealed that 18 host immune-related genes showed significantly up or downregulated post-infection compare to that of pre-infection. Finally, results of the RNA-seq revealed 10 hub DEGs and 7 encoded proteins play essential role to the anti-A. hydrophila infection in American eels. Pathogenicity of A. hydrophila to American eels and RNA-seq of host anti-A. hydrophila infection were firstly reported in this study, shedding new light on our understanding of the A. hydrophila pathogenesis and the host immune response to the A. hydrophila infection strategies in gene transcript.

Galleria mellonella larvae as an alternative model to determine the pathogenicity of avian pathogenic Escherichia coli.

RESEARCH HIGHLIGHTS: Galleria mellonella larvae are a viable model for determining APEC pathogenicity.Larval disease score is the main variable for determining APEC pathogenicity.Response variables should be evaluated up to 24 h post-inoculation.

Interference between Escherichia coli genotypes from the E. coli peritonitis syndrome given simultaneously to productive SPF White Leghorn hens by intratracheal inoculation.

RESEARCH HIGHLIGHTS: All four or only two E. coli genotypes were found in groups of hens given mixes of four genotypes.In contrast, only one genotype was found in individual hens.E. coli genotypes interfere with each other in hens after given as a mix.Interference is likely based on a random process.Broad protection can best be assessed by challenging with single genotypes.

Measurement practices of alanine aminotransferase in children: Temporal changes and etiology for increased values.

Data on alanine aminotransferase (ALT) measurement practices and diagnoses associated with increased values are limited. We evaluated these issues by collecting ALT measurements from 1- to 16-year-old patients investigated in 1992-2018 in a tertiary center. Diagnoses were gathered in 2008-2018. Altogether 145,092 measurements from 28,118 children were taken 42% undergoing repeated testing. Testing increased from 21/1000 to 81/1000 children and the prevalence of elevated values fluctuated between 18% and 26%. An increase was seen especially in emergency care and departments of rheumatology, gastroenterology, hemato-oncology, and psychiatry. Common acute causes associated with elevated ALT were infections (45%), hemato-oncologic conditions (17%), and external reasons (13%), whereas autoimmune diseases (28%), psychiatric conditions (14%), and metabolic-dysfunction associated steatotic liver disease (10%) were common chronic causes. In conclusion, ALT testing increased 3.9-fold while the proportion of increased values remained stable, indicating that increased testing was justified. However, in some departments the testing efficiency was low.

Variation in Pesticide Toxicity in the Western Honey Bee (Apis mellifera) Associated with Consuming Phytochemically Different Monofloral Honeys.

Insecticide toxicity to insect herbivores has long been known to vary across different host plants; this phenomenon has been widely documented in both foliage-feeders and sap-feeders. Species-specific phytochemical content of hostplant tissues is assumed to determine the pattern of induction of insect enzymes that detoxify insecticides, but specific phytochemicals have rarely been linked to host plant-associated variation in pesticide toxicity. Moreover, no studies to date have examined the effects of nectar source identity and phytochemical composition on the toxicity of insecticides to pollinators. In this study, we compared LD50 values for the insecticide bifenthrin, a frequent contaminant of nectar and pollen in agroecosystems, in the western honey bee, Apis mellifera, consuming three phytochemically different monofloral honeys: Nyssa ogeche (tupelo), Robinia pseudoacacia (black locust), and Fagopyrum esculentum (buckwheat). We found that bifenthrin toxicity (LD50) values for honey bees across different honey diets is linked to their species-specific phytochemical content. The profiles of phenolic acids and flavonoids of buckwheat and locust honeys are richer than is the profile of tupelo honey, with buckwheat honey containing the highest total content of phytochemicals and associated with the highest bifenthrin LD50 in honey bees. The vector fitting in the ordination analysis revealed positive correlations between LD50 values and two honey phytochemical richness estimates, Chao1 and Abundance-based Coverage Estimator (ACE). These findings suggest unequal effects among different phytochemicals, consistent with the interpretation that certain compounds, including ones that are rare, may have a more pronounced effect in mitigating pesticide toxicity.

Advantages of a synchrotron light source for fluorescence-detected linear dichroism.

Fluorescence-detected linear dichroism (FD-LD) enables one to collect linear dichroism spectra for oriented fluorophores in the presence of other absorbing species and light scattering. The experiment proceeds by scanning the excitation wavelength and using a filter to collect only emitted photons from the fluorophore. Thus, it has the potential to give data with enhanced selectivity and quality. By using a synchrotron radiation light source and fluorescence-detection, we show data for a range of fluorophores in different orienting environments. Film and flow-oriented FD-LD spectra were collected down to 170 nm. Even for flow-oriented liposomes, we have data collected down to 210 nm. For strongly scattering samples, for example, liposomes, FD-LD has the clear advantage that scattering is absent for the longer wavelength fluorescence photons. The collimated and smaller beam size of the synchrotron radiation also gives rise to sharper and more well-defined features in the spectra.

Linear dichroism reveals the perpendicular orientation of DNA bases in the RecA and Rad51 recombinase filaments: A possible mechanism for the strand exchange reaction.

Linear dichroism spectroscopy is used to investigate the structure of RecA family recombinase filaments (RecA and Rad51 proteins) with DNA for clarifying the molecular mechanism of DNA strand exchange promoted by these proteins and its activation. The measurements show that the recombinases promote the perpendicular base orientation of single-stranded DNA only in the presence of activators, indicating the importance of base orientation in the reaction. We summarize the results and discuss the role of DNA base orientation.

Evaluation of in silico model predictions for mammalian acute oral toxicity and regulatory application in pesticide hazard and risk assessment.

The United States Environmental Protection Agency (USEPA) uses the lethal dose 50% (LD50) value from in vivo rat acute oral toxicity studies for pesticide product label precautionary statements and environmental risk assessment (RA). The Collaborative Acute Toxicity Modeling Suite (CATMoS) is a quantitative structure-activity relationship (QSAR)-based in silico approach to predict rat acute oral toxicity that has the potential to reduce animal use when registering a new pesticide technical grade active ingredient (TGAI). This analysis compared LD50 values predicted by CATMoS to empirical values from in vivo studies for the TGAIs of 177 conventional pesticides. The accuracy and reliability of the model predictions were assessed relative to the empirical data in terms of USEPA acute oral toxicity categories and discrete LD50 values for each chemical. CATMoS was most reliable at placing pesticide TGAIs in acute toxicity categories III (>500-5000 mg/kg) and IV (>5000 mg/kg), with 88% categorical concordance for 165 chemicals with empirical in vivo LD50 values ≥ 500 mg/kg. When considering an LD50 for RA, CATMoS predictions of 2000 mg/kg and higher were found to agree with empirical values from limit tests (i.e., single, high-dose tests) or definitive results over 2000 mg/kg with few exceptions.