Mapping Geological Events and Nitrogen Fixation Evolution Onto the Timetree of the Evolution of Nitrogen-Fixation Genes.

Nitrogen is essential for all organisms, but biological nitrogen fixation (BNF) occurs only in a small fraction of prokaryotes. Previous studies divided nitrogenase-gene-carrying prokaryotes into Groups I to IV and provided evidence that BNF first evolved in bacteria. This study constructed a timetree of the evolution of nitrogen-fixation genes and estimated that archaea evolved BNF much later than bacteria and that nitrogen-fixing cyanobacteria evolved later than 1,900 MYA, considerably younger than the previous estimate of 2,200 MYA. Moreover, Groups III and II/I diverged ∼2,280 MYA, after the Kenorland supercontinent breakup (∼2,500-2,100 MYA) and the Great Oxidation Event (∼2,400-2,100 MYA); Groups III and Vnf/Anf diverged ∼2,086 MYA, after the Yarrabubba impact (∼2,229 MYA); and Groups II and I diverged ∼1,920 MYA, after the Vredefort impact (∼2,023 MYA). In summary, this study provided a timescale of BNF events and discussed the possible effects of geological events on BNF evolution.

Regulators of early maize leaf development inferred from transcriptomes of laser capture microdissection (LCM)-isolated embryonic leaf cells.

The superior photosynthetic efficiency of C4 leaves over C3 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.

Horizontally Transferred DNA in the Genome of the Fungus Pyricularia oryzae is Associated With Repressive Histone Modifications.

Horizontal gene transfer (HGT) is a means of exchanging genetic material asexually. The process by which horizontally transferred genes are domesticated by the host genome is of great interest but is not well understood. In this study, we determined the telomere-to-telomere genome sequence of the wheat-infecting Pyricularia oryzae strain Br48. SNP analysis indicated that the Br48 strain is a hybrid of wheat- and Brachiaria-infecting strains by a sexual or parasexual cross. Comparative genomic analysis identified several megabase-scale "insertions" in the Br48 genome, some of which were possibly gained by HGT-related events from related species, such as P. pennisetigena or P. grisea. Notably, the mega-insertions often contained genes whose phylogeny is not congruent with the species phylogeny. Moreover, some of the genes have a close homolog even in distantly related organisms, such as basidiomycetes or prokaryotes, implying the involvement of multiple HGT events. Interestingly, the levels of the silent epigenetic marks H3K9me3 and H3K27me3 in a genomic region tended to be negatively correlated with the phylogenetic concordance of genes in the same region, suggesting that horizontally transferred DNA is preferentially targeted for epigenetic silencing. Indeed, the putative HGT-derived genes were activated when MoKmt6, the gene responsible for H3K27me3 modification, was deleted. Notably, these genes also tended to be up-regulated during infection, suggesting that they are now under host control and have contributed to establishing a fungal niche. In conclusion, this study suggests that epigenetic modifications have played an important role in the domestication of HGT-derived genes in the P. oryzae genome.

Discovering unknown human and mouse transcription factor binding sites and their characteristics from ChIP-seq data.

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.

Origin and Evolution of Nitrogen Fixation in Prokaryotes.

