Neofunctionalization of Mitochondrial Proteins and Incorporation into Signaling Networks in Plants.

Because of their symbiotic origin, many mitochondrial proteins are well conserved across eukaryotic kingdoms. It is however less obvious how specific lineages have obtained novel nuclear-encoded mitochondrial proteins. Here, we report a case of mitochondrial neofunctionalization in plants. Phylogenetic analysis of genes containing the Domain of Unknown Function 295 (DUF295) revealed that the domain likely originated in Angiosperms. The C-terminal DUF295 domain is usually accompanied by an N-terminal F-box domain, involved in ubiquitin ligation via binding with ASK1/SKP1-type proteins. Due to gene duplication, the gene family has expanded rapidly, with 94 DUF295-related genes in Arabidopsis thaliana alone. Two DUF295 family subgroups have uniquely evolved and quickly expanded within Brassicaceae. One of these subgroups has completely lost the F-box, but instead obtained strongly predicted mitochondrial targeting peptides. We show that several representatives of this DUF295 Organellar group are effectively targeted to plant mitochondria and chloroplasts. Furthermore, many DUF295 Organellar genes are induced by mitochondrial dysfunction, whereas F-Box DUF295 genes are not. In agreement, several Brassicaceae-specific DUF295 Organellar genes were incorporated in the evolutionary much older ANAC017-dependent mitochondrial retrograde signaling pathway. Finally, a representative set of DUF295 T-DNA insertion mutants was created. No obvious aberrant phenotypes during normal growth and mitochondrial dysfunction were observed, most likely due to the large extent of gene duplication and redundancy. Overall, this study provides insight into how novel mitochondrial proteins can be created via "intercompartmental" gene duplication events. Moreover, our analysis shows that these newly evolved genes can then be specifically integrated into relevant, pre-existing coexpression networks.

Evolution of P2A and P5A ATPases: ancient gene duplications and the red algal connection to green plants revisited.

In a search for slowly evolving nuclear genes that may cast light on the deep evolution of plants, we carried out phylogenetic analyses of two well-characterized subfamilies of P-type pumps (P2A and P5A ATPases) from representative branches of the eukaryotic tree of life. Both P-type ATPase genes were duplicated very early in eukaryotic evolution and before the divergence of the present eukaryotic supergroups. Synapomorphies identified in the sequences provide evidence that green plants and red algae are more distantly related than are green plants and eukaryotic supergroups in which secondary or tertiary plastids are common, such as several groups belonging to the clade that includes Stramenopiles, Alveolata, Rhizaria, Cryptophyta and Haptophyta (SAR). We propose that red algae branched off soon after the first photosynthesizing eukaryote had acquired a primary plastid, while in another lineage that led to SAR, the primary plastid was lost but, in some cases, regained as a secondary or tertiary plastid.

Genome Sequencing Reveals the Origin of the Allotetraploid Arabidopsis suecica.

Polyploidy is an example of instantaneous speciation when it involves the formation of a new cytotype that is incompatible with the parental species. Because new polyploid individuals are likely to be rare, establishment of a new species is unlikely unless polyploids are able to reproduce through self-fertilization (selfing), or asexually. Conversely, selfing (or asexuality) makes it possible for polyploid species to originate from a single individual-a bona fide speciation event. The extent to which this happens is not known. Here, we consider the origin of Arabidopsis suecica, a selfing allopolyploid between Arabidopsis thaliana and Arabidopsis arenosa, which has hitherto been considered to be an example of a unique origin. Based on whole-genome re-sequencing of 15 natural A. suecica accessions, we identify ubiquitous shared polymorphism with the parental species, and hence conclusively reject a unique origin in favor of multiple founding individuals. We further estimate that the species originated after the last glacial maximum in Eastern Europe or central Eurasia (rather than Sweden, as the name might suggest). Finally, annotation of the self-incompatibility loci in A. suecica revealed that both loci carry non-functional alleles. The locus inherited from the selfing A. thaliana is fixed for an ancestral non-functional allele, whereas the locus inherited from the outcrossing A. arenosa is fixed for a novel loss-of-function allele. Furthermore, the allele inherited from A. thaliana is predicted to transcriptionally silence the allele inherited from A. arenosa, suggesting that loss of self-incompatibility may have been instantaneous.

