Search

Brain L-serine is critical for neurodevelopment and is thought to be synthesized solely from glucose. In contrast, we found that the influx of L-serine across the blood-brain barrier (BBB) is essential for brain development. We identified the endothelial Slc38a5, previously thought to be a glutamine transporter, as an L-serine transporter expressed at the BBB in early postnatal life. Young Slc38a5 knockout (KO) mice exhibit developmental alterations and a decrease in brain L-serine and D-serine, without changes in serum or liver amino acids. Slc38a5-KO brains exhibit accumulation of neurotoxic deoxysphingolipids, synaptic and mitochondrial abnormalities, and decreased neurogenesis at the dentate gyrus. Slc38a5-KO pups exhibit motor impairments that are affected by the administration of L-serine at concentrations that replenish the serine pool in the brain. Our results highlight a critical role of Slc38a5 in supplying L-serine via the BBB for proper brain development.

The superior photosynthetic efficiency of C 4 leaves over C 3 leaves is owing to their unique Kranz anatomy, in which the vein is surrounded by one layer of bundle sheath (BS) cells and one layer of mesophyll (M) cells. Kranz anatomy development starts from three contiguous ground meristem (GM) cells, but its regulators and underlying molecular mechanism are largely unknown. To identify the regulators, we obtained the transcriptomes of 11 maize embryonic leaf cell types from five stages of pre-Kranz cells starting from median GM cells and six stages of pre-M cells starting from undifferentiated cells. Principal component and clustering analyses of transcriptomic data revealed rapid pre-Kranz cell differentiation in the first two stages but slow differentiation in the last three stages, suggesting early Kranz cell fate determination. In contrast, pre-M cells exhibit a more prolonged transcriptional differentiation process. Differential gene expression and coexpression analyses identified gene coexpression modules, one of which included 3 auxin transporter and 18 transcription factor (TF) genes, including known regulators of Kranz anatomy and/or vascular development. In situ hybridization of 11 TF genes validated their expression in early Kranz development. We determined the binding motifs of 15 TFs, predicted TF target gene relationships among the 18 TF and 3 auxin transporter genes, and validated 67 predictions by electrophoresis mobility shift assay. From these data, we constructed a gene regulatory network for Kranz development. Our study sheds light on the regulation of early maize leaf development and provides candidate leaf development regulators for future study.

Transcription factor binding sites (TFBSs) are essential for gene regulation, but the number of known TFBSs remains limited. We aimed to discover and characterize unknown TFBSs by developing a computational pipeline for analyzing ChIP-seq (chromatin immunoprecipitation followed by sequencing) data. Applying it to the latest ENCODE ChIP-seq data for human and mouse, we found that using the irreproducible discovery rate as a quality-control criterion resulted in many experiments being unnecessarily discarded. By contrast, the number of motif occurrences in ChIP-seq peak regions provides a highly effective criterion, which is reliable even if supported by only one experimental replicate. In total, we obtained 2,058 motifs from 1,089 experiments for 354 human TFs and 163 motifs from 101 experiments for 34 mouse TFs. Among these motifs, 487 have not previously been reported. Mapping the canonical motifs to the human genome reveals a high TFBS density ±2 kb around transcription start sites (TSSs) with a peak at -50 bp. On average, a promoter contains 5.7 TFBSs. However, 70% of TFBSs are in introns (41%) and intergenic regions (29%), whereas only 12% are in promoters (-1 kb to +100 bp from TSSs). Notably, some TFs (e.g., CTCF, JUN, JUNB, and NFE2) have motifs enriched in intergenic regions, including enhancers. We inferred 142 cobinding TF pairs and 186 (including 115 completely) tethered binding TF pairs, indicating frequent interactions between TFs and a higher frequency of tethered binding than cobinding. This study provides a large number of previously undocumented motifs and insights into the biological and genomic features of TFBSs., Competing Interests: The authors declare no competing interest.

Arabidopsis AINTEGUMENTA (ANT), an AP2 transcription factor, is known to control plant growth and floral organogenesis. In this study, our transcriptome analysis and in situ hybridization assays of maize embryonic leaves suggested that maize ANT1 (ZmANT1) regulates vascular development. To better understand ANT1 functions, we determined the binding motif of ZmANT1 and then showed that ZmANT1 binds the promoters of millet SCR1 , GNC , and AN3 , which are key regulators of Kranz anatomy, chloroplast development, and plant growth, respectively. We generated a mutant with a single-codon deletion and two frameshift mutants of the ANT1 ortholog in the C4 millet Setaria viridis by the CRISPR/Cas9 technique. The two frameshift mutants displayed reduced photosynthesis efficiency and growth rate, smaller leaves, and lower grain yields than wild-type (WT) plants. Moreover, their leaves sporadically exhibited distorted Kranz anatomy and vein spacing. Conducting transcriptomic analysis of developing leaves in the WT and the three mutants we identified differentially expressed genes (DEGs) in the two frameshift mutant lines and found many down-regulated DEGs enriched in photosynthesis, heme, tetrapyrrole binding, and antioxidant activity. In addition, we predicted many target genes of ZmANT1 and chose 13 of them to confirm binding of ZmANT1 to their promoters. Based on the above observations, we proposed a model for ANT1 regulation of cell proliferation and leaf growth, vascular and vein development, chloroplast development, and photosynthesis through its target genes. Our study revealed biological roles of ANT1 in several developmental processes beyond its known roles in plant growth and floral organogenesis., Competing Interests: The authors declare no competing interest.

Cellulosomes, which are multienzyme complexes from anaerobic bacteria, are considered nature's finest cellulolytic machinery. Thus, constructing a cellulosome in an industrial yeast has long been a goal pursued by scientists. However, it remains highly challenging due to the size and complexity of cellulosomal genes. Here, we overcame the difficulties by synthesizing the Clostridium thermocellum scaffoldin gene ( CipA ) and the anchoring protein gene ( OlpB ) using advanced synthetic biology techniques. The engineered Kluyveromyces marxianus , a probiotic yeast, secreted a mixture of dockerin-fused fungal cellulases, including an endoglucanase ( TrEgIII ), exoglucanase ( CBHII ), β-glucosidase ( NpaBGS ), and cellulase boosters ( TaLPMO and MtCDH ). The confocal microscopy results confirmed the cell-surface display of OlpB- Sc GPI and fluorescence-activated cell sorting analysis results revealed that almost 81% of yeast cells displayed OlpB- Sc GPI. We have also demonstrated the cellulosome complex formation using purified and crude cellulosomal proteins. Native polyacrylamide gel electrophoresis and mass spectrometric analysis further confirmed the cellulosome complex formation. Our engineered cellulosome can accommodate up to 63 enzymes, whereas the largest engineered cellulosome reported thus far could accommodate only 12 enzymes and was expressed by a plasmid instead of chromosomal integration. Interestingly, CipA 2B9C (with two cellulose binding modules, CBM) released significantly higher quantities of reducing sugars compared with other CipA variants, thus confirming the importance of cohesin numbers and CBM domain on cellulosome complex. The engineered yeast host efficiently degraded cellulosic substrates and released 3.09 g/L and 8.61 g/L of ethanol from avicel and phosphoric acid-swollen cellulose, respectively, which is higher than any previously constructed yeast cellulosome., Competing Interests: The authors declare no competing interest.