The origin of nitrogen fixation is an important issue in evolutionary biology. While nitrogen is required by all living organisms, only a small fraction of bacteria and archaea can fix nitrogen. The prevailing view is that nitrogen fixation first evolved in archaea and was later transferred to bacteria. However, nitrogen-fixing (Nif) bacteria are far larger in number and far more diverse in ecological niches than Nif archaea. We, therefore, propose the bacteria-first hypothesis, which postulates that nitrogen fixation first evolved in bacteria and was later transferred to archaea. As >30,000 prokaryotic genomes have been sequenced, we conduct an in-depth comparison of the two hypotheses. We first identify the six genes involved in nitrogen fixation in all sequenced prokaryotic genomes and then reconstruct phylogenetic trees using the six Nif proteins individually or in combination. In each of these trees, the earliest lineages are bacterial Nif protein sequences and in the oldest clade (group) the archaeal sequences are all nested inside bacterial sequences, suggesting that the Nif proteins first evolved in bacteria. The bacteria-first hypothesis is further supported by the observation that the majority of Nif archaea carry the major bacterial Mo (molybdenum) transporter (ModABC) rather than the archaeal Mo transporter (WtpABC). Moreover, in our phylogeny of all available ModA and WtpA protein sequences, the earliest lineages are bacterial sequences while archaeal sequences are nested inside bacterial sequences. Furthermore, the bacteria-first hypothesis is supported by available isotopic data. In conclusion, our study strongly supports the bacteria-first hypothesis.

Simultaneous detection of miRNA and mRNA at the single-cell level in plant tissues.

Detecting the simultaneous presence of a microRNA (miRNA) and a mRNA in a specific tissue can provide support for the prediction that the miRNA regulates the mRNA. Although two such methods have been developed for mammalian tissues, they have a low signal-noise ratio and/or poor resolution at the single-cell level. To overcome these drawbacks, we develop a method that uses sequence-specific miRNA-locked nucleic acid (LNA) and mRNA-LNA probes. Moreover, it augments the detection signal by rolling circle amplification, achieving a high signal-noise ratio at the single-cell level. Dot signals are counted for determining the expression levels of mRNA and miRNA molecules in specific cells. We show a high sequence specificity of our miRNA-LNA probe, revealing that it can discriminate single-base mismatches. Numerical quantification by our method is tested in transgenic rice lines with different gene expression levels. We conduct several applications. First, the spatial expression profiling of osa-miR156 and OsSPL12 in rice leaves reveals their specific expression in mesophyll cells. Second, studying rice and its mutant lines with our method reveals opposite expression patterns of miRNA and its target mRNA in tissues. Third, the dynamic expression profiles of ZmGRF8 and zma-miR396 during maize leaf development provide evidence that zma-miR396 regulates the preferential spatial expression of ZmGRF8 in bundle sheath cells. Finally, our method can be scaled up to simultaneously detect multiple miRNAs and mRNAs in a tissue. Thus, it is a sensitive and versatile technique for studying miRNA regulation of plant tissue development.

Chicken HOXC8 and HOXC10 genes may play a role in the altered skull morphology associated with the Crest phenotype.

One of the most intriguing traits found in domestic chickens is the Crest phenotype. This trait, characterized by a tuft of elongated feathers sprouted from the head, is found in breeds such as Polish chickens and Silkie chickens. Moreover, some crested chicken breeds also exhibit a protuberance in their anterodorsal skull region. Previous studies have strived to identify the causative factors of this trait. This study aimed to elucidate the role of chicken HOXC8 and HOXC10 in the formation of the Crest phenotype. We explored the effect of ectopic expression of HOXC8 or HOXC10 on the chicken craniofacial morphology using the RCAS retrovirus transformation system. Microcomputed tomography scanning was conducted to measure the 3D structure of the cranial bone of transgenic embryos for geometric morphometric analysis. We found that the ectopic expression of HOXC8 or HOXC10 in chicken heads caused mild morphological changes in the skull compared with the GFP-transgenic control group. Geometric morphometric analysis showed that HOXC8 and HOXC10 transgenic groups expressed a mild upward shape change in the frontal region of the skull compared with the control group, which is similar to what is seen in the crested chicken breeds. In conclusion, this study supports findings in previous studies in which HOX genes play a role in the formation of the altered skull morphology related to the Crest phenotype. It also supports that mutations in HOX genes may contribute to intra- and inter-specific variation in morphological traits in vertebrates.

Maize Golden2-like transcription factors boost rice chloroplast development, photosynthesis, and grain yield.

