Genomic approach to organ-specific growth control

Growth is a fundamental aspect of life cycle of all organisms. Body size varies highly in most animal groups, such as mammals. Moreover, growth of a multicellular organism is not uniform enlargement of size, but different body parts and organs grow to their characteristic sizes at different times. Currently very little is known about the molecular mechanisms governing this organ-specific growth. The genome sequencing projects have provided complete genomic DNA sequences of several species over the past decade. The amount of genomic sequence information, including sequence variants within species, is constantly increasing. Based on the universal genetic code, we can make sense of this sequence information as far as it codes proteins. However, less is known about the molecular mechanisms that control expression of genes, and about the variations in gene expression that underlie many pathological states in humans. This is caused in part by lack of information about the second genetic code that consists of the binding specificities of transcription factors and the combinatorial code by which transcription factor binding sites are assembled to form tissue-specific and/or ligand-regulated enhancer elements. This thesis presents a high-throughput assay for identification of transcription factor binding specificities, which were then used to measure the DNA binding profiles of transcription factors involved in growth control. We developed ‘enhancer element locator’, a computational tool, which can be used to predict functional enhancer elements. A genome-wide prediction of human and mouse enhancer elements generated a large database of enhancer elements. This database can be used to identify target genes of signaling pathways, and to predict activated transcription factors based on changes in gene expression. Predictions validated in transgenic mouse embryos revealed the presence of multiple tissue-specific enhancers in mouse c- and N-Myc genes, which has implications to organ specific growth control and tumor type specificity of oncogenes. Furthermore, we were able to locate a variation in a single nucleotide, which carries a susceptibility to colorectal cancer, to an enhancer element and propose a mechanism by which this SNP might be involved in generation of colorectal cancer.

Mapping of assembled epitopic regions of human chorionic gonadotropin reveals proximity of CTP alpha to the determinant loop beta 93-100

Three overlapping assembled epitopes of beta hCG have been mapped using MAb probes and a single step solid phase radioimmunoassay. These epitopes have been shown to be at receptor binding region comprising of the loop region beta Cys93-Cys100. Importance of disulphide bonds in maintaining integrity of these epitopes is assessed. Two MAbs (INN 58 and INN 22) interact with the beta region as well as the alpha C-terminal peptide, while the other MAb INN 24 interacts with only the beta region. Cross-reactivity pattern with beta hCG and hLH as web as the reported crystal structure of hCG substantiates the epitope identification. The results demonstrate utility of MAbs as probes in investigations on three-dimensional structure of gonadatropins.

Microarrays in Lung Cancer Research: From Comparative Analyses to Verified Findings

Keuhkosyöpä on yleisimpiä syöpätauteja. Se jaetaan kahteen päätyyppiin: pienisoluiseen ja ei-pienisoluiseen keuhkosyöpään. Ei-pienisoluinen keuhkosyöpä jaetaan lisäksi alatyyppeihin, joista suurimmat ovat levyepiteeli-, adeno- ja suurisoluinen karsinooma. Keuhkosyövän tärkein riskitekijä on tupakointi, mutta muutkin työ- ja elinympäristön altisteet, kuten asbesti, voivat johtaa syöpään. Väitöstyössä tutkittiin kahdenlaisten keuhkosyöpäryhmien erityispiirteitä. Työssä kartoitettiin, onko löydettävissä muutoksia, jotka erottavat asbestikeuhkosyövät muista syövistä sekä luuytimeen varhaisessa vaiheessa leviävät keuhkosyövät leviämättömistä syövistä. Tutkimusten ensimmäisessä vaiheessa käytettiin mikrosirupohjaisia menetelmiä, jotka mahdollistavat jopa kaikkien geenien tarkastelun yhden kokeen avulla. Vertailevien mikrosirututkimusten avulla on mahdollista paikantaa geenejä tai kromosomialueita, joiden muutokset erottelevat ryhmät toisistaan. Asbestiin liittyvissä tutkimuksissa paikannettiin kuusi kromosomialuetta, joissa geenien kopiolukumäärän sekä ilmenemistason muutokset erottelivat potilaat altistushistorian mukaan. Riippumattomilla laboratoriomenetelmillä tehtyjen jatkoanalyysien avulla pystyttiin varmistamaan, että 19p-alueen häviämä oli yhteydessä asbestialtistukseen. Työssä osoitettiin myös, että 19p-alueen muutoksia voidaan indusoida altistamalla soluja asbestille in vitro. Tutkimuksessa saatiin lisäksi viitteitä asbestispesifisistä muutoksista signaalinvälitysreiteissä, sillä yhdessä toimivien geenien ilmentymisessä havaittiin eroja asbestille altistuneiden ja altistumattomien välillä. Vertailemalla luuytimeen syövän aikaisessa vaiheessa levinneiden ja leviämättömien keuhkoadenokarsinoomien muutosprofiileita toisiinsa, paikannettiin viisi aluetta, joilla geenien kopiolukumäärä- sekä ilmenemistason muutokset erottelivat ryhmät toisistaan. Jatkoanalyyseissä havaittiin, että 4q-alueen häviämää esiintyi adenokarsinoomien lisäksi levyepiteelikarsinoomiin, jotka olivat levinneet luuytimeen. Myös keuhkosyöpien aivometastaaseissa alue oli toistuvasti hävinnyt. Väitöstyön tutkimukset osoittavat, että vertailevien mikrosiruanalyysien avulla saadaan tietoa syöpäryhmien erityispiirteistä. Työssä saadut tulokset osoittavat, että 19p-alueen muutokset ovat tyypillisiä asbestikeuhkosyöville ja 4q-alueen muutokset luuytimeen aikaisessa vaiheessa leviäville keuhkosyöville.

