Development and maintenance of the neuronal cytoskeleton in aggregated cell cultures of fetal rat telencephalon and influence of elevated K+ concentrations.

Serum-free aggregating cell cultures of fetal rat telencephalon were examined by biochemical and immunocytochemical methods for their development-dependent expression of several cytoskeletal proteins, including the heavy- and medium-sized neurofilament subunits (H-NF and M-NF, respectively); brain spectrin; synapsin I; beta-tubulin; and the microtubule-associated proteins (MAPs) 1, 2, and 5 and tau protein. It was found that with time in culture the levels of most of these cytoskeletal proteins increased greatly, with the exceptions of the particular beta-tubulin form studied, which remained unchanged, and MAP 5, which greatly decreased. Among the neurofilament proteins, expression of M-NF preceded that of H-NF, with the latter being detectable only after approximately 3 weeks in culture. Furthermore, MAP 2 and tau protein showed a development-dependent change in expression from the juvenile toward the adult form. The comparison of these developmental changes in cytoskeletal protein levels with those observed in rat brain tissue revealed that protein expression in aggregate cultures is nearly identical to that in vivo during maturation of the neuronal cytoskeleton. Aggregate cultures deprived of glial cells, i.e., neuron-enriched cultures prepared by treating early cultures with the antimitotic drug cytosine arabinoside, exhibited pronounced deficits in M-NF, H-NF, MAP 2, MAP 1, synapsin I, and brain spectrin, with increased levels of a 145-kDa brain spectrin breakdown product. These adverse effects of glial cell deprivation could be reversed by the maintenance of neuron-enriched cultures at elevated concentrations of KCl (30 mM). This chronic treatment had to be started at an early developmental stage to be effective, a finding suggesting that sustained depolarization by KCl is able to enhance the developmental expression and maturation of the neuronal cytoskeleton.

Rapid Improvement and Maintenance of the Lipid Profile After Roux-en-Y Gastric Bypass (RYGBP)

Background: Dyslipidemia, a major component of the metabolic syndrome and an important cardiovascular risk factor, is one of the commonest comorbidity associated with morbid obesity. The aim of this paper is to show that RYGBP markedly improves dyslipidemia and that this improvement maintains over time. Patients and Methods: Prospectively updated databank for bariatric patients. Patients undergoing RYGBP have yearly blood tests during follow-up. The results for lipids at one to five years were compared with preoperative values. Results: The mean excess BMI loss after one and five years was 77,9 % and 72,3%respectively. After one year, there was a significant reduction of the mean total cholesterol, LDL-cholesterol, total cholesterol/HDL ratio and triglyceride values, which maintained up to five years, and an increase of the HDL fraction, which progressed until five years. The proportion of patients with abnormal values decreased from 24,3 to 6,2% for total cholesterol, from 45,1 to 11,7 %for HDL, from 53,3 to 21,9 for LDL, and from 40,5 to 10 % for triglycerides, with no significant change between three and five years, despite some weight regain. Conclusions: RYGBP rapidly improves all components of dyslipidemia, and thereby reduces the overall cardiovascular risk in operated patients.

The Complete Chromosomal Organization of the Reference Strain of the Leishmania Genome Project, L. major ;Friedlin'.