Chloroplasts are the sites for photosynthesis, and two Golden2-like factors act as transcriptional activators of chloroplast development in rice (Oryza sativa L.) and maize (Zea mays L.). Rice OsGLK1 and OsGLK2 are orthologous to maize ZmGLK1 (ZmG1) and ZmGLK2 (ZmG2), respectively. However, while rice OsGLK1 and OsGLK2 act redundantly to regulate chloroplast development in mesophyll cells, maize ZmG1 and ZmG2 are functionally specialized and expressed in different cell-specific manners. To boost rice chloroplast development and photosynthesis, we generated transgenic rice plants overexpressing ZmG1 and ZmG2, individually or simultaneously, with constitutive promoters (pZmUbi::ZmG1 and p35S::ZmG2) or maize promoters (pZmG1::ZmG1, pZmG2::ZmG2, and pZmG1::ZmG1/pZmG2::ZmG2). Both ZmG1 and ZmG2 genes were highly expressed in transgenic rice leaves. Moreover, ZmG1 and ZmG2 showed coordinated expression in pZmG1::ZmG1/pZmG2::ZmG2 plants. All Golden2-like (GLK) transgenic plants had higher chlorophyll and protein contents, Rubisco activities and photosynthetic rates per unit leaf area in flag leaves. However, the highest grain yields occurred when maize promoters were used; pZmG1::ZmG1, pZmG2::ZmG2, and pZmG1::ZmG1/pZmG2::ZmG2 transgenic plants showed increases in grain yield by 51%, 47%, and 70%, respectively. In contrast, the pZmUbi::ZmG1 plant produced smaller seeds without yield increases. Transcriptome analysis indicated that maize GLKs act as master regulators promoting the expression of both photosynthesis-related and stress-responsive regulatory genes in both rice shoot and root. Thus, by promoting these important functions under the control of their own promoters, maize GLK1 and GLK2 genes together dramatically improved rice photosynthetic performance and productivity. A similar approach can potentially improve the productivity of many other crops.

Genome analysis of SARS-CoV-2 isolates from a population reveals the rapid selective sweep of a haplotype carrying many pre-existing and new mutations.

To understand the mechanism underlying the evolution of SARS-CoV-2 in a population, we sequenced 92 viral genomes from Assam, India. Analysis of these and database sequences revealed a complete selective sweep of a haplotype in Assam carrying 13 pre-existing variants, including a high leap in frequency of a variant on ORF8, which is involved in immune evasion. A comparative study between sequences of same lineage and similar time frames in and outside Assam showed that 10 of the 13 pre-existing variants had a frequency ranging from 96 to 99%, and the remaining 3 had a low frequency outside Assam. Using a phylogenetic approach to infer sequential occurrences of variants we found that the variant Phe120del on ORF8, which had a low frequency (1.75%) outside Assam, is at the base of the phylogenetic tree of variants and became totally fixed (100%) in Assam population. Based on this observation, we inferred that the variant on ORF8 had a selective advantage, so it carried the haplotype to reach the100% frequency. The haplotype also carried 32 pre-existing variants at a frequency from 1.00 to 80.00% outside Assam. Those of these variants that are more closely linked to the S-protein locus, which often carries advantageous mutations and is tightly linked to the ORF8 locus, retained higher frequencies, while the less tightly linked variants showed lower frequencies, likely due to recombination among co- circulating variants in Assam. The ratios of non-synonymous substitutions to synonymous substitutions suggested that some genes such as those coding for the S-protein and non-structural proteins underwent positive selection while others were subject to purifying selection during their evolution in Assam. Furthermore, we observed negative correlation of the Ct value of qRT-PCR of the patients with abundant ORF6 transcripts, suggesting that ORF6 can be used as a marker for estimating viral titer. In conclusion, our in-depth analysis of SARS-CoV-2 genomes in a regional population reveals the mechanism and dynamics of viral evolution.

Lactobacillus plantarum isolated from kefir enhances immune responses and survival of white shrimp (Penaeus vannamei) challenged with Vibrio alginolyticus.