How negative selection, positive selection, and population size contribute to the large variation in nucleotide substitution rates among RNA viruses remains unclear. Here, we studied the ratios of nonsynonymous-to-synonymous substitution rates ( d N / d S ) in protein-coding genes of human RNA and DNA viruses and mammals. Among the 21 RNA viruses studied, 18 showed a genome-average d N / d S from 0.01 to 0.10, indicating that over 90% of nonsynonymous mutations are eliminated by negative selection. Only HIV-1 showed a d N / d S (0.31) higher than that (0.22) in mammalian genes. By comparing the d N / d S values among genes in the same genome and among species or strains, we found that both positive selection and population size play significant roles in the d N / d S variation among genes and species. Indeed, even in flaviviruses and picornaviruses, which showed the lowest ratios among the 21 species studied, positive selection appears to have contributed significantly to d N / d S We found the view that positive selection occurs much more frequently in influenza A subtype H3N2 than subtype H1N1 holds only for the hemagglutinin and neuraminidase genes, but not for other genes. Moreover, we found no support for the view that vector-borne RNA viruses have lower d N / d S ratios than non-vector-borne viruses. In addition, we found a correlation between d N and d S , implying a correlation between d N and the mutation rate. Interestingly, only 2 of the 8 DNA viruses studied showed a d N / d S < 0.10, while 4 showed a d N / d S > 0.22. These observations increase our understanding of the mechanisms of RNA virus evolution., Competing Interests: The authors declare no conflict of interest.

Tetep is a rice cultivar known for broad-spectrum resistance to blast, a devastating fungal disease. The molecular basis for its broad-spectrum resistance is still poorly understood. Is it because Tetep has many more NLR genes than other cultivars? Or does Tetep possess multiple major NLR genes that can individually confer broad-spectrum resistance to blast? Moreover, are there many interacting NLR pairs in the Tetep genome? We sequenced its genome, obtained a high-quality assembly, and annotated 455 nucleotide-binding site leucine-rich repeat ( NLR ) genes. We cloned and tested 219 NLR genes as transgenes in 2 susceptible cultivars using 5 to 12 diversified pathogen strains; in many cases, fewer than 12 strains were successfully cultured for testing. Ninety cloned NLR s showed resistance to 1 or more pathogen strains and each strain was recognized by multiple NLR s. However, few NLR s showed resistance to >6 strains, so multiple NLR s are apparently required for Tetep's broad-spectrum resistance to blast. This was further supported by the pedigree analyses, which suggested a correlation between resistance and the number of Tetep-derived NLR s. In developing a method to identify NLR pairs each of which functions as a unit, we found that >20% of the NLR s in the Tetep and 3 other rice genomes are paired. Finally, we designed an extensive set of molecular markers for rapidly introducing clustered and paired NLR s in the Tetep genome for breeding new resistant cultivars. This study increased our understanding of the genetic basis of broad-spectrum blast resistance in rice., Competing Interests: The authors declare no conflict of interest.

Time-series transcriptomes of a biological process obtained under different conditions are useful for identifying the regulators of the process and their regulatory networks. However, such data are 3D (gene expression, time, and condition), and there is currently no method that can deal with their full complexity. Here, we developed a method that avoids time-point alignment and normalization between conditions. We applied it to analyze time-series transcriptomes of developing maize leaves under light-dark cycles and under total darkness and obtained eight time-ordered gene coexpression networks (TO-GCNs), which can be used to predict upstream regulators of any genes in the GCNs. One of the eight TO-GCNs is light-independent and likely includes all genes involved in the development of Kranz anatomy, which is a structure crucial for the high efficiency of photosynthesis in C 4 plants. Using this TO-GCN, we predicted and experimentally validated a regulatory cascade upstream of SHORTROOT1 , a key Kranz anatomy regulator. Moreover, we applied the method to compare transcriptomes from maize and rice leaf segments and identified regulators of maize C 4 enzyme genes and RUBISCO SMALL SUBUNIT2 Our study provides not only a powerful method but also novel insights into the regulatory networks underlying Kranz anatomy development and C 4 photosynthesis., Competing Interests: The authors declare no conflict of interest.

Repeated artificial selection of a complex trait facilitates the identification of genes underlying the trait, especially if multiple selected descendant lines are available. Here we developed a pedigree-based approach to identify genes underlying the Green Revolution (GR) phenotype. From a pedigree analysis, we selected 30 cultivars including the "miracle rice" IR8, a GR landmark, its ancestors and descendants, and also other related cultivars for identifying high-yield genes. Through sequencing of these genomes, we identified 28 ancestral chromosomal blocks that were maintained in all the high-yield cultivars under study. In these blocks, we identified six genes of known function, including the GR gene sd1 , and 123 loci with genes of unknown function. We randomly selected 57 genes from the 123 loci for knockout or knockdown studies and found that a high proportion of these genes are essential or have phenotypic effects related to rice production. Notably, knockout lines have significant changes in plant height ( P < 0.003), a key GR trait, compared with wild-type lines. Some gene knockouts or knockdowns were especially interesting. For example, knockout of Os10g0555100, a putative glucosyltransferase gene, showed both reduced growth and altered panicle architecture. In addition, we found that in some retained chromosome blocks several GR-related genes were clustered, although they have unrelated sequences, suggesting clustering of genes with similar functions. In conclusion, we have identified many high-yield genes in rice. Our method provides a powerful means to identify genes associated with a specific trait., Competing Interests: The authors declare no conflict of interest.

High vein density, a distinctive trait of C 4 leaves, is central to both C 3 -to-C 4 evolution and conversion of C 3 to C 4 -like crops. We tested the hypothesis that high vein density in C 4 leaves is due to elevated auxin biosynthesis and transport in developing leaves. Up-regulation of genes in auxin biosynthesis pathways and higher auxin content were found in developing C 4 leaves compared with developing C 3 leaves. The same observation held for maize foliar (C 4 ) and husk (C 3 ) leaf primordia. Moreover, auxin content and vein density were increased in loss-of-function mutants of Arabidopsis MYC2 , a suppressor of auxin biosynthesis. Treatment with an auxin biosynthesis inhibitor or an auxin transport inhibitor led to much fewer veins in new leaves. Finally, both Arabidopsis thaliana auxin efflux transporter pin1 and influx transporter lax2 mutants showed reduced vein numbers. Thus, development of high leaf vein density requires elevated auxin biosynthesis and transport., Competing Interests: The authors declare no conflict of interest.

