Origin and diversity of the Sox transcription factor gene family: genome-wide analysis in Fugu rubripes

The SOX family of transcription factors are found throughout the animal kingdom and are important in a variety of developmental contexts. Genome analysis has identified 20 Sox genes in human and mouse, which can be subdivided into 8 groups, based on sequence comparison and intron-exon structure. Most of the SOX groups identified in mammals are represented by a single SOX sequence in invertebrate model organisms, suggesting a duplication and divergence mechanism has operated during vertebrate evolution. We have now analysed the Sox gene complement in the pufferfish, Fugu rubripes, in order to shed further light on the diversity and origins of the Sox gene family. Major differences were found between the Sox family in Fugu and those in humans and mice. In particular, Fugu does not have orthologues of Sry, Sox,15 and Sox30, which appear to be specific to mammals, while Sox19, found in Fugu and zebrafish but absent in mammals, seems to be specific to fishes. Six mammalian Sox genes are represented by two copies each in Fugu, indicating a large-scale gene duplication in the fish lineage. These findings point to recent Sox gene loss, duplication and divergence occurring during the evolution of tetrapod and teleost lineages, and provide further evidence for large-scale segmental or a whole-genome duplication occurring early in the radiation of teleosts. (C) 2004 Elsevier B.V. All rights reserved.

Repressor- and activator-type ethylene response factors functioning in jasmonate signaling and disease resistance identified via a genome-wide screen of Arabidopsis transcription factor gene expression

To identify transcription factors (TFs) involved in jasmonate (JA) signaling and plant defense, we screened 1,534 Arabidopsis (Arabidopsis thaliana) TFs by real-time quantitative reverse transcription-PCR for their altered transcript at 6 h following either methyl JA treatment or inoculation with the incompatible pathogen Alternaria brassicicola. We identified 134 TFs that showed a significant change in expression, including many APETALA2/ethylene response factor (AP2/ERF), MYB, WRKY, and NACTF genes with unknown functions. Twenty TF genes were induced by both the pathogen and methyl JA and these included 10 members of the AP2/ERF TF family, primarily from the B1a and B3 subclusters. Functional analysis of the B1a TF AtERF4 revealed that AtERF4 acts as a novel negative regulator of JA-responsive defense gene expression and resistance to the necrotrophic fungal pathogen Fusarium oxysporum and antagonizes JA inhibition of root elongation. In contrast, functional analysis of the B3 TF AtERF2 showed that AtERF2 is a positive regulator of JA-responsive defense genes and resistance to F. oxysporum and enhances JA inhibition of root elongation. Our results suggest that plants coordinately express multiple repressor-and activator-type AP2/ERFs during pathogen challenge to modulate defense gene expression and disease resistance.

Assumption-free estimation of heritability from genome-wide identity-by-descent sharing between full siblings

The study of continuously varying, quantitative traits is important in evolutionary biology, agriculture, and medicine. Variation in such traits is attributable to many, possibly interacting, genes whose expression may be sensitive to the environment, which makes their dissection into underlying causative factors difficult. An important population parameter for quantitative traits is heritability, the proportion of total variance that is due to genetic factors. Response to artificial and natural selection and the degree of resemblance between relatives are all a function of this parameter. Following the classic paper by R. A. Fisher in 1918, the estimation of additive and dominance genetic variance and heritability in populations is based upon the expected proportion of genes shared between different types of relatives, and explicit, often controversial and untestable models of genetic and non-genetic causes of family resemblance. With genome-wide coverage of genetic markers it is now possible to estimate such parameters solely within families using the actual degree of identity-by-descent sharing between relatives. Using genome scans on 4,401 quasi-independent sib pairs of which 3,375 pairs had phenotypes, we estimated the heritability of height from empirical genome-wide identity-by-descent sharing, which varied from 0.374 to 0.617 (mean 0.498, standard deviation 0.036). The variance in identity-by-descent sharing per chromosome and per genome was consistent with theory. The maximum likelihood estimate of the heritability for height was 0.80 with no evidence for non-genetic causes of sib resemblance, consistent with results from independent twin and family studies but using an entirely separate source of information. Our application shows that it is feasible to estimate genetic variance solely from within- family segregation and provides an independent validation of previously untestable assumptions. Given sufficient data, our new paradigm will allow the estimation of genetic variation for disease susceptibility and quantitative traits that is free from confounding with non-genetic factors and will allow partitioning of genetic variation into additive and non-additive components.

Genome-wide scan of IQ finds significant linkage to a quantitative trait locus on 2q