In recent years, analysis of the genomes of many organisms has received increasing international attention. The bulk of the effort to date has centred on the Human Genome Project and analysis of model organisms such as yeast, Drosophila and Caenorhabditis elegans. More recently, the revolution in genome sequencing and gene identification has begun to impact on infectious disease organisms. Initially, much of the effort was concentrated on prokaryotes, but small eukaryotic genomes, including the protozoan parasites Plasmodium, Toxoplasma and trypanosomatids (Leishmania, Trypanosoma brucei and T. cruzi), as well as some multicellular organisms, such as Brugia and Schistosoma, are benefiting from the technological advances of the genome era. These advances promise a radical new approach to the development of novel diagnostic tools, chemotherapeutic targets and vaccines for infectious disease organisms, as well as to the more detailed analysis of cell biology and function.Several networks or consortia linking laboratories around the world have been established to support these parasite genome projects[1] (for more information, see http://www.ebi.ac.uk/ parasites/paratable.html). Five of these networks were supported by an initiative launched in 1994 by the Specific Programme for Research and Tropical Diseases (TDR) of the WHO[2, 3, 4, 5, 6]. The Leishmania Genome Network (LGN) is one of these[3]. Its activities are reported at http://www.ebi.ac.uk/parasites/leish.html, and its current aim is to map and sequence the genome of Leishmania by the year 2002. All the mapping, hybridization and sequence data are also publicly available from LeishDB, an AceDB-based genome database (http://www.ebi.ac.uk/parasites/LGN/leissssoft.html).

The genome of the leaf-cutting ant Acromyrmex echinatior suggests key adaptations to advanced social life and fungus farming.

We present a high-quality (>100× depth) Illumina genome sequence of the leaf-cutting ant Acromyrmex echinatior, a model species for symbiosis and reproductive conflict studies. We compare this genome with three previously sequenced genomes of ants from different subfamilies and focus our analyses on aspects of the genome likely to be associated with known evolutionary changes. The first is the specialized fungal diet of A. echinatior, where we find gene loss in the ant's arginine synthesis pathway, loss of detoxification genes, and expansion of a group of peptidase proteins. One of these is a unique ant-derived contribution to the fecal fluid, which otherwise consists of "garden manuring" fungal enzymes that are unaffected by ant digestion. The second is multiple mating of queens and ejaculate competition, which may be associated with a greatly expanded nardilysin-like peptidase gene family. The third is sex determination, where we could identify only a single homolog of the feminizer gene. As other ants and the honeybee have duplications of this gene, we hypothesize that this may partly explain the frequent production of diploid male larvae in A. echinatior. The fourth is the evolution of eusociality, where we find a highly conserved ant-specific profile of neuropeptide genes that may be related to caste determination. These first analyses of the A. echinatior genome indicate that considerable genetic changes are likely to have accompanied the transition from hunter-gathering to agricultural food production 50 million years ago, and the transition from single to multiple queen mating 10 million years ago.

Genome sequencing of the plant pathogen Taphrina deformans, the causal agent of peach leaf curl.

Taphrina deformans is a fungus responsible for peach leaf curl, an important plant disease. It is phylogenetically assigned to the Taphrinomycotina subphylum, which includes the fission yeast and the mammalian pathogens of the genus Pneumocystis. We describe here the genome of T. deformans in the light of its dual plant-saprophytic/plant-parasitic lifestyle. The 13.3-Mb genome contains few identifiable repeated elements (ca. 1.5%) and a relatively high GC content (49.5%). A total of 5,735 protein-coding genes were identified, among which 83% share similarities with other fungi. Adaptation to the plant host seems reflected in the genome, since the genome carries genes involved in plant cell wall degradation (e.g., cellulases and cutinases), secondary metabolism, the hallmark glyoxylate cycle, detoxification, and sterol biosynthesis, as well as genes involved in the biosynthesis of plant hormones. Genes involved in lipid metabolism may play a role in its virulence. Several locus candidates for putative MAT cassettes and sex-related genes akin to those of Schizosaccharomyces pombe were identified. A mating-type-switching mechanism similar to that found in ascomycetous yeasts could be in effect. Taken together, the findings are consistent with the alternate saprophytic and parasitic-pathogenic lifestyles of T. deformans. IMPORTANCE: Peach leaf curl is an important plant disease which causes significant losses of fruit production. We report here the genome sequence of the causative agent of the disease, the fungus Taphrina deformans. The genome carries characteristic genes that are important for the plant infection process. These include (i) proteases that allow degradation of the plant tissues; (ii) secondary metabolites which are products favoring interaction of the fungus with the environment, including the host; (iii) hormones that are responsible for the symptom of severely distorted leaves on the host; and (iv) drug detoxification enzymes that confer resistance to fungicides. The availability of the genome allows the design of new drug targets as well as the elaboration of specific management strategies to fight the disease.

The complete genome sequence of the gram-positive bacterium Bacillus subtilis.

Bacillus subtilis is the best-characterized member of the Gram-positive bacteria. Its genome of 4,214,810 base pairs comprises 4,100 protein-coding genes. Of these protein-coding genes, 53% are represented once, while a quarter of the genome corresponds to several gene families that have been greatly expanded by gene duplication, the largest family containing 77 putative ATP-binding transport proteins. In addition, a large proportion of the genetic capacity is devoted to the utilization of a variety of carbon sources, including many plant-derived molecules. The identification of five signal peptidase genes, as well as several genes for components of the secretion apparatus, is important given the capacity of Bacillus strains to secrete large amounts of industrially important enzymes. Many of the genes are involved in the synthesis of secondary metabolites, including antibiotics, that are more typically associated with Streptomyces species. The genome contains at least ten prophages or remnants of prophages, indicating that bacteriophage infection has played an important evolutionary role in horizontal gene transfer, in particular in the propagation of bacterial pathogenesis.

Genome sequence of the Drosophila melanogaster male-killing Spiroplasma strain MSRO endosymbiont.

Spiroplasmas are helical and motile members of a cell wall-less eubacterial group called Mollicutes. Although all spiroplasmas are associated with arthropods, they exhibit great diversity with respect to both their modes of transmission and their effects on their hosts; ranging from horizontally transmitted pathogens and commensals to endosymbionts that are transmitted transovarially (i.e., from mother to offspring). Here we provide the first genome sequence, along with proteomic validation, of an endosymbiotic inherited Spiroplasma bacterium, the Spiroplasma poulsonii MSRO strain harbored by Drosophila melanogaster. Comparison of the genome content of S. poulsonii with that of horizontally transmitted spiroplasmas indicates that S. poulsonii has lost many metabolic pathways and transporters, demonstrating a high level of interdependence with its insect host. Consistent with genome analysis, experimental studies showed that S. poulsonii metabolizes glucose but not trehalose. Notably, trehalose is more abundant than glucose in Drosophila hemolymph, and the inability to metabolize trehalose may prevent S. poulsonii from overproliferating. Our study identifies putative virulence genes, notably, those for a chitinase, the H2O2-producing glycerol-3-phosphate oxidase, and enzymes involved in the synthesis of the eukaryote-toxic lipid cardiolipin. S. poulsonii also expresses on the cell membrane one functional adhesion-related protein and two divergent spiralin proteins that have been implicated in insect cell invasion in other spiroplasmas. These lipoproteins may be involved in the colonization of the Drosophila germ line, ensuring S. poulsonii vertical transmission. The S. poulsonii genome is a valuable resource to explore the mechanisms of male killing and symbiont-mediated protection, two cardinal features of many facultative endosymbionts. IMPORTANCE: Most insect species, including important disease vectors and crop pests, harbor vertically transmitted endosymbiotic bacteria. These endosymbionts play key roles in their hosts' fitness, including protecting them against natural enemies and manipulating their reproduction in ways that increase the frequency of symbiont infection. Little is known about the molecular mechanisms that underlie these processes. Here, we provide the first genome draft of a vertically transmitted male-killing Spiroplasma bacterium, the S. poulsonii MSRO strain harbored by D. melanogaster. Analysis of the S. poulsonii genome was complemented by proteomics and ex vivo metabolic experiments. Our results indicate that S. poulsonii has reduced metabolic capabilities and expresses divergent membrane lipoproteins and potential virulence factors that likely participate in Spiroplasma-host interactions. This work fills a gap in our knowledge of insect endosymbionts and provides tools with which to decipher the interaction between Spiroplasma bacteria and their well-characterized host D. melanogaster, which is emerging as a model of endosymbiosis.

Enzymic, phylogenetic, and structural characterization of the unusual papain-like protease domain of Plasmodium falciparum SERA5.

Serine repeat antigen 5 (SERA5) is an abundant antigen of the human malaria parasite Plasmodium falciparum and is the most strongly expressed member of the nine-gene SERA family. It appears to be essential for the maintenance of the erythrocytic cycle, unlike a number of other members of this family, and has been implicated in parasite egress and/or erythrocyte invasion. All SERA proteins possess a central domain that has homology to papain except in the case of SERA5 (and some other SERAs), where the active site cysteine has been replaced with a serine. To investigate if this domain retains catalytic activity, we expressed, purified, and refolded a recombinant form of the SERA5 enzyme domain. This protein possessed chymotrypsin-like proteolytic activity as it processed substrates downstream of aromatic residues, and its activity was reversed by the serine protease inhibitor 3,4-diisocoumarin. Although all Plasmodium SERA enzyme domain sequences share considerable homology, phylogenetic studies revealed two distinct clusters across the genus, separated according to whether they possess an active site serine or cysteine. All Plasmodia appear to have at least one member of each group. Consistent with separate biological roles for members of these two clusters, molecular modeling studies revealed that SERA5 and SERA6 enzyme domains have dramatically different surface properties, although both have a characteristic papain-like fold, catalytic cleft, and an appropriately positioned catalytic triad. This study provides impetus for the examination of SERA5 as a target for antimalarial drug design.

Analysis of the regulation of Nodal function by SPC-mediated cleavage during early mammalian development

RÉSUMÉ Après implantation dans l'utérus, le foetus de mammifère est composé de trois populations différentes de cellules: l'epiblast, l'ectoderme extraembryonnaire et l'endoderme viscéral. Pendant la gastrulation, les cellules de l'epiblast donnent naissance aux trois lignées germinales: l'ectoderme, le mésoderme et l'endodermes. Les lignées germinales produisent par la suite les différents tissus et organes du corps embryonnaire et adulte. Les cellules de l'ectoderme extraembryonnaire donnent par la suite le composant foetal du placenta qui est essentiel à la survie de l'embryon dans l'utérus. L'épiblast et l'ectoderme extraembryonnaire sont entourés par l'endoderme viscéral et forment une structure connue sous le nom de bouton embryonnaire. L'endoderme viscéral joue un rôle important dans l'embryogenèse car il comporte une sous-population de cellules appelées l'endoderme viscéral antérieur dont les signaux influencent l'épiblast adjacent et déterminent le futur axe antéro-postérieur de l'embryon. La protéine de signalisation Nodal de la famille des TGFß est essentielle dans l'épiblast pour spécifier le mésendoderme, l'endoderme viscéral antérieur, ainsi que pour maintenir les cellules souche de l'ectoderme extraembryonnaire. Ainsi, dans les embryons mutants pour Nodal, aucun axe antéro-postérieur n'est établi, les lignées germinales ne sont pas spécifiés et le placenta ne se développe pas. Au niveau moléculaire, comme pour les protéines de la famille des TGFß, Nodal est initialement synthétisée sous forme de précurseur avant d'être clivée de façon endoproteolytique par des protéanes sécrétées, les proprotéines convertases du type subtilisin (SPC), qui suppriment la partie inhibitrice N-terminale du pro peptide. Dans ce contexte, le projet de ma thèse a été d'analyser l'influence des SPC sur la fonction de Nodal en employant une combinaison d'approches génétiques et biochimiques. Premièrement, nous avons constaté que le clivage du précurseur par les protéases active Nodal, mais en même temps augmente son turn-over et diminue la portée de son action. Deuxièmement, dans l'embryon, il apparaît que Nodal est activé par l'action combinée de Furin et de PACE4, deux protéases sécrétées qui sont spécifiquement exprimées dans les cellules de l'ectoderme extraembryonnaire, donc adjacentes au domaine d'expression de Nodal. De manière similaire aux mutants de Nodal, les embryons mutants pour les deux protéases ne forment pas d'endoderme viscéral antérieur et ne gastrulent pas. Cependant, certains gènes cible de Nodal restent exprimés, suggérant que toutes les activités de Nodal ne sont pas dépendent du clivage par les SPCs. En effet, la génération et l'analyse de mutants portant un allèle knock-in qui code pour une forme mutante de Nodal résistante aux SPC, ont montré que ces mutants ont les caractères phénotypique des mutants de Nodal seulement de façon partielle. La formation de mésoderme est partiellement induite, et de façon remarquable, la forme de Nodal résistante aux SPC est capable d'agir à une distance de sa source, maintenant l'expression de ses propres protéases et d'autres gènes essentiels pour la spécification de l'ectoderme extraembryonnaire. Ensemble, ces résultats prouvent que par leur action directe les protéases extraembryonnaire modulent la signalisation de Nodal pendant le développement mammifère précoce. SUMMARY : Early after implantation in the uterus, the mammalian conceptus is composed of three different cell populations: the epiblast, the extraembryonic ectoderm and the visceral endoderm. During gastrulation, epiblast cells give rise to the three embryonic germ layers: the ectoderm, the mesoderm and the endoderm. These germ layers then generate the different tissues and organs of the embryonic and adult bodies. In parallel, extraembryonic ectoderm cells give rise to the fetal component of the placenta, which is essential for the survival of the embryo in the uterus. Both the epiblast and extraembryonic ectoderm are surrounded by the visceral endoderm to form a structure known as the egg cylinder. The visceral endoderm plays an important role as it harbours a subpopulation of cells called the anterior visceral endoderm, from which signals influence the adjacent epiblast and determine the future antero-posterior embryonic axis. The TGFß-related signalling protein Nodal is required within the epiblast to specify the mesoderm, the endoderm,the anterior visceral endoderm and is also essential to maintain stem cells in the extraembryonic ectoderm. Thus, in Nodal null conceptuses, no antero-posterior axis is established, the germ layers are not specified and the placenta does not develop. At the molecular level, Nodal, like related proteins of the TGFß family, is initially synthesized as a precursor and undergoes endoproteolytic cleavage by secreted proteases of the subtilisin-like proprotein convertases (SPC) to remove an inhibitory N-terminal pro peptide. In the embryo, Nodal is activated by the combined action of Furin and PACE4, two secreted SPCs that are specifically expressed in cells of the extraembryonic ectoderm, thus adjacent to the Nodal expression domain. Similar to Nodal null .embryos, mutant embryos lacking both these proteases fail to specify the anterior visceral endoderm and to undergo gastrulation. However, these mutants still express a subset of Nodal target genes, suggesting that part of Nodal activity is independent on cleavage by SPCs. Indeed, by generating and analyzing mutants with a knock-in allele that encodes an SPC-resistant mutant form of Nodal, I could show that they retain a subset of Nodal activities. Mesoderm formation is partially induced, but most remarkably, SPC-resistant Nodal form is able to act at a distance from its source, maintaining the expression of its proteases and of other genes essential for maintenance of the extraembryonic ectoderm.

Nucleotide sequence of the 5' noncoding region and part of the gag gene of mouse mammary tumor virus; identification of the 5' splicing site for subgenomic mRNAs.

We have determined the sequence of the first 1371 nucleotides at the 5' end of the genome of mouse mammary tumor virus using molecularly cloned proviral DNA of the GR virus strain. The most likely initiation codon used for the gag gene of mouse mammary tumor virus is the first one, located 312 nucleotides from the 5' end of the viral RNA. The 5' splicing site for the subgenomic mRNA's is located approximately 288 nucleotides downstream from the 5' end of the viral RNA. From the DNA sequence the amino acid sequence of the N-terminal half of the gag precursor protein, including p10 and p21, was deduced (353 amino acids).

Gene organization of the mating type regions in the ectomycorrhizal fungus Laccaria bicolor reveals distinct evolution between the two mating type loci

In natural conditions, basidiomycete ectomycorrhizal fungi such as Laccaria bicolor are typically in the dikaryotic state when forming symbioses with trees, meaning that two genetically different individuals have to fuse or 'mate'. Nevertheless, nothing is known about the molecular mechanisms of mating in these ecologically important fungi. Here, advantage was taken of the first sequenced genome of the ectomycorrhizal fungus, Laccaria bicolor, to determine the genes that govern the establishment of cell-type identity and orchestrate mating. The L. bicolor mating type loci were identified through genomic screening. The evolutionary history of the genomic regions that contained them was determined by genome-wide comparison of L. bicolor sequences with those of known tetrapolar and bipolar basidiomycete species, and by phylogenetic reconstruction of gene family history. It is shown that the genes of the two mating type loci, A and B, are conserved across the Agaricales, but they are contained in regions of the genome with different evolutionary histories. The A locus is in a region where the gene order is under strong selection across the Agaricales. By contrast, the B locus is in a region where the gene order is likely under a low selection pressure but where gene duplication, translocation and transposon insertion are frequent.

Mapping of the genes coding for the two major vaccinia virus core polypeptides.

We have mapped the genes coding for two major structural polypeptides of the vaccinia virus core by hybrid selection and transcriptional mapping. First, RNA was selected by hybridization to restriction fragments of the vaccinia virus genome, translated in vitro and the products were immunoprecipitated with antibodies against the two polypeptides. This approach allowed us to map the genes to the left hand end of the largest Hind III restriction fragment of 50 kilobase pairs. Second, transcriptional mapping of this region of the genome revealed the presence of the two expected RNAs. Both RNAs are transcribed from the leftward reading strand and the 5'-ends of the genes are separated by about 7.5 kilobase pairs of DNA. Thus, two genes encoding structural polypeptides with a similar location in the vaccinia virus particle are clustered at approximately 105 kilobase pairs from the left hand end of the 180 kilobase pair vaccinia virus genome.

A taxonomical re-evaluation of the Valais chromosome race of the common shrew Sorex araneus (Insectivora: Soricidae)

The common shrew Sorex araneus Linnaeus, 1758 is subject to intense chromosomal polymorphism. About 65 chromosome races are presently known. One of these chromosome races (the Valais race) is karyologically, morphologically, biochemically, and genetically clearly distinct from all other chromosome races of the species. Recent studies of hybrid zones between the Valais race and other chromosome races in the Swiss and French Alps add further strong evidence for the specific taxonomic status of the Valais race. Chromosomes and diagnostic protein markers reveal sharp frequency clines and strong heterozygote deficits. In one hybrid zone, the maintenance of the strong genetic differentiation of the hybridizing taxa was confirmed by a study with autosomal microsatellites indicating minimal gene flow. A microsatellite marker on the Y-chromosome showed complete absence of male mediated gene flow suggesting hybrid male sterility. To clarify the taxonomic status of this taxon, additional analyses were conducted. A morphometric analysis of the mandible indicated the Valais race is morphologically as distinct from neighbouring chromosome races of S. araneus as from other related Sorex species. In a phylogeny based on complete mitochondrial DNA cytochrome b gene sequences, the Valais race clearly appears as the sister taxon to all other races of S. araneus. Therefore, the chromosome race Valais of S. araneus herein is elevated to specific status and the name Sorex antinorii Bonaparte, 1840 is applied.

N-Glycans on secretory component: mediators of the interaction between secretory IgA and gram-positive commensals sustaining intestinal homeostasis.

Human beings live in symbiosis with billions of microorganisms colonizing mucosal surfaces. The understanding of the mechanisms underlying this fine-tuned intestinal balance has made significant processes during the last decades. We have recently demonstrated that the interaction of SIgA with Gram-positive bacteria is essentially based on Fab-independent, glycan-mediated recognition. Results obtained using mouse hybridoma- and colostrum-derived secretory IgA (SIgA) consistently show that N-glycans present on secretory component (SC) play a crucial role in the process. Natural coating may involve specific Gram-positive cell wall components, which may explain selective recognition at the molecular level. More widely, the existence of these complexes is involved in the modulation of intestinal epithelial cell (IEC) responses in vitro and the formation of intestinal biofilms. Thus, SIgA may act as one of the pillars in homeostatic maintenance of the microbiota in the gut, adding yet another facet to its multiple roles in the mucosal environment.

The complex SNP and CNV genetic architecture of the increased risk of congenital heart defects in Down syndrome.

Congenital heart defect (CHD) occurs in 40% of Down syndrome (DS) cases. While carrying three copies of chromosome 21 increases the risk for CHD, trisomy 21 itself is not sufficient to cause CHD. Thus, additional genetic variation and/or environmental factors could contribute to the CHD risk. Here we report genomic variations that in concert with trisomy 21, determine the risk for CHD in DS. This case-control GWAS includes 187 DS with CHD (AVSD = 69, ASD = 53, VSD = 65) as cases, and 151 DS without CHD as controls. Chromosome 21-specific association studies revealed rs2832616 and rs1943950 as CHD risk alleles (adjusted genotypic P-values <0.05). These signals were confirmed in a replication cohort of 92 DS-CHD cases and 80 DS-without CHD (nominal P-value 0.0022). Furthermore, CNV analyses using a customized chromosome 21 aCGH of 135K probes in 55 DS-AVSD and 53 DS-without CHD revealed three CNV regions associated with AVSD risk (FDR ≤ 0.05). Two of these regions that are located within the previously identified CHD region on chromosome 21 were further confirmed in a replication study of 49 DS-AVSD and 45 DS- without CHD (FDR ≤ 0.05). One of these CNVs maps near the RIPK4 gene, and the second includes the ZBTB21 (previously ZNF295) gene, highlighting the potential role of these genes in the pathogenesis of CHD in DS. We propose that the genetic architecture of the CHD risk of DS is complex and includes trisomy 21, and SNP and CNV variations in chromosome 21. In addition, a yet-unidentified genetic variation in the rest of the genome may contribute to this complex genetic architecture.