Freshwater protists do not go with the flow: population structure in Gonyostomum semen independent of connectivity among lakes.

Many recent studies have found genetically differentiated populations in microorganisms despite potentially high dispersal. We designed a study to specifically examine the importance of physical dispersal barriers, i.e. geographic distance and lack of hydrological connectivity, in restricting gene flow and enhancing divergence in limnic microorganisms. We focused on the nuisance microalga Gonyostomum semen, which has recently expanded in Northern Europe and differentiated into genetically distinct populations. G. semen was sampled from six lakes distributed in two adjacent watersheds, which thereby comprised, both connected and non-connected lakes. The individual isolates were genotyped by amplified fragment length polymorphism. Several lake populations were differentiated from each other, but connectivity within watersheds could not explain the observed population genetic pattern. However, isolation by distance was moderate and might limit the gene flow among distant populations. In addition, we found low, but significant linkage disequilibrium, which indicates regular sexual recombination in this species, despite its high degree of asexual reproduction. Therefore, we conclude that the genetic properties of microalgae with occasional sexual reproduction essentially mirror regularly recombining species. Furthermore, the data indicated bottlenecks supporting the hypothesized recent range expansion of this species.

Genetic variation in KLK2 and KLK3 is associated with concentrations of hK2 and PSA in serum and seminal plasma in young men.

BACKGROUND: Genetic variants in KLK2 and KLK3 have been associated with increased serum concentrations of their encoded proteins, human kallikrein-related peptidase 2 (hK2) and prostate-specific antigen (PSA), and with prostate cancer in older men. Low PSA concentrations in seminal plasma (SP) have been associated with low sperm motility. To evaluate whether KLK2 and KLK3 genetic variants affect physiological prostatic secretion, we studied the association of SNPs with hK2 and PSA concentrations in SP and serum of young, healthy men. METHODS: Leukocyte DNA was extracted from 303 male military conscripts (median age 18.1 years). Nine SNPs across KLK2-KLK3 were genotyped. We measured PSA and hK2 in SP and serum using immunofluorometric assays. The association of genotype frequencies with hK2 and PSA concentrations was tested with the Kruskal-Wallis test. RESULTS: Four KLK2 SNPs (rs198972, rs198977, rs198978, and rs80050017) were strongly associated with hK2 concentrations in SP and serum, with individuals homozygous for the major alleles having 3- to 7-fold higher concentrations than the intermediate concentrations found in other homozygotes and heterozygotes (all P < 0.001). Three of these SNPs were significantly associated with percentage of free PSA (%fPSA) in serum (all P < 0.007). Three KLK3 SNPs showed associations with PSA in SP, and the rs1058205 SNP was associated with total PSA in serum (P = 0.001) and %fPSA (P = 0.015). CONCLUSIONS: Associations observed in young, healthy men between the SP and serum concentrations of hK2 and PSA and several genetic variants in KLK2 and KLK3 could be useful to refine models of PSA cutoff values in prostate cancer testing.

Investigating highly replicated asthma genes as candidate genes for allergic rhinitis.

BACKGROUND: Asthma genetics has been extensively studied and many genes have been associated with the development or severity of this disease. In contrast, the genetic basis of allergic rhinitis (AR) has not been evaluated as extensively. It is well known that asthma is closely related with AR since a large proportion of individuals with asthma also present symptoms of AR, and patients with AR have a 5-6 fold increased risk of developing asthma. Thus, the relevance of asthma candidate genes as predisposing factors for AR is worth investigating. The present study was designed to investigate if SNPs in highly replicated asthma genes are associated with the occurrence of AR. METHODS: A total of 192 SNPs from 21 asthma candidate genes reported to be associated with asthma in 6 or more unrelated studies were genotyped in a Swedish population with 246 AR patients and 431 controls. Genotypes for 429 SNPs from the same set of genes were also extracted from a Singapore Chinese genome-wide dataset which consisted of 456 AR cases and 486 controls. All SNPs were subsequently analyzed for association with AR and their influence on allergic sensitization to common allergens. RESULTS: A limited number of potential associations were observed and the overall pattern of P-values corresponds well to the expectations in the absence of an effect. However, in the tests of allele effects in the Chinese population the number of significant P-values exceeds the expectations. The strongest signals were found for SNPs in NPSR1 and CTLA4. In these genes, a total of nine SNPs showed P-values <0.001 with corresponding Q-values <0.05. In the NPSR1 gene some P-values were lower than the Bonferroni correction level. Reanalysis after elimination of all patients with asthmatic symptoms excluded asthma as a confounding factor in our results. Weaker indications were found for IL13 and GSTP1 with respect to sensitization to birch pollen in the Swedish population. CONCLUSIONS: Genetic variation in the majority of the highly replicated asthma genes were not associated to AR in our populations which suggest that asthma and AR could have less in common than previously anticipated. However, NPSR1 and CTLA4 can be genetic links between AR and asthma and associations of polymorphisms in NPSR1 with AR have not been reported previously.