GPRA and the asthma locus on chromosome 7p14-p15

The basis of this work was the identification of a genomic region on chromosome 7p14-p15 that strongly associated with asthma and high serum total immunoglobulin E in a Finnish founder population from Kainuu. Using a hierarchical genotyping approach the linkage region was narrowed down until an evolutionary collectively inherited 133-kb haplotype block was discovered. The results were confirmed in two independent data sets: Asthma families from Quebec and allergy families from North-Karelia. In all the three cohorts studied, single nucleotide polymorphisms tagging seven common gene variants (haplotypes) were identified. Over half of the asthma patients carried three evolutionary closely related susceptibility haplotypes as opposed to approximately one third of the healthy controls. The risk effects of the gene variants varied from 1.4 to 2.5. In the disease-associated region, there was one protein-coding gene named GPRA (G Protein-coupled Receptor for Asthma susceptibility also known as NPSR1) which displayed extensive alternative splicing. Only the two isoforms with distinct intracellular tail sequences, GPRA-A and -B, encoded a full-length G protein-coupled receptor with seven transmembrane regions. Using various techniques, we showed that GPRA is expressed in multiple mucosal surfaces including epithelial cells throughout the respiratory tract. GPRA-A has additional expression in respiratory smooth muscle cells. However, in bronchial biopsies with unknown haplotypes, GPRA-B was upregulated in airways of all patient samples in contrast to the lack of expression in controls. Further support for GPRA as a common mediator of inflammation was obtained from a mouse model of ovalbumin-induced inflammation, where metacholine-induced airway hyperresponsiveness correlated with elevated GPRA mRNA levels in the lung and increased GPRA immunostaining in pulmonary macrophages. A novel GPRA agonist, Neuropeptide S (NPS), stimulated phagocytosis of Esterichia coli bacteria in a mouse macrophage cell line indicating a role for GPRA in the removal of inhaled allergens. The suggested GPRA functions prompted us to study, whether GPRA haplotypes associate with respiratory distress syndrome (RDS) and bronchopulmonary dysplasia (BPD) in infants sharing clinical symptoms with asthma. According to the results, near-term RDS and asthma may also share the same susceptibility and protective GPRA haplotypes. As in asthma, GPRA-B isoform expression was induced in bronchial smooth muscle cells in RDS and BPD suggesting a role for GPRA in bronchial hyperresponsiveness. In conclusion, the results of the present study suggest that the dysregulation of the GPRA/NPS pathway may not only be limited to the individuals carrying the risk variants of the gene but is also involved in the regulation of immune functions of asthma.

Human parvovirus B19: tissue persistence and prevalence of prototypic and new variants