Lactobacillus plantarum is known for its probiotics benefit to host, although the effects vary among strains. This study conducted a feeding experiment of three Lactobacillus strains, MRS8, MRS18 and MRS20, which were isolated from kefir and incorporated into the diets of shrimp to evaluate the effects of non-specific immunity, immune-related gene expression, and disease resistance of white shrimp (Penaeus vannamei) against Vibrio alginolyticus. To prepare the experimental feed groups, the basic feed was mixed with different concentrations of L. plantarum strains MRS8, MRS18, and MRS 20, which were incorporated at 0 CFU (control), 1 × 106 CFU (groups 8-6, 18-6, and 20-6), and 1 × 109 CFU (groups 8-9, 18-9, and 20-9) per gram of diet for an in vivo assay. During the rearing period for 28 days of feeding each group, immune responses, namely the total hemocyte count (THC), phagocytic rate (PR), phenoloxidase activity, and respiratory burst were examined on days 0, 1, 4, 7, 14, and 28. The results showed that groups 20-6, 18-9 and 20-9 improved THC, and groups 18-9 and 20-9 improved phenoloxidase activity and respiratory burst as well. The expression of immunity-related genes was also examined. Group 8-9 increased the expression of LGBP, penaeidin 2 (PEN2) and CP, group 18-9 increased the expression of proPO1, ALF, Lysozyme, penaeidin 3 (PEN3) and SOD, and group 20-9 increased the expression of LGBP, ALF, crustin, PEN2, PEN3, penaeidin 4 (PEN4) and CP (p < 0.05). Groups 18-6, 18-9, 2-6, and 20-9 were further used in the challenge test. After feeding for 7 days and 14 days, Vibrio alginolyticus was injected into white shrimp and observed the shrimp survival for 168 h. The results showed that compared to the control, all groups improved the survival rate. Especially, feeding group 18-9 for 14 days improved the survival rate of white shrimp (p < 0.05). After the challenge test for 14 days, the midgut DNA of survival white shrimps was extracted to analyze the colonization of L. plantarum. Among the groups, (6.61 ± 3.58) × 105 CFU/pre shrimp of L. plantarum in feeding group 18-9 and (5.86 ± 2.27) × 105 CFU/pre shrimp in group 20-9 were evaluated by qPCR. Taken together, group 18-9 had the best effects on the non-specific immunity, the immune-related gene expression, and the disease resistance, which might be due to the benefit of the probiotic colonization.

Constructing a yeast to express the largest cellulosome complex on the cell surface.

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-ScGPI and fluorescence-activated cell sorting analysis results revealed that almost 81% of yeast cells displayed OlpB-ScGPI. 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.

Maize ANT1 modulates vascular development, chloroplast development, photosynthesis, and plant growth.

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.

Behavioral and brain- transcriptomic synchronization between the two opponents of a fighting pair of the fish Betta splendens.

Conspecific male animals fight for resources such as food and mating opportunities but typically stop fighting after assessing their relative fighting abilities to avoid serious injuries. Physiologically, how the fighting behavior is controlled remains unknown. Using the fighting fish Betta splendens, we studied behavioral and brain-transcriptomic changes during the fight between the two opponents. At the behavioral level, surface-breathing, and biting/striking occurred only during intervals between mouth-locking. Eventually, the behaviors of the two opponents became synchronized, with each pair showing a unique behavioral pattern. At the physiological level, we examined the expression patterns of 23,306 brain transcripts using RNA-sequencing data from brains of fighting pairs after a 20-min (D20) and a 60-min (D60) fight. The two opponents in each D60 fighting pair showed a strong gene expression correlation, whereas those in D20 fighting pairs showed a weak correlation. Moreover, each fighting pair in the D60 group showed pair-specific gene expression patterns in a grade of membership analysis (GoM) and were grouped as a pair in the heatmap clustering. The observed pair-specific individualization in brain-transcriptomic synchronization (PIBS) suggested that this synchronization provides a physiological basis for the behavioral synchronization. An analysis using the synchronized genes in fighting pairs of the D60 group found genes enriched for ion transport, synaptic function, and learning and memory. Brain-transcriptomic synchronization could be a general phenomenon and may provide a new cornerstone with which to investigate coordinating and sustaining social interactions between two interacting partners of vertebrates.