A genome-wide linkage scan of 795 microsatellite markers (761 autosomal, 34 X chromosome) was performed on Multidimensional Aptitude Battery subtests and verbal, performance and full scale scores, the WAIS-R Digit Symbol subtest, and two word-recognition tests (Schonell Graded Word Reading Test, Cambridge Contextual Reading Test) highly predictive of IQ. The sample included 361 families comprising 2-5 siblings who ranged in age from 15.7 to 22.2 years; genotype, but not phenotype, data were available for 81% of parents. A variance components analysis which controlled for age and sex effects showed significant linkage for the Cambridge reading test and performance IQ to the same region on chromosome 2, with respective LOD scores of 4.15 and 3.68. Suggestive linkage (LOD score > 2.2) for various measures was further supported on chromosomes 6, 7, 11, 14, 21 and 22. Where location of linkage peaks converged for IQ subtests within the same scale, the overall scale score provided increased evidence for linkage to that region over any individual subtest. Association studies of candidate genes, particularly those involved in neural transmission and development, will be directed to genes located under the linkage peaks identified in this study.

Genome-wide analysis of mammalian promoter architecture and evolution

Mammalian promoters can be separated into two classes, conserved TATA box-enriched promoters, which initiate at a welldefined site, and more plastic, broad and evolvable CpG-rich promoters. We have sequenced tags corresponding to several hundred thousand transcription start sites (TSSs) in the mouse and human genomes, allowing precise analysis of the sequence architecture and evolution of distinct promoter classes. Different tissues and families of genes differentially use distinct types of promoters. Our tagging methods allow quantitative analysis of promoter usage in different tissues and show that differentially regulated alternative TSSs are a common feature in protein-coding genes and commonly generate alternative N termini. Among the TSSs, we identified new start sites associated with the majority of exons and with 3' UTRs. These data permit genome-scale identification of tissue-specific promoters and analysis of the cis-acting elements associated with them.

Genome-wide linkage scan for loci influencing plasma triglycerides

Background: Plasma triglyceride concentration is known to be a significant risk factor for cardiovascular disease (CVD). Previous studies have found that the level of triglycerides is strongly influenced by genetic factors. Methods: To identify quantitative trait loci influencing triglycerides, we conducted a genome-wide linkage scan on data from 485 Australian adult dizygotic twin pairs. Prior to linkage analysis, triglyceride values were adjusted for the effects of covariates including age, sex, time since last meal, time of blood collection (CT) and time to plasma separation. Results: The heritability estimate for ln(triglyceride) adjusted for all above fixed effects was 0.49. The highest multipoint LOD score observed was 2.94 (genome-wide p=0.049) on chromosome 7 (at 65cM). This 7p region contains several candidate genes. Two other regions with suggestive multipoint LOD scores were also identified on chromosome 4 (LOD score=2.26 at 62cM) and chromosome X (LOD score=2.01 at 81cM). Conclusions: The linkage peaks found represent newly identified regions for more detailed study, in particular the significant linkage observed on chromosome 7p13. \ (c) 2006 Elsevier B.V. All rights reserved.

Major quantitative trait locus for eosinophil count is located on chromosome 2q

Background: Eosinophils are granulocytic white blood cells implicated in asthma and atopic disease. The degree of eosinophilia in the blood of patients with asthma correlates with the severity of asthmatic symptoms. Quantitative trait loci (QTL) linkage analysis of eosinophil count may be a more powerful strategy of mapping genes involved in asthma than linkage analysis using affected relative pairs. 1 Objective: To identify QTLs responsible for variation in eosinophil count in adolescent twins. Methods: We measured eosinophil count longitudinally in 738 pairs of twins at 12, 14, and 16 years of age. We typed 757 highly polymorphic microsatellite markers at an average spacing of similar to5 centimorgans across the genome. We then used multipoint variance components linkage analysis to test for linkage between marker loci and eosinophil concentrations at each age across the genome. Results: We found highly significant linkage on chromosome 2q33 in 12-year-old twins (logarithm of the odds = 4.6; P = .000002) and suggestive evidence of linkage in the same region in 14-year-olds (logarithm of the odds = 1.0; P = .016). We also found suggestive evidence of linkage at other areas of the genome, including regions on chromosomes 2, 3, 4, 8, 9, 11, 12, 17, 20, and 22. Conclusion: A QTL for eosinophil count is present on chromosome 2q33. This QTL might represent a gene involved in asthma pathophysiology.

Lack of association between genetic markers on chromosome 16q22-Q24 and type 1 diabetes in Russian affected families

Aim To evaluate whether the T1D susceptibility locus on chromosome 16q contributes to the genetic susceptibility to T1D in Russian patients. Method Thirteen microsatellite markers, spanning a 47-centimorgan genomic region on 16q22-q24 were evaluated for linkage to T1D in 98 Russian multiplex families. Multipoint logarithm of odds (LOD) ratio (MLS) and nonparametric LOD (NPL) values were computed for each marker, using GENEHUNTER 2.1 software. Four microsatellites (D16S422, D16S504, D16S3037, and D16S3098) and 6 biallelic markers in 2 positional candidate genes, ICSBP1 and NQO1, were additionally tested for association with T1D in 114 simplex families, using transmission disequilibrium test (TDT). Results A peak of linkage (MLS = 1.35, NPL = 0.91) was shown for marker D16S750, but this was not significant (P = 0.18). The subsequent linkage analysis in the subset of 46 multiplex families carrying a common risk HLA-DR4 haplotype increased peak MLS and NPL values to 1.77 and 1.22, respectively, but showed no significant linkage (P = 0.11) to T1D in the 16q22-q24 genomic region. TDT analysis failed to find significant association between these markers and disease, even after the conditioning for the predisposing HLA-DR4 haplotype. Conclusion Our results did not support the evidence for the susceptibility locus to T1D on chromosome 16q22-24 in the Russian family data set. The lack of association could reflect genetic heterogeneity of type 1 diabetes in diverse ethnic groups.