Human parvovirus B19 is a minute ssDNA virus causing a wide variety of diseases, including erythema infectiosum, arthropathy, anemias, and fetal death. After primary infection, genomic DNA of B19 has been shown to persist in solid tissues of not only symptomatic but also of constitutionally healthy, immunocompetent individuals. In this thesis, the viral DNA was shown to persist as an apparently intact molecule of full length, and without persistence-specific mutations. Thus, although the mere presence of B19 DNA in tissue can not be used as a diagnostic criterion, a possible role in the pathogenesis of diseases e.g. through mRNA or protein production can not be excluded. The molecular mechanism, the host-cell type and the possible clinical significance of B19 DNA tissue persistence are yet to be elucidated. In the beginning of this work, the B19 genomic sequence was considered highly conserved. However, new variants were found: V9 was detected in 1998 in France, in serum of a child with aplastic crisis. This variant differed from the prototypic B19 sequences by ~10 %. In 2002 we found, persisting in skin of constitutionally healthy humans, DNA of another novel B19 variant, LaLi. Genetically this variant differed from both the prototypic sequences and the variant V9 also by ~10%. Simultaneously, B19 isolates with DNA sequences similar to LaLi were introduced by two other groups, in the USA and France. Based on phylogeny, a classification scheme based on three genotypes (B19 types 1-3) was proposed. Although the B19 virus is mainly transmitted via the respiratory route, blood and plasma-derived products contaminated with high levels of B19 DNA have also been shown to be infectious. The European Pharmacopoeia stipulates that, in Europe, from the beginning of 2004, plasma pools for manufacture must contain less than 104 IU/ml of B19 DNA. Quantitative PCR screening is therefore a prerequisite for restriction of the B19 DNA load and obtaining of safe plasma products. Due to the DNA sequence variation among the three B19 genotypes, however, B19 PCR methods might fail to detect the new variants. We therefore examined the suitability of the two commercially available quantitative B19 PCR tests, LightCycler-Parvovirus B19 quantification kit (Roche Diagnostics) and RealArt Parvo B19 LC PCR (Artus), for detection, quantification and differentiation of the three B19 types known, including B19 types 2 and 3. The former method was highly sensitive for detection of the B19 prototype but was not suitable for detection of types 2 and 3. The latter method detected and differentiated all three B19 virus types. However, one of the two type-3 strains was detected at a lower sensitivity. Then, we assessed the prevalence of the three B19 virus types among Finnish blood donors, by screening pooled plasma samples derived from >140 000 blood-donor units: none of the pools contained detectable levels of B19 virus types 2 or 3. According to the results of other groups, B19 type 2 was absent also among Danish blood-donors, and extremely rare among symptomatic European patients. B19 type 3 has been encountered endemically in Ghana and (apparently) in Brazil, and sporadical cases have been detected in France and the UK. We next examined the biological characteristics of these virus types. The p6 promoter regions of virus types 1-3 were cloned in front of a reporter gene, the constructs were transfected into different cell lines, and the promoter activities were measured. As a result, we found that the activities of the three p6 promoters, although differing in sequence by >20%, were of equal strength, and most active in B19-permissive cells. Furthermore, the infectivity of the three B19 types was examined in two B19-permissive cell lines. RT-PCR revealed synthesis of spliced B19 mRNAs, and immunofluorescence verified the production of NS1 and VP proteins in the infected cells. These experiments suggested similar host-cell tropism and showed that the three virus types are strains of the same species, i.e. human parvovirus B19. Last but not least, the sera from subjects infected in the past either with B19 type 1 or type 2 (as evidenced by tissue persistence of the respective DNAs), revealed in VP1/2- and VP2-EIAs a 100 % cross-reactivity between virus types 1 and 2. These results, together with similar studies by others, indicate that the three B19 genotypes constitute a single serotype.

Molecular Alterations in Asbestos-Related Lung Cancer

