High-Density SNP Screening of the Major Histocompatibility Complex in Systemic Lupus Erythematosus Demonstrates Strong Evidence for Independent Susceptibility Regions

A substantial genetic contribution to systemic lupus erythematosus (SLE) risk is conferred by major histocompatibility complex (MHC) gene(s) on chromosome 6p21. Previous studies in SLE have lacked statistical power and genetic resolution to fully define MHC influences. We characterized 1,610 Caucasian SLE cases and 1,470 parents for 1,974 MHC SNPs, the highly polymorphic HLA-DRB1 locus, and a panel of ancestry informative markers. Single-marker analyses revealed strong signals for SNPs within several MHC regions, as well as with HLA-DRB1 (global p = 9.99×10⁻¹⁶). The most strongly associated DRB1 alleles were: *0301 (odds ratio, OR = 2.21, p = 2.53×10⁻¹²), *1401 (OR = 0.50, p = 0.0002), and *1501 (OR = 1.39, p = 0.0032). The MHC region SNP demonstrating the strongest evidence of association with SLE was rs3117103, with OR = 2.44 and p = 2.80×10⁻¹³. Conditional haplotype and stepwise logistic regression analyses identified strong evidence for association between SLE and the extended class I, class I, class III, class II, and the extended class II MHC regions. Sequential removal of SLE–associated DRB1 haplotypes revealed independent effects due to variation within OR2H2 (extended class I, rs362521, p = 0.006), CREBL1 (class III, rs8283, p = 0.01), and DQB2 (class II, rs7769979, p = 0.003, and rs10947345, p = 0.0004). Further, conditional haplotype analyses demonstrated that variation within MICB (class I, rs3828903, p = 0.006) also contributes to SLE risk independent of HLA-DRB1*0301. Our results for the first time delineate with high resolution several MHC regions with independent contributions to SLE risk. We provide a list of candidate variants based on biologic and functional considerations that may be causally related to SLE risk and warrant further investigation.     

Systematic Analysis of Circadian Genes in a Population-Based Sample Reveals Association of TIMELESS with Depression and Sleep Disturbance

Disturbances in the circadian pacemaker system are commonly found in individuals with depression and sleep-related problems. We hypothesized that some of the canonical circadian clock genes would be associated with depression accompanied by signs of disturbed sleep, early morning awakening, or daytime fatigue. We tested this hypothesis in a population-based sample of the Health 2000 dataset from Finland, including 384 depressed individuals and 1270 controls, all with detailed information on sleep and daytime vigilance, and analyzed this set of individuals with regard to 113 single-nucleotide polymorphisms of 18 genes of the circadian system. We found significant association between TIMELESS variants and depression with fatigue (D+FAT+) (rs7486220: pointwise P = 0.000099, OR = 1.66; corrected empirical P for the model of D+FAT+  = 0.0056; haplotype ‘C-A-A-C’ of rs2291739-rs2291738-rs7486220-rs1082214: P = 0.0000075, OR = 1.72) in females, and association to depression with early morning awakening (D+EMA+) (rs1082214: pointwise P = 0.0009, OR = 2.70; corrected empirical P = 0.0374 for the model D+EMA+; haplotype ‘G-T’ of rs7486220 and rs1082214: P = 0.0001, OR = 3.01) in males. There was significant interaction of gender and TIMELESS (for example with rs1082214, P = 0.000023 to D+EMA+ and P = 0.005 to D+FAT+). We obtained supported evidence for involvement of TIMELESS in sleeping problems in an independent set of control individuals with seasonal changes in mood, sleep duration, energy level and social activity in females (P = 0.036, ® = 0.123 for rs1082214) and with early morning awakening or fatigue in males (P = 0.038 and P = 0.0016, respectively, for rs1082214). There was also some evidence of interaction between TIMELESS and PER1 in females to D+FAT+ as well as between TIMELESS and ARNTL, RORA or NR1D1 in males to D+EMA+. These findings support a connection between circadian genes and gender-dependent depression and defective sleep regulation.     

Positive selection is the main driving force for evolution of citrus canker-causing Xanthomonas

