Computational Methods for Locating and Analyzing Conserved Gene Regulatory DNA Elements

This thesis presents methods for locating and analyzing cis-regulatory DNA elements involved with the regulation of gene expression in multicellular organisms. The regulation of gene expression is carried out by the combined effort of several transcription factor proteins collectively binding the DNA on the cis-regulatory elements. Only sparse knowledge of the 'genetic code' of these elements exists today. An automatic tool for discovery of putative cis-regulatory elements could help their experimental analysis, which would result in a more detailed view of the cis-regulatory element structure and function. We have developed a computational model for the evolutionary conservation of cis-regulatory elements. The elements are modeled as evolutionarily conserved clusters of sequence-specific transcription factor binding sites. We give an efficient dynamic programming algorithm that locates the putative cis-regulatory elements and scores them according to the conservation model. A notable proportion of the high-scoring DNA sequences show transcriptional enhancer activity in transgenic mouse embryos. The conservation model includes four parameters whose optimal values are estimated with simulated annealing. With good parameter values the model discriminates well between the DNA sequences with evolutionarily conserved cis-regulatory elements and the DNA sequences that have evolved neutrally. In further inquiry, the set of highest scoring putative cis-regulatory elements were found to be sensitive to small variations in the parameter values. The statistical significance of the putative cis-regulatory elements is estimated with the Two Component Extreme Value Distribution. The p-values grade the conservation of the cis-regulatory elements above the neutral expectation. The parameter values for the distribution are estimated by simulating the neutral DNA evolution. The conservation of the transcription factor binding sites can be used in the upstream analysis of regulatory interactions. This approach may provide mechanistic insight to the transcription level data from, e.g., microarray experiments. Here we give a method to predict shared transcriptional regulators for a set of co-expressed genes. The EEL (Enhancer Element Locator) software implements the method for locating putative cis-regulatory elements. The software facilitates both interactive use and distributed batch processing. We have used it to analyze the non-coding regions around all human genes with respect to the orthologous regions in various other species including mouse. The data from these genome-wide analyzes is stored in a relational database which is used in the publicly available web services for upstream analysis and visualization of the putative cis-regulatory elements in the human genome.

Rhizoctonia -Scots pine interactions: detection, impact on seedling performance and host defence gene response

Rhizoctonia spp. are ubiquitous soil inhabiting fungi that enter into pathogenic or symbiotic associations with plants. In general Rhizoctonia spp. are regarded as plant pathogenic fungi and many cause root rot and other plant diseases which results in considerable economic losses both in agriculture and forestry. Many Rhizoctonia strains enter into symbiotic mycorrhizal associations with orchids and some hypovirulent strains are promising biocontrol candidates in preventing host plant infection by pathogenic Rhizoctonia strains. This work focuses on uni- and binucleate Rhizoctonia (respectively UNR and BNR) strains belonging to the teleomorphic genus Ceratobasidium, but multinucleate Rhizoctonia (MNR) belonging to teleomorphic genus Thanatephorus and ectomycorrhizal fungal species, such as Suillus bovinus, were also included in DNA probe development work. Strain specific probes were developed to target rDNA ITS (internal transcribed spacer) sequences (ITS1, 5.8S and ITS2) and applied in Southern dot blot and liquid hybridization assays. Liquid hybridization was more sensitive and the size of the hybridized PCR products could be detected simultaneously, but the advantage in Southern hybridization was that sample DNA could be used without additional PCR amplification. The impacts of four Finnish BNR Ceratorhiza sp. strains 251, 266, 268 and 269 were investigated on Scot pine (Pinus sylvestris) seedling growth, and the infection biology and infection levels were microscopically examined following tryphan blue staining of infected roots. All BNR strains enhanced early seedling growth and affected the root architecture, while the infection levels remained low. The fungal infection was restricted to the outer cortical regions of long roots and typical monilioid cells detected with strain 268. The interactions of pathogenic UNR Ceratobasidium bicorne strain 1983-111/1N, and endophytic BNR Ceratorhiza sp. strain 268 were studied in single or dual inoculated Scots pine roots. The fungal infection levels and host defence-gene activity of nine transcripts [phenylalanine ammonia lyase (pal1), silbene synthase (STS), chalcone synthase (CHS), short-root specific peroxidase (Psyp1), antimicrobial peptide gene (Sp-AMP), rapidly elicited defence-related gene (PsACRE), germin-like protein (PsGER1), CuZn- superoxide dismutase (SOD), and dehydrin-like protein (dhy-like)] were measured from differentially treated and un-treated control roots by quantitative real time PCR (qRT-PCR). The infection level of pathogenic UNR was restricted in BNR- pre-inoculated Scots pine roots, while UNR was more competitive in simultaneous dual infection. The STS transcript was highly up-regulated in all treated roots, while CHS, pal1, and Psyp1 transcripts were more moderately activated. No significant activity of Sp-AMP, PsACRE, PsGER1, SOD, or dhy-like transcripts were detected compared to control roots. The integrated experiments presented, provide tools to assist in the future detection of these fungi in the environment and to understand the host infection biology and defence, and relationships between these interacting fungi in roots and soils. This study further confirms the complexity of the Rhizoctonia group both phylogenetically and in their infection biology and plant host specificity. The knowledge obtained could be applied in integrated forestry nursery management programmes.