CD38 catalyzes the synthesis of the Ca 2+ messenger, cyclic ADP-ribose (cADPR). It is generally considered to be a type II protein with the catalytic domain facing outside. How it can catalyze the synthesis of intracellular cADPR that targets the endoplasmic Ca 2+ stores has not been resolved. We have proposed that CD38 can also exist in an opposite type III orientation with its catalytic domain facing the cytosol. Here, we developed a method using specific nanobodies to immunotarget two different epitopes simultaneously on the catalytic domain of the type III CD38 and firmly established that it is naturally occurring in human multiple myeloma cells. Because type III CD38 is topologically amenable to cytosolic regulation, we used yeast-two-hybrid screening to identify cytosolic Ca 2+ and integrin-binding protein 1 (CIB1), as its interacting partner. The results from immunoprecipitation, ELISA, and bimolecular fluorescence complementation confirmed that CIB1 binds specifically to the catalytic domain of CD38, in vivo and in vitro. Mutational studies established that the N terminus of CIB1 is the interacting domain. Using shRNA to knock down and Cas9/guide RNA to knock out CIB1, a direct correlation between the cellular cADPR and CIB1 levels was demonstrated. The results indicate that the type III CD38 is functionally active in producing cellular cADPR and that the activity is specifically modulated through interaction with cytosolic CIB1., Competing Interests: The authors declare no conflict of interest.

Accelerated losses of biodiversity are a hallmark of the current era. Large declines of population size have been widely observed and currently 22,176 species are threatened by extinction. The time at which a threatened species began rapid population decline (RPD) and the rate of RPD provide important clues about the driving forces of population decline and anticipated extinction time. However, these parameters remain unknown for the vast majority of threatened species. Here we analyzed the genetic diversity data of nuclear and mitochondrial loci of 2,764 vertebrate species and found that the mean genetic diversity is lower in threatened species than in related nonthreatened species. Our coalescence-based modeling suggests that in many threatened species the RPD began ∼123 y ago (a 95% confidence interval of 20-260 y). This estimated date coincides with widespread industrialization and a profound change in global living ecosystems over the past two centuries. On average the population size declined by ∼25% every 10 y in a threatened species, and the population size was reduced to ∼5% of its ancestral size. Moreover, the ancestral size of threatened species was, on average, ∼22% smaller than that of nonthreatened species. Because the time period of RPD is short, the cumulative effect of RPD on genetic diversity is still not strong, so that the smaller ancestral size of threatened species may be the major cause of their reduced genetic diversity; RPD explains 24.1-37.5% of the difference in genetic diversity between threatened and nonthreatened species., Competing Interests: The authors declare no conflict of interest.

Many secretory tissues release Zn(II) ions along with other molecules in response to external stimuli. Here we demonstrate that secretion of Zn(II) ions from normal, healthy prostate tissue is stimulated by glucose in fasted mice and that release of Zn(II) can be monitored by MRI. An ∼50% increase in water proton signal enhancement is observed in T1-weighted images of the healthy mouse prostate after infusion of a Gd-based Zn(II) sensor and an i.p. bolus of glucose. Release of Zn(II) from intracellular stores was validated in human epithelial prostate cells in vitro and in surgically exposed prostate tissue in vivo using a Zn(II)-sensitive fluorescent probe known to bind to the extracellular surface of cells. Given the known differences in intracellular Zn(II) stores in healthy versus malignant prostate tissues, the Zn(II) sensor was then evaluated in a transgenic adenocarcinoma of the mouse prostate (TRAMP) model in vivo. The agent proved successful in detecting small malignant lesions as early as 11 wk of age, making this noninvasive MR imaging method potentially useful for identifying prostate cancer in situations where it may be difficult to detect using current multiparametric MRI protocols., Competing Interests: A.D.S., W.-H.L., and C.P. have a financial interest in VitalQuan, LLC.

Avian integumentary organs include feathers, scales, claws, and beaks. They cover the body surface and play various functions to help adapt birds to diverse environments. These keratinized structures are mainly composed of corneous materials made of α-keratins, which exist in all vertebrates, and β-keratins, which only exist in birds and reptiles. Here, members of the keratin gene families were used to study how gene family evolution contributes to novelty and adaptation, focusing on tissue morphogenesis. Using chicken as a model, we applied RNA-seq and in situ hybridization to map α- and β-keratin genes in various skin appendages at embryonic developmental stages. The data demonstrate that temporal and spatial α- and β-keratin expression is involved in establishing the diversity of skin appendage phenotypes. Embryonic feathers express a higher proportion of β-keratin genes than other skin regions. In feather filament morphogenesis, β-keratins show intricate complexity in diverse substructures of feather branches. To explore functional interactions, we used a retrovirus transgenic system to ectopically express mutant α- or antisense β-keratin forms. α- and β-keratins show mutual dependence and mutations in either keratin type results in disrupted keratin networks and failure to form proper feather branches. Our data suggest that combinations of α- and β-keratin genes contribute to the morphological and structural diversity of different avian skin appendages, with feather-β-keratins conferring more possible composites in building intrafeather architecture complexity, setting up a platform of morphological evolution of functional forms in feathers.

The prevailing view that the evolution of cells in a tumor is driven by Darwinian selection has never been rigorously tested. Because selection greatly affects the level of intratumor genetic diversity, it is important to assess whether intratumor evolution follows the Darwinian or the non-Darwinian mode of evolution. To provide the statistical power, many regions in a single tumor need to be sampled and analyzed much more extensively than has been attempted in previous intratumor studies. Here, from a hepatocellular carcinoma (HCC) tumor, we evaluated multiregional samples from the tumor, using either whole-exome sequencing (WES) (n = 23 samples) or genotyping (n = 286) under both the infinite-site and infinite-allele models of population genetics. In addition to the many single-nucleotide variations (SNVs) present in all samples, there were 35 "polymorphic" SNVs among samples. High genetic diversity was evident as the 23 WES samples defined 20 unique cell clones. With all 286 samples genotyped, clonal diversity agreed well with the non-Darwinian model with no evidence of positive Darwinian selection. Under the non-Darwinian model, MALL (the number of coding region mutations in the entire tumor) was estimated to be greater than 100 million in this tumor. DNA sequences reveal local diversities in small patches of cells and validate the estimation. In contrast, the genetic diversity under a Darwinian model would generally be orders of magnitude smaller. Because the level of genetic diversity will have implications on therapeutic resistance, non-Darwinian evolution should be heeded in cancer treatments even for microscopic tumors.

Maize is a major crop and a model plant for studying C4 photosynthesis and leaf development. However, a genomewide regulatory network of leaf development is not yet available. This knowledge is useful for developing C3 crops to perform C4 photosynthesis for enhanced yields. Here, using 22 transcriptomes of developing maize leaves from dry seeds to 192 h post imbibition, we studied gene up- and down-regulation and functional transition during leaf development and inferred sets of strongly coexpressed genes. More significantly, we developed a method to predict transcription factor binding sites (TFBSs) and their cognate transcription factors (TFs) using genomic sequence and transcriptomic data. The method requires not only evolutionary conservation of candidate TFBSs and sets of strongly coexpressed genes but also that the genes in a gene set share the same Gene Ontology term so that they are involved in the same biological function. In addition, we developed another method to predict maize TF-TFBS pairs using known TF-TFBS pairs in Arabidopsis or rice. From these efforts, we predicted 1,340 novel TFBSs and 253 new TF-TFBS pairs in the maize genome, far exceeding the 30 TF-TFBS pairs currently known in maize. In most cases studied by both methods, the two methods gave similar predictions. In vitro tests of 12 predicted TF-TFBS interactions showed that our methods perform well. Our study has significantly expanded our knowledge on the regulatory network involved in maize leaf development.