The TAF5L gene on chromosome 1q42 is associated with type 1 diabetes in Russian affected patients

Type 1 diabetes (T1D) is a multifactorial autoimmune disease, with strong genetic component. Several susceptibility loci contribute to genetic predisposition to T1D. One of these loci have been mapped to chromosome 1q42 in UK and US joined affected family data sets but needs to be replicated in other populations. In this study, we evaluated sixteen microsatellites located on 1q42 for linkage with T1D in 97 Russian affected sibling pairs. A 2.7-cm region of suggestive linkage to T1D between markers D1S1644 and D1S225 was found by multipoint linkage analysis. The peak of linkage was shown for D1S2847 (P = 0.0005). Transmission disequilibrium test showed significant undertransmission of the 156-bp allele of D1S2847 from parents to diabetic children (28 transmissions vs. 68 nontransmissions, P = 0.043) in Russian affected families. A preferential transmission from parents to diabetic offspring was also shown for the T(-25) and T1362 alleles of the C/T(-25) and C/T1362 dimorphisms, both located at the TAF5L gene, which is situated 103 kb from D1S2847. Together with the A/C744 TAF5L SNP, these markers share common T(-25)/A744/T1362 and C(-25)/C744/T1362 haplotypes associated with higher and lower risk of diabetes (Odds Ratio = 2.15 and 0.62, respectively). Our results suggest that the TAF5L gene, encoding TAF5L-like RNA polymerase II p300/CBP associated factor (PCAF)-associated factor, could represent the susceptibility gene for T1D on chromosome 1q42 in Russian affected patients.

Clusters of internally primed transcripts reveal novel long noncoding RNAs

Non- protein- coding RNAs ( ncRNAs) are increasingly being recognized as having important regulatory roles. Although much recent attention has focused on tiny 22- to 25- nucleotide microRNAs, several functional ncRNAs are orders of magnitude larger in size. Examples of such macro ncRNAs include Xist and Air, which in mouse are 18 and 108 kilobases ( Kb), respectively. We surveyed the 102,801 FANTOM3 mouse cDNA clones and found that Air and Xist were present not as single, full- length transcripts but as a cluster of multiple, shorter cDNAs, which were unspliced, had little coding potential, and were most likely primed from internal adenine- rich regions within longer parental transcripts. We therefore conducted a genome- wide search for regional clusters of such cDNAs to find novel macro ncRNA candidates. Sixty- six regions were identified, each of which mapped outside known protein- coding loci and which had a mean length of 92 Kb. We detected several known long ncRNAs within these regions, supporting the basic rationale of our approach. In silico analysis showed that many regions had evidence of imprinting and/ or antisense transcription. These regions were significantly associated with microRNAs and transcripts from the central nervous system. We selected eight novel regions for experimental validation by northern blot and RT- PCR and found that the majority represent previously unrecognized noncoding transcripts that are at least 10 Kb in size and predominantly localized in the nucleus. Taken together, the data not only identify multiple new ncRNAs but also suggest the existence of many more macro ncRNAs like Xist and Air.

A genome scan for eye color in 502 twin families: Most variation is due to a QTL on chromosome 15q

We have rated eye color on a 3-point scale (1=blue/grey, 2=hazel/green, 3=brown) in 502 twin families and carried out a 5-10 cM genome scan (400-757 markers). We analyzed eye color as a threshold trait and performed multipoint sib pair linkage analysis using variance components analysis in Mx. A lod of 19.2 was found at the marker D15S1002, less than 1 cM from OCA2, which has been previously implicated in eye color variation. We estimate that 74% of variance in eye color liability is due to this QTL and a further 18% due to polygenic effects. However, a large shoulder on this peak suggests that other loci affecting eye color may be telomeric of OCA2 and inflating the QTL estimate. No other peaks reached genome-wide significance, although lods >2 were seen on 5p and 14q and lods >1 were additionally seen on chromosomes 2, 3, 6, 7, 8, 9, 17 and 18. Most of these secondary peaks were reduced or eliminated when we repeated the scan as a two locus analysis with the 15q linkage included, although this does not necessarily exclude them as false positives. We also estimated the interaction between the 15q QTL and the other marker locus but there was only minor evidence for additive x additive epistasis. Elaborating the analysis to the full two-locus model including non-additive main effects and interactions did not strengthen the evidence for epistasis. We conclude that most variation in eye color in Europeans is due to polymorphism in OCA2 but that there may be modifiers at several other loci.