Toll-like receptor gene polymorphisms are associated with allergic rhinitis: a case control study.

BACKGROUND: The Toll-like receptor proteins are important in host defense and initiation of the innate and adaptive immune responses. A number of studies have identified associations between genetic variation in the Toll-like receptor genes and allergic disorders such as asthma and allergic rhinitis. The present study aim to search for genetic variation associated with allergic rhinitis in the Toll-like receptor genes. METHODS: A first association analysis genotyped 73 SNPs in 182 cases and 378 controls from a Swedish population. Based on these results an additional 24 SNPs were analyzed in one Swedish population with 352 cases and 709 controls and one Chinese population with 948 cases and 580 controls. RESULTS: The first association analysis identified 4 allergic rhinitis-associated SNPs in the TLR7-TLR8 gene region. Subsequent analysis of 24 SNPs from this region identified 7 and 5 significant SNPs from the Swedish and Chinese populations, respectively. The corresponding risk-associated haplotypes are significant after Bonferroni correction and are the most common haplotypes in both populations. The associations are primarily detected in females in the Swedish population, whereas it is seen in males in the Chinese population. Further independent support for the involvement of this region in allergic rhinitis was obtained from quantitative skin prick test data generated in both populations. CONCLUSIONS: Haplotypes in the TLR7-TLR8 gene region were associated with allergic rhinitis in one Swedish and one Chinese population. Since this region has earlier been associated with asthma and allergic rhinitis in a Danish linkage study this speaks strongly in favour of this region being truly involved in the development of this disease.

Variation in the VWF gene in Swedish patients with type 1 von Willebrand Disease.

The spectrum of mutations in the von Willebrand factor (VWF) gene in a Swedish type 1 von Willebrand disease (VWD) population was investigated. To gain more knowledge about the dynamics of VWD mutations, the data were analyzed from a population genetics perspective. The VWF gene was resequenced in 54 Swedish patients diagnosed with type 1 VWD. Fifty-five variable sites were located in exons, 10 in the promoter and 38 in introns. The spectrum of mutations was similar to a European study, but included 10 new candidate mutations. The synonymous sites were evenly distributed along the coding sequence, whereas nonsynonymous sites were located into three clusters. Overall, 44% of patients had no mutations or candidate mutations and no promoter haplotype was significantly associated with disease. In 11 patients (20%), more than one mutation or candidate mutation was detected. The allelic identity for the putative disease-causing mutations was approximately 0.1, compatible with an overall disease frequency of 1%. VWF sequences for exon 28 from eight monkey species were compared with the variable positions found in our patients. Positions classified as mutations were overrepresented among sites that were fixed in all eight monkey species. No general increase of the mutation rate was found for the pseudogene region.

Natural variation in autumn expression is the major adaptive determinant distinguishing Arabidopsis FLC haplotypes.

In Arabidopsis thaliana, winter is registered during vernalization through the temperature-dependent repression and epigenetic silencing of floral repressor FLOWERING LOCUS C (FLC). Natural Arabidopsis accessions show considerable variation in vernalization. However, which aspect of the FLC repression mechanism is most important for adaptation to different environments is unclear. By analysing FLC dynamics in natural variants and mutants throughout winter in three field sites, we find that autumnal FLC expression, rather than epigenetic silencing, is the major variable conferred by the distinct Arabidopsis FLChaplotypes. This variation influences flowering responses of Arabidopsis accessions resulting in an interplay between promotion and delay of flowering in different climates to balance survival and, through a post-vernalization effect, reproductive output. These data reveal how expression variation through non-coding cis variation at FLC has enabled Arabidopsis accessions to adapt to different climatic conditions and year-on-year fluctuations.