Development of tissue culture and virus-induced gene silencing for Calendula officinalis

Calendula officinalis is grown widely as an ornamental plant across Europe. It belongs to the large. Asteraceae family. In this study, the aim was to explore the possibilities to use Calendula officinalis as a new model organism for flower development and secondary mechanism studies in Asteraceae. Tissue culture of Calendula officinalis was established using nine different cultivars. Murashige & Skoog (MS) medium with four different combinations of plant growth regulators were tested. Of all these combinations, the medium containing 1mg/l BAP, 0.1 mg/l IAA, and 1mg/l Zeatin achieved highest frequency of adventitious shoot regeneration from hypocotyl and cotyledon explants. Virus-induced gene silencing is a recent developed genetic tool for charactering the gene functions in plants, and extends the range of host plants that are not accessible for Agrobacterium transformation. Here, tobacco rattle virus (TRV)-based VIGS technique was tested in calendula (cv. Single Orange). We used TRV carrying Gerbera hybrid phytoene desaturase (PDS) gene fragment to induce PDS silencing in calendula. Vacuum infiltration and syringe infiltration methods both resulted in photo-bleaching phenotypes in leaves, bracts and petals. Loss-of-function phenotypes occurred on calendula 13 days post-infiltration. In conclusion, the data indicates that calendula explants can be regenerated through tissue culture which is a prerequisite for development of stable transformation methods. However, further optimization is still needed to improve the frequency. In addition, VIGS was applied to silence PDS marker gene expression indicating that this method has potential for gene functional studies in future.

Improving oncolytic adenoviral therapies for gastrointestinal cancers and tumor initiating cells

Although the treatment of most cancers has improved steadily, only few metastatic solid tumors can be cured. Despite responses, refractory clones often emerge and the disease becomes refractory to available treatment modalities. Furthermore, resistance factors are shared between different treatment regimens and therefore loss of response typically occurs rapidly, and there is a tendency for cross-resistance between agents. Therefore, new agents with novel mechanisms of action and lacking cross-resistance to currently available approaches are needed. Modified oncolytic adenoviruses, featuring cancer-celective cell lysis and spread, constitute an interesting drug platform towards the goals of tumor specificity and the implementation of potent multimodal treatment regimens. In this work, we demonstrate the applicability of capsid-modified, transcriptionally targeted oncolytic adenoviruses in targeting gastric, pancreatic and breast cancer. A variety of capsid modified adenoviruses were tested for transductional specificity first in gastric and pancreatic cancer cells and patient tissues and then in mice. Then, oncolytic viruses featuring the same capsid modifications were tested to confirm that successful transductional targeting translates into enhanced oncolytic potential. Capsid modified oncolytic viruses also prolonged the survival of tumor bearing orthotopic models of gastric and pancreatic cancer. Taken together, oncolytic adenoviral gene therapy could be a potent drug for gastric and pancreatic cancer, and its specificity, potency and safety can be modulated by means of capsid modification. We also characterized a new intraperitoneal virus delivery method in benefit for the persistence of gene delivery to intraperitoneal gastric and pancreatic cancer tumors. With a silica implant a steady and sustained virus release to the vicinity of the tumor improved the survival of the orthotopic tumor bearing mice. Furthermore, silica gel-based virus delivery lowered the toxicity mediating proimflammatory cytokine response and production of total and anti-adenovirus neutralizing antibodies (NAbs). On the other hand, silica shielded the virus against pre-excisting NAbs, resulting in a more favourable biodistribution in the preimmunized mice. The silica implant might therefore be of interest in treating intraperitoneally disseminated disease. Cancer stem cells are thought to be resistant to conventional cancer drugs and might play an important role in cancer relapse and the formation of metastasis. Therefore, we examined if transcriptionally modified oncolytic adenoviruses are able to kill these cells. Complete eradication of CD44+CD24-/low putative breast cancer stem cells was seen in vitro, and significant antitumor activity was detected in CD44+CD24-/low –derived tumor bearing mice. Thus, genetically engineered oncolytic adenoviruses have potential in destroying cancer initiating cells, which may have relevance for the elimination of cancer stem cells in humans.