Understanding the evolutionary history and potential of bacterial pathogens is critical to prevent the emergence of new infectious bacterial diseases. Xanthomonas axonopodis subsp. citri (Xac) (synonym X. citri subsp. citri), which causes citrus canker, is one of the hardest-fought plant bacterial pathogens in US history. Here, we sequenced 21 Xac strains (14 XacA, 3 XacA* and 4 XacA^w) with different host ranges from North America and Asia and conducted comparative genomic and evolutionary analyses. Our analyses suggest that acquisition of beneficial genes and loss of detrimental genes most likely allowed XacA to infect a broader range of hosts as compared with XacA^w and XacA*. Recombination was found to have occurred frequently on the relative ancient branches, but rarely on the young branches of the clonal genealogy. The ratio of recombination/mutation ρ/θ was 0.0790±0.0005, implying that the Xac population was clonal in structure. Positive selection has affected 14% (395 out of 2822) of core genes of the citrus canker-causing Xanthomonas. The genes affected are enriched in ‘carbohydrate transport and metabolism'' and ‘DNA replication, recombination and repair'' genes (P<0.05). Many genes related to virulence, especially genes involved in the type III secretion system and effectors, are affected by positive selection, further highlighting the contribution of positive selection to the evolution of citrus canker-causing Xanthomonas. Our results suggest that both metabolism and virulence genes provide advantages to endow XacA with higher virulence and a wider host range. Our analysis advances our understanding of the genomic basis of specialization by positive selection in bacterial evolution.     

Genetically Dependent ERBB3 Expression Modulates Antigen Presenting Cell Function and Type 1 Diabetes Risk

Type 1 diabetes (T1D) is an autoimmune disease resulting from the complex interaction between multiple susceptibility genes, environmental factors and the immune system. Over 40 T1D susceptibility regions have been suggested by recent genome-wide association studies; however, the specific genes and their role in the disease remain elusive. The objective of this study is to identify the susceptibility gene(s) in the 12q13 region and investigate the functional link to the disease pathogenesis. A total of 19 SNPs in the 12q13 region were analyzed by the TaqMan assay for 1,434 T1D patients and 1,865 controls. Thirteen of the SNPs are associated with T1D (best p = 4×10⁻¹¹), thus providing confirmatory evidence for at least one susceptibility gene in this region. To identify candidate genes, expression of six genes in the region was analyzed by real-time RT-PCR for PBMCs from 192 T1D patients and 192 controls. SNP genotypes in the 12q13 region are the main factors that determine ERBB3 mRNA levels in PBMCs. The protective genotypes for T1D are associated with higher ERBB3 mRNA level (p<10⁻¹⁰). Furthermore, ERBB3 protein is expressed on the surface of CD11c⁺ cells (dendritic cells and monocytes) in peripheral blood after stimulation with LPS, polyI:C or CpG. Subjects with protective genotypes have significantly higher percentages of ERBB3⁺ monocytes and dendritic cells (p = 1.1×10⁻⁹); and the percentages of ERBB3⁺ cells positively correlate with the ability of APC to stimulate T cell proliferation (R² = 0.90, p<0.0001). Our results indicate that ERBB3 plays a critical role in determining APC function and potentially T1D pathogenesis.     

Evolutionary Genomics Reveals Lineage-Specific Gene Loss and Rapid Evolution of a Sperm-Specific Ion Channel Complex: CatSpers and CatSperβ

The mammalian CatSper ion channel family consists of four sperm-specific voltage-gated Ca²⁺ channels that are crucial for sperm hyperactivation and male fertility. All four CatSper subunits are believed to assemble into a heteromultimeric channel complex, together with an auxiliary subunit, CatSperβ. Here, we report a comprehensive comparative genomics study and evolutionary analysis of CatSpers and CatSperβ, with important correlation to physiological significance of molecular evolution of the CatSper channel complex. The development of the CatSper channel complex with four CatSpers and CatSperβ originated as early as primitive metazoans such as the Cnidarian Nematostella vectensis. Comparative genomics revealed extensive lineage-specific gene loss of all four CatSpers and CatSperβ through metazoan evolution, especially in vertebrates. The CatSper channel complex underwent rapid evolution and functional divergence, while distinct evolutionary constraints appear to have acted on different domains and specific sites of the four CatSper genes. These results reveal unique evolutionary characteristics of sperm-specific Ca²⁺ channels and their adaptation to sperm biology through metazoan evolution.     