Antrodia cinnamomea, a polyporus mushroom of Taiwan, has long been used as a remedy for cancer, hypertension, and hangover, with an annual market of over $100 million (US) in Taiwan. We obtained a 32.15-Mb genome draft containing 9,254 genes. Genome ontology enrichment and pathway analyses shed light on sexual development and the biosynthesis of sesquiterpenoids, triterpenoids, ergostanes, antroquinonol, and antrocamphin. We identified genes differentially expressed between mycelium and fruiting body and 242 proteins in the mevalonate pathway, terpenoid pathways, cytochrome P450s, and polyketide synthases, which may contribute to the production of medicinal secondary metabolites. Genes of secondary metabolite biosynthetic pathways showed expression enrichment for tissue-specific compounds, including 14-α-demethylase (CYP51F1) in fruiting body for converting lanostane to ergostane triterpenoids, coenzymes Q (COQ) for antroquinonol biosynthesis in mycelium, and polyketide synthase for antrocamphin biosynthesis in fruiting body. Our data will be useful for developing a strategy to increase the production of useful metabolites.

Assembling a large genome using next generation sequencing reads requires large computer memory and a long execution time. To reduce these requirements, we propose an extension-based assembler, called JR-Assembler, where J and R stand for "jumping" extension and read "remapping." First, it uses the read count to select good quality reads as seeds. Second, it extends each seed by a whole-read extension process, which expedites the extension process and can jump over short repeats. Third, it uses a dynamic back trimming process to avoid extension termination due to sequencing errors. Fourth, it remaps reads to each assembled sequence, and if an assembly error occurs by the presence of a repeat, it breaks the contig at the repeat boundaries. Fifth, it applies a less stringent extension criterion to connect low-coverage regions. Finally, it merges contigs by unused reads. An extensive comparison of JR-Assembler with current assemblers using datasets from small, medium, and large genomes shows that JR-Assembler achieves a better or comparable overall assembly quality and requires lower memory use and less central processing unit time, especially for large genomes. Finally, a simulation study shows that JR-Assembler achieves a superior performance on memory use and central processing unit time than most current assemblers when the read length is 150 bp or longer, indicating that the advantages of JR-Assembler over current assemblers will increase as the read length increases with advances in next generation sequencing technology.

DNA damage triggers cell cycle arrest to provide a time window for DNA repair. Failure of arrest could lead to genomic instability and tumorigenesis. DNA damage-induced G1 arrest is generally achieved by the accumulation of Cyclin-dependent kinase inhibitor 1 (p21). However, p21 is degraded and does not play a role in UV-induced G1 arrest. The mechanism of UV-induced G1 arrest thus remains elusive. Here, we have identified a critical role for CUE domain-containing protein 2 (CUEDC2) in this process. CUEDC2 binds to and inhibits anaphase-promoting complex/cyclosome-Cdh1 (APC/C(Cdh1)), a critical ubiquitin ligase in G1 phase, thereby stabilizing Cyclin A and promoting G1-S transition. In response to UV irradiation, CUEDC2 undergoes ERK1/2-dependent phosphorylation and ubiquitin-dependent degradation, leading to APC/C(Cdh1)-mediated Cyclin A destruction, Cyclin-dependent kinase 2 inactivation, and G1 arrest. A nonphosphorylatable CUEDC2 mutant is resistant to UV-induced degradation. Expression of this stable mutant effectively overrides UV-induced G1-S block. These results establish CUEDC2 as an APC/C(Cdh1) inhibitor and indicate that regulated CUEDC2 degradation is critical for UV-induced G1 arrest.

Our anatomical analysis revealed that a dry maize seed contains four to five embryonic leaves at different developmental stages. Rudimentary kranz structure (KS) is apparent in the first leaf with a substantial density, but its density decreases toward younger leaves. Upon imbibition, leaf expansion occurs rapidly with new KSs initiated from the palisade-like ground meristem cells in the middle of the leaf. In parallel to the anatomical analysis, we obtained the time course transcriptomes for the embryonic leaves in dry and imbibed seeds every 6 h up to hour 72. Over this time course, the embryonic leaves exhibit transcripts of 30,255 genes at a level that can be regarded as "expressed." In dry seeds, ∼25,500 genes are expressed, showing functional enrichment in transcription, RNA processing, protein synthesis, primary metabolic pathways, and calcium transport. During the 72-h time course, ∼13,900 genes, including 590 transcription factor genes, are differentially expressed. Indeed, by 30 h postimbibition, ∼2,200 genes expressed in dry seeds are already down-regulated, and ∼2,000 are up-regulated. Moreover, the top 1% expressed genes at 54 h or later are very different from those before 30 h, reflecting important developmental and physiological transitions. Interestingly, clusters of genes involved in hormone metabolism, signaling, and responses are differentially expressed at various time points and TF gene expression is also modular and stage specific. Our dataset provides an opportunity for hypothesizing the timing of regulatory actions, particularly in the context of KS development.

Protein structure and function are closely related, especially in functional surfaces, which are local spatial regions that perform the biological functions. Also, protein structures tend to evolve more slowly than amino acid sequences. We have therefore developed a method to classify proteins using the structures of functional surfaces; we call it protein surface classification (PSC). PSC may reflect functional relationships among proteins and may detect evolutionary relationships among highly divergent sequences. We focused on the surfaces of ligand-bound regions because they represent well-defined structures. Specifically, we used structural attributes to measure similarities between binding surfaces and constructed a PSC library of ~2,000 binding surface types from the bound forms. Using flavin mononucleotide-binding proteins and glycosidases as examples, we show how the evolutionary position of an uncharacterized protein can be defined and its function inferred from the characterized members of the same surface subtype. We found that proteins with the same enzyme nomenclature may be divided into subtypes and that two proteins in the same CATH (Class, Architecture, Topology, Homologous superfamily) fold may belong to two different surface types. In conclusion, our approach complements the sequence-based and fold-domain classifications and has the advantage of associating the shape of a protein with its biological function. As an expandable library, PSC provides a resource of spatial patterns for studying the evolution of protein structure and function.

Current methods of monitoring insulin secretion lack the required spatial and temporal resolution to adequately map the dynamics of exocytosis of native insulin granules in intact cell populations in three dimensions. Exploiting the fact that insulin granules contain a high level of Zn(2+), and that Zn(2+) is coreleased with insulin during secretion, we have developed a fluorescent, cell surface-targeted zinc indicator for monitoring induced exocytotic release (ZIMIR). ZIMIR displayed a robust fluorescence enhancement on Zn(2+) chelation and bound Zn(2+) with high selectivity against Ca(2+) and Mg(2+). When added to cultured β cells or intact pancreatic islets at low micromolar concentrations, ZIMIR labeled cells rapidly, noninvasively, and stably, and it reliably reported changes in Zn(2+) concentration near the sites of granule fusion with high sensitivity that correlated well with membrane capacitance measurement. Fluorescence imaging of ZIMIR-labeled β cells followed the dynamics of exocytotic activity at subcellular resolution, even when using simple epifluorescence microscopy, and located the chief sites of insulin release to intercellular junctions. Moreover, ZIMIR imaging of intact rat islets revealed that Zn(2+)/insulin release occurred largely in small groups of adjacent β cells, with each forming a "secretory unit." Concurrent imaging of ZIMIR and Fura-2 showed that the amplitude of cytosolic Ca(2+) elevation did not necessarily correlate with insulin secretion activity, suggesting that events downstream of Ca(2+) signaling underlie the cell-cell heterogeneity in insulin release. In addition to studying stimulation-secretion coupling in cells with Zn(2+)-containing granules, ZIMIR may find applications in β-cell engineering and screening for molecules regulating insulin secretion on high-throughput platforms.