Molecular Pathways of Disturbed Sleep and Depression: Studies on Adenosine and Gene Expression Patterns

Background: Adenosine is a potent sleep-promoting substance, and one of its targets is the basal forebrain. Fairly little is known about its mechanism of action in the basal forebrain and about the receptor subtype mediating its regulating effects on sleep homeostasis. Homeostatic deficiency might be one of the causes of the profoundly disturbed sleep pattern in major depressive disorder, which could explain the reduced amounts of delta-activity-rich stages 3 and 4. Since major depression has a relatively high heritability, and on the other hand adenosine regulates sleep homeostasis and might also be involved in mood modulation, adenosine-related genes should be considered for their possible contribution to a predisposition for depression and disturbed sleep in humans. Depression is a complex disorder likely involving the abnormal functioning of several genes. Novel target genes which could serve as the possible common substrates for depression and comorbid disturbed sleep should be identified. In this way specific brain areas related to sleep regulation should be studied by using animal model of depression which represents more homogenous phenotype as compared to humans. It is also important to study these brain areas during the development of depressive-like features to understand how early changes could facilitate pathophysiological changes in depression. Aims and methods: We aimed to find out whether, in the basal forebrain, adenosine induces recovery non-rapid eye movement (NREM) sleep after prolonged waking through the A1 or/and A2A receptor subtype. A1 and A2A receptor antagonists were perfused into the rat basal forebrain during 3 h of sleep deprivation, and the amount of NREM sleep and delta power during recovery NREM sleep were analyzed. We then explored whether polymorphisms in genes related to the metabolism, transport and signaling of adenosine could predispose to depression accompanied by signs of disturbed sleep. DNA from 1423 individuals representative of the Finnish population and including controls and cases with depression, depression accompanied by early morning awakenings and depression accompanied by fatigue, was used in the study to investigate the possible association between polymorphisms from adenosine-related genes and cases. Finally to find common molecular substrates of depression and disturbed sleep, gene expression changes were investigated in specific brain areas in the rat clomipramine model of depression. We focused on the basal forebrain of 3-week old clomipramine-treated rats which develop depressive-like symptoms later in adulthood and on the hypothalamus of adult female clomipramine-treated rats. Results: Blocking of the A1 receptor during sleep deprivation resulted in a reduction of the recovery NREM sleep amount and delta power, whereas A2A receptor antagonism had no effect. Polymorphisms in adenosine-related genes SLC29A3 (equilibrative nucleoside transporter type 3) in women and SLC28A1 (concentrative nucleoside transporter type 1) in men associated with depression alone as well as when accompanied by early morning awakenings and fatigue. In Study III the basal forebrain of postnatal rats treated with clomipramine displayed disturbances in gamma-aminobutyric acid (GABA) receptor type A signaling, in synaptic transmission and possible epigenetic changes. CREB1 was identified as a common transcription denominator which also mediates epigenetic regulation. In the hypothalamus the major changes included the expression of genes in GABA-A receptor pathway, K+ channel-related, glutamatergic and mitochondrial genes, as well as an overexpression of genes related to RNA and mRNA processing. Conclusions: Adenosine plays an important role in sleep homeostasis by promoting recovery NREM sleep via the A1 receptor subtype in the basal forebrain. Also adenosine levels might contribute to the risk of depression with disturbed sleep, since the genes encoding nucleoside transporters showed the strongest associations with depression alone and when accompanied by signs of disturbed sleep in both women and men. Sleep and mood abnormalities in major depressive disorder could be a consequence of multiple changes at the transcriptional level, GABA-A receptor signaling and synaptic transmission in sleep-related basal forebrain and the hypothalamus.