Five Quantitative Trait Loci Control Radiation-Induced Adenoma Multiplicity in Mom1R ApcMin/+ Mice

Ionising radiation is a carcinogen capable of inducing tumours, including colorectal cancer, in both humans and animals. By backcrossing a recombinant line of Apc^Min/+ mice to the inbred BALB/c mouse strain, which is unusually sensitive to radiation–induced tumour development, we obtained panels of 2Gy-irradiated and sham-irradiated N2 Apc^Min/+ mice for genotyping with a genome-wide panel of microsatellites at ∼15 cM density and phenotyping by counting adenomas in the small intestine. Interval and composite interval mapping along with permutation testing identified five significant susceptibility quantitative trait loci (QTLs) responsible for radiation induced tumour multiplicity in the small intestine. These were defined as Mom (Modifier of Min) radiation-induced polyposis (Mrip1-5) on chromosome 2 (log of odds, LOD 2.8, p = 0.0003), two regions within chromosome 5 (LOD 5.2, p<0.00001, 6.2, p<0.00001) and two regions within chromosome 16 respectively (LOD 4.1, p  = 4×10⁻⁵, 4.8, p<0.00001). Suggestive QTLs were found for sham-irradiated mice on chromosomes 3, 6 and 13 (LOD 1.7, 1.5 and 2.0 respectively; p<0.005). Genes containing BALB/c specific non-synonymous polymorphisms were identified within Mrip regions and prediction programming used to locate potentially functional polymorphisms. Our study locates the QTL regions responsible for increased radiation-induced intestinal tumorigenesis in Apc^Min/+ mice and identifies candidate genes with predicted functional polymorphisms that are involved in spindle checkpoint and chromosomal stability (Bub1b, Casc5, and Bub1), DNA repair (Recc1 and Prkdc) or inflammation (Duox2, Itgb2l and Cxcl5). Our study demonstrates use of in silico analysis in candidate gene identification as a way of reducing large-scale backcross breeding programmes.     

Increased Frequency and Compromised Function of T Regulatory Cells in Systemic Sclerosis (SSc) Is Related to a Diminished CD69 and TGFβ Expression

Background

Regulatory T cells (Tregs) are essential in the control of tolerance. Evidence implicates Tregs in human autoimmune conditions. Here we investigated their role in systemic sclerosis (SSc).

Methods/Principal Findings

Patients were subdivided as having limited cutaneous SSc (lcSSc, n = 20) or diffuse cutaneous SSc (dcSSc, n = 48). Further subdivision was made between early dcSSc (n = 24) and late dcSSc (n = 24) based upon the duration of disease. 26 controls were studied for comparison. CD3+ cells were isolated using FACS and subsequently studied for the expression of CD4, CD8, CD25, FoxP3, CD127, CD62L, GITR, CD69 using flow cytometry. T cell suppression assays were performed using sorted CD4CD25^highCD127^- and CD4CD25^lowCD127^high and CD3⁺ cells. Suppressive function was correlated with CD69 surface expression and TGFβ secretion/expression. The frequency of CD4⁺CD25⁺ and CD25^highFoxP3^highCD127^neg T cells was highly increased in all SSc subgroups. Although the expression of CD25 and GITR was comparable between groups, expression of CD62L and CD69 was dramatically lower in SSc patients, which correlated with a diminished suppressive function. Co-incubation of Tregs from healthy donors with plasma from SSc patients fully abrogated suppressive activity. Activation of Tregs from healthy donors or SSc patients with PHA significantly up regulated CD69 expression that could be inhibited by SSc plasma.

Conclusions/Significance

These results indicate that soluble factors in SSc plasma inhibit Treg function specifically that is associated with altered Treg CD69 and TGFβ expression. These data suggest that a defective Treg function may underlie the immune dysfunction in systemic sclerosis.     

A Major Histocompatibility Class I Locus Contributes to Multiple Sclerosis Susceptibility Independently from HLA-DRB1*15:01

Background