Whole-Genome Duplication Facilitated the Evolution of C4 Photosynthesis in Gynandropsis gynandra.

In higher plants, whole-genome duplication (WGD) is thought to facilitate the evolution of C4 photosynthesis from C3 photosynthesis. To understand this issue, we used new and existing leaf-development transcriptomes to construct two coding sequence databases for C4Gynandropsis gynandra and C3Tarenaya hassleriana, which shared a WGD before their divergence. We compared duplicated genes in the two species and found that the WGD contributed to four aspects of the evolution of C4 photosynthesis in G. gynandra. First, G. gynandra has retained the duplicates of ALAAT (alanine aminotransferase) and GOGAT (glutamine oxoglutarate aminotransferase) for nitrogen recycling to establish a photorespiratory CO2 pump in bundle sheath (BS) cells for increasing photosynthesis efficiency, suggesting that G. gynandra experienced a C3-C4 intermediate stage during the C4 evolution. Second, G. gynandra has retained almost all known vein-development-related paralogous genes derived from the WGD event, likely contributing to the high vein complexity of G. gynandra. Third, the WGD facilitated the evolution of C4 enzyme genes and their recruitment into the C4 pathway. Fourth, several genes encoding photosystem I proteins were derived from the WGD and are upregulated in G. gynandra, likely enabling the NADH dehydrogenase-like complex to produce extra ATPs for the C4 CO2 concentration mechanism. Thus, the WGD apparently played an enabler role in the evolution of C4 photosynthesis in G. gynandra. Importantly, an ALAAT duplicate became highly expressed in BS cells in G. gynandra, facilitating nitrogen recycling and transition to the C4 cycle. This study revealed how WDG may facilitate C4 photosynthesis evolution.

Identifying Primate ACE2 Variants That Confer Resistance to SARS-CoV-2.

SARS-CoV-2 infects humans through the binding of viral S-protein (spike protein) to human angiotensin I converting enzyme 2 (ACE2). The structure of the ACE2-S-protein complex has been deciphered and we focused on the 27 ACE2 residues that bind to S-protein. From human sequence databases, we identified nine ACE2 variants at ACE2-S-protein binding sites. We used both experimental assays and protein structure analysis to evaluate the effect of each variant on the binding affinity of ACE2 to S-protein. We found one variant causing complete binding disruption, two and three variants, respectively, strongly and mildly reducing the binding affinity, and two variants strongly enhancing the binding affinity. We then collected the ACE2 gene sequences from 57 nonhuman primates. Among the 6 apes and 20 Old World monkeys (OWMs) studied, we found no new variants. In contrast, all 11 New World monkeys (NWMs) studied share four variants each causing a strong reduction in binding affinity, the Philippine tarsier also possesses three such variants, and 18 of the 19 prosimian species studied share one variant causing a strong reduction in binding affinity. Moreover, one OWM and three prosimian variants increased binding affinity by >50%. Based on these findings, we proposed that the common ancestor of primates was strongly resistant to and that of NWMs was completely resistant to SARS-CoV-2 and so is the Philippine tarsier, whereas apes and OWMs, like most humans, are susceptible. This study increases our understanding of the differences in susceptibility to SARS-CoV-2 infection among primates.