Vascular growth factors and progenitor cells in cardiac allograft arteriosclerosis

Heart transplantation is the only therapeutic modality for many end-stage heart diseases but poor long-term survival remains a challenging problem. This is mainly due to the development of cardiac allograft arteriosclerosis (TxCAD) that is an accelerated form of coronary artery disease. Both traditional cardiovascular and transplantation-related risk factors for TxCAD have been identified but options for therapy are limited. TxCAD involves dysfunction of cardiac allograft vascular cells. Activated endothelial cells (EC) regulate allograft inflammation and secrete smooth muscle cell (SMC) growth factors. In turn, SMC and their progenitors invade the intima of the injured vessels and occlude the affected coronary arteries. Different vascular growth factors have to be delicately regulated in normal vascular development. In the present study, experimental heterotopic transplantation models were used to study the role of angiogenic and pro-inflammatory vascular endothelial growth factor (VEGF), EC growth factor angiopoietin (Ang), and SMC mitogen platelet-derived growth factor (PDGF) in the development of TxCAD. Pharmacological and gene transfer approaches were used to target these growth factors and to assess their therapeutic potential. This study shows that alloimmune response in heart transplants upregulates VEGF expression, and induces allograft angiogenesis that involves donor-derived primitive EC. Intracoronary adenoviral VEGF gene transfer increased macrophage infiltration, intimal angiogenesis and TxCAD. VEGF inhibition with PTK787 decreased allograft inflammation and TxCAD, and simultaneous PDGF inhibition with imatinib further decreased TxCAD. Specific inhibition of two VEGF-receptors (VEGFR) decreased allograft inflammation and TxCAD, and VEGFR-2 inhibition normalized the density of primitive and mature capillaries in the allografts. Adenovirus-mediated transient Ang1 expression in the allograft had anti-inflammatory and anti-arteriosclerotic effects. Adeno-associated virus (AAV)-mediated prolonged Ang1 or Ang2 expression had similar anti-inflammatory effects. However, AAV-Ang1 activated allograft SMC whereas AAV-Ang2 had no effects on SMC activation and decreased the development of TxCAD. These studies indicate an interplay of inflammation, angiogenesis and arteriosclerosis in cardiac allografts, and show that vascular growth factors are important regulators in the process. Also, VEGF inhibition, PDGF inhibition and angiopoietin therapy with clinically-relevant pharmacological agents or novel gene therapy approaches may counteract vascular dysfunction in cardiac allografts, and have beneficial effects on the survival of heart transplant patients in the future.

Associations of gene polymorphisms and nutrition with calcium homeostasis and bone mineral density : Studies on skeletal nutrigenetics

Heredity explains a major part of the variation in calcium homeostasis and bone strength, and the susceptibility to osteoporosis is polygenetically regulated. Bone phenotype results from the interplay between lifestyle and genes, and several nutritional factors modulate bone health throughout life. Thus, nutrigenetics examining the genetic variation in nutrient intake and homeostatic control is an important research area in the etiology of osteoporosis. Despite continuing progress in the search for candidate genes for osteoporosis, the results thus far have been inconclusive. The main objective of this thesis was to investigate the associations of lactase, vitamin D receptor (VDR), calcium sensing receptor (CaSR) and parathyroid hormone (PTH) gene polymorphisms and lifestyle factors and their interactions with bone health in Finns at varying stages of the skeletal life span. Markers of calcium homeostasis and bone remodelling were measured from blood and urine samples. Bone strength was measured at peripheral and central bone sites. Lifestyle factors were assessed with questionnaires and interviews. Genetic lactase non-persistence (the C/C-13910 genotype) was associated with lower consumption of milk from childhood, predisposing females in particular to inadequate calcium intake. Consumption of low-lactose milk and milk products was shown to decrease the risk for inadequate calcium intake. In young adulthood, bone loss was more common in males than in females. Males with the lactase C/C-13910 genotype may be more susceptible to bone loss than males with the other lactase genotypes, although calcium intake predicts changes in bone mass more than the lactase genotype. The BsmI and FokI polymorphisms of the VDR gene were associated with bone mass in growing adolescents, but the associations weakened with age. In young adults, the A986S polymorphism of the calcium sensing receptor gene was associated with serum ionized calcium concentrations, and the BstBI polymorphism of the parathyroid gene was related to bone strength. The FokI polymorphism and sodium intake showed an interaction effect on urinary calcium excretion. A novel gene-gene interaction between the VDR FokI and PTH BstBI gene polymorphisms was found in the regulation of PTH secretion and urinary calcium excretion. Further research should be carried out with more number of Finns at varying stages of the skeletal life span and more detailed measurements of bone strength. Research should concern mechanisms by which genetic variants affect calcium homeostasis and bone strength, and the role of diet-gene and gene-gene interactions in the pathogenesis of osteoporosis.