In Northern European descended populations, genetic susceptibility for multiple sclerosis (MS) is associated with alleles of the human leukocyte antigen (HLA) Class II gene DRB1. Whether other major histocompatibility complex (MHC) genes contribute to MS susceptibility is controversial.

Methodology/Principal Findings

A case control analysis was performed using 958 single nucleotide polymorphisms (SNPs) spanning the MHC assayed in two independent datasets. The discovery dataset consisted of 1,018 cases and 1,795 controls and the replication dataset was composed of 1,343 cases and 1,379 controls. The most significantly MS-associated SNP in the discovery dataset was rs3135391, a Class II SNP known to tag the HLA-DRB1*15:01 allele, the primary MS susceptibility allele in the MHC (O.R. = 3.04, p<1×10⁻⁷⁸). To control for the effects of the HLA-DRB1*15:01 haplotype, case control analysis was performed adjusting for this HLA-DRB1*15:01 tagging SNP. After correction for multiple comparisons (false discovery rate = .05) 52 SNPs in the Class I, II and III regions were significantly associated with MS susceptibility in both datasets using the Cochran Armitage trend test. The discovery and replication datasets were merged and subjects carrying the HLA-DRB1*15:01 tagging SNP were excluded. Association tests showed that 48 of the 52 replicated SNPs retained significant associations with MS susceptibility independently of the HLA-DRB1*15:01 as defined by the tagging SNP. 20 Class I SNPs were associated with MS susceptibility with p-values ≤1×10⁻⁸. The most significantly associated SNP was rs4959039, a SNP in the downstream un-translated region of the non-classical HLA-G gene (Odds ratio 1.59, 95% CI 1.40, 1.81, p = 8.45×10⁻¹³) and is in linkage disequilibrium with several nearby SNPs. Logistic regression modeling showed that this SNP''s contribution to MS susceptibility was independent of the Class II and Class III SNPs identified in this screen.

Conclusions

A MHC Class I locus contributes to MS susceptibility independently of the HLA-DRB1*15:01 haplotype.     

Prospectively Isolated Cancer-Associated CD10+ Fibroblasts Have Stronger Interactions with CD133+ Colon Cancer Cells than with CD133? Cancer Cells

Although CD133 has been reported to be a promising colon cancer stem cell marker, the biological functions of CD133⁺ colon cancer cells remain controversial. In the present study, we investigated the biological differences between CD133⁺ and CD133⁻ colon cancer cells, with a particular focus on their interactions with cancer-associated fibroblasts, especially CD10⁺ fibroblasts. We used 19 primary colon cancer tissues, 30 primary cultures of fibroblasts derived from colon cancer tissues and 6 colon cancer cell lines. We isolated CD133⁺ and CD133⁻ subpopulations from the colon cancer tissues and cultured cells. In vitro analyses revealed that the two populations showed similar biological behaviors in their proliferation and chemosensitivity. In vivo analyses revealed that CD133⁺ cells showed significantly greater tumor growth than CD133⁻ cells (P = 0.007). Moreover, in cocultures with primary fibroblasts derived from colon cancer tissues, CD133⁺ cells exhibited significantly more invasive behaviors than CD133⁻ cells (P<0.001), especially in cocultures with CD10⁺ fibroblasts (P<0.0001). Further in vivo analyses revealed that CD10⁺ fibroblasts enhanced the tumor growth of CD133⁺ cells significantly more than CD10⁻ fibroblasts (P<0.05). These data demonstrate that the in vitro invasive properties and in vivo tumor growth of CD133⁺ colon cancer cells are enhanced in the presence of specific cancer-associated fibroblasts, CD10⁺ fibroblasts, suggesting that the interactions between these specific cell populations have important roles in cancer progression. Therefore, these specific interactions may be promising targets for new colon cancer therapies.     

Genome-Wide Analysis of DNA Methylation Differences in Muscle and Fat from Monozygotic Twins Discordant for Type 2 Diabetes

Background

Monozygotic twins discordant for type 2 diabetes constitute an ideal model to study environmental contributions to type 2 diabetic traits. We aimed to examine whether global DNA methylation differences exist in major glucose metabolic tissues from these twins.

Methodology/Principal Findings