High Diversity of Microcystin Chemotypes within a Summer Bloom of the Cyanobacterium Microcystis botrys.

The fresh-water cyanobacterium Microcystis is known to form blooms world-wide, and is often responsible for the production of microcystins found in lake water. Microcystins are non-ribosomal peptides with toxic effects, e.g. on vertebrates, but their function remains largely unresolved. Moreover, not all strains produce microcystins, and many different microcystin variants have been described. Here we explored the diversity of microcystin variants within Microcystis botrys, a common bloom-former in Sweden. We isolated a total of 130 strains through the duration of a bloom in eutrophic Lake Vomb, and analyzed their microcystin profiles with tandem mass spectrometry (LC-MS/MS). We found that microcystin producing (28.5%) and non-producing (71.5%) M. botrys strains, co-existed throughout the bloom. However, microcystin producing strains were more prevalent towards the end of the sampling period. Overall, 26 unique M. botrys chemotypes were identified, and while some chemotypes re-occurred, others were found only once. The M. botrys chemotypes showed considerable variation both in terms of number of microcystin variants, as well as in what combinations the variants occurred. To our knowledge, this is the first report on microcystin chemotype variation and dynamics in M. botrys. In addition, our study verifies the co-existence of microcystin and non-microcystin producing strains, and we propose that environmental conditions may be implicated in determining their composition.

RNAi as a Tool to Study Virulence in the Pathogenic Yeast Candida glabrata.

The yeast Candida glabrata is a major opportunistic pathogen causing mucosal and systemic infections in humans. Systemic infections caused by this yeast have high mortality rates and are difficult to treat due to this yeast's intrinsic and frequently adapting antifungal resistance. To understand and treat C. glabrata infections, it is essential to investigate the molecular basis of C. glabrata virulence and resistance. We established an RNA interference (RNAi) system in C. glabrata by expressing the Dicer and Argonaute genes from Saccharomyces castellii (a budding yeast with natural RNAi). Our experiments with reporter genes and putative virulence genes showed that the introduction of RNAi resulted in 30 and 70% gene-knockdown for the construct-types antisense and hairpin, respectively. The resulting C. glabrata RNAi strain was used for the screening of a gene library for new virulence-related genes. Phenotypic profiling with a high-resolution quantification of growth identified genes involved in the maintenance of cell integrity, antifungal drugs, and ROS resistance. The genes identified by this approach are promising targets for the treatment of C. glabrata infections.

The effect of pedigree structure on detection of deletions and other null alleles.

Deletions and other null alleles for genetic markers can be detected as a special case of non-Mendelian inheritance, ie when a parent and a child appear to be homozygous for different alleles. The probability to detect a deletion for a fixed overall number of investigated individuals was calculated for biallelic and multiallelic markers with varying allele frequencies. To determine the effect of increasing the number of parents and grandparents, the probability for this event was derived for a parent and one child, a trio, a trio with one grandparent and a trio with two grandparents. The results for biallelic markers show that for a fixed total number of individuals, a sample of trios with two grandparents is always more efficient than the other family types, despite a lower total number of founder chromosomes in the sample. For multiallelic markers the outcome varies. The effect of adding additional children to a nuclear family was also investigated. For nuclear families, the optimal number of children is two or three, depending on the allele frequencies. It is shown that adding children is more efficient than adding grandparents.

Comprehensive evaluation of genetic variation in S100A7 suggests an association with the occurrence of allergic rhinitis.