Stress and developmental responses in Arabidopsis thaliana: regulation through the transcription factor-interacting protein RCD1

Plants constantly face adverse environmental conditions, such as drought or extreme temperatures that threaten their survival. They demonstrate astonishing metabolic flexibility in overcoming these challenges and one of the key responses to stresses is changes in gene expression leading to alterations in cellular functions. This is brought about by an intricate network of transcription factors and associated regulatory proteins. Protein-protein interactions and post-translational modifications are important steps in this control system along with carefully regulated degradation of signaling proteins. This work concentrates on the RADICAL-INDUCED CELL DEATH1 (RCD1) protein which is an important regulator of abiotic stress-related and developmental responses in Arabidopsis thaliana. Plants lacking this protein function display pleiotropic phenotypes including sensitivity to apoplastic reactive oxygen species (ROS) and salt, ultraviolet B (UV-B) and paraquat tolerance, early flowering and senescence. Additionally, the mutant plants overproduce nitric oxide, have alterations in their responses to several plant hormones and perturbations in gene expression profiles. The RCD1 gene is transcriptionally unresponsive to environmental signals and the regulation of the protein function is likely to happen post-translationally. RCD1 belongs to a small protein family and, together with its closest homolog SRO1, contains three distinguishable domains: In the N-terminus, there is a WWE domain followed by a poly(ADP-ribose) polymerase-like domain which, despite sequence conservation, does not seem to be functional. The C-terminus of RCD1 contains a novel domain called RST. It is present in RCD1-like proteins throughout the plant kingdom and is able to mediate physical interactions with multiple transcription factors. In conclusion, RCD1 is a key point of signal integration that links ROS-mediated cues to transcriptional regulation by yet unidentified means, which are likely to include post-translational mechanisms. The identification of RCD1-interacting transcription factors, most of whose functions are still unknown, opens new avenues for studies on plant stress as well as developmental responses.

Prognostic molecular factors and algorithms in diffuse large B-cell lymphoma

Diffuse large B-cell lymphoma (DLBCL) is the most common of the non-Hodgkin lymphomas. As DLBCL is characterized by heterogeneous clinical and biological features, its prognosis varies. To date, the International Prognostic Index has been the strongest predictor of outcome for DLBCL patients. However, no biological characters of the disease are taken into account. Gene expression profiling studies have identified two major cell-of-origin phenotypes in DLBCL with different prognoses, the favourable germinal centre B-cell-like (GCB) and the unfavourable activated B-cell-like (ABC) phenotypes. However, results of the prognostic impact of the immunohistochemically defined GCB and non-GCB distinction are controversial. Furthermore, since the addition of the CD20 antibody rituximab to chemotherapy has been established as the standard treatment of DLBCL, all molecular markers need to be evaluated in the post-rituximab era. In this study, we aimed to evaluate the predictive value of immunohistochemically defined cell-of-origin classification in DLBCL patients. The GCB and non-GCB phenotypes were defined according to the Hans algorithm (CD10, BCL6 and MUM1/IRF4) among 90 immunochemotherapy- and 104 chemotherapy-treated DLBCL patients. In the chemotherapy group, we observed a significant difference in survival between GCB and non-GCB patients, with a good and a poor prognosis, respectively. However, in the rituximab group, no prognostic value of the GCB phenotype was observed. Likewise, among 29 high-risk de novo DLBCL patients receiving high-dose chemotherapy and autologous stem cell transplantation, the survival of non-GCB patients was improved, but no difference in outcome was seen between GCB and non-GCB subgroups. Since the results suggested that the Hans algorithm was not applicable in immunochemotherapy-treated DLBCL patients, we aimed to further focus on algorithms based on ABC markers. We examined the modified activated B-cell-like algorithm based (MUM1/IRF4 and FOXP1), as well as a previously reported Muris algorithm (BCL2, CD10 and MUM1/IRF4) among 88 DLBCL patients uniformly treated with immunochemotherapy. Both algorithms distinguished the unfavourable ABC-like subgroup with a significantly inferior failure-free survival relative to the GCB-like DLBCL patients. Similarly, the results of the individual predictive molecular markers transcription factor FOXP1 and anti-apoptotic protein BCL2 have been inconsistent and should be assessed in immunochemotherapy-treated DLBCL patients. The markers were evaluated in a cohort of 117 patients treated with rituximab and chemotherapy. FOXP1 expression could not distinguish between patients, with favourable and those with poor outcomes. In contrast, BCL2-negative DLBCL patients had significantly superior survival relative to BCL2-positive patients. Our results indicate that the immunohistochemically defined cell-of-origin classification in DLBCL has a prognostic impact in the immunochemotherapy era, when the identifying algorithms are based on ABC-associated markers. We also propose that BCL2 negativity is predictive of a favourable outcome. Further investigational efforts are, however, warranted to identify the molecular features of DLBCL that could enable individualized cancer therapy in routine patient care.