Skeletal muscle (n = 11 pairs) and subcutaneous adipose tissue (n = 5 pairs) biopsies were collected from 53–80 year-old monozygotic twin pairs discordant for type 2 diabetes. DNA methylation was measured by microarrays at 26,850 cytosine-guanine dinucleotide (CpG) sites in the promoters of 14,279 genes. Bisulfite sequencing was applied to validate array data and to quantify methylation of intergenic repetitive DNA sequences. The overall intra-pair variation in DNA methylation was large in repetitive (LINE1, D4Z4 and NBL2) regions compared to gene promoters (standard deviation of intra-pair differences: 10% points vs. 4% points, P<0.001). Increased variation of LINE1 sequence methylation was associated with more phenotypic dissimilarity measured as body mass index (r = 0.77, P = 0.007) and 2-hour plasma glucose (r = 0.66, P = 0.03) whereas the variation in promoter methylation did not associate with phenotypic differences. Validated methylation changes were identified in the promoters of known type 2 diabetes-related genes, including PPARGC1A in muscle (13.9±6.2% vs. 9.0±4.5%, P = 0.03) and HNF4A in adipose tissue (75.2±3.8% vs. 70.5±3.7%, P<0.001) which had increased methylation in type 2 diabetic individuals. A hypothesis-free genome-wide exploration of differential methylation without correction for multiple testing identified 789 and 1,458 CpG sites in skeletal muscle and adipose tissue, respectively. These methylation changes only reached some percentage points, and few sites passed correction for multiple testing.

Conclusions/Significance

Our study suggests that likely acquired DNA methylation changes in skeletal muscle or adipose tissue gene promoters are quantitatively small between type 2 diabetic and non-diabetic twins. The importance of methylation changes in candidate genes such as PPARGC1A and HNF4A should be examined further by replication in larger samples.     

Properties of an Inwardly Rectifying ATP-sensitive K+ Channel in the Basolateral Membrane of Renal Proximal Tubule

The potassium conductance of the basolateral membrane (BLM) of proximal tubule cells is a critical regulator of transport since it is the major determinant of the negative cell membrane potential and is necessary for pump-leak coupling to the Na⁺,K⁺-ATPase pump. Despite this pivotal physiological role, the properties of this conductance have been incompletely characterized, in part due to difficulty gaining access to the BLM. We have investigated the properties of this BLM K⁺ conductance in dissociated, polarized Ambystoma proximal tubule cells. Nearly all seals made on Ambystoma cells contained inward rectifier K⁺ channels (γ_{slope, in} = 24.5 ± 0.6 pS, γ_{chord, out} = 3.7 ± 0.4 pS). The rectification is mediated in part by internal Mg²⁺. The open probability of the channel increases modestly with hyperpolarization. The inward conducting properties are described by a saturating binding–unbinding model. The channel conducts Tl⁺ and K⁺, but there is no significant conductance for Na⁺, Rb⁺, Cs⁺, Li⁺, NH₄⁺, or Cl⁻. The channel is inhibited by barium and the sulfonylurea agent glibenclamide, but not by tetraethylammonium. Channel rundown typically occurs in the absence of ATP, but cytosolic addition of 0.2 mM ATP (or any hydrolyzable nucleoside triphosphate) sustains channel activity indefinitely. Phosphorylation processes alone fail to sustain channel activity. Higher doses of ATP (or other nucleoside triphosphates) reversibly inhibit the channel. The K⁺ channel opener diazoxide opens the channel in the presence of 0.2 mM ATP, but does not alleviate the inhibition of millimolar doses of ATP. We conclude that this K⁺ channel is the major ATP-sensitive basolateral K⁺ conductance in the proximal tubule.     

Complex Adaptations Can Drive the Evolution of the Capacitor [PSI+], Even with Realistic Rates of Yeast Sex