Many human RNA viruses show extraordinarily stringent selective constraints on protein evolution.

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 (dN/dS) in protein-coding genes of human RNA and DNA viruses and mammals. Among the 21 RNA viruses studied, 18 showed a genome-average dN/dS from 0.01 to 0.10, indicating that over 90% of nonsynonymous mutations are eliminated by negative selection. Only HIV-1 showed a dN/dS (0.31) higher than that (0.22) in mammalian genes. By comparing the dN/dS 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 dN/dS 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 dN/dS 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 dN/dS ratios than non-vector-borne viruses. In addition, we found a correlation between dN and dS, implying a correlation between dN and the mutation rate. Interestingly, only 2 of the 8 DNA viruses studied showed a dN/dS < 0.10, while 4 showed a dN/dS > 0.22. These observations increase our understanding of the mechanisms of RNA virus evolution.

Large-scale identification and functional analysis of NLR genes in blast resistance in the Tetep rice genome sequence.

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 NLRs showed resistance to 1 or more pathogen strains and each strain was recognized by multiple NLRs. However, few NLRs showed resistance to >6 strains, so multiple NLRs 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 NLRs. In developing a method to identify NLR pairs each of which functions as a unit, we found that >20% of the NLRs 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 NLRs 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.

Comparative transcriptomics method to infer gene coexpression networks and its applications to maize and rice leaf transcriptomes.

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 C4 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 C4 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 C4 photosynthesis.

Chromosomal-level genome assembly of the semi-dwarf rice Taichung Native 1, an initiator of Green Revolution.

Here we report the 409.5 Mb chromosome-level assembly of the first bred semi-dwarf rice, the Taichung Native 1 (TN1), which served as the template for the development of the Green Revolution (GR) cultivar IR8 "miracle rice". We sequenced the TN1 genome utilizing multiple platforms and produced PacBio long reads, Illumina paired-end reads, Illumina mate-pair reads and 10x Genomics linked reads. We used a hybrid approach to assemble the 226× coverage of sequences by a combination of de novo and reference-guided approaches. The assembled TN1 genome has an N50 scaffold size of 33.1 Mb with the longest measuring 45.5 Mb. We annotated 37,526 genes, in which 24,102 (64.23%) were assigned Blast2GO annotations. The genome has 4672 or 95.4% complete BUSCOs and a repeat content of 51.52%. We developed our own method of creating a GR pangenome using the orthologous relationships of the proteins of TN1, IR8, MH63 and IR64, identifying 16,999 core orthologue groups of Green Revolution. From the pangenome, we identified a set of shared and unique gene ontology terms for the accessory clusters, characterizing TN1, IR8, MH63 and IR64. This TN1 genome assembly and GR pangenome will be a resource for new genomic discoveries about Green Revolution, and for improving the disease and insect resistances and the yield of rice.

Rice transcription factor GAMYB modulates bHLH142 and is homeostatically regulated by TDR during anther tapetal and pollen development.

GIBBERELLIN MYB GENE (GAMYB), UNDEVELOPED TAPETUM1 (UDT1), TDR INTERACTING PROTEIN2 (TIP2/bHLH142), TAPETUM DEGENERATION RETARDATION (TDR), and ETERNAL TAPETUM 1/DELAYED TAPETUM DEGENERATION (EAT1/DTD) are important transcription factors that play a crucial role during pollen development in rice. This study demonstrates that bHLH142 acts downstream of UDT1 and GAMYB and works as a 'hub' in these two pollen pathways. We show that GAMYB modulates bHLH142 expression through specific binding to the MYB motif of the bHLH142 promoter during the early stage of pollen development, while TDR acts as a transcriptional repressor of the GAMYB modulation of bHLH142 by binding to the E-box close to the MYB motif on the promoter. Altered expression of these transcription factors highlights that a tight, precise, and coordinated regulation among them is essential for normal pollen development. Most notably, we show that the regulatory pathways of GAMYB and UDT1 rely on bHLH142 in a direct and indirect manner, respectively, and function in different tissues with distinct biological roles during pollen development. This study advances our understanding of the molecular mechanisms of rice pollen development.