Characterization of LRRTM and NGR gene families: Expression and functions

Some leucine-rich repeat (LRR) -containing membrane proteins are known regulators of neuronal growth and synapse formation. In this work I characterize two gene families encoding neuronal LRR membrane proteins, namely the LRRTM (leucine-rich repeat, transmembrane neuronal) and NGR (Nogo-66 receptor) families. I studied LRRTM and NGR family member's mRNA tissue distribution by RT-PCR and by in situ hybridization. Subcellular localization of LRRTM1 protein was studied in neurons and in non-neuronal cells. I discovered that LRRTM and NGR family mRNAs are predominantly expressed in the nervous system, and that each gene possesses a specific expression pattern. I also established that LRRTM and NGR family mRNAs are expressed by neurons, and not by glial cells. Within neurons, LRRTM1 protein is not transported to the plasma membrane; rather it localizes to endoplasmic reticulum. Nogo-A (RTN4), MAG, and OMgp are myelin-associated proteins that bind to NgR1 to limit axonal regeneration after central nervous system injury. To better understand the functions of NgR2 and NgR3, and to explore the possible redundancy in the signaling of myelin inhibitors of neurite growth, I mapped the interactions between NgR family and the known and candidate NgR1 ligands. I identified high-affinity interactions between RTN2-66, RTN3-66 and NgR1. I also demonstrate that Rtn3 mRNA is expressed in the same glial cell population of mouse spinal cord white matter as Nogo-A mRNA, and thus it could have a role in myelin inhibition of axonal growth. To understand how NgR1 interacts with multiple structurally divergent ligands, I aimed first to map in more detail the nature of Nogo-A:NgR1 interactions, and then to systematically map the binding sites of multiple myelin ligands in NgR1 by using a library of NgR1 expression constructs encoding proteins with one or multiple surface residues mutated to alanine. My analysis of the Nogo-A:NgR1 -interactions revealed a novel interaction site between the proteins, suggesting a trivalent Nogo-A:NgR1-interaction. Our analysis also defined a central binding region on the concave side of NgR1's LRR domain that is required for the binding of all known ligands, and a surrounding region critical for binding MAG and OMgp. To better understand the biological role of LRRTMs, I generated Lrrtm1 and Lrrtm3 knock out mice. I show here that reporter genes expressed from the targeted loci can be used for maping the neuronal connections of Lrrtm1 and Lrrtm3 expressing neurons in finer detail. With regard to LRRTM1's role in humans, we found a strong association between a 70 kb-spanning haplotype in the proposed promoter region of LRRTM1 gene and two possibly related phenotypes: left-handedness and schizophrenia. Interestingly, the responsible haplotype was linked to phenotypic variability only when paternally inherited. In summary, I identified two families of neuronal receptor-like proteins, and mapped their expression and certain protein-protein interactions. The identification of a central binding region in NgR1 shared by multiple ligands may facilitate the design and development of small molecule therapeutics blocking binding of all NgR1 ligands. Additionally, the genetic association data suggests that allelic variation upstream of LRRTM1 may play a role in the development of left-right brain asymmetry in humans. Lrrtm1 and Lrrtm3 knock out mice developed as a part of this study will likely be useful for schizophrenia and Alzheimer s disease research.