The [PSI⁺] prion may enhance evolvability by revealing previously cryptic genetic variation, but it is unclear whether such evolvability properties could be favored by natural selection. Sex inhibits the evolution of other putative evolvability mechanisms, such as mutator alleles. This paper explores whether sex also prevents natural selection from favoring modifier alleles that facilitate [PSI⁺] formation. Sex may permit the spread of “cheater” alleles that acquire the benefits of [PSI⁺] through mating without incurring the cost of producing [PSI⁺] at times when it is not adaptive. Using recent quantitative estimates of the frequency of sex in Saccharomyces paradoxus, we calculate that natural selection for evolvability can drive the evolution of the [PSI⁺] system, so long as yeast populations occasionally require complex adaptations involving synergistic epistasis between two loci. If adaptations are always simple and require substitution at only a single locus, then the [PSI⁺] system is not favored by natural selection. Obligate sex might inhibit the evolution of [PSI⁺]-like systems in other species.     

Heat-shock promoters: targets for evolution by P transposable elements in Drosophila

Transposable elements are potent agents of genomic change during evolution, but require access to chromatin for insertion—and not all genes provide equivalent access. To test whether the regulatory features of heat-shock genes render their proximal promoters especially susceptible to the insertion of transposable elements in nature, we conducted an unbiased screen of the proximal promoters of 18 heat-shock genes in 48 natural populations of Drosophila. More than 200 distinctive transposable elements had inserted into these promoters; greater than 96% are P elements. By contrast, few or no P element insertions segregate in natural populations in a “negative control” set of proximal promoters lacking the distinctive regulatory features of heat-shock genes. P element transpositions into these same genes during laboratory mutagenesis recapitulate these findings. The natural P element insertions cluster in specific sites in the promoters, with up to eight populations exhibiting P element insertions at the same position; laboratory insertions are into similar sites. By contrast, a “positive control” set of promoters resembling heat-shock promoters in regulatory features harbors few P element insertions in nature, but many insertions after experimental transposition in the laboratory. We conclude that the distinctive regulatory features that typify heat-shock genes (in Drosophila) are especially prone to mutagenesis via P elements in nature. Thus in nature, P elements create significant and distinctive variation in heat-shock genes, upon which evolutionary processes may act.     

Apoptosis of CD4+CD25high T Cells in Type 1 Diabetes May Be Partially Mediated by IL-2 Deprivation

Background

Type 1 diabetes (T1D) is a T-cell mediated autoimmune disease targeting the insulin-producing pancreatic β cells. Naturally occurring FOXP3⁺CD4⁺CD25^high regulatory T cells (T_regs) play an important role in dominant tolerance, suppressing autoreactive CD4⁺ effector T cell activity. Previously, in both recent-onset T1D patients and β cell antibody-positive at-risk individuals, we observed increased apoptosis and decreased function of polyclonal T_regs in the periphery. Our objective here was to elucidate the genes and signaling pathways triggering apoptosis in T_regs from T1D subjects.

Principal Findings

Gene expression profiles of unstimulated T_regs from recent-onset T1D (n = 12) and healthy control subjects (n = 15) were generated. Statistical analysis was performed using a Bayesian approach that is highly efficient in determining differentially expressed genes with low number of replicate samples in each of the two phenotypic groups. Microarray analysis showed that several cytokine/chemokine receptor genes, HLA genes, GIMAP family genes and cell adhesion genes were downregulated in T_regs from T1D subjects, relative to control subjects. Several downstream target genes of the AKT and p53 pathways were also upregulated in T1D subjects, relative to controls. Further, expression signatures and increased apoptosis in T_regs from T1D subjects partially mirrored the response of healthy T_regs under conditions of IL-2 deprivation. CD4⁺ effector T-cells from T1D subjects showed a marked reduction in IL-2 secretion. This could indicate that prior to and during the onset of disease, T_regs in T1D may be caught up in a relatively deficient cytokine milieu.

Conclusions

In summary, expression signatures in T_regs from T1D subjects reflect a cellular response that leads to increased sensitivity to apoptosis, partially due to cytokine deprivation. Further characterization of these signaling cascades should enable the detection of genes that can be targeted for restoring T_reg function in subjects predisposed to T1D.     

Complex Adaptations Can Drive the Evolution of the Capacitor [PSI +], Even with Realistic Rates of Yeast Sex