Since the sensational discovery of a living coelacanth off the east coast of South Africa, the geographic distribution of viable coelacanth populations has been a subject of debate. In the past, the coelacanths off the African mainland were thought to be strays from the Comoros because most coelacanths captured were caught in the waters surrounding the Comoros archipelagos. However, in recent years, a large number of coelacanths were captured off the coast of Tanzania, including nine living specimens observed in a remotely operated vehicles survey. Thus, it is possible that there is a reproducing population inhabiting waters off the Tanzania coast. We have sequenced the complete mitochondrial genomes of 21 Tanzanian and 2 Comoran coelacanths and analyzed these sequences together with two additional full mitochondrial genomes and 47 d-loop sequences from the literature. We found that the coelacanth population off the northern Tanzanian coast is genetically differentiated from those of the southern Tanzania coast and the Comoros, whereas no significant genetic differentiation occurs between the latter two localities. The differentiation between the northern and southern Tanzanian coast populations is consistent with the hypothesis that the existence of northward-flowing ocean current along the Tanzanian coast may reduce or prevent gene flow from the northern to the southern population. Finally, we estimated that the population localized to the southern Tanzanian coast and the Comoros diverged from other coelacanths at least 200,000 y ago. These results indicate that the coelacanths off the northern Tanzania coast are not strays but a genetically distinct group. Our study provides important information for the conservation of this threatened "living fossil."

Protein binding site residues, especially catalytic residues, play a central role in protein function. Because more than 99% of the ∼ 12 million protein sequences in the nonredundant protein database have no structural information, it is desirable to develop methods to predict the binding site residues of a protein from its primary sequence. This task is highly challenging, because the binding site residues constitute only a small portion of a protein. However, the binding site residues of a protein are clustered in its functional pocket(s), and their spatial patterns tend to be conserved in evolution. To take advantage of these evolutionary and structural principles, we constructed a database of ∼ 50,000 templates (called the pocket-containing segment database), each of which includes not only a sequence segment that contains a functional pocket but also the structural attributes of the pocket. To use this database, we designed a template-matching technique, termed residue-matching profiling, and established a criterion for selecting templates for a query sequence. Finally, we developed a probabilistic model for assigning spatial scores to matched residues between the template and query sequence in local alignments using a set of selected scoring matrices and for computing the binding likelihood of each matched residue in the query sequence. From the likelihoods, one can predict the binding site residues in the query sequence. An automated computational pipeline was developed for our method. A performance evaluation shows that our method achieves a 70% precision in predicting binding site residues at 60% sensitivity.

A balance between gene expression stability and evolvability is essential for the long-term maintenance of a living system. In this paper, we studied whether the genetic and epigenetic properties of the promoter affect gene expression variability. We hypothesized that upstream distance and orientation (head-to-head or head-to-tail) are important for the promoter architecture and gene expression variability. We found that in budding yeast genes with a short upstream distance tend to have low gene expression variability, and their promoter is flanked by strongly positioned nucleosomes and tends to have low nucleosome occupancy. These observations suggest that in vivo positioning of the flanking nucleosomes facilitates stable nucleosome depletion at the core promoter region and enhances gene expression stability. Head-to-head genes have, on average, lower gene expression variability, greater nucleosome depletion at the core promoter region, and more strongly positioned nucleosomes that flank the core promoter than do head-to-tail genes. These observations hold for diverse eukaryotes. In complex organisms such as mammals, only a small fraction of head-to-tail genes have retained a short upstream distance, probably because the promoter may not be flanked by a strongly positioned nucleosome on the upstream side.

It is commonly assumed but not proven that microRNAs (miRNAs) and their targets coevolve. Under this assumption, miRNAs and targets from different species may interact adversely, resulting in reduced fitness. However, the strength of the adverse interactions may not be detectable because even outright deletions of miRNAs often manifest only subtle fitness effects. We tested and measured the strength of heterospecific interactions by carrying out transgenic experiments across Drosophila species by overexpressing the miR310s cluster of Drosophila melanogaster (Dm310s) and Drosophila pseudoobscura (Dp310s) in D. melanogaster. Flies overexpressing the heterospecific Dp310s are only one-third as viable as those overexpressing the conspecific Dm310s. The viability effect is easily detectable in comparison to the effect of the deletion of miR310s. The number of genes significantly misexpressed under the influence of Dp310s is 3-10 times greater than under Dm310s. Importantly, the numbers of predicted targets are similar between them. Expression analysis of the predicted target genes suggests that miRNAs may sometimes function to buffer fluctuations in the transcriptome output. After the buffering function has evolved, heterospecific combinations may cause adverse effects.

MicroRNAs (miRNAs) play an important role in posttranscriptional regulation of genes. We developed a method to predict human miRNAs without requiring cross-species conservation. We first identified lowly/moderately expressed tissue-selective genes using EST data and then identified overrepresented motifs of seven nucleotides in the 3' UTRs of these genes. Using these motifs as potential target sites of miRNAs, we recovered more than two-thirds of the known human miRNAs. We then used those motifs that did not match any known human miRNA seed region to infer novel miRNAs. We predicted 36 new human miRNA genes with 44 mature forms and 4 novel alternative mature forms of 2 known miRNA genes when a stringent criterion was used and many more novel miRNAs when a less stringent criterion was used. We tested the expression of 11 predicted miRNAs in three human cell lines and found 5 of them expressed in all three cell lines and 1 expressed in one cell line. We selected 2 of them, P-2 and P-27-5p, to do functional validation, using their mimics and inhibitors and using both luciferase assay and Western blotting. These experiments provided strong evidence that both P-2 and P-27-5p are novel miRNAs and that CREB3L3, which encodes cAMP-responsive element binding protein 3-like 3, is a target gene of P-2, whereas LAMB3, which encodes laminin beta3, is a target gene of P-27-5p.