Background: Asbestos is a well known cancer-causing mineral fibre, which has a synergistic effect on lung cancer risk in combination with tobacco smoking. Several in vitro and in vivo experiments have demonstrated that asbestos can evoke chromosomal damage and cause alterations as well as gene expression changes. Lung tumours, in general, have very complex karyotypes with several recurrently gained and lost chromosomal regions and this has made it difficult to identify specific molecular changes related primarily to asbestos exposure. The main aim of these studies has been to characterize asbestos-related lung cancer at a molecular level. Methods: Samples from asbestos-exposed and non-exposed lung cancer patients were studied using array comparative genomic hybridization (aCGH) and fluorescent in situ hybridization (FISH) to detect copy number alterations (CNA) as well as microsatellite analysis to detect allelic imbalance (AI). In addition, asbestos-exposed cell lines were studied using gene expression microarrays. Results: Eighteen chromosomal regions showing differential copy number in the lung tumours of asbestos-exposed patients compared to those of non-exposed patients were identified. The most significant differences were detected at 2p21-p16.3, 5q35.3, 9q33.3-q34.11, 9q34.13-q34.3, 11p15.5, 14q11.2 and 19p13.1-p13.3 (p<0.005). The alterations at 2p and 9q were validated and characterized in detail using AI and FISH analysis in a larger study population. Furthermore, in vitro studies were performed to examine the early gene expression changes induced by asbestos in three different lung cell lines. The results revealed specific asbestos-associated gene expression profiles and biological processes as well as chromosomal regions enriched with genes believed to contribute to the common asbestos-related responses in the cell lines. Interestingly, the most significant region enriched with asbestos-response genes was identified at 2p22, close to the previously identified region showing asbestos-related CNA in lung tumours. Additionally, in this thesis, the dysregulated biological processes (Gene Ontology terms) detected in the cell line experiment were compared to dysregulated processes identified in patient samples in a later study (Ruosaari et al., 2008a). Commonly affected processes such as those related to protein ubiquitination, ion transport and surprisingly sensory perception of smell were identified. Conclusions: The identification of specific CNA and dysregulated biological processes shed some light on the underlying genes acting as mediators in asbestos-related lung carcinogenesis. It is postulated that the combination of several asbestos-specific molecular alterations could be used to develop a diagnostic method for the identification of asbestos-related lung cancer.

Association of FOXE1 Polyalanine Repeat Region with Papillary Thyroid Cancer

CONTEXT: Polyalanine tract variations in transcription factors have been identified for a wide spectrum of developmental disorders. The thyroid transcription factor forkhead factor E1 (FOXE1) contains a polymorphic polyalanine tract with 12-22 alanines. Single-nucleotide polymorphisms (SNP) close to this locus are associated with papillary thyroid cancer (PTC), and a strong linkage disequilibrium block extends across this region. OBJECTIVE: The objective of the study was to assess whether the FOXE1 polyalanine repeat region was associated with PTC and to assess the effect of polyalanine repeat region variants on protein expression, DNA binding, and transcriptional function on FOXE1-responsive promoters. DESIGN: This was a case-control study. SETTING: The study was conducted at a tertiary referral hospital. PATIENTS AND METHODS: The FOXE1 polyalanine repeat region and tag SNP were genotyped in 70 PTC, with a replication in a further 92 PTC, and compared with genotypes in 5767 healthy controls (including 5667 samples from the Wellcome Trust Case Control Consortium). In vitro studies were performed to examine the protein expression, DNA binding, and transcriptional function for FOXE1 variants of different polyalanine tract lengths. RESULTS: All the genotyped SNP were in tight linkage disequilibrium, including the FOXE1 polyalanine repeat region. We confirmed the strong association of rs1867277 with PTC (overall P = 1 × 10(-7), odds ratio 1.84, confidence interval 1.31-2.57). rs1867277 was in tight linkage disequilibrium with the FOXE1 polyalanine repeat region (r(2) = 0.95). FOXE1(16Ala) was associated with PTC with an odds ratio of 2.23 (confidence interval 1.42-3.50; P = 0.0005). Functional studies in vitro showed that FOXE1(16Ala) was transcriptionally impaired compared with FOXE1(14Ala), which was not due to differences in protein expression or DNA binding. CONCLUSIONS: We have confirmed the previous association of FOXE1 with PTC. Our data suggest that the coding polyalanine expansion in FOXE1 may be responsible for the observed association between FOXE1 and PTC.

Relative stability of DNA as a generic criterion for promoter prediction: whole genome annotation of microbial genomes with varying nucleotide base composition