The [PSI⁺] prion may enhance evolvability by revealing previously cryptic genetic variation, but it is unclear whether such evolvability properties could be favored by natural selection. Sex inhibits the evolution of other putative evolvability mechanisms, such as mutator alleles. This paper explores whether sex also prevents natural selection from favoring modifier alleles that facilitate [PSI⁺] formation. Sex may permit the spread of “cheater” alleles that acquire the benefits of [PSI⁺] through mating without incurring the cost of producing [PSI⁺] at times when it is not adaptive. Using recent quantitative estimates of the frequency of sex in Saccharomyces paradoxus, we calculate that natural selection for evolvability can drive the evolution of the [PSI⁺] system, so long as yeast populations occasionally require complex adaptations involving synergistic epistasis between two loci. If adaptations are always simple and require substitution at only a single locus, then the [PSI⁺] system is not favored by natural selection. Obligate sex might inhibit the evolution of [PSI⁺]-like systems in other species.     

Loss of the IR region in conifer plastomes: Changes in the selection pressure and substitution rate of protein‐coding genes

Plastid genomes (plastomes) have a quadripartite structure, but some species have drastically reduced or lost inverted repeat (IR) regions. IR regions are important for genome stability and the evolution rate. In the evolutionary process of gymnosperms, the typical IRs of conifers were lost, possibly affecting the evolutionary rate and selection pressure of genomic protein‐coding genes. In this study, we selected 78 gymnosperm species (51 genera, 13 families) for evolutionary analysis. The selection pressure analysis results showed that negative selection effects were detected in all 50 common genes. Among them, six genes in conifers had higher ω values than non‐conifers, and 12 genes had lower ω values. The evolutionary rate analysis results showed that 9 of 50 common genes differed between conifers and non‐conifers. It is more obvious that in non‐conifers, the rates of psbA (trst, trsv, ratio, dN, dS, and ω) were 2.6‐ to 3.1‐fold of conifers. In conifers, trsv, ratio, dN, dS, and ω of ycf2 were 1.2‐ to 3.6‐fold of non‐conifers. In addition, the evolution rate of ycf2 in the IR was significantly reduced. psbA is undergoing dynamic change, with an abnormally high evolution rate as a small portion of it enters the IR region. Although conifers have lost the typical IR regions, we detected no change in the substitution rate or selection pressure of most protein‐coding genes due to gene function, plant habitat, or newly acquired IRs.     

Systemic Simvastatin Rescues Retinal Ganglion Cells from Optic Nerve Injury Possibly through Suppression of Astroglial NF-κB Activation

Neuroinflammation is involved in the death of retinal ganglion cells (RGCs) after optic nerve injury. The purpose of this study was to determine whether systemic simvastatin can suppress neuroinflammation in the optic nerve and rescue RGCs after the optic nerve is crushed. Simvastatin or its vehicle was given through an osmotic minipump beginning one week prior to the crushing. Immunohistochemistry and real-time PCR were used to determine the degree of neuroinflammation on day 3 after the crushing. The density of RGCs was determined in Tuj-1 stained retinal flat mounts on day 7. The effect of simvastain on the TNF-α-induced NF-κB activation was determined in cultured optic nerve astrocytes. On day 3, CD68-positive cells, most likely microglia/macrophages, were accumulated at the crushed site. Phosphorylated NF-κB was detected in some astrocytes at the border of the lesion where the immunoreactivity to MCP-1 was intensified. There was an increase in the mRNA levels of the CD68 (11.4-fold), MCP-1 (22.6-fold), ET-1 (2.3-fold), GFAP (1.6-fold), TNF-α (7.0-fold), and iNOS (14.8-fold) genes on day 3. Systemic simvastatin significantly reduced these changes. The mean ± SD number of RGCs was 1816.3±232.6/mm² (n = 6) in the sham controls which was significantly reduced to 831.4±202.5/mm² (n = 9) on day 7 after the optic nerve was crushed. This reduction was significantly suppressed to 1169.2±201.3/mm² (P = 0.01, Scheffe; n = 9) after systemic simvastatin. Simvastatin (1.0 µM) significantly reduced the TNF-α-induced NF-κB activation in cultured optic nerve astrocytes. We conclude that systemic simvastatin can reduce the death of RGCs induced by crushing the optic nerve possibly by suppressing astroglial NF-κB activation.