A flexible summary statistics-based colocalization method with application to the mucin cystic fibrosis lung disease modifier locus

Mucus obstruction is a central feature in the cystic fibrosis (CF) airways. A genome-wide association study (GWAS) of lung disease by the CF Gene Modifier Consortium (CFGMC) identified a significant locus containing two mucin genes, MUC20 and MUC4. Expression quantitative trait locus (eQTL) analysis using human nasal epithelia (HNE) from 94 CF-affected Canadians in the CFGMC demonstrated MUC4 eQTLs that mirrored the lung association pattern in the region, suggesting that MUC4 expression may mediate CF lung disease. Complications arose, however, with colocalization testing using existing methods: the locus is complex and the associated SNPs span a 0.2 Mb region with high linkage disequilibrium (LD) and evidence of allelic heterogeneity. We previously developed the Simple Sum (SS), a powerful colocalization test in regions with allelic heterogeneity, but SS assumed eQTLs to be present to achieve type I error control. Here we propose a two-stage SS (SS2) colocalization test that avoids a priori eQTL assumptions, accounts for multiple hypothesis testing and the composite null hypothesis, and enables meta-analysis. We compare SS2 to published approaches through simulation and demonstrate type I error control for all settings with the greatest power in the presence of high LD and allelic heterogeneity. Applying SS2 to the MUC20/MUC4 CF lung disease locus with eQTLs from CF HNE revealed significant colocalization with MUC4 (p = 1.31 × 10^?5) rather than with MUC20. The SS2 is a powerful method to inform the responsible gene(s) at a locus and guide future functional studies. SS2 has been implemented in the application LocusFocus.     

沿海养殖池塘对红树林生态系统影响的时空变化监测与分析