Accurate pre-mRNA splicing requires primary splicing signals, including the splice sites, a polypyrimidine tract, and a branch site, other splicing-regulatory elements (SREs). The SREs include exonic splicing enhancers (ESEs), exonic splicing silencers (ESSs), intronic splicing enhancers (ISEs), and intronic splicing silencers (ISSs), which are typically located near the splice sites. However, it is unclear to what extent splicing-driven selective pressure constrains exonic and intronic sequences, especially those distant from the splice sites. Here, we studied the distribution of SREs in human genes in terms of DNA strand-asymmetry patterns. Under a neutral evolution model, each mononucleotide or oligonucleotide should have a symmetric (Chargaff's second parity rule), or weakly asymmetric yet uniform, distribution throughout a pre-mRNA transcript. However, we found that large sets of unbiased, experimentally determined SREs show a distinct strand-asymmetry pattern that is inconsistent with the neutral evolution model, and reflects their functional roles in splicing. ESEs are selected in exons and depleted in introns and vice versa for ESSs. Surprisingly, this trend extends into deep intronic sequences, accounting for one third of the genome. Selection is detectable even at the mononucleotide level, so that the asymmetric base compositions of exons and introns are predictive of ESEs and ESSs. We developed a method that effectively predicts SREs based on strand asymmetry, expanding the current catalog of SREs. Our results suggest that human genes have been optimized for exon and intron discrimination through an RNA landscape shaped during evolution.

A gapped transcription factor-binding site (TFBS) contains one or more highly degenerate positions. Discovering gapped motifs is difficult, because allowing highly degenerate positions in a motif greatly enlarges the search space and complicates the discovery process. Here, we propose a method for discovering TFBSs, especially gapped motifs. We use ChIP-chip data to judge the binding strength of a TF to a putative target promoter and use orthologous sequences from related species to judge the degree of evolutionary conservation of a predicted TFBS. Candidate motifs are constructed by growing compact motif blocks and by concatenating two candidate blocks, allowing 0-15 degenerate positions in between. The resultant patterns are statistically evaluated for their ability to distinguish between target and nontarget genes. Then, a position-based ranking procedure is proposed to enhance the signals of true motifs by collecting position concurrences. Empirical tests on 32 known yeast TFBSs show that the method is highly accurate in identifying gapped motifs, outperforming current methods, and it also works well on ungapped motifs. Predictions on additional 54 TFs successfully discover 11 gapped and 38 ungapped motifs supported by literature. Our method achieves high sensitivity and specificity for predicting experimentally verified TFBSs.

The HA1 domain of HA, the major antigenic protein of influenza A viruses, contains all of the antigenic sites of HA and is under continual immune-driven selection. To resolve controversies on whether only a few or many residue sites of HA1 have undergone positive selection, whether positive selection at HA1 is continual or punctuated, and whether antigenic change is punctuated, we introduce an approach to analyze 2,248 HA1 sequences collected from 1968 to 2005. We identify 95 substitutions at 63 sites from 1968 to 2005 and show that each substitution occurred very rapidly. The rapid substitution and the fact that 57 of the 63 sites are antigenic sites indicate that hitchhiking plays a minor role and that most of these sites, many more than previously found, have undergone positive selection. Strikingly, 88 of the 95 substitutions occurred in groups, and multiple mutations at antigenic sites sped up the fixation process. Our results suggest that positive selection has been ongoing most of the time, not sporadic, and that multiple mutations at antigenic sites cumulatively enhance antigenic drift, indicating that antigenic change is less punctuated than recently proposed.

MicroRNAs (miRNAs) function as endogenous translational repressors of protein-coding genes in animals by binding to target sites in the 3' UTRs of mRNAs. Because a single nucleotide change in the sequence of a target site can affect miRNA regulation, naturally occurring SNPs in target sites are candidates for functional variation that may be of interest for biomedical applications and evolutionary studies. However, little is known to date about variation among humans at miRNAs and their target sites. In this study, we analyzed publicly available SNP data in context with miRNAs and their target sites throughout the human genome, and we found a relatively low level of variation in functional regions of miRNAs, but an appreciable level of variation at target sites. Approximately 400 SNPs were found at experimentally verified target sites or predicted target sites that are otherwise evolutionarily conserved across mammals. Moreover, approximately 250 SNPs potentially create novel target sites for miRNAs in humans. If some variants have functional effects, they might confer phenotypic differences among humans. Although the majority of these SNPs appear to be evolving under neutrality, interestingly, some of these SNPs are found at relatively high population frequencies even in experimentally verified targets, and a few variants are associated with atypically long-range haplotypes that may have been subject to recent positive selection.

Many Saccharomyces cerevisiae duplicate genes that were derived from an ancient whole-genome duplication (WGD) unexpectedly show a small synonymous divergence (K(S)), a higher sequence similarity to each other than to orthologues in Saccharomyces bayanus, or slow evolution compared with the orthologue in Kluyveromyces waltii, a non-WGD species. This decelerated evolution was attributed to gene conversion between duplicates. Using approximately 300 WGD gene pairs in four species and their orthologues in non-WGD species, we show that codon-usage bias and protein-sequence conservation are two important causes for decelerated evolution of duplicate genes, whereas gene conversion is effective only in the presence of strong codon-usage bias or protein-sequence conservation. Furthermore, we find that change in mutation pattern or in tDNA copy number changed codon-usage bias and increased the K(S) distance between K. waltii and S. cerevisiae. Intriguingly, some proteins showed fast evolution before the radiation of WGD species but little or no sequence divergence between orthologues and paralogues thereafter, indicating that functional conservation after the radiation may also be responsible for decelerated evolution in duplicates.

The question of how duplicate genes are retained in a population remains controversial. The duplication-degeneration-complementation model, which involves no positive selection, stipulates a higher retention rate of duplicate genes in a small population than in a large one. This model has been accepted by many evolutionists. However, we found considerably more retentions and fewer losses of duplicate genes in the mouse genome than in the human genome, although the population size of rodents is in general larger than that of primates. Indeed, in nearly every interval of synonymous divergence between duplicate genes, the number of gene retentions in mouse is larger than that in human. Our findings suggest a more important role of positive selection in duplicate retention than duplication-degeneration-complementation. In addition, certain functional categories show a higher tendency of lineage-specific expansion than expected, suggesting lineage-specific selection or functional bias in retained duplicates.

It is firmly believed that ancestral primates were nocturnal, with nocturnality having been maintained in most prosimian lineages. Under this traditional view, the opsin genes in all nocturnal prosimians should have undergone similar degrees of functional relaxation and accumulated similar extents of deleterious mutations. This expectation is rejected by the short-wavelength (S) opsin gene sequences from 14 representative prosimians. We found severe defects of the S opsin gene only in lorisiforms, but no defect in five nocturnal and two diurnal lemur species and only minor defects in two tarsiers and two nocturnal lemurs. Further, the nonsynonymous-to-synonymous rate ratio of the S opsin gene is highest in the lorisiforms and varies among the other prosimian branches, indicating different time periods of functional relaxation among lineages. These observations suggest that the ancestral primates were diurnal or cathemeral and that nocturnality has evolved several times in the prosimians, first in the lorisiforms but much later in other lineages. This view is further supported by the distribution pattern of the middle-wavelength (M) and long-wavelength (L) opsin genes among prosimians.

Knowing transcription factors (TFs) involved in the yeast cell cycle is helpful for understanding the regulation of yeast cell cycle genes. We therefore developed two methods for predicting (i) individual cell cycle TFs and (ii) synergistic TF pairs. The essential idea is that genes regulated by a cell cycle TF should have higher (lower, if it is a repressor) expression levels than genes not regulated by it during one or more phases of the cell cycle. This idea can also be used to identify synergistic interactions of TFs. Applying our methods to chromatin immunoprecipitation data and microarray data, we predict 50 cell cycle TFs and 80 synergistic TF pairs, including most known cell cycle TFs and synergistic TF pairs. Using these and published results, we describe the behaviors of 50 known or inferred cell cycle TFs in each cell cycle phase in terms of activation/repression and potential positive/negative interactions between TFs. In addition to the cell cycle, our methods are also applicable to other functions.

Although the evolutionary significance of gene duplication has long been recognized, it remains unclear what determines gene duplicability. We find protein complexity to be an important determinant because the proportion of unduplicated genes (P) increases with the number of subunits in a protein. However, P is high (>or=65%) for both monomers and multimers in yeast, but

To study the evolution of the yeast protein interaction network, we first classified yeast proteins by their evolutionary histories into isotemporal categories, then analyzed the interaction tendencies within and between the categories, and finally reconstructed the main growth path. We found that two proteins tend to interact with each other if they are in the same or similar categories, but tended to avoid each other otherwise, and that network evolution mirrors the universal tree of life. These observations suggest synergistic selection during network evolution and provide insights into the hierarchical modularity of cellular networks.

Human DNA sequence variation data are useful for studying the origin, evolution, and demographic history of modern humans and the mechanisms of maintenance of genetic variability in human populations, and for detecting linkage association of disease. Here, we report worldwide variation data from a approximately 10-kilobase noncoding autosomal region. We identified 75 variant sites in 64 humans (128 sequences) and 463 variant sites among the human, chimpanzee, and orangutan sequences. Statistical tests suggested that the region is selectively neutral. The average nucleotide diversity (pi) across the region was 0.088% among all of the human sequences obtained, 0.085% among African sequences, and 0.082% among non-African sequences, supporting the view of a low nucleotide diversity ( approximately 0.1%) in humans. The comparable pi value in non-Africans to that in Africans indicates no severe bottleneck during the evolution of modern non-Africans; however, the possibility of a mild bottleneck cannot be excluded because non-Africans showed considerably fewer variants than Africans. The present and two previous large data sets all show a strong excess of low frequency variants in comparison to that expected from an equilibrium population, indicating a relatively recent population expansion. The mutation rate was estimated to be 1.15 x 10(-9) per nucleotide per year. Estimates of the long-term effective population size N(e) by various statistical methods were similar to those in other studies. The age of the most recent common ancestor was estimated to be approximately 1.29 million years ago among all of the sequences obtained and approximately 634,000 years ago among the non-African sequences, providing the first evidence from a noncoding autosomal region for ancient human histories, even among non-Africans.

The Pax-6 gene encodes a transcription factor containing both a paired and a homeodomain and is highly conserved among Metazoa. In both vertebrates and invertebrates, Pax-6 is required for eye morphogenesis, development of parts of the central nervous system, and, in some phyla, for the development of olfactory sense organs. Ectopic expression of Pax-6 from insects, mammals, cephalopods, and ascidians induces ectopic eyes in Drosophila, suggesting that Pax-6 may be a universal master control gene for eye morphogenesis. Platyhelminthes are an ancient phylum, originating from the base of spiralian protostomes, that bear primitive eyes, consisting of a group of rhabdomeric photoreceptor cells enclosed in a cup of pigment cells. The analysis of Pax-6 and its expression pattern should provide insights into the ancestral function of Pax-6 in eye morphogenesis. We have identified the Pax-6 gene of the planarian Dugesia(G)tigrina (Platyhelminthes; Turbellaria; Tricladida). This gene shares significant sequence identity and conserved genomic organization with Pax-6 proteins from other phyla. Phylogenetic analysis indicates that it clusters with the other Pax-6 genes, but in the most basal position. DtPax-6 is expressed as a single transcript in both regenerating and fully grown eyes, and electron microscopy studies show strong expression in the perykarion of both photoreceptor and pigment cells. Very low levels of expression also are detectable in other body regions. Because a bona fide Pax-6 homolog so far has not been detected in diploblastic animals, we speculate that Pax-6 may be typical for triploblasts and that the appearance of additional Pax genes may have coincided with increasingly complex body plans.

It is known that the squirrel monkey, marmoset, and other related New World (NW) monkeys possess three high-frequency alleles at the single X-linked photopigment locus, and that the spectral sensitivity peaks of these alleles are within those delimited by the human red and green pigment genes. The three alleles in the squirrel monkey and marmoset have been sequenced previously. In this study, the three alleles were found and sequenced in the saki monkey, capuchin, and tamarin. Although the capuchin and tamarin belong to the same family as the squirrel monkey and marmoset, the saki monkey belongs to a different family and is one of the species that is most divergent from the squirrel monkey and marmoset, suggesting the presence of the triallelic system in many NW monkeys. The nucleotide sequences of these alleles from the five species studied indicate that gene conversion occurs frequently and has partially or completely homogenized intronic and exonic regions of the alleles in each species, making it appear that a triallelic system arose independently in each of the five species studied. Nevertheless, a detailed analysis suggests that the triallelic system arose only once in the NW monkey lineage, from a middle wavelength (green) opsin gene, and that the amino acid differences at functionally critical sites among alleles have been maintained by natural selection in NW monkeys for >20 million years. Moreover, the two X-linked opsin genes of howler monkeys (a NW monkey genus) were evidently derived from the incorporation of a middle (green) and a long wavelength (red) allele into one chromosome; these two genes together with the (autosomal) blue opsin gene would immediately enable even a male monkey to have trichromatic vision.

Estimation of evolutionary distances has always been a major issue in the study of molecular evolution because evolutionary distances are required for estimating the rate of evolution in a gene, the divergence dates between genes or organisms, and the relationships among genes or organisms. Other closely related issues are the estimation of the pattern of nucleotide substitution, the estimation of the degree of rate variation among sites in a DNA sequence, and statistical testing of the molecular clock hypothesis. Mathematical treatments of these problems are considerably simplified by the assumption of a stationary process in which the nucleotide compositions of the sequences under study have remained approximately constant over time, and there now exist fairly extensive studies of stationary models of nucleotide substitution, although some problems remain to be solved. Nonstationary models are much more complex, but significant progress has been recently made by the development of the paralinear and LogDet distances. This paper reviews recent studies on the above issues and reports results on correcting the estimation bias of evolutionary distances, the estimation of the pattern of nucleotide substitution, and the estimation of rate variation among the sites in a sequence.

Pax proteins are a family of transcription factors with a highly conserved paired domain; many members also contain a paired-type homeodomain and/or an octapeptide. Nine mammalian Pax genes are known and classified into four subgroups: Pax-1/9, Pax-2/5/8, Pax-3/7, and Pax-4/6. Most of these genes are involved in nervous system development. In particular, Pax-6 is a key regulator that controls eye development in vertebrates and Drosophila. Although the Pax-4/6 subgroup seems to be more closely related to Pax-2/5/8 than to Pax-3/7 or Pax-1/9, its evolutionary origin is unknown. We therefore searched for a Pax-6 homolog and related genes in Cnidaria, which is the lowest phylum of animals that possess a nervous system and eyes. A sea nettle (a jellyfish) genomic library was constructed and two pax genes (Pax-A and -B) were isolated and partially sequenced. Surprisingly, unlike most known Pax genes, the paired box in these two genes contains no intron. In addition, the complete cDNA sequences of hydra Pax-A and -B were obtained. Hydra Pax-B contains both the homeodomain and the octapeptide, whereas hydra Pax-A contains neither. DNA binding assays showed that sea nettle Pax-A and -B and hydra Pax-A paired domains bound to a Pax-5/6 site and a Pax-5 site, although hydra Pax-B paired domain bound neither. An alignment of all available paired domain sequences revealed two highly conserved regions, which cover the DNA binding contact positions. Phylogenetic analysis showed that Pax-A and especially Pax-B were more closely related to Pax-2/5/8 and Pax-4/6 than to Pax-1/9 or Pax-3/7 and that the Pax genes can be classified into two supergroups: Pax-A/Pax-B/Pax-2/5/8/4/6 and Pax-1/9/3/7. From this analysis and the gene structure, we propose that modern Pax-4/6 and Pax-2/5/8 genes evolved from an ancestral gene similar to cnidarian Pax-B, having both the homeodomain and the octapeptide.

As additivity is a very useful property for a distance measure, a general additive distance is proposed under the stationary time-reversible (SR) model of nucleotide substitution or, more generally, under the stationary, time-reversible, and rate variable (SRV) model, which allows rate variation among nucleotide sites. A method for estimating the mean distance and the sampling variance is developed. In addition, a method is developed for estimating the variance-covariance matrix of distances, which is useful for the statistical test of phylogenies and molecular clocks. Computer simulation shows (i) if the sequences are longer than, say, 1000 bp, the SR method is preferable to simpler methods; (ii) the SR method is robust against deviations from time-reversibility; (iii) when the rate varies among sites, the SRV method is much better than the SR method because the distance is seriously underestimated by the SR method; and (iv) our method for estimating the sampling variance is accurate for sequences longer than 500 bp. Finally, a test is constructed for testing whether DNA evolution follows a general Markovian model.

In humans and rodents the male-to-female ratio of mutation rate (alpha m) has been suggested to be extremely large, so that the process of nucleotide substitution is almost completely male-driven. However, our sequence data from the last intron of the X chromosome-linked (Zfx) and Y chromosome-linked (Zfy) zinc finger protein genes suggest that alpha m is only approximately 2 in rodents with a 95% confidence interval from 1 to 3. Moreover, from published data on oogenesis and spermatogenesis we estimate the male-to-female ratio of the number of germ cell divisions per generation to be approximately 2 in rodents, confirming our estimate of alpha m and suggesting that errors in DNA replication are the primary source of mutation. As the estimated alpha m for rodents is only one-third of our previous estimate of approximately 6 for higher primates, there appear to be generation-time effects--i.e., alpha m decreases with decreasing generation time.

Phylogenetic analysis of DNA sequences from primates, rodents, lagomorphs, artiodactyls, carnivores, and birds strongly suggests that the order Rodentia is an outgroup to the other four mammalian orders and that Artiodactyla and Carnivora belong to a superordinal clade. Further, there is strong evidence against the Glires concept, which unites Lagomorpha and Rodentia. The radiation among Lagomorpha, Primates, and Artiodactyla--Carnivora is very bush-like, but there is some evidence that Lagomorpha has branched off first. Thus, the branching sequence for these five orders of mammals seems to be Rodentia, Lagomorpha, Primates, Artiodactyla, and Carnivora. The branching date for Rodentia could be as early as 100 million years ago. The rate of nucleotide substitution in the rodent lineage is shown to be at least 1.5 times higher than those in the other four mammalian lineages.

Comparison of plant mitochondrial (mt), chloroplast (cp) and nuclear (n) DNA sequences shows that the silent substitution rate in mtDNA is less than one-third that in cpDNA, which in turn evolves only half as fast as plant nDNA. The slower rate in mtDNA than in cpDNA is probably due to a lower mutation rate. Silent substitution rates in plant and mammalian mtDNAs differ by one or two orders of magnitude, whereas the rates in nDNAs may be similar. In cpDNA, the rate of substitution both at synonymous sites and in noncoding sequences in the inverted repeat is greatly reduced in comparison to single-copy sequences. The rate of cpDNA evolution appears to have slowed in some dicot lineages following the monocot/dicot split, and the slowdown is more conspicuous at nonsynonymous sites than at synonymous sites.

When the coding regions of 11 genes from rodents (mouse or rat) and man are compared with those from another mammalian species (usually bovine), it is found that rodents evolve significantly faster than man. The ratio of the number of nucleotide substitutions in the rodent lineage to that in the human lineage since their divergence is 2.0 for synonymous substitutions and 1.3 for nonsynonymous substitutions. Rodents also evolve faster in the 5' and 3' untranslated regions of five different mRNAs; the ratios are 2.6 and 3.1, respectively. The numbers of nucleotide substitutions between members of the beta-globin gene family that were duplicated before the man-mouse split are also higher in mouse than in man. The difference is, again, greater for synonymous substitutions than for nonsynonymous substitutions. This tendency is more consistent with the neutralist view of molecular evolution than with the selectionist view. A simple explanation for the higher rates in rodents is that rodents have shorter generation times and, thus, higher mutation rates. The implication of our findings for the study of molecular phylogeny is discussed.

The divergence between monocots and dicots represents a major event in higher plant evolution, yet the date of its occurrence remains unknown because of the scarcity of relevant fossils. We have estimated this date by reconstructing phylogenetic trees from chloroplast DNA sequences, using two independent approaches: the rate of synonymous nucleotide substitution was calibrated from the divergence of maize, wheat, and rice, whereas the rate of nonsynonymous substitution was calibrated from the divergence of angiosperms and bryophytes. Both methods lead to an estimate of the monocot-dicot divergence at 200 million years (Myr) ago (with an uncertainty of about 40 Myr). This estimate is also supported by analyses of the nuclear genes encoding large and small subunit ribosomal RNAs. These results imply that the angiosperm lineage emerged in Jurassic-Triassic time, which considerably predates its appearance in the fossil record (approximately 120 Myr ago). We estimate the divergence between cycads and angiosperms to be approximately 340 Myr, which can be taken as an upper bound for the age of angiosperms.

Formulas are developed for the distribution of allele frequencies and the mean and variance of heterozygosity under mutation and selection pressures. In large populations, even slight selection drastically changes the shape of the distribution of allele frequencies and reduces heterozygosity. On the other hand, the number of rare alleles in a sample is much less affected by selection. Under genic selection, heterozygosity may decrease with increasing population size. As a test statistic, the variance of heterozygosity can be used to detect the presence of selection, though it is not efficient when selection is very slight.

A simple method is proposed for constructing phylogenetic trees from distance matrices. The procedure for constructing tree topologies is similar to that of the unweighted pair-group method (UPG method) but makes corrections for unequal rates of evolution among lineages. The procedure for estimating branch lengths is the same as that of the Fitch and Margoliash method (F-M method) except that it allows no negative branch lengths. The performance of the present procedure for the construction of tree topologies is compared with that of the UPG method, the F-M method, Farris' method, and the modified Farris method by using Tateno's simulation outputs for nucleotide sequence divergence and his results for the performances of the latter four methods [Tateno, Y. (1978) Dissertation (Univ. Texas, Houston, TX). In this limited comparison, the present method performs considerably better than the UPG method and the F-M method and about equally well as the last two methods. The present method appears to be preferable to the UPG method for analysis of data from populations that have not differentiated much. Indeed, an application of the present method to gene frequency data from some Amerindian populations gives a tree topology far more reasonable than that obtained by the UPG method.