BACKGROUND: S100A7 is a calcium-binding protein with chemotactic and antimicrobial properties. S100A7 protein levels are decreased in nasal lavage fluid from individuals with ongoing allergic rhinitis, suggesting a role for S100A7 in allergic airway inflammation. The aims of this study were to describe genetic variation in S100A7 and search for associations between this variation and allergic rhinitis. METHODS: Peripheral blood was collected from 184 atopic patients with a history of pollen-induced allergic rhinitis and 378 non-atopic individuals, all of Swedish origin. DNA was extracted and the S100A7 gene was resequenced in a subset of 47 randomly selected atopic individuals. Nine polymorphisms were genotyped in 184 atopic and 378 non-atopic individuals and subsequently investigated for associations with allergic rhinitis as well as skin prick test results. Haplotypes were estimated and compared in the two groups. RESULTS: Thirteen polymorphisms were identified in S100A7, of which 7 were previously undescribed. rs3014837 (G/C), which gives rise to an Asp --> Glu amino acid shift, had significantly increased minor allele frequency in atopic individuals. The major haplotype, containing the major allele at all sites, was more common in non-atopic individuals, while the haplotype containing the minor allele at rs3014837 was equally more common among the atopic individuals. Additionally, heterozygotes at this site had significantly higher scores in skin prick tests for 9 out of 11 tested allergens, compared to homozygotes. CONCLUSION: This is the first study describing genetic variation, associated with allergy, in S100A7. The results indicate that rs3014837 is linked to allergic rhinitis in our Swedish population and render S100A7 a strong candidate for further investigations regarding its role in allergic inflammation.

Temperature Sensing Is Distributed throughout the Regulatory Network that Controls FLC Epigenetic Silencing in Vernalization.

Many organisms need to respond to complex, noisy environmental signals for developmental decision making. Here, we dissect how Arabidopsis plants integrate widely fluctuating field temperatures over month-long timescales to progressively upregulate VERNALIZATION INSENSITIVE3 (VIN3) and silence FLOWERING LOCUS C (FLC), aligning flowering with spring. We develop a mathematical model for vernalization that operates on multiple timescales-long term (month), short term (day), and current (hour)-and is constrained by experimental data. Our analysis demonstrates that temperature sensing is not localized to specific nodes within the FLC network. Instead, temperature sensing is broadly distributed, with each thermosensory process responding to specific features of the plants' history of exposure to warm and cold. The model accurately predicts FLC silencing in new field data, allowing us to forecast FLC expression in changing climates. We suggest that distributed thermosensing may be a general property of thermoresponsive regulatory networks in complex natural environments.

Absence of warmth permits epigenetic memory of winter in Arabidopsis.

Plants integrate widely fluctuating temperatures to monitor seasonal progression. Here, we investigate the temperature signals in field conditions that result in vernalisation, the mechanism by which flowering is aligned with spring. We find that multiple, distinct aspects of the temperature profile contribute to vernalisation. In autumn, transient cold temperatures promote transcriptional shutdown of Arabidopsis FLOWERING LOCUS C (FLC), independently of factors conferring epigenetic memory. As winter continues, expression of VERNALIZATION INSENSITIVE3 (VIN3), a factor needed for epigenetic silencing, is upregulated by at least two independent thermosensory processes. One integrates long-term cold temperatures, while the other requires the absence of daily temperatures above 15 °C. The lack of spikes of high temperature, not just prolonged cold, is thus the major driver for vernalisation. Monitoring of peak daily temperature is an effective mechanism to judge seasonal progression, but is likely to have deleterious consequences for vernalisation as the climate becomes more variable.

AS3MT-mediated tolerance to arsenic evolved by multiple independent horizontal gene transfers from bacteria to eukaryotes.

Organisms have evolved the ability to tolerate toxic substances in their environments, often by producing metabolic enzymes that efficiently detoxify the toxicant. Inorganic arsenic is one of the most toxic and carcinogenic substances in the environment, but many organisms, including humans, metabolise inorganic arsenic to less toxic metabolites. This multistep process produces mono-, di-, and trimethylated arsenic metabolites, which the organism excretes. In humans, arsenite methyltransferase (AS3MT) appears to be the main metabolic enzyme that methylates arsenic. In this study, we examined the evolutionary origin of AS3MT and assessed the ability of different genotypes to produce methylated arsenic metabolites. Phylogenetic analysis suggests that multiple, independent horizontal gene transfers between different bacteria, and from bacteria to eukaryotes, increased tolerance to environmental arsenic during evolution. These findings are supported by the observation that genetic variation in AS3MT correlates with the capacity to methylate arsenic. Adaptation to arsenic thus serves as a model for how organisms evolve to survive under toxic conditions.