随着沿海人类活动的日益加剧,其对红树林生态系统健康和可持续发展的影响也逐渐凸显,实现红树林周边典型人类活动时空动态变化监测对红树林生态系统的保护与修复意义重大。该研究基于Landsat多时相遥感数据和GoogleEarth Engine平台,通过面向对象的机器学习方法,融入水体季节波动信息作为分类特征,获取了1990、2000、2010和2020年4个不同时期中国沿海红树林分布省区(包括广东、福建、浙江、台湾、广西及海南) 30 m分辨率的养殖池塘空间格局及其变化特征,并进一步解析养殖池塘对红树林生态系统的影响。研究结果表明:(1)研究区域内4个时间节点的沿海养殖池塘面积总量分别为2 963、5 200、5 377及4 805 km²,呈先增加后减少的趋势,于2010–2020年间达到峰值。沿海养殖池塘面积变化趋势和达峰时间存在明显区域差异性,其主要原因是红树林保护政策、养殖池塘规范管理和阶段性经济目标的区域差异化。(2)我国沿海养殖池塘集中分布在21°–24°N区域(广东和广西),与红树林沿纬度的分布格局呈错峰分布。其中,红树林与沿海养殖池塘集中分布区(21°–...     

CRBP-III:lacZ expression pattern reveals a novel heterogeneity of vascular endothelial cells

Vascular endothelial cells are structurally and functionally heterogeneous. However, the molecular basis of this heterogeneity remains poorly defined. We used subtractive and differential screening to identify genes that exhibit heterogeneous expression patterns among vascular endothelial cells. One such gene is cellular retinol binding protein III (CRBP-III/Rbp7). Analysis of the lacZ knockin line for this gene (CRBP-III:lacZ) revealed a novel organ-specific vascular endothelial expression pattern. LacZ was expressed in vascular endothelial cells in heart, skeletal muscle, adipose tissues, thymus, and salivary gland. However, it was not detected in other tissues such as brain, liver, and lung. Furthermore, the expression within each organ was primarily restricted to small capillary endothelial cells, but could not be detected in larger vessels. This organ-specific vascular endothelial expression of CRPB:lacZ is relatively resistant to the changes of organ microenvironment. However, the level of expression can be modified by vitamin A deficiency. Therefore, our results provide novel molecular evidence for the heterogeneity of vascular endothelial cells.     

The influence of history and contemporary stream hydrology on the evolution of genetic diversity within species: an examination of microsatellite DNA variation in bull trout,Salvelinus confluentus (Pisces: Salmonidae)

An understanding of the relative roles of historical and contemporary factors in structuring genetic variation is a fundamental, but understudied aspect of geographic variation. We examined geographic variation in microsatellite DNA allele frequencies in bull trout (Salvelinus confluentus, Salmonidae) to test hypotheses concerning the relative roles of postglacial dispersal (historical) and current landscape features (contemporary) in structuring genetic variability and population differentiation. Bull trout exhibit relatively low intrapopulation microsatellite variation (average of 1.9 alleles per locus, average He = 0.24), but high levels of interpopulation divergence (F(ST) = 0.39). We found evidence of historical influences on microsatellite variation in the form of a decrease in the number of alleles and heterozygosities in populations on the periphery of the range relative to populations closer to putative glacial refugia. In addition, one region of British Columbia that was colonized later during deglaciation and by more indirect watershed connections showed less developed and more variable patterns of isolation by distance than a similar region colonized earlier and more directly from refugia. Current spatial and drainage interconnectedness among sites and the presence of migration barriers (falls and cascades) within individual streams were found to be important contemporary factors influencing historical patterns of genetic variability and interpopulation divergence. Our work illustrates the limited utility of equilibrium models to delineate population structure and patterns of genetic diversity in recently founded populations or those inhabiting highly heterogeneous environments, and it highlights the need for approaches incorporating a landscape context for population divergence. Substantial microsatellite DNA divergence among bull trout populations may also signal divergence in traits important to population persistence in specific environments.     

Structure and Conformational Heterogeneity of a Weak Toxin from the Cobra Naja kaouthiaVenom

Resonances in the two-dimensional ¹H NMR spectra of a weak toxin (WTX) from the venom of cobra Naja kaouthiafor all 65 amino acid residues were assigned. The amino acid sequence of WTX, determined by the sequential assignment of spin systems, was found to be similar to that of the CM-9a toxin from the N. kaouthiavenom. Unlike CM-9a, WTX contains an additional Trp36 residue; Lys50 and Tyr52 are interchanged; and there is a Thr residue in place of Arg2. For some residues of WTX, the presence of two components of approximately equal intensities in the spectra was shown, which is explained by the conformational heterogeneity of the polypeptide owing to the cis–transisomerization of the peptide bond Arg32–Pro33. The data (contacts of the nuclear Overhauser effect, constants of spin–spin coupling of protons, and rates of exchange of amide protons for deuterium of the solvent) made it possible to determine the secondary structure of two forms of WTX, which is characterized by the presence of two antiparallel -sheets, one of which consists of two strands (regions 1–5 and 13–17) and the other, of three strands (regions 23–28, 38–43, and 55–59).     

兔左室壁三层心肌细胞的分离及动作电位、钙和钾电流分布的异质性

探讨兔左室壁三层心肌单个细胞的分离方法以及电生理特征,实验以胶原酶按二步消化法分离兔心肌细胞,其中用剃须刀分离左室游离壁内,中,外三层心肌,采用全细胞膜片钳记录AP和离子电流,结果显示：（1）中层细胞上的动作电位时程明显长于内膜下心肌和外膜下心肌,且存在显著的1相切迹和2相驼峰;（2）中层细胞的Ica,L和Lto较内,外膜下的大,IK,s相反,可见三层心肌细胞上多种电流存在显著差异。     

Accelerated molecular evolution in halophilic crustaceans

In contrast to the stable ionic composition of the oceans, inland waters show striking diversity, possessing salt concentrations varying from I mM to 5 M. Although species diversity is highest in fresh water, some lineages have colonized hypersaline environments where they encounter elevated levels of both ultraviolet (UV) radiation and osmotic stress. This study compares rates of evolution in halophilic and freshwater taxa for two groups of microcrustaceans, anostracans and daphniids, from Australia and North America. The results establish that halophilic species show consistent rate acceleration, involving elevated levels of both insertion/deletion events and of nucleotide substitutions. The elevated pace of molecular evolution does not appear to be linked to selection or to other agents that are known to influence the supply rate of mutations, such as UV exposure, generation length, or shifts in metabolic rate. However, variance in ionic strength, which is known to have potent effects on DNA-protein interactions as well as on the structural properties of DNA and proteins, might account for the lowered fidelity of DNA replication in life from hypersaline settings. Regardless of its cause, the consistent rate acceleration in halophiles suggests that past efforts to employ sequence divergences to date events, such as the age of asexual lineages in Artemia, have resulted in serious overestimates. More generally, the results indicate that coordinated shifts in rates of molecular evolution may occur in lineages exposed to extreme environmental conditions.     

Some general points in estimating heterogeneity variance with the DerSimonian-Laird estimator

In this paper we consider estimating heterogeneity variance with the DerSimonian-Laird (DSL) estimator as typically used in meta-analysis. In its general form the DSL estimator requires inverse population-averaged study-specific variances as weights, in which case the estimator is unbiased. It has become common practice, however, to use estimates of the study-specific variances instead of their population-averaged versions. This can lead to considerable bias. Simulations illustrate these findings.     

A geostatistical analysis of small-scale spatial variability in bacterial abundance and community structure in salt marsh creek bank sediments

Small-scale variations in bacterial abundance and community structure were examined in salt marsh sediments from Virginia's eastern shore. Samples were collected at 5 cm intervals (horizontally) along a 50 cm elevation gradient, over a 215 cm horizontal transect. For each sample, bacterial abundance was determined using acridine orange direct counts and community structure was analyzed using randomly amplified polymorphic DNA fingerprinting of whole-community DNA extracts. A geostatistical analysis was used to determine the degree of spatial autocorrelation among the samples, for each variable and each direction (horizontal and vertical). The proportion of variance in bacterial abundance that could be accounted for by the spatial model was quite high (vertical: 60%, horizontal: 73%); significant autocorrelation was found among samples separated by 25 cm in the vertical direction and up to 115 cm horizontally. In contrast, most of the variability in community structure was not accounted for by simply considering the spatial separation of samples (vertical: 11%, horizontal: 22%), and must reflect variability from other parameters (e.g., variation at other spatial scales, experimental error, or environmental heterogeneity). Microbial community patch size based upon overall similarity in community structure varied between 17 cm (vertical) and 35 cm (horizontal). Overall, variability due to horizontal position (distance from the creek bank) was much smaller than that due to vertical position (elevation) for both community properties assayed. This suggests that processes more correlated with elevation (e.g., drainage and redox potential) vary at a smaller scale (therefore producing smaller patch sizes) than processes controlled by distance from the creek bank.     

Rates of gene rearrangement and nucleotide substitution are correlated in the mitochondrial genomes of insects

A number of studies indicated that lineages of animals with high rates of mitochondrial (mt) gene rearrangement might have high rates of mt nucleotide substitution. We chose the hemipteroid assemblage and the Insecta to test the idea that rates of mt gene rearrangement and mt nucleotide substitution are correlated. For this purpose, we sequenced the mt genome of a lepidopsocid from the Psocoptera, the only order of hemipteroid insects for which an entire mtDNA sequence is not available. The mt genome of this lepidopsocid is circular, 16,924 bp long, and contains 37 genes and a putative control region; seven tRNA genes and a protein-coding gene in this genome have changed positions relative to the ancestral arrangement of mt genes of insects. We then compared the relative rates of nucleotide substitution among species from each of the four orders of hemipteroid insects and among the 20 insects whose mt genomes have been sequenced entirely. All comparisons among the hemipteroid insects showed that species with higher rates of gene rearrangement also had significantly higher rates of nucleotide substitution statistically than did species with lower rates of gene rearrangement. In comparisons among the 20 insects, where the mt genomes of the two species differed by more than five breakpoints, the more rearranged species always had a significantly higher rate of nucleotide substitution than the less rearranged species. However, in comparisons where the mt genomes of two species differed by five or less breakpoints, the more rearranged species did not always have a significantly higher rate of nucleotide substitution than the less rearranged species. We tested the statistical significance of the correlation between the rates of mt gene rearrangement and mt nucleotide substitution with nine pairs of insects that were phylogenetically independent from one another. We found that the correlation was positive and statistically significant (R2 = 0.73, P = 0.01; Rs = 0.67, P < 0.05). We propose that increased rates of nucleotide substitution may lead to increased rates of gene rearrangement in the mt genomes of insects.