The rapid increase in genome sequence information has necessitated the annotation of their functional elements, particularly those occurring in the non-coding regions, in the genomic context. Promoter region is the key regulatory region, which enables the gene to be transcribed or repressed, but it is difficult to determine experimentally. Hence an in silico identification of promoters is crucial in order to guide experimental work and to pin point the key region that controls the transcription initiation of a gene. In this analysis, we demonstrate that while the promoter regions are in general less stable than the flanking regions, their average free energy varies depending on the GC composition of the flanking genomic sequence. We have therefore obtained a set of free energy threshold values, for genomic DNA with varying GC content and used them as generic criteria for predicting promoter regions in several microbial genomes, using an in-house developed tool `PromPredict'. On applying it to predict promoter regions corresponding to the 1144 and 612 experimentally validated TSSs in E. coli (50.8% GC) and B. subtilis (43.5% GC) sensitivity of 99% and 95% and precision values of 58% and 60%, respectively, were achieved. For the limited data set of 81 TSSs available for M. tuberculosis (65.6% GC) a sensitivity of 100% and precision of 49% was obtained.

Designing Dense Radar Networks for Weather Surveillance in Tropical Regions.

Design considerations are presented for a dense weather radar network to support multiple services including aviation. Conflicts, tradeoffs and optimization issues in the context of operation in a tropical region are brought out. The upcoming Indian radar network is used as a case study. Algorithms for data mosaicing are briefly outlined.

Spatial and vertical heterogeneities in aerosol properties over oceanic regions around India: Implications for radiative forcing

The influence of atmospheric aerosols on Earth's radiation budget and hence climate, though well recognized and extensively investigated in recent years, remains largely uncertain mainly because of the large spatio-temporal heterogeneity and the lack of data with adequate resolution. To characterize this diversity, a major multi-platform field campaign ICARB (Integrated Campaign for Aerosols, gases and Radiation Budget) was carried out during the pre-monsoon period of 2006 over the Indian landmass and surrounding oceans, which was the biggest such campaign ever conducted over this region. Based on the extensive and concurrent measurements of the optical and physical properties of atmospheric aerosols during ICARB, the spatial distribution of aerosol radiative forcing was estimated over the entire Bay of Bengal (BoB), northern Indian Ocean and Arabian Sea (AS) as well as large spatial variations within these regions. Besides being considerably lower than the mean values reported earlier for this region, our studies have revealed large differences in the forcing components between the BoB and the AS. While the regionally averaged aerosol-induced atmospheric forcing efficiency was 31 +/- 6 W m(-2) tau(-1) for the BoB, it was only similar to 18 +/- 7 W m(-2) tau(-1) for the AS. Airborne measurements revealed the presence of strong, elevated aerosol layers even over the oceans, leading to vertical structures in the atmospheric forcing, resulting in significant warming in the lower troposphere. These observations suggest serious climate implications and raise issues ranging from the impact of aerosols on vertical thermal structure of the atmospheric and hence cloud formation processes to monsoon circulation.

Functional diversity of human protein kinase splice variants marks significant expansion of human kinome

Background: Protein kinases are involved in diverse spectrum of cellular processes. Availability of draft version of the human genomic data in the year 2001 enabled recognition of repertoire of protein kinases. However, over the years the human genomic data is being refined and the current release of human genomic data has helped us to recognize a larger repertoire of over 900 human protein kinases represented mainly by splice variants. Results: Many of these identified protein kinases are alternatively spliced products. Interestingly, some of the human kinase splice variants appear to be significantly diverged in terms of their functional properties as represented by incorporation or absence of one or more domains. Many sets of protein kinase splice variants have substantially different domain organization and in a few sets of splice variants kinase domains belong to different subfamilies of kinases suggesting potential participation in different signal transduction pathways. Conclusions: Addition or deletion of a domain between splice variants of multi-domain kinases appears to be a means of generating differences in the functional features of otherwise similar kinases. It is intriguing that marked sequence diversity within the catalytic regions of some of the splice variant kinases result in kinases belonging to different subfamilies. These human kinase splice variants with different functions might contribute to diversity of eukaryotic cellular signaling.

ATM, RAD50 and p53 in breast cancer

Breast cancer is the most commonly occurring cancer among women, and its incidence is increasing worldwide. Positive family history is a well established risk factor for breast cancer, and it is suggested that the proportion of breast cancer that can be attributed to genetic factors may be as high as 30%. However, all the currently known breast cancer susceptibility genes are estimated to account for 20-30% of familial breast cancer, and only 5% of the total breast cancer incidence. It is thus likely that there are still other breast cancer susceptibility genes to be found. Cellular responses to DNA damage are crucial for maintaining genomic integrity and preventing the development of cancer. The genes operating in DNA damage response signaling network are thus good candidates for breast cancer susceptibility genes. The aim of this study was to evaluate the role of three DNA damage response associated genes, ATM, RAD50, and p53, in breast cancer. ATM, a gene causative for ataxia telangiectasia (A-T), has long been a strong candidate for a breast cancer susceptibility gene because of its function as a key DNA damage signal transducer. We analyzed the prevalence of known Finnish A-T related ATM mutations in large series of familial and unselected breast cancer cases from different geographical regions in Finland. Of the seven A-T related mutations, two were observed in the studied familial breast cancer patients. Additionally, a third mutation previously associated with breast cancer susceptibility was also detected. These founder mutations may be responsible for excess familial breast cancer regionally in Northern and Central Finland, but in Southern Finland our results suggest only a minor effect, if any, of any ATM genetic variants on familial breast cancer. We also screened the entire coding region of the ATM gene in 47 familial breast cancer patients from Southern Finland, and evaluated the identified variants in additional cases and controls. All the identified variants were too rare to significantly contribute to breast cancer susceptibility. However, the role of ATM in cancer development and progression was supported by the results of the immunohistochemical studies of ATM expression, as reduced ATM expression in breast carcinomas was found to correlate with tumor differentiation and hormone receptor status. Aberrant ATM expression was also a feature shared by the BRCA1/2 and the difficult-to-treat ER/PR/ERBB2-triple-negative breast carcinomas. From the clinical point of view, identification of phenotypic and genetic similarities between the BRCA1/2 and the triple-negative breast tumors could have an implication in designing novel targeted therapies to which both of these classes of breast cancer might be exceptionally sensitive. Mutations of another plausible breast cancer susceptibility gene, RAD50, were found to be very rare, and RAD50 can only be making a minor contribution to familial breast cancer predisposition in UK and Southern Finland. The Finnish founder mutation RAD50 687delT seems to be a null allele and may carry a small increased risk of breast cancer. RAD50 is not acting as a classical tumor suppressor gene, but it is possible that RAD50 haploinsufficiency is contributing to cancer. In addition to relatively rare breast cancer susceptibility alleles, common polymorphisms may also be associated with increased breast cancer risk. Furthermore, these polymorphisms may have an impact on the progression and outcome of the disease. Our results suggest no effect of the common p53 R72P polymorphism on familial breast cancer risk or breast cancer risk in the population, but R72P seems to be associated with histopathologic features of the tumors and survival of the patients; 72P homozygous genotype was an independent prognostic factor among the unselected breast cancer patients, with a two-fold increased risk of death. These results present important novel findings also with clinical significance, as codon 72 genotype could be a useful additional prognostic marker in breast cancer, especially among the subgroup of patients with wild-type p53 in their tumors.

Structure of Mycobacterium tuberculosis RuvA, a protein involved in recombination

The process of recombinational repair is crucial for maintaining genomic integrity and generating biological diversity. In association with RuvB and RuvC, RuvA plays a central role in processing and resolving Holliday junctions, which are a critical intermediate in homologous recombination. Here, the cloning, purification and structure determination of the RuvA protein from Mycobacterium tuberculosis (MtRuvA) are reported. Analysis of the structure and comparison with other known RuvA proteins reveal an octameric state with conserved subunit-subunit interaction surfaces, indicating the requirement of octamer formation for biological activity. A detailed analysis of plasticity in the RuvA molecules has led to insights into the invariant and variable regions, thus providing a framework for understanding regional flexibility in various aspects of RuvA function.

A theoretical model for the magnetic helicity of solar active regions

Active regions on the solar surface are known to possess magnetic helicity, which is predominantly negative in the northern hemisphere and positive in the southern hemisphere. Choudhuri et al. [Choudhuri, A.R. On the connection between mean field dynamo theory and flux tubes. Solar Phys. 215, 31–55, 2003] proposed that the magnetic helicity arises due to the wrapping up of the poloidal field of the convection zone around rising flux tubes which form active regions. Choudhuri [Choudhuri, A.R., Chatterjee, P., Nandy, D. Helicity of solar active regions from a dynamo model. ApJ 615, L57–L60, 2004] used this idea to calculate magnetic helicity from their solar dynamo model. Apart from getting broad agreements with observational data, they also predict that the hemispheric helicity rule may be violated at the beginning of a solar cycle. Chatterjee et al. [Chatterjee, P., Choudhuri, A.R., Petrovay, K. Development of twist in an emerging magnetic flux tube by poloidal field accretion. A&A 449, 781–789, 2006] study the penetration of the wrapped poloidal field into the rising flux tube due to turbulent diffusion using a simple 1-d model. They find that the extent of penetration of the wrapped field will depend on how weak the magnetic field inside the rising flux tube becomes before its emergence. They conclude that more detailed observational data will throw light on the physical conditions of flux tubes just before their emergence to the photosphere.