Multivariate analysis of chromosomal imbalances in breast cancer delineates cytogenetic pathways and reveals complex relationships among imbalances.

More than 550 breast adenocarcinomas with clonal chromosomal abnormalities have been reported. Although the aberration pattern is clearly nonrandom, no specific primary or secondary karyotypic abnormality has been identified, and furthermore the chronological order in which the aberrations appear during disease progression is not well known. The high degree of karyotypic complexity in epithelial tumors such as breast cancer is one reason why our understanding of the sequential order of cytogenetic evolution is unclear. To overcome some of these difficulties, we have used several statistical methods that allow identification and interpretation of karyotypic pathways. These methods were applied on 538 breast cancer karyotypes. The distribution of the number of imbalances/tumor showed a monomodal appearance, indicating that one single mode of karyotypic evolution is operating in this tumor type. We show that there exists a temporal order with respect to the appearance of chromosomal imbalances. The imbalances +1pq, 1q-, 3p-, and +7 appear earlier than expected from random events, and two cytogenetic pathways, one initiated by +1q and followed by 11q- and -22, the other initiated by either 3p- or 1q- and followed by 1p-, 3q-, and 6q-, can be discerned. We also show that +7 and +8q behave independently of the other imbalances and cannot, by simple means, be incorporated in the identified pathway scheme. Although the cytogenetic pathways are well separated at earlier stages, they later converge and include a common set of late imbalances.

Ovarian carcinoma develops through multiple modes of chromosomal evolution.

Ovarian carcinoma has the highest mortality of all of the gynecologic cancers. The chromosomal changes in this tumor type are highly complex, and the karyotypes typically show severe aneuploidy. Despite the abundance of cytogenetic information, with approximately 400 published karyotypes, very little is known about the mode of karyotypic evolution and the possible presence of cytogenetic pathways related to tumor development. In the present investigation we used 387 ovarian carcinoma karyotypes to identify the most frequent genomic imbalances. Tumor cases were then classified with respect to the presence or absence of these imbalances and statistically analyzed to assess the order of appearance of chromosomal imbalances, as well as possible karyotypic pathways and cytogenetic subtypes. We establish the temporal order by which the different imbalances occur and show that at least two cytogenetic pathways exist, one characterized by +7, +8q, and +12, and one by 6q- and 1q-. We show that ovarian carcinomas develop through at least three phases of karyotypic evolution. At the early stages, Phase I, the karyotypic evolution seems to proceed though step-wise acquisition of changes. The transition to Phase II showed signs of an increased chromosomal instability, most probably caused by extensive telomere crisis and the onset of breakage-fusion-bridge cycles. This process was linked to the presence of imbalances characteristic for the 6q-/1q- pathway. The transition to Phase III involved triploidization and was also linked to the presence of the 6q-/1q- pathway.

Dissecting karyotypic patterns in colorectal tumors: two distinct but overlapping pathways in the adenoma-carcinoma transition.

More than 500 colorectal tumors with clonal chromosomal abnormalities have been reported. Although the pattern of aberrations is nonrandom, no specific primary or secondary karyotypic abnormality has been identified. Also, the chronological order in which the aberrations appear during disease progression is not well known. One reason why our understanding of the cytogenetic evolution is unclear is the high degree of karyotypic complexity seen in these tumors. To overcome some of these difficulties we have previously used several statistical methods that allow identification and interpretation of karyotypic pathways as well as establishment of a temporal order of appearance of the imbalances. These methods were applied on 531 colorectal tumor karyotypes. By using a resampling strategy, 1p-, +7, 7q-, and +12p were identified as early events. Two major and two minor cytogenetic pathways were identified by means of principal component analysis. The two major pathways were initiated with 1p- and +7, and the minor pathways were initiated with +12p and 7q-. The +7/+12p tumors were found to be hyperdiploid, whereas those with 1p-/7q- were pseudodiploid. We also show that the adenoma-carcinoma transition in the 1p- pathway is strongly linked to karyoytypic evolution, whereas the +7 pathway is not, and that the cytogenetic pathways are separated at both early and late stages.